Return to r_mod.c CVS log

Up to [Research Unix] / researchv10no / libF77

Annotation of researchv10no/libF77/r_mod.c, revision 1.1

1.1     ! root        1: #include "f2c.h"
        !             2: 
        !             3: #ifdef KR_headers
        !             4: #ifdef IEEE_drem
        !             5: double drem();
        !             6: #else
        !             7: double floor();
        !             8: #endif
        !             9: double r_mod(x,y) real *x, *y;
        !            10: #else
        !            11: #ifdef IEEE_drem
        !            12: double drem(double, double);
        !            13: #else
        !            14: #undef abs
        !            15: #include "math.h"
        !            16: #endif
        !            17: double r_mod(real *x, real *y)
        !            18: #endif
        !            19: {
        !            20: #ifdef IEEE_drem
        !            21:        double xa, ya, z;
        !            22:        if ((ya = *y) < 0.)
        !            23:                ya = -ya;
        !            24:        z = drem(xa = *x, ya);
        !            25:        if (xa > 0) {
        !            26:                if (z < 0)
        !            27:                        z += ya;
        !            28:                }
        !            29:        else if (z > 0)
        !            30:                z -= ya;
        !            31:        return z;
        !            32: #else
        !            33:        double quotient;
        !            34:        if( (quotient = (double)*x / *y) >= 0)
        !            35:                quotient = floor(quotient);
        !            36:        else
        !            37:                quotient = -floor(-quotient);
        !            38:        return(*x - (*y) * quotient );
        !            39: #endif
        !            40: }

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