Return to r_tan.c CVS log

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

Annotation of 43BSDReno/lib/libF77/r_tan.c, revision 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_tan.c     5.4     2/25/90
        !             7:  */
        !             8: 
        !             9: #ifndef tahoe
        !            10: float r_tan(x)
        !            11: float *x;
        !            12: {
        !            13: double tan();
        !            14: return( tan(*x) );
        !            15: }
        !            16: 
        !            17: #else tahoe
        !            18: 
        !            19: /*
        !            20:        SINGLE PRECISION floating point tangent
        !            21: 
        !            22:        sin/cos is used after argument reduction to [0,pi/4] range.
        !            23:        since x is in this range, tan(x) is in [0,1] range and
        !            24:        no overflow can occur here.
        !            25: */
        !            26: 
        !            27: #include <errno.h>
        !            28: 
        !            29: int    errno;
        !            30: static double invpi = 1.27323954473516268;  /* 4/pi */
        !            31: 
        !            32: #ifndef __GNUC__
        !            33: fortran float sin(), cos();
        !            34: #else
        !            35: #define        sin(x) \
        !            36: ({ \
        !            37:        float __result; \
        !            38:        asm("ldd %1; cvdf; sinf; stf %0" : "=rm" (__result) : "rm" (x)); \
        !            39:        __result; \
        !            40: })
        !            41: #define        cos(x) \
        !            42: ({ \
        !            43:        float __result; \
        !            44:        asm("ldd %1; cvdf; cosf; stf %0" : "=rm" (__result) : "rm" (x)); \
        !            45:        __result; \
        !            46: })
        !            47: #endif
        !            48: 
        !            49: float
        !            50: r_tan(parg)
        !            51: float *parg;
        !            52: {
        !            53:        double arg;
        !            54:        double modf();
        !            55:        float flmax_();
        !            56:        double temp, e, x, xsq;
        !            57:        int sign;
        !            58:        int flag, i;
        !            59: 
        !            60:        arg = *parg;
        !            61:        flag = 0;
        !            62:        sign = 1.;
        !            63:        if(arg < 0.){           /* tan(-arg) = -tan(arg) */
        !            64:                arg = -arg;
        !            65:                sign = -1.;
        !            66:        }
        !            67:        arg = arg*invpi;   /*overflow?*/
        !            68:        x = modf(arg,&e);
        !            69:        i = e;
        !            70:        switch(i%4) {
        !            71:        case 1:                 /* 2nd octant: tan(x) = 1/tan(1-x) */
        !            72:                x = 1. - x;
        !            73:                flag = 1;
        !            74:                break;
        !            75: 
        !            76:        case 2:                 /* 3rd octant: tan(x) = -1/tan(x) */
        !            77:                sign = - sign;
        !            78:                flag = 1;
        !            79:                break;
        !            80: 
        !            81:        case 3:                 /* 4th octant: tan(x) = -tan(1-x) */
        !            82:                x = 1. - x;
        !            83:                sign = - sign;
        !            84:                break;
        !            85: 
        !            86:        case 0:                 /* 1st octant */
        !            87:                break;
        !            88:        }
        !            89:        x = x/invpi;
        !            90: 
        !            91:        temp = sin(x)/cos(x);
        !            92: 
        !            93:        if(flag == 1) {
        !            94:                if(temp == 0.) {        /* check for singular "point" */
        !            95:                        errno = ERANGE;
        !            96:                        if (sign>0)
        !            97:                                return(flmax_());
        !            98:                        return(-flmax_());
        !            99:                }
        !           100:                temp = 1./temp;
        !           101:        }
        !           102:        return(sign*temp);
        !           103: }
        !           104: 
        !           105: #endif tahoe

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