Mercurial > hg > octave-nkf

diff liboctave/floatCHOL.cc @ 10350:12884915a8e4

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
merge MArray classes & improve Array interface
author Jaroslav Hajek <highegg@gmail.com>
date Sat, 23 Jan 2010 21:41:03 +0100
parents 07ebe522dac2
children 8a5e980da6aa
line wrap: on
line diff
--- a/liboctave/floatCHOL.cc
+++ b/liboctave/floatCHOL.cc
@@ -123,9 +123,9 @@
       octave_idx_type spocon_info = 0;
 
       // Now calculate the condition number for non-singular matrix.
-      Array<float> z (3*n);
+      Array<float> z (3*n, 1);
       float *pz = z.fortran_vec ();
-      Array<octave_idx_type> iz (n);
+      Array<octave_idx_type> iz (n, 1);
       octave_idx_type *piz = iz.fortran_vec ();
       F77_XFCN (spocon, SPOCON, (F77_CONST_CHAR_ARG2 ("U", 1), n, h,
                                  n, anorm, xrcond, pz, piz, spocon_info
@@ -416,7 +416,7 @@
   else
     {
       FloatMatrix a = chol_mat.transpose () * chol_mat;
-      Array<octave_idx_type> p (n);
+      Array<octave_idx_type> p (n, 1);
       for (octave_idx_type k = 0; k < n; k++) p(k) = k;
       if (i < j)
         {