Changeset 774 in Sophya for trunk/SophyaLib/NTools/linfit.cc


Ignore:
Timestamp:
Mar 10, 2000, 4:24:53 PM (26 years ago)
Author:
ansari
Message:

Adaptation aux modifs de TMatrixRC<T> - Reza 10/3/2000

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SophyaLib/NTools/linfit.cc

    r540 r774  
    11#include "machdefs.h"
    22#include "peida.h"
     3#include "sopemtx.h"
    34#include "linfit.h"
    45
     
    3839  for (int j=0; j<nf; j++)
    3940    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 */
    4143   
    4244  c = fx * y;
    4345
    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 */
    4547
    4648  double xi2 = 0;
     
    105107  }
    106108
    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 */
    108110
    109111  for (int l=0; l<nf; l++) {
Note: See TracChangeset for help on using the changeset viewer.