Changeset 1207 in Sophya for trunk/SophyaPI/PIext/nomtmatvecadapter.cc

Timestamp:

Sep 29, 2000, 4:26:15 PM (25 years ago)

Author:

ercodmgr

Message:

• refonte de piaffiting avec les fits adapteurs
• addapteur de fit pour les objets existants
• entree des Image<T> dans les fits (retour !!!)

cmv 29/9/00

File:

• trunk/SophyaPI/PIext/nomtmatvecadapter.cc (modified) (4 diffs)

trunk/SophyaPI/PIext/nomtmatvecadapter.cc

@@ -7 +7 @@
 
 #include "tvector.h"
+#include "objfitter.h"
 #include "nomtmatvecadapter.h"
 #include "piscdrawwdg.h"
@@ -119 +120 @@
 }
 
+
+/* --Methode-- */
+template <class T> 
+GeneralFitData* NOMAdapter_TMatrix<T>::GetGeneralFitData(bool& adel
+ ,GeneralFitData::FitErrType errtype,double errscale,double errmin
+ ,int i1,int i2,int j1,int j2)
+{
+adel = false;
+if(!mMtx) return(NULL);
+
+int nx,ny;
+TVector<T>* vec = dynamic_cast<TVector<T> *>(mMtx);
+
+if(vec) {  // C'est un TVector
+  nx = vec->NElts();
+  ny = 1;
+} else {
+  nx = mMtx->NRows();
+  ny = mMtx->NCol();
+  if(nx<=0 || ny<=0) return(NULL);
+}
+
+i1 = (i1<0||i1>=nx)? 0: i1;
+i2 = (i2<0||i2>=nx||i2<i1)? nx-1: i2;
+j1 = (j1<0||j1>=ny)? 0: j1;
+j2 = (j2<0||j2>=ny||j2<j1)? ny-1: j2;
+
+GeneralFitData* mGData = NULL;
+
+if(vec) {  // C'est un TVector
+  mGData = new GeneralFitData(1,(i2-i1+1),0);
+  adel = true;
+  for(int i=i1;i<=i2;i++) {
+    double x = (double) i;
+    double f = (*vec)(i);
+    double e = 1.;
+    e = GeneralFitData::ComputeError(f,e,errtype,errscale,errmin);
+    mGData->AddData1(x,f,e);
+  }
+} else {
+  mGData = new GeneralFitData(2,(i2-i1+1)*(j2-j1+1),0);
+  adel = true;
+  for(int i=i1;i<=i2;i++) for(int j=j1;j<=j2;j++) {
+    double x = (double) i;
+    double y = (double) j;
+    double f = (*mMtx)(i,j);
+    double e = 1.;
+    e = GeneralFitData::ComputeError(f,e,errtype,errscale,errmin);
+    mGData->AddData2(x,y,f,e);
+  }
+}
+
+return mGData;
+}
+
+template <class T> 
+AnyDataObj* NOMAdapter_TMatrix<T>::FitResidusObj(GeneralFit& mfit)
+{
+TVector<T>* vec = dynamic_cast<TVector<T> *>(mMtx);
+if(vec) {
+  TVector<T>* v = new TVector<T>(ObjectFitter::FitResidus(*vec,mfit),true);
+  return v;
+} else {
+  TMatrix<T>* m = new TMatrix<T>(ObjectFitter::FitResidus(*mMtx,mfit),true);
+  return m;
+}
+}
+
+template <class T> 
+AnyDataObj* NOMAdapter_TMatrix<T>::FitFunctionObj(GeneralFit& mfit)
+{
+TVector<T>* vec = dynamic_cast<TVector<T> *>(mMtx);
+if(vec) {
+  TVector<T>* v = new TVector<T>(ObjectFitter::FitFunction(*vec,mfit),true);
+  return v;
+} else {
+  TMatrix<T>* m = new TMatrix<T>(ObjectFitter::FitFunction(*mMtx,mfit),true);
+  return m;
+}
+}
 
 
@@ -226 +307 @@
 #pragma define_template NOMAdapter_TMatrix<r_4>
 #pragma define_template NOMAdapter_TMatrix<r_8>
-#pragma define_template NOMAdapter_TMatrix< complex<r_4> >
-#pragma define_template NOMAdapter_TMatrix< complex<r_8> >
+//#pragma define_template NOMAdapter_TMatrix< complex<r_4> >
+//#pragma define_template NOMAdapter_TMatrix< complex<r_8> >
 //#pragma define_template NTupInt_TMatrix<uint_2>
 //#pragma define_template NTupInt_TMatrix<int_2>
@@ -243 +324 @@
 template class NOMAdapter_TMatrix<r_4>;
 template class NOMAdapter_TMatrix<r_8>;
-template class NOMAdapter_TMatrix< complex<r_4> >;
-template class NOMAdapter_TMatrix< complex<r_8> >;
+//template class NOMAdapter_TMatrix< complex<r_4> >;
+//template class NOMAdapter_TMatrix< complex<r_8> >;
 // template class NTupInt_TMatrix<uint_2>;
 // template class NTupInt_TMatrix<int_2>;