BSD 4_3_Tahoe release

[unix-history] / usr / src / usr.lib / libF77 / r_mod.c

diff --git a/usr/src/usr.lib/libF77/r_mod.c b/usr/src/usr.lib/libF77/r_mod.c
index d82a721..55408da 100644 (file)
--- a/usr/src/usr.lib/libF77/r_mod.c
+++ b/usr/src/usr.lib/libF77/r_mod.c
@@ -1,8 +1,15 @@
 /*
 /*

- *     "@(#)r_mod.c    1.1"
+ * Copyright (c) 1980 Regents of the University of California.
+ * All rights reserved.  The Berkeley software License Agreement
+ * specifies the terms and conditions for redistribution.
+ *
+ *     @(#)r_mod.c     5.4     11/3/86

  */
 
  */
 

-double r_mod(x,y)
+#ifndef tahoe
+float flt_retval;
+
+float r_mod(x,y)

 float *x, *y;
 {
 double floor(), quotient = *x / *y;
 float *x, *y;
 {
 double floor(), quotient = *x / *y;


@@ -10,5 +17,27 @@ if (quotient >= 0.0)
        quotient = floor(quotient);
 else
        quotient = -floor(-quotient);
        quotient = floor(quotient);
 else
        quotient = -floor(-quotient);

+flt_retval = *x - (*y) * quotient ;
+return(flt_retval);
+}
+
+#else tahoe
+
+/*   THIS IS BASED ON THE TAHOE REPR. FOR FLOATING POINT */
+#include <tahoemath/FP.h>
+
+double r_mod(x,y)
+float *x, *y;
+{
+double floor(), quotient = *x / *y;
+if (quotient >= 0.0)
+       quotient = floor(quotient);
+else {
+       *(unsigned long *)&quotient ^= SIGN_BIT;
+       quotient = floor(quotient);
+       if (quotient != 0)
+               *(unsigned long *)&quotient ^= SIGN_BIT;
+       }

 return(*x - (*y) * quotient );
 }
 return(*x - (*y) * quotient );
 }

+#endif tahoe