Return to psin.c CVS log

Up to [Research Unix] / researchv10no / cmd / pico

Annotation of researchv10no/cmd/pico/psin.c, revision 1.1

1.1     ! root        1: static
        !             2: short costab[91]={
        !             3:        10000,  9998,   9993,   9986,   9975,
        !             4:        9961,   9945,   9925,   9902,   9876,
        !             5:        9848,   9816,   9781,   9743,   9702,
        !             6:        9659,   9612,   9563,   9510,   9455,
        !             7:        9396,   9335,   9271,   9205,   9135,
        !             8:        9063,   8987,   8910,   8829,   8746,
        !             9:        8660,   8571,   8480,   8386,   8290,
        !            10:        8191,   8090,   7986,   7880,   7771,
        !            11:        7660,   7547,   7431,   7313,   7193,
        !            12:        7071,   6946,   6819,   6691,   6560,
        !            13:        6427,   6293,   6156,   6018,   5877,
        !            14:        5735,   5591,   5446,   5299,   5150,
        !            15:        5000,   4848,   4694,   4539,   4383,
        !            16:        4226,   4067,   3907,   3746,   3583,
        !            17:        3420,   3255,   3090,   2923,   2756,
        !            18:        2588,   2419,   2249,   2079,   1908,
        !            19:        1736,   1564,   1391,   1218,   1045,
        !            20:        871,    697,    523,    348,    174,
        !            21:        0,
        !            22: };
        !            23: 
        !            24: pcos(x)
        !            25:        register x;
        !            26: {
        !            27:        x %= 360;
        !            28:        while(x<0)
        !            29:                x+=360;
        !            30:        if(x<=180)
        !            31:                return(x<90? costab[x] : -costab[180-x]);
        !            32:        return(x<180+90? -costab[x-180] : costab[360-x]);
        !            33: }
        !            34: 
        !            35: psin(x)
        !            36:        register x;
        !            37: {
        !            38:        return(pcos(x-90));
        !            39: }
        !            40: 
        !            41: Psin(x)
        !            42:        register x;
        !            43: {
        !            44:        return(Pcos(x-9000));
        !            45: }
        !            46: 
        !            47: Pcos(x)
        !            48:        register x;
        !            49: {
        !            50:        register int a = x/100;
        !            51:        register int b = x%100;
        !            52:        register int c = pcos(a);
        !            53:        register int d = pcos(a+1);
        !            54: 
        !            55:        a = c + (b*(d-c))/100;
        !            56:        return a;
        !            57: }

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