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

Timestamp:

May 9, 2011, 12:10:57 PM (14 years ago)

Author:

ansari

Message:

encore debug/correction de Four2DRespRatio, Reza 08/05/2011

File:

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

trunk/Cosmo/RadioBeam/mdish.cc

@@ -168 +168 @@
 //---------------------------------------------------------------
 Four2DRespRatio::Four2DRespRatio(Four2DResponse& a, Four2DResponse& b, double maxratio)
-  : Four2DResponse(0, a.D(), a.D()), a_(a), b_(b), maxratio_(maxratio)
+  : Four2DResponse(0, a.D(), a.D()), a_(a), b_(b), maxratio_(maxratio), zerothr_(.5/maxratio)
 {
 }
@@ -176 +176 @@
   double ra = a_.Value(kx,ky);
   double rb = b_.Value(kx,ky);
-  if (ra<rb) {
-    if (rb>1.e-19)  return(ra/rb);  
-    else return 0.;
-  }
-  double rval=maxratio_;
-  if (rb>1.e-19) rval=ra/rb;
+  if ((ra<zerothr_)||(rb<zerothr_))  return 0.;
+  double rval=ra/rb;
   if (rval<maxratio_) return rval;
   return maxratio_;