Return to oper.c CVS log

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

Annotation of researchv10no/cmd/cfront/libC/complex/oper.c, revision 1.1

1.1     ! root        1: 
        !             2: #include "complex.h"
        !             3: #include "const.h"
        !             4: 
        !             5: complex operator*(complex a1,complex a2) 
        !             6: {
        !             7:        return complex(a1.re*a2.re-a1.im*a2.im, a1.re*a2.im+a1.im*a2.re);
        !             8: }
        !             9: 
        !            10: 
        !            11: complex operator/(complex a1, complex a2)
        !            12: {
        !            13:        register double r = a2.re;      /* (r,i) */
        !            14:        register double i = a2.im;
        !            15:        register double ti;             /* (tr,ti) */
        !            16:        register double tr;
        !            17: 
        !            18:        tr = ABS(r);
        !            19:        ti = ABS(i);
        !            20: 
        !            21:        if (tr <= ti) {
        !            22:                ti = r/i;
        !            23:                tr = i * (1 + ti*ti);
        !            24:                r = a1.re;
        !            25:                i = a1.im;
        !            26:        }
        !            27:        else {
        !            28:                ti = -i/r;
        !            29:                tr = r * (1 + ti*ti);
        !            30:                r = -a1.im;
        !            31:                i = a1.re;
        !            32:        }
        !            33: 
        !            34:        return complex( (r*ti + i)/tr, (i*ti - r)/tr );
        !            35: }
        !            36: 
        !            37: void complex::operator*=(complex a)
        !            38: {
        !            39:        register double r = re*a.re - im*a.im;
        !            40:        register double i = re*a.im + im*a.re;
        !            41:        re = r;
        !            42:        im = i;
        !            43: }
        !            44: 
        !            45: void complex::operator/=(complex a)
        !            46: {
        !            47:        complex quot, temp1, temp2;
        !            48: 
        !            49:        if ( (temp2.re = a.re) < 0 ) temp2.re = -temp2.re;
        !            50:        if ( (temp2.im = a.im) < 0 ) temp2.im = -temp2.im;
        !            51:        if ( temp2.re <= temp2.im) {
        !            52:                temp2.im = a.re/a.im;
        !            53:                temp2.re = a.im * (1 + temp2.im*temp2.im);
        !            54:                temp1 = *this;
        !            55:        }
        !            56:        else {
        !            57:                temp2.im = -a.im/a.re;
        !            58:                temp2.re = a.re * (1 + temp2.im*temp2.im);
        !            59:                temp1.re = -im;
        !            60:                temp1.im = re;
        !            61:        }
        !            62:        re = (temp1.re * temp2.im + temp1.im) / temp2.re;
        !            63:        im = (temp1.im * temp2.im - temp1.re) / temp2.re;
        !            64: }

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