Changeset 774 in Sophya for trunk/SophyaLib/NTools/linfit.cc
- Timestamp:
- Mar 10, 2000, 4:24:53 PM (26 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaLib/NTools/linfit.cc
r540 r774 1 1 #include "machdefs.h" 2 2 #include "peida.h" 3 #include "sopemtx.h" 3 4 #include "linfit.h" 4 5 … … 38 39 for (int j=0; j<nf; j++) 39 40 for (int k=j; k<nf; k++) 40 a(j,k) = a(k,j) = fx.Row(j) * fx.Row(k); 41 a(j,k) = a(k,j) = TMatrixRC<r_8>::Row(const_cast<Matrix &>(fx), j) 42 * TMatrixRC<r_8>::Row(const_cast<Matrix &>(fx), k); /* $CHECK$ Reza 10/3/2000 */ 41 43 42 44 c = fx * y; 43 45 44 if ( Matrix::GausPiv(a,c) == 0) THROW(singMatxErr);46 if (SimpleMatrixOperation<r_8>::GausPiv(a,c) == 0) THROW(singMatxErr); /* $CHECK$ Reza 10/3/2000 */ 45 47 46 48 double xi2 = 0; … … 105 107 } 106 108 107 if ( Matrix::GausPiv(a,d) == 0) THROW(singMatxErr);109 if (SimpleMatrixOperation<r_8>::GausPiv(a,d) == 0) THROW(singMatxErr); /* $CHECK$ Reza 10/3/2000 */ 108 110 109 111 for (int l=0; l<nf; l++) {
Note:
See TracChangeset
for help on using the changeset viewer.