Return to r_mod.c CVS log

Up to [CSRG BSD Unix] / 43BSDReno / lib / libF77

Annotation of 43BSDReno/lib/libF77/r_mod.c, revision 1.1.1.1

1.1       root        1: /*
                      2:  * Copyright (c) 1980 Regents of the University of California.
                      3:  * All rights reserved.  The Berkeley software License Agreement
                      4:  * specifies the terms and conditions for redistribution.
                      5:  *
                      6:  *     @(#)r_mod.c     5.4     11/3/86
                      7:  */
                      8: 
                      9: #ifndef tahoe
                     10: float flt_retval;
                     11: 
                     12: float r_mod(x,y)
                     13: float *x, *y;
                     14: {
                     15: double floor(), quotient = *x / *y;
                     16: if (quotient >= 0.0)
                     17:        quotient = floor(quotient);
                     18: else
                     19:        quotient = -floor(-quotient);
                     20: flt_retval = *x - (*y) * quotient ;
                     21: return(flt_retval);
                     22: }
                     23: 
                     24: #else tahoe
                     25: 
                     26: /*   THIS IS BASED ON THE TAHOE REPR. FOR FLOATING POINT */
                     27: #include <tahoemath/FP.h>
                     28: 
                     29: double r_mod(x,y)
                     30: float *x, *y;
                     31: {
                     32: double floor(), quotient = *x / *y;
                     33: if (quotient >= 0.0)
                     34:        quotient = floor(quotient);
                     35: else {
                     36:        *(unsigned long *)&quotient ^= SIGN_BIT;
                     37:        quotient = floor(quotient);
                     38:        if (quotient != 0)
                     39:                *(unsigned long *)&quotient ^= SIGN_BIT;
                     40:        }
                     41: return(*x - (*y) * quotient );
                     42: }
                     43: #endif tahoe

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.