Changeset 1545 in Sophya for trunk/SophyaProg/PrgMap


Ignore:
Timestamp:
Jun 18, 2001, 6:07:04 PM (24 years ago)
Author:
ansari
Message:

compil/portage SGI-CC , Reza 18/6/2001

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SophyaProg/PrgMap/extrapsph.cc

    r1528 r1545  
    9595cout<<"...Filling Reducted Spheres"<<endl;
    9696uint_4 n=0;
    97 for(int_4 i=0;i<sphin.NbPixels();i++) {
     97int_4 i, ir;
     98for(i=0;i<sphin.NbPixels();i++) {
    9899  if(sphinw(i)<=0.) continue;
    99100  double theta,phi;
    100101  sphin.PixThetaPhi(i,theta,phi);
    101   int_4 ir = sphred.PixIndexSph(theta,phi);
     102  ir = sphred.PixIndexSph(theta,phi);
    102103  sphred(ir)  += sphin(i)*sphinw(i);
    103104  sphredw(ir) += sphinw(i);
     
    107108cout<<"...Computing Reducted Spheres"<<endl;
    108109n=0;
    109 for(int_4 ir=0;ir<sphred.NbPixels();ir++)
     110for(ir=0;ir<sphred.NbPixels();ir++)
    110111  if(sphredw(ir) > 0.) {sphred(ir) /= sphredw(ir); n++;}
    111112cout<<"      Reducted Sphere: Number of pixels filled "<<n<<endl;
     
    114115cout<<"...Filling hole for Output Spheres"<<endl;
    115116n=0;
    116 for(int_4 i=0;i<sphin.NbPixels();i++) {
     117for(i=0;i<sphin.NbPixels();i++) {
    117118  if(sphinw(i)>0.) continue;
    118119  double theta,phi;
Note: See TracChangeset for help on using the changeset viewer.