Mercurial > hg > octave-lyh

diff liboctave/oct-norm.cc @ 14520:c4ed0fdf2b62

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix initialization in Higham's method for norm estimate (bug #36031). * oct-norm.cc (higham): Initialize mu to 1.
author Marco Caliari <marco.caliari@univr.it>
date Tue, 03 Apr 2012 20:38:05 -0700
parents 72c96de7a403
children
line wrap: on
line diff
--- a/liboctave/oct-norm.cc
+++ b/liboctave/oct-norm.cc
@@ -415,7 +415,7 @@
   // the OSE part
   VectorT y(m.rows (), 1, 0), z(m.rows (), 1);
   typedef typename VectorT::element_type RR;
-  RR lambda = 0, mu = 0;
+  RR lambda = 0, mu = 1;
   for (octave_idx_type k = 0; k < m.columns (); k++)
     {
       octave_quit ();