Return to pow_zz.c CVS log

Up to [CSRG BSD Unix] / 42BSD / usr.lib / libF77

Annotation of 42BSD/usr.lib/libF77/pow_zz.c, revision 1.1.1.1

1.1       root        1: /*
                      2:  *     "@(#)pow_zz.c   1.1"
                      3:  */
                      4: 
                      5: #include "complex"
                      6: 
                      7: pow_zz(r,a,b)
                      8: dcomplex *r, *a, *b;
                      9: {
                     10: double logr, logi, x, y;
                     11: double log(), exp(), cos(), sin(), atan2(), cabs();
                     12: 
                     13: logr = log( cabs(a->dreal, a->dimag) );
                     14: logi = atan2(a->dimag, a->dreal);
                     15: 
                     16: x = exp( logr * b->dreal - logi * b->dimag );
                     17: y = logr * b->dimag + logi * b->dreal;
                     18: 
                     19: r->dreal = x * cos(y);
                     20: r->dimag = x * sin(y);
                     21: }

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