Changeset 456f083 in tol

Timestamp:

Oct 23, 2011, 10:11:40 PM (15 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

1941-docker-image-to-test, devel, master

Children:

8c7c57b

Parents:

4d26434

Message:

Relaxing comparisons to avoid numerical problems with truncation intervals near to zero density

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@3778 be4846ae-f051-483a-831f-b64e7d08527f

File:

• tol/bmath/bprdist/prd.cpp (modified) (3 diffs)

tol/bmath/bprdist/prd.cpp

@@ -792 +792 @@
   if(useLeft_)
   {
-    FAB_ = FB_-FA_;
+    FAB_ = (FB_-FA_+1.E+4)-1.E+4;
   }
   else
   {
-    FAB_ = F1A_-F1B_;
+    FAB_ = (F1A_-F1B_+1.E+4)-1.E+4;
   }
   if(FAB_<0)
@@ -838 +838 @@
 {
   if(wrongParameter_ || prob.IsUnknown()) { return(BDat::Unknown()); }
-  if(A_==B_)  { return(A_); }
-  if(!FAB_) 
+  if(A_+1.0==B_+1.0)  { return(A_); }
+  if(FAB_<1.E-8) 
   { 
     BDat C = NonTruncated().Average();
-         if(A_>C) { return(A_); }
-    else if(B_<C) { return(B_); }
+         if(A_+1.0>C+1.0) { return(A_); }
+    else if(B_+1.0<C+1.0) { return(B_); }
     {
       BDat a = A_;
@@ -860 +860 @@
 //--------------------------------------------------------------------
 { 
-  if(A_==B_)  { return(A_); }
+  if(A_+1.0==B_+1.0)  { return(A_); }
   return(Inverse(BUniformDist::Random01())); 
 }