Mercurial > hg > octave-nkf
changeset 4330:9f86c2055b58
[project @ 2003-02-18 22:31:55 by jwe]
author | jwe |
---|---|
date | Tue, 18 Feb 2003 22:31:55 +0000 |
parents | d53c33d93440 |
children | d3278845f764 |
files | liboctave/CMatrix.cc liboctave/dMatrix.cc |
diffstat | 2 files changed, 28 insertions(+), 28 deletions(-) [+] |
line wrap: on
line diff
--- a/liboctave/CMatrix.cc +++ b/liboctave/CMatrix.cc @@ -960,9 +960,9 @@ Complex *tmp_data = retval.fortran_vec (); Array<Complex> z(1); - int lwork = 1; - - // Query the optimum work array size + int lwork = -1; + + // Query the optimum work array size. F77_XFCN (zgetri, ZGETRI, (nc, tmp_data, nr, pipvt, z.fortran_vec (), lwork, info)); @@ -981,7 +981,7 @@ info = 0; - /* Calculate the norm of the matrix, for later use */ + // Calculate the norm of the matrix, for later use. double anorm; if (calc_cond) anorm = retval.abs().sum().row(0).max(); @@ -992,13 +992,13 @@ (*current_liboctave_error_handler) ("unrecoverable error in zgetrf"); else { - /* Throw-away extra info LAPACK gives so as to not change output */ + // Throw-away extra info LAPACK gives so as to not change output. rcond = 0.; if ( info != 0) info = -1; else if (calc_cond) { - /* Now calculate the condition number for non-singular matrix */ + // Now calculate the condition number for non-singular matrix. char job = '1'; Array<double> rz (2 * nc); double *prz = rz.fortran_vec (); @@ -1436,7 +1436,7 @@ info = 0; - /* Calculate the norm of the matrix, for later use */ + // Calculate the norm of the matrix, for later use. double anorm = 0; if (calc_cond) anorm = atmp.abs().sum().row(0).max(); @@ -1447,7 +1447,7 @@ (*current_liboctave_error_handler) ("unrecoverable error in zgetrf"); else { - /* Throw-away extra info LAPACK gives so as to not change output */ + // Throw-away extra info LAPACK gives so as to not change output. rcond = 0.; if ( info != 0) { @@ -1458,7 +1458,7 @@ { if (calc_cond) { - /* Now calc the condition number for non-singular matrix */ + // Now calc the condition number for non-singular matrix. char job = '1'; Array<Complex> z (2*nr); Complex *pz = z.fortran_vec (); @@ -1583,7 +1583,7 @@ Array<double> rz (2 * nc); double *prz = rz.fortran_vec (); - /* Calculate the norm of the matrix, for later use */ + // Calculate the norm of the matrix, for later use. double anorm = atmp.abs().sum().row(0).max(); F77_XFCN (zgetrf, ZGETRF, (nr, nr, tmp_data, nr, pipvt, info)); @@ -1592,7 +1592,7 @@ (*current_liboctave_error_handler) ("unrecoverable error in zgetrf"); else { - /* Throw-away extra info LAPACK gives so as to not change output */ + // Throw-away extra info LAPACK gives so as to not change output. rcond = 0.; if ( info != 0) { @@ -1607,7 +1607,7 @@ } else { - /* Now calculate the condition number for non-singular matrix */ + // Now calculate the condition number for non-singular matrix. char job = '1'; F77_XFCN (zgecon, ZGECON, ( &job, nc, tmp_data, nr, anorm, rcond, pz, prz, info)); @@ -1732,7 +1732,7 @@ Array<double> rz (2 * nc); double *prz = rz.fortran_vec (); - /* Calculate the norm of the matrix, for later use */ + // Calculate the norm of the matrix, for later use. double anorm = atmp.abs().sum().row(0).max(); F77_XFCN (zgetrf, ZGETRF, (nr, nr, tmp_data, nr, pipvt, info)); @@ -1741,7 +1741,7 @@ (*current_liboctave_error_handler) ("unrecoverable error in zgetrf"); else { - /* Throw-away extra info LAPACK gives so as to not change output */ + // Throw-away extra info LAPACK gives so as to not change output. rcond = 0.; if ( info != 0) { @@ -1756,7 +1756,7 @@ } else { - /* Now calculate the condition number for non-singular matrix */ + // Now calculate the condition number for non-singular matrix. char job = '1'; F77_XFCN (zgecon, ZGECON, ( &job, nc, tmp_data, nr, anorm, rcond, pz, prz, info));
--- a/liboctave/dMatrix.cc +++ b/liboctave/dMatrix.cc @@ -630,7 +630,7 @@ Array<double> z(1); int lwork = -1; - // Query the optimum work array size + // Query the optimum work array size. F77_XFCN (dgetri, DGETRI, (nc, tmp_data, nr, pipvt, z.fortran_vec (), lwork, info)); @@ -648,7 +648,7 @@ info = 0; - /* Calculate the norm of the matrix, for later use */ + // Calculate the norm of the matrix, for later use. double anorm = 0; if (calc_cond) anorm = retval.abs().sum().row(0).max(); @@ -659,13 +659,13 @@ (*current_liboctave_error_handler) ("unrecoverable error in dgetrf"); else { - /* Throw-away extra info LAPACK gives so as to not change output */ + // Throw-away extra info LAPACK gives so as to not change output. rcond = 0.; if ( info != 0) info = -1; else if (calc_cond) { - /* Now calculate the condition number for non-singular matrix */ + // Now calculate the condition number for non-singular matrix. char job = '1'; Array<int> iz (nc); int *piz = iz.fortran_vec (); @@ -1101,7 +1101,7 @@ info = 0; - /* Calculate the norm of the matrix, for later use */ + // Calculate the norm of the matrix, for later use. double anorm = 0; if (calc_cond) anorm = atmp.abs().sum().row(0).max(); @@ -1112,7 +1112,7 @@ (*current_liboctave_error_handler) ("unrecoverable error in dgetrf"); else { - /* Throw-away extra info LAPACK gives so as to not change output */ + // Throw-away extra info LAPACK gives so as to not change output. rcond = 0.; if ( info != 0) { @@ -1123,7 +1123,7 @@ { if (calc_cond) { - /* Now calc the condition number for non-singular matrix */ + // Now calc the condition number for non-singular matrix. char job = '1'; Array<double> z (4 * nc); double *pz = z.fortran_vec (); @@ -1219,7 +1219,7 @@ Array<int> iz (nc); int *piz = iz.fortran_vec (); - /* Calculate the norm of the matrix, for later use */ + // Calculate the norm of the matrix, for later use. double anorm = atmp.abs().sum().row(0).max(); F77_XFCN (dgetrf, DGETRF, (nr, nr, tmp_data, nr, pipvt, info)); @@ -1228,7 +1228,7 @@ (*current_liboctave_error_handler) ("unrecoverable error in dgetrf"); else { - /* Throw-away extra info LAPACK gives so as to not change output */ + // Throw-away extra info LAPACK gives so as to not change output. rcond = 0.; if ( info != 0) { @@ -1243,7 +1243,7 @@ } else { - /* Now calculate the condition number for non-singular matrix */ + // Now calculate the condition number for non-singular matrix. char job = '1'; F77_XFCN (dgecon, DGECON, ( &job, nc, tmp_data, nr, anorm, rcond, pz, piz, info)); @@ -1366,7 +1366,7 @@ Array<int> iz (nc); int *piz = iz.fortran_vec (); - /* Calculate the norm of the matrix, for later use */ + // Calculate the norm of the matrix, for later use. double anorm = atmp.abs().sum().row(0).max(); F77_XFCN (dgetrf, DGETRF, (nr, nr, tmp_data, nr, pipvt, info)); @@ -1375,7 +1375,7 @@ (*current_liboctave_error_handler) ("unrecoverable error in dgetrf"); else { - /* Throw-away extra info LAPACK gives so as to not change output */ + // Throw-away extra info LAPACK gives so as to not change output. rcond = 0.; if ( info > 0) { @@ -1390,7 +1390,7 @@ } else { - /* Now calculate the condition number for non-singular matrix */ + // Now calculate the condition number for non-singular matrix. char job = '1'; F77_XFCN (dgecon, DGECON, ( &job, nc, tmp_data, nr, anorm, rcond, pz, piz, info));