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

Timestamp:

Feb 18, 2010, 10:45:29 AM (16 years ago)

Author:

cmv

Message:

correction pour petits cubes (demande assez de pys pour zref_), cmv 18/02/2010

File:

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

trunk/Cosmo/SimLSS/genefluct3d.cc

@@ -307 +307 @@
  // TO BE FIXED TO BE FIXED TO BE FIXED TO BE FIXED TO BE FIXED TO BE FIXED
  if(loscom_min_<=1.e-50)
-   for(int i=0;i<50;i++)
+   for(int i=0;i<10;i++)
      cout<<"ATTENTION TOUTES LES PARTIES DU CODE NE MARCHENT PAS POUR UN OBSERVATEUR DANS LE CUBE"<<endl;
  // TO BE FIXED TO BE FIXED TO BE FIXED TO BE FIXED TO BE FIXED TO BE FIXED
@@ -337 +337 @@
  // Be shure to have one dlc <loscom_min and one >loscom_max
  if(zinc<=0.) zinc = 0.01;
- for(double z=0.; ; z+=zinc) {
-   double dlc = cosmo_->Dloscom(z);
-   if(dlc<loscom_min_) {zred_.resize(0); loscom_.resize(0);}
-   zred_.push_back(z);
-   loscom_.push_back(dlc);
-   z += zinc;
-   if(dlc>loscom_max_) break; // on sort apres avoir stoque un dlc>dlcmax
+ double zmin = 0., dlcmin=0.;
+ while(1) {
+   if(lp_>0)
+     cout<<"...Filling zred starting at zmin="<<zmin<<" with zinc="<<zinc
+         <<", loscom_min-max=["<<loscom_min_<<","<<loscom_max_<<"]"<<endl;
+   zred_.resize(0); loscom_.resize(0);
+   for(double z=zmin; ; z+=zinc) {
+     double dlc = cosmo_->Dloscom(z);
+     if(dlc<loscom_min_) {
+       dlcmin = dlc;
+       zmin = z;
+       zred_.resize(0); loscom_.resize(0);
+     }
+     zred_.push_back(z);
+     loscom_.push_back(dlc);
+     z += zinc;
+     if(dlc>loscom_max_) {
+       if(lp_>0)
+         cout<<"  Min: z="<<zmin<<" dlc="<<dlcmin<<", Max: z="<<z<<" dlc="<<dlc<<endl;
+       break; // on sort apres avoir stoque un dlc>dlcmax
+     }
+   }
+   if(zred_.size()>=10) break;
+   zinc /= 10.;
+   cout<<"  not enough points ("<<zred_.size()<<") for zref, retry with zinc="<<zinc<<endl;
  }