Changeset 1528 in Sophya for trunk/SophyaProg/PrgMap/extrapsph.cc

Timestamp:

Jun 14, 2001, 5:16:20 PM (24 years ago)

Author:

cmv

Message:

optimisation cmv 14/6/01

File:

• trunk/SophyaProg/PrgMap/extrapsph.cc (modified) (10 diffs)

trunk/SophyaProg/PrgMap/extrapsph.cc

@@ -1 +1 @@
 // Pour boucher les trous d'une sphere HEALPIX en clusterisant dans une sphere
 // plus petite             cmv 13/6/01
+// extrapsph -r 2 sph143k05.fits sph143k05_e.fits
+//                sphout.fits sphoutw.fits sphred.fits sphredw.fits
 #include "machdefs.h"
 #include <unistd.h>
@@ -17 +19 @@
     <<"                     sphout.fits [sphout_w.fits"
     <<" sphred.fits sphred_w.fits]"<<endl
-    <<" -r reduc : reduction factor for nlat in clustering must be 2^n (def=2)"<<endl
-    <<" sphin.fits : input sphere"<<endl
-    <<" sphin_w.fits : input sphere filling weight"<<endl
-    <<" sphout.fits  : output sphere"<<endl
-    <<" sphout_w.fits  : output sphere filling weight"<<endl
-    <<" sphred.fits  : output reducted sphere"<<endl
-    <<" sphred_w.fits  : output reducted sphere filling weight"<<endl;
+    <<" -r reduc : reduction factor for clustering (must be 2^n, def=2)"<<endl
+    <<" sphin.fits    : input sphere"<<endl
+    <<" sphin_w.fits  : input sphere filling weight"<<endl
+    <<" sphout.fits   : output sphere"<<endl
+    <<" sphout_w.fits : output sphere filling weight"<<endl
+    <<" sphred.fits   : output reducted sphere"<<endl
+    <<" sphred_w.fits : output reducted sphere filling weight"<<endl;
 }
 
@@ -52 +54 @@
 if(optind+5<narg) fsphredw = arg[optind+5];
 
-cout<<"Sphere Input : "<<fsphin<<endl
-    <<"Weight Sphere Input : "<<fsphinw<<endl
-    <<"Sphere Output : "<<fsphout<<endl
-    <<"Weight Sphere Output : "<<fsphoutw<<endl
-    <<"Reducted Sphere Output : "<<fsphred<<endl
-    <<"Reducted Weight Sphere Output : "<<fsphredw<<endl
+cout<<"Input Sphere : "<<fsphin<<endl
+    <<"Input Weight Sphere : "<<fsphinw<<endl
+    <<"Output Sphere : "<<fsphout<<endl
+    <<"Output Weight Sphere : "<<fsphoutw<<endl
+    <<"Output Reducted Sphere : "<<fsphred<<endl
+    <<"Output Reducted Weight Sphere : "<<fsphredw<<endl
     <<"Reduction : "<<red<<endl;
 
@@ -65 +67 @@
 int nlat = sphin.SizeIndex();
 int nlatr = nlat/red;
-cout<<"Opening Sphere Input :"<<endl
+cout<<"Opening Input Sphere :"<<endl
     <<"          Type of map : "<<sphin.TypeOfMap()<<endl
     <<"     Number of pixels : "<<sphin.NbPixels()<<endl
@@ -72 +74 @@
 SphereHEALPix<r_8> sphinw;
 {FitsInFile sfits(fsphinw); sfits >> sphinw;}
-cout<<"Opening Weight Sphere Input :"<<endl
+cout<<"Opening Input Weight Sphere :"<<endl
     <<"          Type of map : "<<sphinw.TypeOfMap()<<endl
     <<"     Number of pixels : "<<sphinw.NbPixels()<<endl
@@ -81 +83 @@
 }
 
-// Creation des spheres output
-cout<<"Creating Spheres for output"<<endl;
-SphereHEALPix<r_8> sphout(sphin,false);
-SphereHEALPix<r_8> sphoutw(sphinw,false);
-
 // Creation des spheres reduites
 cout<<"Creating Reducted Spheres : nlatr = "<<nlatr<<endl;
@@ -92 +89 @@
 SphereHEALPix<r_8> sphredw(nlatr);
   sphredw.SetPixels(0.);
-cout<<"Creating Reducted Sphere :"<<endl
-    <<"     Number of pixels : "<<sphred.NbPixels()<<endl
+cout<<"     Number of pixels : "<<sphred.NbPixels()<<endl
     <<"                 Nlat : "<<sphred.SizeIndex()<<endl;
 
 // Filling reducted spheres
 cout<<"...Filling Reducted Spheres"<<endl;
+uint_4 n=0;
 for(int_4 i=0;i<sphin.NbPixels();i++) {
   if(sphinw(i)<=0.) continue;
@@ -105 +102 @@
   sphred(ir)  += sphin(i)*sphinw(i);
   sphredw(ir) += sphinw(i);
+  n++;
 }
+cout<<"      Input Sphere: Number of filled pixels "<<n<<endl;
 cout<<"...Computing Reducted Spheres"<<endl;
-uint_4 n=0;
+n=0;
 for(int_4 ir=0;ir<sphred.NbPixels();ir++)
   if(sphredw(ir) > 0.) {sphred(ir) /= sphredw(ir); n++;}
-cout<<"      Number of filled pixels "<<n<<endl;
+cout<<"      Reducted Sphere: Number of pixels filled "<<n<<endl;
 
 // Filling hole for Output Spheres
 cout<<"...Filling hole for Output Spheres"<<endl;
 n=0;
-for(int_4 i=0;i<sphout.NbPixels();i++) {
-  if(sphoutw(i)>0.) continue;
+for(int_4 i=0;i<sphin.NbPixels();i++) {
+  if(sphinw(i)>0.) continue;
   double theta,phi;
-  sphout.PixThetaPhi(i,theta,phi);
+  sphin.PixThetaPhi(i,theta,phi);
   int_4 ir = sphred.PixIndexSph(theta,phi);
   if(sphredw(ir)<=0.) continue;
-  sphout(i)  = sphred(ir);
-  sphoutw(i) = -sphredw(ir);
+  sphin(i)  = sphred(ir);
+  sphinw(i) = -sphredw(ir);
   n++;
 }
-cout<<"      Number of filled pixels "<<n<<endl;
+cout<<"      Output Sphere: Number of pixels filled "<<n<<endl;
 
 // Ecriture des spheres
@@ -132 +131 @@
 FitsOutFile sfits(fsphout);
 cout<<"Writing Output Sphere Fits file"<<endl;
-sfits<<sphout;
+sfits<<sphin;
 }
 if(fsphoutw) {
@@ -138 +137 @@
 FitsOutFile sfits(fsphoutw);
 cout<<"Writing Output Sphere Weight Fits file"<<endl;
-sfits<<sphoutw;
+sfits<<sphinw;
 }
 if(fsphred) {