Changeset 3991 in Sophya

Timestamp:

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

Author:

ansari

Message:

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

Location:

trunk/Cosmo/RadioBeam

Files:

• applobe.cc (modified) (1 diff)
• mdish.cc (modified) (2 diffs)
• mdish.h (modified) (1 diff)

trunk/Cosmo/RadioBeam/applobe.cc

@@ -90 +90 @@
     else fgoptarg=false;
   }
-  if (narg < 3) {
+  if (narg < 4) {
     cout << " applobe/error arguments , applobe -h for help " << endl;
     return 2;

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_;

trunk/Cosmo/RadioBeam/mdish.h

@@ -95 +95 @@
   Four2DResponse& a_;
   Four2DResponse& b_;
-  double maxratio_;
+  double maxratio_, zerothr_;
 };