Return to sin.c CVS log

Up to [MW Coherent from dump] / coherent / b / lib / libm

Annotation of coherent/b/lib/libm/sin.c, revision 1.1

1.1     ! root        1: /*
        !             2:  * libm/sin.c
        !             3:  * Evaluate the sine function.
        !             4:  */
        !             5: 
        !             6: #include <math.h>
        !             7: 
        !             8: #if    EMU87
        !             9: #include "emumath.h"
        !            10: #endif
        !            11: 
        !            12: #define        THRESHOLD       1.111111111111111111e-3
        !            13: 
        !            14: /*
        !            15:  * When |x| is very small, use the series
        !            16:  *     x - x**3 / 6 + x**5 / 120 - ...
        !            17:  * to get an accurate result while avoiding loss of precision.
        !            18:  * The THRESHOLD value above is determined empirically.
        !            19:  */
        !            20: static readonly double sintab[] ={
        !            21:         0.000000000000000000000000,
        !            22:         1.000000000000000000000000,
        !            23:         0.000000000000000000000000,
        !            24:        -0.166666666666666666666666,
        !            25:         0.000000000000000000000000,
        !            26:         0.008333333333333333333333
        !            27: };
        !            28: 
        !            29: double
        !            30: sin(x) double x;
        !            31: {
        !            32:        if (fabs(x) < THRESHOLD)
        !            33:                return _pol(x, sintab, 6);
        !            34:        return cos(PI/2.0 - x);
        !            35: }
        !            36: 
        !            37: /* end of sin.c */

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