# HG changeset patch
# User Jaroslav Hajek <highegg@gmail.com>
# Date 1280226718 -7200
# Node ID 3d89d262f5d4734fdad3ab8049d8a72f49b06ea8
# Parent  23d2378512a02d8ce17558645584f7aa0204397c
slight simplification in rsf2csf code

diff -r 23d2378512a0 -r 3d89d262f5d4 libcruft/lapack-xtra/crsf2csf.f
--- a/libcruft/lapack-xtra/crsf2csf.f	Tue Jul 27 11:26:43 2010 +0200
+++ b/libcruft/lapack-xtra/crsf2csf.f	Tue Jul 27 12:31:58 2010 +0200
@@ -35,8 +35,8 @@
 c 2x2 block, form Givens rotation [c, i*s; i*s, c] 
            x = t(j,j)
            z = t(j,j+1)
-           c(j) = sqrt(abs(z/(y-z)))
-           s(j) = sign(sqrt(abs(y/(y-z))),z)
+           c(j) = sqrt(z/(z-y))
+           s(j) = sign(sqrt(-y/(z-y)),z)
 c apply new rotation to t(j:j+1,j)
            call crcrot1(2,t(j,j),c(j),s(j))
 c apply all rotations to t(1:j+1,j+1)
diff -r 23d2378512a0 -r 3d89d262f5d4 libcruft/lapack-xtra/zrsf2csf.f
--- a/libcruft/lapack-xtra/zrsf2csf.f	Tue Jul 27 11:26:43 2010 +0200
+++ b/libcruft/lapack-xtra/zrsf2csf.f	Tue Jul 27 12:31:58 2010 +0200
@@ -35,8 +35,8 @@
 c 2x2 block, form Givens rotation [c, i*s; i*s, c] 
            x = t(j,j)
            z = t(j,j+1)
-           c(j) = sqrt(abs(z/(y-z)))
-           s(j) = sign(sqrt(abs(y/(y-z))),z)
+           c(j) = sqrt(z/(z-y))
+           s(j) = sign(sqrt(-y/(z-y)),z)
 c apply new rotation to t(j:j+1,j)
            call zrcrot1(2,t(j,j),c(j),s(j))
 c apply all rotations to t(1:j+1,j+1)