Changeset 3789 in Sophya for trunk/Cosmo/RadioBeam/mdish.cc

Timestamp:

Jun 27, 2010, 4:57:36 PM (15 years ago)

Author:

ansari

Message:

Corrections diverses, Reza 27/06/2010

File:

• trunk/Cosmo/RadioBeam/mdish.cc (modified) (3 diffs)

trunk/Cosmo/RadioBeam/mdish.cc

@@ -9 +9 @@
 //--------------------------------------------------
 // Constructor
-Four2DResponse::Four2DResponse(int typ, double dx, double dy)
+Four2DResponse::Four2DResponse(int typ, double dx, double dy, double lambda)
   : typ_(typ), dx_((dx>1.e-3)?dx:1.), dy_((dy>1.e-3)?dy:1.)
 {
-  setLambdaRef();
-  setLambda();
+  setLambdaRef(lambda);
+  setLambda(lambda);
 }
 
@@ -78 +78 @@
 // -- Four2DRespRatio : rapport de la reponse entre deux objets Four2DResponse
 //---------------------------------------------------------------
-Four2DRespRatio::Four2DRespRatio(Four2DResponse& a, Four2DResponse& b)
-  : Four2DResponse(0, a.D(), a.D()), a_(a), b_(b) 
+Four2DRespRatio::Four2DRespRatio(Four2DResponse& a, Four2DResponse& b, double divzthr)
+  : Four2DResponse(0, a.D(), a.D()), a_(a), b_(b), divzthr_(divzthr)
 {
 }
@@ -87 +87 @@
   double ra = a_.Value(kx,ky);
   double rb = b_.Value(kx,ky);
-  if (rb<1.e-19) rb = 1.e-19;
+  if (rb<divzthr_) {
+    if (ra<rb)  return 0.;
+    else rb=divzthr_;
+  }
   return (ra/rb);
 }