Mercurial > hg > octave-nkf

diff liboctave/floatHESS.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/floatHESS.cc
+++ b/liboctave/floatHESS.cc
@@ -80,17 +80,17 @@
   hess_mat = a;
   float *h = hess_mat.fortran_vec ();
 
-  Array<float> scale (n);
+  Array<float> scale (n, 1);
   float *pscale = scale.fortran_vec ();
 
   F77_XFCN (sgebal, SGEBAL, (F77_CONST_CHAR_ARG2 (&job, 1),
                              n, h, n, ilo, ihi, pscale, info
                              F77_CHAR_ARG_LEN (1)));
 
-  Array<float> tau (n-1);
+  Array<float> tau (n-1, 1);
   float *ptau = tau.fortran_vec ();
 
-  Array<float> work (lwork);
+  Array<float> work (lwork, 1);
   float *pwork = work.fortran_vec ();
 
   F77_XFCN (sgehrd, SGEHRD, (n, ilo, ihi, h, n, ptau, pwork,