Return to cos.c CVS log

Up to [Research Unix] / researchv9 / cmd / cfront / libC / complex

Annotation of researchv9/cmd/cfront/libC/complex/cos.c, revision 1.1.1.1

1.1       root        1: 
                      2: #include "complex.h"
                      3: #include "const.h"
                      4: 
                      5: complex cos(complex z)
                      6: /*
                      7:        The cosine of z: cos(z)=cosh(iz).
                      8: */
                      9: {
                     10:        complex y(-z.im, z.re); 
                     11:        return cosh(y);
                     12: }
                     13: 
                     14: 
                     15: 
                     16: complex cosh(complex z)
                     17: /*
                     18:        The complex hyperbolic cosine of z
                     19: */
                     20: {
                     21:        double  cosh_x, sinh_x, cos_y, sin_y;
                     22: 
                     23: #define        COSH_GOOD       1e7
                     24: #define        COSH_HUGE       1e38
                     25: 
                     26:        if (z.re > MAX_EXPONENT) {
                     27:                complex_error(C_COSH_RE, z.re);
                     28:                cosh_x = sinh_x = COSH_HUGE;
                     29:        }
                     30:        else if (z.re < MIN_EXPONENT) {
                     31:                complex_error(C_COSH_RE, z.re);
                     32:                cosh_x = COSH_HUGE;
                     33:                sinh_x = -COSH_HUGE;
                     34:        }
                     35:        else {
                     36:                double pos_exp = exp(z.re); 
                     37:                double neg_exp = 1/pos_exp;
                     38:                cosh_x = (pos_exp + neg_exp)/2;
                     39:                sinh_x = (pos_exp - neg_exp)/2;
                     40:        }
                     41: 
                     42:        if (ABS(z.im) > COSH_GOOD) {
                     43:                complex_error(C_COSH_IM, z.im);
                     44:                cos_y = sin_y = 0;
                     45:        }
                     46:        else {
                     47:                cos_y = cos(z.im);
                     48:                sin_y = sin(z.im);
                     49:        }
                     50: 
                     51:        return complex(cos_y*cosh_x, sin_y*sinh_x);
                     52: }

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