Changeset 833 in Sophya for trunk/SophyaLib/Samba/sphericaltransformserver.cc

Timestamp:

Apr 7, 2000, 12:28:03 PM (25 years ago)

Author:

ansari

Message:

Compil SGI-CC (Pb for(int i; i redeclare)) Reza 7/4/2000

File:

• trunk/SophyaLib/Samba/sphericaltransformserver.cc (modified) (16 diffs)

trunk/SophyaLib/Samba/sphericaltransformserver.cc

@@ -80 +80 @@
       LambdaLMBuilder lb(theta,nlmax,nmmax);
       //  somme sur m de 0 a l'infini
-      for (int m = 0; m <= nmmax; m++)
+      int m;
+      for (m = 0; m <= nmmax; m++)
         {
           //  somme sur l de m a l'infini
@@ -104 +105 @@
     //        obtains the negative m of b(m,theta) (= complex conjugate)
 
-      for (int m=1;m<=nmmax;m++)
+      for (m=1;m<=nmmax;m++)
         {
           //compiler doesn't have conj()
@@ -169 +170 @@
 
   for (int kk=0; kk<bw.NElts(); kk++) bw(kk)=(T)0.;
-  for (int m=-b_m.Mmax();m<=-1;m++)
+  int m;
+  for (m=-b_m.Mmax();m<=-1;m++)
     {
       int maux=m;
@@ -177 +179 @@
       bw(iw) += b_m(m) * complex<T>( (T)cos(aux),(T)sin(aux) )  ;
     }
-  for (int m=0;m<=b_m.Mmax();m++)
+  for (m=0;m<=b_m.Mmax();m++)
     {
       //      int iw=((m % nph) +nph) % nph; //between 0 and nph = m'
@@ -248 +250 @@
 
   for (int kk=0; kk<bw.NElts(); kk++) bw(kk)=(T)0.;
-  for (int m=-b_m.Mmax();m<=-1;m++)
+  int m;
+  for (m=-b_m.Mmax();m<=-1;m++)
     {
       int maux=m;
@@ -256 +259 @@
       bw(iw) += b_m(m) * complex<T>( (T)cos(aux),(T)sin(aux) )  ;
     }
-  for (int m=0;m<=b_m.Mmax();m++)
+  for (m=0;m<=b_m.Mmax();m++)
     {
       //      int iw=((m % nph) +nph) % nph; //between 0 and nph = m'
@@ -365 +368 @@
 
   int im_max=min(nph,nmmax+1);
-  for (int i=0;i< dataout.NElts();i++) dataout(i)=complex<T>((T)0.,(T)0.);
-  for (int i=0;i<im_max;i++) dataout(i)=transformedData(i);
+  int i;
+  for (i=0;i< dataout.NElts();i++) dataout(i)=complex<T>((T)0.,(T)0.);
+  for (i=0;i<im_max;i++) dataout(i)=transformedData(i);
 
 
@@ -373 +377 @@
   //  }
   for (int kk=nph; kk<dataout.NElts(); kk++) dataout(kk)=dataout(kk%nph);
-  for (int i = 0;i <dataout.NElts();i++){
+  for (i = 0;i <dataout.NElts();i++){
     dataout(i)*= (complex<T>)(complex<double>(cos(-i*phi0),sin(-i*phi0)));
   }
@@ -415 +419 @@
   int sizeOfTransformToGet = min(transformedData.NElts(),nmmax+1);
   //  int im_max=min(transformedData.NElts()-1,nmmax);
-  for (int i=0;i<sizeOfTransformToGet;i++) dataout(i)=transformedData(i);
+  int i;
+  for (i=0;i<sizeOfTransformToGet;i++) dataout(i)=transformedData(i);
 
 
@@ -422 +427 @@
   if (sizeOfTransformToGet == (transformedData.NElts()))
     {
-      for (int i=transformedData.NElts(); i<min(nph,dataout.NElts()); i++)
+      for (i=transformedData.NElts(); i<min(nph,dataout.NElts()); i++)
         {
 
@@ -434 +439 @@
         }
     }
-  for (int i = 0;i <dataout.NElts();i++){
+  for (i = 0;i <dataout.NElts();i++){
     dataout(i)*= (complex<T>)(complex<double>(cos(-i*phi0),sin(-i*phi0)));
   }
@@ -640 +645 @@
       phasem(i)=0; 
     }
-
-  for(int kk=0; kk<nph; kk++) datain(kk)=complex<T>(dataq(kk),datau(kk));
+  int kk;
+  for(kk=0; kk<nph; kk++) datain(kk)=complex<T>(dataq(kk),datau(kk));
 
   phasep = CFromFourierAnalysis(nmmax,datain,phi0); 
 
-  for(int kk=0; kk<nph; kk++) datain(kk)=complex<T>(dataq(kk),-datau(kk));
+  for(kk=0; kk<nph; kk++) datain(kk)=complex<T>(dataq(kk),-datau(kk));
   phasem = CFromFourierAnalysis(nmmax,datain,phi0); 
   LambdaPMBuilder lpmb(theta,nlmax,nmmax);
@@ -695 +700 @@
       //      LambdaPMBuilder lpmb(theta,nlmax,nmmax);
       r_8 sqr2inv=1/Rac2;
-      for (int m = 0; m <= nmmax; m++)
+      int m;
+      for (m = 0; m <= nmmax; m++)
         {
           r_8 lambda_w=0.;
@@ -718 +724 @@
         }
       //        obtains the negative m of b(m,theta) (= complex conjugate)
-      for (int m=1;m<=nmmax;m++)
+      for (m=1;m<=nmmax;m++)
         {
           b_m_theta_q(-m) = conj(b_m_theta_q(m));
@@ -764 +770 @@
 
       LambdaPMBuilder lpmb(theta,nlmax,nmmax);
-      for (int m = 0; m <= nmmax; m++)
+      int m;
+      for (m = 0; m <= nmmax; m++)
         {
           r_8 lambda_p=0.;
@@ -784 +791 @@
       
       //        obtains the negative m of b(m,theta) (= complex conjugate)
-      for (int m=1;m<=nmmax;m++)
+      for (m=1;m<=nmmax;m++)
         {
           b_m_theta_p(-m) = conj(b_m_theta_m(m));