Mercurial > hg > octave-lyh

diff scripts/linear-algebra/housh.m @ 3426:f8dde1807dee

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2000-01-13 08:40:00 by jwe]
author jwe
date Thu, 13 Jan 2000 08:40:53 +0000
parents f7e4a95916f2
children e098ebb77023
line wrap: on
line diff
--- a/scripts/linear-algebra/housh.m
+++ b/scripts/linear-algebra/housh.m
@@ -28,7 +28,7 @@
   #   beta: If beta = 0, then no reflection need be applied (zer set to 0)
   #   housv: householder vector
   # mar 6,1987 : rev dec 17,1988
-  #		rev sep 19,1991 (blas)
+  #             rev sep 19,1991 (blas)
   # translated from FORTRAN Aug 1995
   # A. S. Hodel
 
@@ -43,10 +43,10 @@
   else
     housv = x;
     m = max(abs(housv));
-    if (m ~= 0.0) 
+    if (m ~= 0.0)
       housv = housv/m;
       alpha = norm(housv);
-      if (alpha > z) 
+      if (alpha > z)
         beta = 1.0/(alpha*(alpha+abs(housv(j))));
         sg = sign(housv(j));
         if( sg == 0)