Return to jacstep CVS log

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

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

1.1     ! root        1:         d = (a(q,q)-a(p,p))*0.5/a(p,q);
        !             2:         t = 1/(abs(d)+sqrt(d*d+1));
        !             3:         if d < 0, t = -t; end;
        !             4:         c = 1/sqrt(1+t*t);  s = t*c;
        !             5:         R = eye(n); r(p,p)=c; r(q,q)=c; r(p,q)=s; r(q,p)=-s;
        !             6:         X = X*R;
        !             7:         A = R'*A*R;

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