1 SUBROUTINE dlaev2( A, B, C, RT1, RT2, CS1, SN1 )
9 DOUBLE PRECISION A, B, C, CS1, RT1, RT2, SN1
68 parameter( one = 1.0d0 )
70 parameter( two = 2.0d0 )
72 parameter( zero = 0.0d0 )
74 parameter( half = 0.5d0 )
78 DOUBLE PRECISION AB, ACMN, ACMX, ACS, ADF, CS, CT, DF, RT, SM,
93 IF( abs( a ).GT.abs( c ) )
THEN
101 rt = adf*sqrt( one+( ab / adf )**2 )
102 ELSE IF( adf.LT.ab )
THEN
103 rt = ab*sqrt( one+( adf / ab )**2 )
110 IF( sm.LT.zero )
THEN
118 rt2 = ( acmx / rt1 )*acmn - ( b / rt1 )*b
119 ELSE IF( sm.GT.zero )
THEN
127 rt2 = ( acmx / rt1 )*acmn - ( b / rt1 )*b
139 IF( df.GE.zero )
THEN
149 sn1 = one / sqrt( one+ct*ct )
152 IF( ab.EQ.zero )
THEN
157 cs1 = one / sqrt( one+tn*tn )
161 IF( sgn1.EQ.sgn2 )
THEN
subroutine dlaev2(A, B, C, RT1, RT2, CS1, SN1)