Changeset f4f6628 in tol

Timestamp:

Nov 24, 2011, 10:33:16 AM (15 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

35aac13

Parents:

523e2d1

Message:

Refs #1342
Forcing initial interior point for main linear block

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

File:

• tol/btol/vmatrix_type/vmat_stats.cpp (modified) (4 diffs)

tol/btol/vmatrix_type/vmat_stats.cpp

@@ -715 +715 @@
       double lowMrg = lower;
       double uppMrg = upper;
-      if(lower>=K)
+      if(lowMrg>K)
       { 
+        if(uppMrg>=2*lowMrg) { uppMrg=2*lowMrg; }
         //Si todo el dominio está por encima de K tomamos como nuevos límites
         //tentativos los puntos con densidad 1.1 y 110 veces las del límite inferior.
         //De esta forma el nuevo límite inferior es por construcción 100 veces más 
         //probable que el superior 
-        double lowMrg_ = sqrt(pow(lower,2) + 2.0 * log(  1.1) );
-        double uppMrg_ = sqrt(pow(upper,2) + 2.0 * log(110.0) );
+        double lowMrg_ = sqrt(pow(lowMrg,2) + 2.0 * log(  1.1) );
+        double uppMrg_ = sqrt(pow(uppMrg,2) + 2.0 * log(110.0) );
         //Si los límites tentativos se salen del dominio hay que meterlos dentro 
         if(lowMrg_>=upper)
         {
-          lowMrg = 0.99*lower+0.01*upper;
-          uppMrg = 0.99*upper+0.01*lower;
+          lowMrg = 0.99*lowMrg+0.01*uppMrg;
+          uppMrg = 0.99*uppMrg+0.01*lowMrg;
         }
         else 
@@ -734 +735 @@
           if(uppMrg_>=upper)
           {
-            uppMrg = 0.99*upper+0.01*lowMrg;
+            uppMrg = 0.99*uppMrg+0.01*lowMrg;
           }
           else
@@ -756 +757 @@
         }
       }
-      else if(upper<=-K)
+      else if(upper<-K)
       {
+        if(lowMrg<=2*uppMrg) { lowMrg=2*uppMrg; }
         //Si todo el dominio está por debajo de -K tomamos como nuevos límites
         //tentativos los puntos con densidad 1.1 y 110 veces las del límite superior.
         //De esta forma el nuevo límite superior es por construcción 100 veces más 
         //probable que el inferior 
-        double uppMrg_ = -sqrt(pow(upper,2) + 2.0 * log(  1.1) );
-        double lowMrg_ = -sqrt(pow(lower,2) + 2.0 * log(110.0) );
+        double uppMrg_ = -sqrt(pow(uppMrg,2) + 2.0 * log(  1.1) );
+        double lowMrg_ = -sqrt(pow(lowMrg,2) + 2.0 * log(110.0) );
         //Si los límites tentativos se salen del dominio hay que meterlos dentro 
         if(uppMrg_<=lower)
         {
-          lowMrg = 0.99*lower+0.01*upper;
-          uppMrg = 0.99*upper+0.01*lower;
+          lowMrg = 0.99*lowMrg+0.01*uppMrg;
+          uppMrg = 0.99*uppMrg+0.01*lowMrg;
         }
         else 
@@ -775 +777 @@
           if(lowMrg_<=lower)
           {
-            lowMrg = 0.99*lower+0.01*uppMrg;
+            lowMrg = 0.99*lowMrg+0.01*uppMrg;
           }
           else