Changeset 2421 in Sophya for trunk/SophyaLib/TArray/sopemtx.cc

Timestamp:

Aug 8, 2003, 3:07:06 PM (22 years ago)

Author:

ansari

Message:

• Modifications methodes TMatrix::NRows() ::NCols() --> Peuvent etre

appelees sur des matrices non allouees.

• Declaration const pour TMatrix::Transpose() et ::Rearrange(..)
• remplacement THROW par throw ds Linfitter avec message explicite

Reza 08/08/2003

File:

• trunk/SophyaLib/TArray/sopemtx.cc (modified) (6 diffs)

trunk/SophyaLib/TArray/sopemtx.cc

@@ -619 +619 @@
 {
 uint_4 n = x.NElts();
-if (n != y.NElts()) THROW(sizeMismatchErr);
+if (n != y.NElts()) 
+  throw SzMismatchError("LinFitter<T>::LinFit(x,y,nf,f,c)/Error x.NElts() <> y.Nelts() ");
  
 TMatrix<T> fx(nf,n);
@@ -643 +644 @@
 {
 uint_4 n = y.NElts();
-if (n != fx.NCol()) THROW(sizeMismatchErr);
+if (n != fx.NCol()) 
+  throw SzMismatchError("LinFitter<T>::LinFit(fx, y, c)/Error y.NElts() <> fx.Nelts() ");
 
 uint_4 nf = fx.NRows();
@@ -656 +658 @@
 c = fx * y;
 
-if(SimpleMatrixOperation<T>::GausPiv(a,c)==(T) 0) THROW(singMatxErr);
+if(SimpleMatrixOperation<T>::GausPiv(a,c)==(T) 0) 
+  throw SingMatrixExc("LinFitter<T>::LinFit(fx, y, c) - Non invertible matrix (by GausPiv())");
 
 r_8 xi2 = 0., ax;
@@ -691 +694 @@
 {
 uint_4 n = x.NElts();
-if (n != y.NElts()) THROW(sizeMismatchErr);
+if (n != y.NElts()) 
+  throw SzMismatchError("LinFitter<T>::LinFit(x,y,errY2,nf,f,c,errC)/Error x.NElts() <> y.Nelts() ");
 
 TMatrix<T> fx(nf, n);
@@ -720 +724 @@
 {
 uint_4 n = y.NElts();
-if( n != errY2.NElts()) THROW(sizeMismatchErr);
-if( n != fx.NCol()) THROW(sizeMismatchErr);
-
+if( n != errY2.NElts()) 
+  throw SzMismatchError("LinFitter<T>::LinFit(fx,y,errY2,c,errC)/Error y.NElts() <> errY2.Nelts() ");
+if( n != fx.NCol()) 
+  throw SzMismatchError("LinFitter<T>::LinFit(fx,y,errY2,c,errC)/Error y.NElts() <> fx.NCols() ");
 uint_4 nf = fx.NRows();
 
@@ -748 +753 @@
 }
 
-if(SimpleMatrixOperation<T>::GausPiv(a,d)==(T) 0) THROW(singMatxErr);
+if(SimpleMatrixOperation<T>::GausPiv(a,d)==(T) 0) 
+  throw SingMatrixExc("LinFitter<T>::LinFit(...ErrY2...) - Non invertible matrix (by GausPiv())");
+
 
 for(uint_4 l=0; l<nf; l++) {