Changeset 3520 in Sophya for trunk/SophyaPI/PIext/pitvmaad.cc

Timestamp:

Sep 11, 2008, 2:45:01 PM (17 years ago)

Author:

ansari

Message:

Ajout/codage methode P2DArrayAdapter::MeanVal() pour les classes adapter heritant de P2DArrayAdapter, Reza 11/09/2008

File:

• trunk/SophyaPI/PIext/pitvmaad.cc (modified) (2 diffs)

trunk/SophyaPI/PIext/pitvmaad.cc

@@ -71 +71 @@
 
 /* --Methode-- */
+template <class T> 
+double POTMatrixAdapter<T>::MeanVal(int ix1, int ix2, int jy1, int jy2)
+{
+int ec;
+if (ix1>ix2) { ec=ix1; ix1=ix2; ix2=ec; }
+if (jy1>jy2) { ec=jy1; jy1=jy2; jy2=ec; }    
+double ss = 0.;
+for(int j=jy1; j<=jy2; j++) 
+  for(int i=ix1; i<=ix2; i++)  ss += (double)((*mMtx)(j, i));
+ss /= (double)((jy2-jy1+1)*(ix2-ix1+1));
+return ss;  
+}   
+
+
+/* --Methode-- */
 DECL_TEMP_SPEC  /* equivalent a template <> , pour SGI-CC en particulier */
 double POTMatrixAdapter< complex<float> >::Value(int ix, int iy)
@@ -89 +104 @@
 return(sqrt(re*re+im*im));
 }
+
+/* --Methode-- */
+DECL_TEMP_SPEC  /* equivalent a template <> , pour SGI-CC en particulier */
+double POTMatrixAdapter< complex<float> >::MeanVal(int ix1, int ix2, int jy1, int jy2)
+{
+int ec;
+if (ix1>ix2) { ec=ix1; ix1=ix2; ix2=ec; }
+if (jy1>jy2) { ec=jy1; jy1=jy2; jy2=ec; }    
+complex<float> ss = 0.;
+for(int j=jy1; j<=jy2; j++) 
+  for(int i=ix1; i<=ix2; i++)  ss += (*mMtx)(j, i);
+ss /= (double)((jy2-jy1+1)*(ix2-ix1+1));
+return sqrt(ss.real()*ss.real()+ss.imag()*ss.imag());  
+}   
+
+/* --Methode-- */
+DECL_TEMP_SPEC  /* equivalent a template <> , pour SGI-CC en particulier */
+double POTMatrixAdapter< complex<double> >::MeanVal(int ix1, int ix2, int jy1, int jy2)
+{
+int ec;
+if (ix1>ix2) { ec=ix1; ix1=ix2; ix2=ec; }
+if (jy1>jy2) { ec=jy1; jy1=jy2; jy2=ec; }    
+complex<double> ss = 0.;
+for(int j=jy1; j<=jy2; j++) 
+  for(int i=ix1; i<=ix2; i++)  ss += (*mMtx)(j, i);
+ss /= (double)((jy2-jy1+1)*(ix2-ix1+1));
+return sqrt(ss.real()*ss.real()+ss.imag()*ss.imag());  
+}   
 
 #ifdef __CXX_PRAGMA_TEMPLATES__