Changeset 3507 in Sophya for trunk/SophyaProg/Tests

Timestamp:

Aug 8, 2008, 12:53:04 PM (17 years ago)

Author:

ansari

Message:

modif prog test pour appel a SphereECP::ExtractPartial() , Reza 08/08/2008

File:

• trunk/SophyaProg/Tests/tsphereecp.cc (modified) (6 diffs)

trunk/SophyaProg/Tests/tsphereecp.cc

@@ -10 +10 @@
 #include <typeinfo>
 #include "timing.h"
+#include "resusage.h"
+
 #include "sphereecp.h"
 #include "fiosphereecp.h"
@@ -21 +23 @@
   try {
     SophyaInit();
+    ResourceUsage ru;
     
     InitTim();
@@ -31 +34 @@
     ecp2 += 4.;
     cout << " ecp2 = ecp+4. = \n " << ecp2 ;
+    
     {
       cout << " Writing ecp ecp2 to POutPersist sphecp.ppf " << endl;
@@ -45 +49 @@
       cout << recp2;
     }
+    
+    // On cree une carte complete et on la remplit 
+    SphereECP<r_4> ecpf(200);
+    for(int_4 kk=0; kk<ecpf.NbPixels(); kk++) { 
+      double tet, phi;
+      ecpf.PixThetaPhi(kk, tet, phi);
+      ecpf(kk) = sin(tet)*cos(phi);
+    }
+    // Extraction d'une carte partielle de ecpf
+    SphereECP<r_4> ecpex = ecpf.ExtractPartial(M_PI/2.-M_PI/12., M_PI/2.+M_PI/12, 0.23, 0.23+M_PI/2.5);
+    cout << " ecpf (full) and ecpex (extracted partiel) created ... " << endl;
+    cout << ecpf << ecpex ;
+    {
+      cout << " Writing ecpf ecpex to POutPersist sphecpfp.ppf " << endl;
+      POutPersist po("sphecpfp.ppf");
+      po << PPFNameTag("ecpf") << ecpf << PPFNameTag("ecpex") << ecpex;
+    }
+         
     // Test avec transforme Ylm 
     int lmax = 92;
@@ -59 +81 @@
     cout << " Computing Cl from map " << endl;
     Vector clout = ylmserver.DecomposeToCl(map, lmax,  0.);
+    
     SphereECP<r_8> mappar(0.25*M_PI, 0.75*M_PI, 128, 0.25*M_PI, 1.25*M_PI, 256);
     cout << " Generating partial map from Cl " << endl;
@@ -65 +88 @@
     Vector cloutpar = ylmserver.DecomposeToCl(mappar, lmax,  0.);
 
+    cout << " Extracting partial map and Computing Cl from mapex " << endl;
+    SphereECP<r_8> mapex = map.ExtractPartial(0.25*M_PI, 0.75*M_PI, 0.25*M_PI, 1.25*M_PI);
+    Vector clmex = ylmserver.DecomposeToCl(mapex, lmax,  0.);
+
     cout << " Writing clin, map, clout to ecpylm.ppf" << endl;
     POutPersist so("ecpylm.ppf");
     so << PPFNameTag("clin") << clin <<  PPFNameTag("map") << map
        << PPFNameTag("clout") << clout ;
+    cout << " Writing mapex, clmex to ecpylm.ppf" << endl;
+    so << PPFNameTag("clmex") << clmex <<  PPFNameTag("mapex") << mapex;
     cout << " Writing mappar, cloutpar to ecpylm.ppf" << endl;
     so << PPFNameTag("mappar") << mappar
        << PPFNameTag("cloutpar") << cloutpar ;
+    
+    cout << " --> End of tsphereecp computations \n" << ru; 
+   
   }
   catch (PThrowable & exc) {