Return to cdiv CVS log

Up to [Research Unix] / researchv10dc / cmd / matlab / exec

Annotation of researchv10dc/cmd/matlab/exec/cdiv, revision 1.1

1.1     ! root        1:    ar = real(a); ai = imag(a); br = real(b); bi = imag(b);
        !             2:    p = bi/br;
        !             3:    t = (ai - p*ar)/(br + p*bi); 
        !             4:    cr = p*t + ar/br; 
        !             5:    ci = t; 
        !             6:    p2 = br/bi; 
        !             7:    t2 = (ai + p2*ar)/(bi + p2*br);
        !             8:    ci2 = p2*t2 - ar/bi;
        !             9:    cr2 = t2;
        !            10:    s = abs(br) + abs(bi);
        !            11:    ars = ar/s;
        !            12:    ais = ai/s;
        !            13:    brs = br/s;
        !            14:    bis = bi/s;
        !            15:    s = brs**2 + bis**2;
        !            16:    cr3 = (ars*brs + ais*bis)/s;
        !            17:    ci3 = (ais*brs - ars*bis)/s;
        !            18:    <cr ci; cr2 ci2; cr3 ci3>

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