--- a/EigenLib3.mod
+++ b/EigenLib3.mod
@@ -781,24 +781,24 @@
(* GOTO enter1; *)
ELSE
IF (ABS(s) > MachEps) THEN Ca := c / Root1; Sa := s / Root1; END;
- Cot2X:= d / Root1;
- CotX := Cot2X + (sig*sqrt(1.0 + Cot2X*Cot2X));
- Sx := sig / sqrt(1.0 + CotX*CotX);
- Cx := Sx*CotX;
- eta := (Er*Br + Bi*Ei) / Root1;
- tse := (Br*Bi - Er*Ei) / Root1;
- Te := sig*( - Root1*De + tse*d) / Root2;
- Tee := (d*De + Root1*tse) / Root2;
- nD := Root2*Root2 + Tee*Tee;
- Tee := Hj*Cx*Sx;
- Cos2a:= Ca*Ca - Sa*Sa;
- Sin2a:= 2.0*Ca*Sa;
- tem := Hr*Cos2a + Hi*Sin2a;
- tep := Hi*Cos2a - Hr*Sin2a;
- Hr := Cx*Cx*Hr - Sx*Sx*tem - Ca*Tee;
- Hi := Cx*Cx*Hi + Sx*Sx*tep - Sa*Tee;
- b := isw*Te*Ca + eta*Sa;
- e := eta*Ca - isw*Te*Sa;
+ Cot2X := d / Root1;
+ CotX := Cot2X + (sig*sqrt(1.0 + Cot2X*Cot2X));
+ Sx := sig / sqrt(1.0 + CotX*CotX);
+ Cx := Sx*CotX;
+ eta := (Er*Br + Bi*Ei) / Root1;
+ tse := (Br*Bi - Er*Ei) / Root1;
+ Te := sig*( - Root1*De + tse*d) / Root2;
+ Tee := (d*De + Root1*tse) / Root2;
+ nD := Root2*Root2 + Tee*Tee;
+ Tee := Hj*Cx*Sx;
+ Cos2a := Ca*Ca - Sa*Sa;
+ Sin2a := 2.0*Ca*Sa;
+ tem := Hr*Cos2a + Hi*Sin2a;
+ tep := Hi*Cos2a - Hr*Sin2a;
+ Hr := Cx*Cx*Hr - Sx*Sx*tem - Ca*Tee;
+ Hi := Cx*Cx*Hi + Sx*Sx*tep - Sa*Tee;
+ b := isw*Te*Ca + eta*Sa;
+ e := eta*Ca - isw*Te*Sa;
END; (* IF (Root1 < MachEps) *)
(*enter1:*) s := Hr - sig*Root2*e;
c := Hi - sig*Root2*b;
@@ -807,26 +807,26 @@
Cb:=1.0; Ch:=1.0; Sb:=0.0; Sh:=0.0;
(* GOTO trans; *)
ELSE
- Cb := - c / Root;
- Sb := s / Root;
- Tee := Cb*b - e*Sb;
- nC := Tee*Tee;
- TanH:= Root / (g + 2.0*(nC + nD));
- Ch := 1.0 / sqrt(1.0 - TanH*TanH);
- Sh := Ch*TanH;
+ Cb := - c / Root;
+ Sb := s / Root;
+ Tee := Cb*b - e*Sb;
+ nC := Tee*Tee;
+ TanH := Root / (g + 2.0*(nC + nD));
+ Ch := 1.0 / sqrt(1.0 - TanH*TanH);
+ Sh := Ch*TanH;
END; (* IF *)
(*trans:*) tem := Sx*Sh*(Sa*Cb - Sb*Ca);
- im := - Sx*Sh*(Ca*Cb + Sa*Sb);
- c1 := CMPLX(Cx*Ch - tem,im);
- c2 := CMPLX(Cx*Ch + tem,im);
+ im := - Sx*Sh*(Ca*Cb + Sa*Sb);
+ c1 := CMPLX(Cx*Ch - tem,im);
+ c2 := CMPLX(Cx*Ch + tem,im);
tep := Sx*Ch*Ca;
tem := Cx*Sh*Sb;
- s1re := tep - tem;
- s2re := - tep - tem;
+ s1re := tep - tem;
+ s2re := - tep - tem;
tep := Sx*Ch*Sa;
tem := Cx*Sh*Cb;
- s1im := tep + tem;
- s2im := tep - tem;
+ s1im := tep + tem;
+ s2im := tep - tem;
s1 := CMPLX(s1re,s1im);
s2 := CMPLX(s2re,s2im);
tem := LongComplexMath.abs(s1); (* M.R.*)