Mercurial > hg > octave-nkf

comparison liboctave/dbleHESS.cc @ 11570:57632dea2446

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
attempt better backward compatibility for Array constructors
author John W. Eaton <jwe@octave.org>
date Wed, 19 Jan 2011 17:55:56 -0500
parents fd0a3ac60b0e
children 72c96de7a403
comparison
equal deleted inserted replaced
11569:7e9a111cae20 11570:57632dea2446
82 octave_idx_type ihi; 82 octave_idx_type ihi;
83 83
84 hess_mat = a; 84 hess_mat = a;
85 double *h = hess_mat.fortran_vec (); 85 double *h = hess_mat.fortran_vec ();
86 86
87 Array<double> scale (n, 1); 87 Array<double> scale (dim_vector (n, 1));
88 double *pscale = scale.fortran_vec (); 88 double *pscale = scale.fortran_vec ();
89 89
90 F77_XFCN (dgebal, DGEBAL, (F77_CONST_CHAR_ARG2 (&job, 1), 90 F77_XFCN (dgebal, DGEBAL, (F77_CONST_CHAR_ARG2 (&job, 1),
91 n, h, n, ilo, ihi, pscale, info 91 n, h, n, ilo, ihi, pscale, info
92 F77_CHAR_ARG_LEN (1))); 92 F77_CHAR_ARG_LEN (1)));
93 93
94 Array<double> tau (n-1, 1); 94 Array<double> tau (dim_vector (n-1, 1));
95 double *ptau = tau.fortran_vec (); 95 double *ptau = tau.fortran_vec ();
96 96
97 Array<double> work (lwork, 1); 97 Array<double> work (dim_vector (lwork, 1));
98 double *pwork = work.fortran_vec (); 98 double *pwork = work.fortran_vec ();
99 99
100 F77_XFCN (dgehrd, DGEHRD, (n, ilo, ihi, h, n, ptau, pwork, 100 F77_XFCN (dgehrd, DGEHRD, (n, ilo, ihi, h, n, ptau, pwork,
101 lwork, info)); 101 lwork, info));
102 102