Mercurial > hg > octave-nkf
diff liboctave/dMatrix.cc @ 532:2ca256b77602
[project @ 1994-07-20 19:56:55 by jwe]
| author | jwe |
|---|---|
| date | Wed, 20 Jul 1994 19:56:55 +0000 |
| parents | 8e9dcc406332 |
| children | 883197c5ad75 |
line wrap: on
line diff
--- a/liboctave/dMatrix.cc +++ b/liboctave/dMatrix.cc @@ -566,7 +566,7 @@ } Matrix -Matrix::inverse (int& info, volatile double& rcond) const +Matrix::inverse (int& info, double& rcond) const { int nr = rows (); int nc = cols (); @@ -585,7 +585,8 @@ F77_FCN (dgeco) (tmp_data, &nr, &nc, ipvt, &rcond, z); - if (rcond + 1.0 == 1.0) + volatile double tmp_rcond = rcond; + if (tmp_rcond + 1.0 == 1.0) { info = -1; copy (tmp_data, data (), len); // Restore matrix contents. @@ -685,7 +686,7 @@ } DET -Matrix::determinant (int& info, volatile double& rcond) const +Matrix::determinant (int& info, double& rcond) const { DET retval; @@ -709,7 +710,8 @@ F77_FCN (dgeco) (tmp_data, &nr, &nr, ipvt, &rcond, z); - if (rcond + 1.0 == 1.0) + volatile double tmp_rcond = rcond; + if (tmp_rcond + 1.0 == 1.0) { info = -1; retval = DET (); @@ -746,7 +748,7 @@ } Matrix -Matrix::solve (const Matrix& b, int& info, volatile double& rcond) const +Matrix::solve (const Matrix& b, int& info, double& rcond) const { Matrix retval; @@ -767,7 +769,8 @@ F77_FCN (dgeco) (tmp_data, &nr, &nr, ipvt, &rcond, z); - if (rcond + 1.0 == 1.0) + volatile double tmp_rcond = rcond; + if (tmp_rcond + 1.0 == 1.0) { info = -2; } @@ -827,7 +830,7 @@ } ColumnVector -Matrix::solve (const ColumnVector& b, int& info, volatile double& rcond) const +Matrix::solve (const ColumnVector& b, int& info, double& rcond) const { ColumnVector retval; @@ -848,7 +851,8 @@ F77_FCN (dgeco) (tmp_data, &nr, &nr, ipvt, &rcond, z); - if (rcond + 1.0 == 1.0) + volatile double tmp_rcond = rcond; + if (tmp_rcond + 1.0 == 1.0) { info = -2; }