Mercurial > hg > octave-jordi

diff src/DLD-FUNCTIONS/qz.cc @ 4629:ed92d574b55b

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2003-11-17 20:43:13 by jwe]
author jwe
date Mon, 17 Nov 2003 20:43:13 +0000 (2003-11-17)
parents 30ba814d6700
children bdbee5282954
line wrap: on
line diff
--- a/src/DLD-FUNCTIONS/qz.cc
+++ b/src/DLD-FUNCTIONS/qz.cc
@@ -107,7 +107,7 @@
 			     F77_CHAR_ARG_LEN_DECL);
 
   F77_RET_T
-  F77_FUNC (dlag2, DLAG2) (double* A, const int& LDA, double* B,
+  F77_FUNC (dlag2, DLAG2) (const double* A, const int& LDA, const double* B,
 			   const int& LDB, const double& SAFMIN,
 			   double& SCALE1, double& SCALE2,
 			   double& WR1, double& WR2, double& WI);
@@ -760,8 +760,10 @@
 		  // fortran_vec instead of &aa(jj,jj) here.
 
 		  double scale1, scale2, wr1, wr2, wi;
+		  const double *aa_ptr = aa.data () + jj*nn+jj;
+		  const double *bb_ptr = bb.data () + jj*nn+jj;
 		  F77_XFCN (dlag2, DLAG2,
-			    (&aa(jj,jj), nn, &bb(jj,jj), nn, safmin,
+			    (aa_ptr, nn, bb_ptr, nn, safmin,
 			     scale1, scale2, wr1, wr2, wi));
 
 #ifdef DEBUG_EIG