Changeset 522 in Sophya

Timestamp:

Oct 26, 1999, 7:15:37 PM (26 years ago)

Author:

ansari

Message:

Fix a dumb bug

File:

• trunk/SophyaLib/Samba/ana2fast.cc (modified) (2 diffs)

trunk/SophyaLib/Samba/ana2fast.cc

@@ -105 +105 @@
     //cout << "fft:"<<phas_np[0]<<" "<<phas_sp[0]<<" "<<phas_nm[0]<<" "<<phas_sm[0]<<endl;
     for (int m = 0; m <= nmmax; m++){
-      //      cout << phas_np[m]<<" "<<phas_sp[m]<<" "<<phas_nm[m]<<" "<<phas_sm[m]<<endl;
+      cout << phas_np[m]<<" "<<phas_sp[m]<<" "<<phas_nm[m]<<" "<<phas_sm[m]<<endl;
       a2lmp[m][m]+=l2b.lam2lmp(m,m)*phas_np[m]+l2b.lam2lmp(m,m,-1)*phas_sp[m];
       a2lmm[m][m]+=l2b.lam2lmm(m,m)*phas_nm[m]+l2b.lam2lmm(m,m,-1)*phas_sm[m];
@@ -127 +127 @@
   float domega=(4.*M_PI)/(12.*nsmax*nsmax);
   for (int m = 0; m <= nmmax; m++){
+      a2lme[m][m]=-(a2lmp[m][m]+a2lmm[m][m])/2.*static_cast<double>(domega);
+      a2lmb[m][m]=im*(a2lmp[m][m]-a2lmm[m][m])/2.*static_cast<double>(domega);
     for (int l = m+1; l<= nlmax; l++){
-      a2lme[l][m]=-(a2lmp[l][m]+a2lmm[l][m])/2.*(double)domega;
-      a2lmb[l][m]=im*(a2lmp[l][m]-a2lmm[l][m])/2.*(double)domega;
+      a2lme[l][m]=-(a2lmp[l][m]+a2lmm[l][m])/2.*static_cast<double>(domega);
+      a2lmb[l][m]=im*(a2lmp[l][m]-a2lmm[l][m])/2.*static_cast<double>(domega);
     }
   }