Changeset 1545 in Sophya for trunk/SophyaProg/PrgMap/extrapsph.cc
- Timestamp:
- Jun 18, 2001, 6:07:04 PM (24 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaProg/PrgMap/extrapsph.cc
r1528 r1545 95 95 cout<<"...Filling Reducted Spheres"<<endl; 96 96 uint_4 n=0; 97 for(int_4 i=0;i<sphin.NbPixels();i++) { 97 int_4 i, ir; 98 for(i=0;i<sphin.NbPixels();i++) { 98 99 if(sphinw(i)<=0.) continue; 99 100 double theta,phi; 100 101 sphin.PixThetaPhi(i,theta,phi); 101 i nt_4 ir = sphred.PixIndexSph(theta,phi);102 ir = sphred.PixIndexSph(theta,phi); 102 103 sphred(ir) += sphin(i)*sphinw(i); 103 104 sphredw(ir) += sphinw(i); … … 107 108 cout<<"...Computing Reducted Spheres"<<endl; 108 109 n=0; 109 for(i nt_4 ir=0;ir<sphred.NbPixels();ir++)110 for(ir=0;ir<sphred.NbPixels();ir++) 110 111 if(sphredw(ir) > 0.) {sphred(ir) /= sphredw(ir); n++;} 111 112 cout<<" Reducted Sphere: Number of pixels filled "<<n<<endl; … … 114 115 cout<<"...Filling hole for Output Spheres"<<endl; 115 116 n=0; 116 for(i nt_4 i=0;i<sphin.NbPixels();i++) {117 for(i=0;i<sphin.NbPixels();i++) { 117 118 if(sphinw(i)>0.) continue; 118 119 double theta,phi;
Note:
See TracChangeset
for help on using the changeset viewer.