Changeset 3783 in Sophya for trunk/Cosmo/RadioBeam/specpk.cc

Timestamp:

Jun 16, 2010, 12:23:38 AM (15 years ago)

Author:

ansari

Message:

Ajout des programmes calcpk.cc calcpk2.cc syncube.cc tjyk.cc (voir fichier README) Reza 15/06/2010

File:

• trunk/Cosmo/RadioBeam/specpk.cc (modified) (1 diff)

trunk/Cosmo/RadioBeam/specpk.cc

@@ -225 +225 @@
   double kxx, kyy, kzz;
   // sa_size_t is large integer type  
-  for(sa_size_t kz=0; kz<fourAmp.SizeZ(); kz++) {
+  // We ignore 0th term in all frequency directions ...
+  for(sa_size_t kz=1; kz<fourAmp.SizeZ(); kz++) {
     kzz =  (kz > fourAmp.SizeZ()/2) ? (double)(fourAmp.SizeZ()-kz)*dkz_ : (double)kz*dkz_; 
-    for(sa_size_t ky=0; ky<fourAmp.SizeY(); ky++) {
+    for(sa_size_t ky=1; ky<fourAmp.SizeY(); ky++) {
       kyy =  (ky > fourAmp.SizeY()/2) ? (double)(fourAmp.SizeY()-ky)*dky_ : (double)ky*dky_; 
-      for(sa_size_t kx=0; kx<fourAmp.SizeX(); kx++) {
+      for(sa_size_t kx=1; kx<fourAmp.SizeX(); kx++) {  // ignore the 0th coefficient (constant term)
         double kxx=(double)kx*dkx_;
         complex<TF> za = fourAmp(kx, ky, kz);