Changeset 2635 in Sophya for trunk/SophyaLib/Samba

Timestamp:

Nov 19, 2004, 12:24:22 PM (21 years ago)

Author:

cmv

Message:

amelioration de la normalisation et protection cmv 19/11/04

File:

• trunk/SophyaLib/Samba/harmspher.cc (modified) (1 diff)

trunk/SophyaLib/Samba/harmspher.cc

@@ -145 +145 @@
  if(expos!=0.) {
    for(int i=0;i<=lmax-m;i++) {
-     double e = log10(fabs(xlm[i]));
-     if(lp>1) cout<<"     xlm["<<i<<"]="<<xlm[i]<<" (e="<<e<<")";
-     if(expos+e<EXPO_MINI) xlm[i] = 0.;
-     else if(expos>EXPO_MINI) xlm[i] *= pow(10.,expos);
-     else xlm[i] = pow(10.,expos+e);
+     double fxml=fabs(xlm[i]);
+     if(lp>1) cout<<"     xlm["<<i<<"]="<<xlm[i];
+     if(fxml>0.) {
+       double e = log10(fxml);
+       if(lp>1) cout<<" (e="<<e<<")";
+       if(expos+e<EXPO_MINI) xlm[i] = 0.;
+       else if(expos>EXPO_MINI) xlm[i] *= pow(10.,expos);
+       else xlm[i] = pow(10.,expos+e);
+     }
      if(lp>1) cout<<" ---> "<<xlm[i]<<endl;
    }