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

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.