Mercurial > octave-nkf

changeset 10823:3d89d262f5d4

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
slight simplification in rsf2csf code
author Jaroslav Hajek <highegg@gmail.com>
date Tue, 27 Jul 2010 12:31:58 +0200
parents 23d2378512a0
children 1e6664326d32
files libcruft/lapack-xtra/crsf2csf.f libcruft/lapack-xtra/zrsf2csf.f
diffstat 2 files changed, 4 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- 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)
--- 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)