Changeset 833 in Sophya
- Timestamp:
- Apr 7, 2000, 12:28:03 PM (25 years ago)
- Location:
- trunk
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaExt/FitsIOServer/fitsioserver.cc
r721 r833 328 328 long repeat; 329 329 int ntest= 0; 330 for(int ii = 0; ii < ncols; ii++) 330 int ii; 331 for(ii = 0; ii < ncols; ii++) 331 332 { 332 333 int DTYPE; … … 353 354 char **clname; 354 355 clname= new char*[tfields.size()]; 355 for(i nt ii = 0; ii < tfields.size(); ii++)356 for(ii = 0; ii < tfields.size(); ii++) 356 357 { 357 358 // ostringstream ne passe pas sur tous les compilos .. Reza 23/12/99 … … 366 367 // create a NTuple 367 368 NTuple nt0(tfields.size(),clname); 368 for(i nt ii = 0; ii < tfields.size(); ii++)369 for(ii = 0; ii < tfields.size(); ii++) 369 370 delete [] clname[ii]; 370 371 delete [] clname; -
trunk/SophyaLib/Samba/alm.cc
r729 r833 29 29 // --- smoothes the initial power spectrum --- 30 30 TVector<T> cl=clin; 31 for (int l=0;l<n_l;l++) 31 int l; 32 for (l=0;l<n_l;l++) 32 33 { 33 34 r_8 gauss=exp(-l*(l+1.)*sig_smooth*sig_smooth); … … 39 40 r_8 hsqrt2 = 1.0 / Rac2; 40 41 41 for ( intl=0;l<n_l;l++)42 for (l=0;l<n_l;l++) 42 43 { 43 44 T rms=sqrt(cl(l)); -
trunk/SophyaLib/Samba/sphericaltransformserver.cc
r746 r833 80 80 LambdaLMBuilder lb(theta,nlmax,nmmax); 81 81 // somme sur m de 0 a l'infini 82 for (int m = 0; m <= nmmax; m++) 82 int m; 83 for (m = 0; m <= nmmax; m++) 83 84 { 84 85 // somme sur l de m a l'infini … … 104 105 // obtains the negative m of b(m,theta) (= complex conjugate) 105 106 106 for ( intm=1;m<=nmmax;m++)107 for (m=1;m<=nmmax;m++) 107 108 { 108 109 //compiler doesn't have conj() … … 169 170 170 171 for (int kk=0; kk<bw.NElts(); kk++) bw(kk)=(T)0.; 171 for (int m=-b_m.Mmax();m<=-1;m++) 172 int m; 173 for (m=-b_m.Mmax();m<=-1;m++) 172 174 { 173 175 int maux=m; … … 177 179 bw(iw) += b_m(m) * complex<T>( (T)cos(aux),(T)sin(aux) ) ; 178 180 } 179 for ( intm=0;m<=b_m.Mmax();m++)181 for (m=0;m<=b_m.Mmax();m++) 180 182 { 181 183 // int iw=((m % nph) +nph) % nph; //between 0 and nph = m' … … 248 250 249 251 for (int kk=0; kk<bw.NElts(); kk++) bw(kk)=(T)0.; 250 for (int m=-b_m.Mmax();m<=-1;m++) 252 int m; 253 for (m=-b_m.Mmax();m<=-1;m++) 251 254 { 252 255 int maux=m; … … 256 259 bw(iw) += b_m(m) * complex<T>( (T)cos(aux),(T)sin(aux) ) ; 257 260 } 258 for ( intm=0;m<=b_m.Mmax();m++)261 for (m=0;m<=b_m.Mmax();m++) 259 262 { 260 263 // int iw=((m % nph) +nph) % nph; //between 0 and nph = m' … … 365 368 366 369 int im_max=min(nph,nmmax+1); 367 for (int i=0;i< dataout.NElts();i++) dataout(i)=complex<T>((T)0.,(T)0.); 368 for (int i=0;i<im_max;i++) dataout(i)=transformedData(i); 370 int i; 371 for (i=0;i< dataout.NElts();i++) dataout(i)=complex<T>((T)0.,(T)0.); 372 for (i=0;i<im_max;i++) dataout(i)=transformedData(i); 369 373 370 374 … … 373 377 // } 374 378 for (int kk=nph; kk<dataout.NElts(); kk++) dataout(kk)=dataout(kk%nph); 375 for (i nt i= 0;i <dataout.NElts();i++){379 for (i = 0;i <dataout.NElts();i++){ 376 380 dataout(i)*= (complex<T>)(complex<double>(cos(-i*phi0),sin(-i*phi0))); 377 381 } … … 415 419 int sizeOfTransformToGet = min(transformedData.NElts(),nmmax+1); 416 420 // int im_max=min(transformedData.NElts()-1,nmmax); 417 for (int i=0;i<sizeOfTransformToGet;i++) dataout(i)=transformedData(i); 421 int i; 422 for (i=0;i<sizeOfTransformToGet;i++) dataout(i)=transformedData(i); 418 423 419 424 … … 422 427 if (sizeOfTransformToGet == (transformedData.NElts())) 423 428 { 424 for (i nt i=transformedData.NElts(); i<min(nph,dataout.NElts()); i++)429 for (i=transformedData.NElts(); i<min(nph,dataout.NElts()); i++) 425 430 { 426 431 … … 434 439 } 435 440 } 436 for (i nt i= 0;i <dataout.NElts();i++){441 for (i = 0;i <dataout.NElts();i++){ 437 442 dataout(i)*= (complex<T>)(complex<double>(cos(-i*phi0),sin(-i*phi0))); 438 443 } … … 640 645 phasem(i)=0; 641 646 } 642 643 for( intkk=0; kk<nph; kk++) datain(kk)=complex<T>(dataq(kk),datau(kk));647 int kk; 648 for(kk=0; kk<nph; kk++) datain(kk)=complex<T>(dataq(kk),datau(kk)); 644 649 645 650 phasep = CFromFourierAnalysis(nmmax,datain,phi0); 646 651 647 for( intkk=0; kk<nph; kk++) datain(kk)=complex<T>(dataq(kk),-datau(kk));652 for(kk=0; kk<nph; kk++) datain(kk)=complex<T>(dataq(kk),-datau(kk)); 648 653 phasem = CFromFourierAnalysis(nmmax,datain,phi0); 649 654 LambdaPMBuilder lpmb(theta,nlmax,nmmax); … … 695 700 // LambdaPMBuilder lpmb(theta,nlmax,nmmax); 696 701 r_8 sqr2inv=1/Rac2; 697 for (int m = 0; m <= nmmax; m++) 702 int m; 703 for (m = 0; m <= nmmax; m++) 698 704 { 699 705 r_8 lambda_w=0.; … … 718 724 } 719 725 // obtains the negative m of b(m,theta) (= complex conjugate) 720 for ( intm=1;m<=nmmax;m++)726 for (m=1;m<=nmmax;m++) 721 727 { 722 728 b_m_theta_q(-m) = conj(b_m_theta_q(m)); … … 764 770 765 771 LambdaPMBuilder lpmb(theta,nlmax,nmmax); 766 for (int m = 0; m <= nmmax; m++) 772 int m; 773 for (m = 0; m <= nmmax; m++) 767 774 { 768 775 r_8 lambda_p=0.; … … 784 791 785 792 // obtains the negative m of b(m,theta) (= complex conjugate) 786 for ( intm=1;m<=nmmax;m++)793 for (m=1;m<=nmmax;m++) 787 794 { 788 795 b_m_theta_p(-m) = conj(b_m_theta_m(m)); -
trunk/SophyaProg/PMixer/easyTest.cc
r761 r833 57 57 Vector vecOfNu(maxVal); 58 58 Vector vecOfFDeNu(maxVal); 59 for (int i=0; i<maxVal; i++) 59 int i; 60 for (i=0; i<maxVal; i++) 60 61 { 61 62 vecOfNu(i) = (double)i+230.; … … 79 80 Vector vecOfNuRSV(SRVmaxVal); 80 81 Vector vecOfFDeNuRSV(SRVmaxVal); 81 for (i nt i=0; i<SRVmaxVal; i++)82 for (i=0; i<SRVmaxVal; i++) 82 83 { 83 84 vecOfNuRSV(i) = (double)i+130.;
Note:
See TracChangeset
for help on using the changeset viewer.