Changeset 4061 in Sophya for trunk/Cosmo/SimLSS/genefluct3d.cc

Timestamp:

Apr 17, 2012, 11:21:23 AM (13 years ago)

Author:

cmv

Message:

cmv, 17/04/2012

File:

• trunk/Cosmo/SimLSS/genefluct3d.cc (modified) (3 diffs)

trunk/Cosmo/SimLSS/genefluct3d.cc

@@ -810 +810 @@
    for(long jj=0;jj<2;jj++) { // selon j
      long j=0;
-     if(jj==1) {if( Ny_%2!=0) continue; else j = Ny_/2;}
+     if(jj==1) {if(Ny_%2!=0) continue; else j = Ny_/2;}
      for(long i=1;i<(Nx_+1)/2;i++) {
        int_8 ip = IndexC(i,j,k);
@@ -820 +820 @@
    for(long ii=0;ii<2;ii++) {
      long i=0;
-     if(ii==1) {if( Nx_%2!=0) continue; else i = Nx_/2;}
+     if(ii==1) {if(Nx_%2!=0) continue; else i = Nx_/2;}
      for(long j=1;j<(Ny_+1)/2;j++) {
        int_8 ip = IndexC(i,j,k);
@@ -884 +884 @@
 // Gestion divergence en 0: sin(y)/y = 1 - y^2/6*(1-y^2/20)
 //                          avec y = k_x*dx/2
+//--- Ex a 1D: TF(k) = sinc(k D/2) = sin(k D/2) / (k D/2)
+// Le premier zero du sinc, cad k1 tq TF(k1)=0 est a: k1 D/2 = Pi -> k1 = 2Pi/D
+// Avec un echantillonage de pas D, on a: k_nyq = Pi/D = k1/2
+//   donc pour la frequence de Nyquist: TF(k_nyq) = sinc(Pi/2) = 2/Pi ~= 0.6366
+//---
 {
  if(lp_>0) cout<<"--- FilterByPixel ---"<<endl;