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

Timestamp:

Jun 25, 2010, 12:00:30 PM (15 years ago)

Author:

ansari

Message:

Ajout classes / programmes de calcul d'effet de lobe sur les radio-sources, Reza 25/06/2010

File:

• trunk/Cosmo/RadioBeam/specpk.cc (modified) (2 diffs)

trunk/Cosmo/RadioBeam/specpk.cc

@@ -151 +151 @@
   // fourAmp represent 3-D fourier transform of a real input array. 
   // The second half of the array along Y and Z contain negative frequencies
-  double kxx, kyy, kzz;
+  double kxx, kyy, kzz, rep, amp;
   // sa_size_t is large integer type  
   for(sa_size_t kz=0; kz<fourAmp.SizeZ(); kz++) {
@@ -158 +158 @@
       kyy =  (ky>fourAmp.SizeY()/2) ? -(double)(fourAmp.SizeY()-ky)*dky_ : (double)ky*dky_; 
       for(sa_size_t kx=0; kx<fourAmp.SizeX(); kx++) {
-        double kxx=(double)kx*dkx_;
-        double rep = resp(kxx, kyy);
+        kxx=(double)kx*dkx_;
+        rep = resp(kxx, kyy);
         if (crmask&&(kz==0))  mask(ky,kx)=((rep<1.e-8)?9.e9:(1./rep));
         if (rep<1.e-8)  fourAmp(kx, ky, kz) = complex<TF>(9.e9,0.);
         else {
-          double amp = 1./sqrt(rep)/sqrt(2.);
+          amp = 1./sqrt(rep)/sqrt(2.);
           fourAmp(kx, ky, kz) = complex<TF>(rg_.Gaussian(amp), rg_.Gaussian(amp));   
         }