Changeset 1270 for applications/robust/robustlib.h

Timestamp:

01/07/11 18:05:28 (13 years ago)

Author:

sindj

Message:

Opravena chyba s chybejici podminkou, zbyva doresit numericke problemy. JS

Files:

• applications/robust/robustlib.h (modified) (5 diffs)

applications/robust/robustlib.h

@@ -19 +19 @@
 using namespace itpp;
 
-const double max_range = 99999999.0;//numeric_limits<double>::max()/10e-10;
+const double max_range = 1000.0;//numeric_limits<double>::max()/10e-10;
 
 enum actions {MERGE, SPLIT};
@@ -726 +726 @@
                                                         else
                                                         {
+                                                                ((toprow*)current_parent)->condition_order++;
+
                                                                 if(current_parent->negativechildren.size()>0)
                                                                 {
@@ -736 +738 @@
                                                                         current_parent->set_state(1, SPLIT);
 
-                                                                        ((toprow*)current_parent)->condition+=toadd;
+                                                                        ((toprow*)current_parent)->condition+=toadd;                                                                    
                                                                 }
                                                                 else
@@ -966 +968 @@
                                         else
                                         {
-                                                toprow* neutral_toprow = new toprow();                                          
+                                                toprow* neutral_toprow = new toprow();
+
+                                                neutral_toprow->condition       = zeros(number_of_parameters+1);
+                                                neutral_toprow->condition_order = ((toprow*)current_polyhedron)->condition_order+1;
 
                                                 new_totally_neutral_child = neutral_toprow;
@@ -984 +989 @@
                                         }
 
-                                        cout << ((toprow*)current_polyhedron)->condition << endl << toadd << endl;
+                                        // cout << ((toprow*)current_polyhedron)->condition << endl << toadd << endl;
 
                                         toprow* positive_poly = new toprow(((toprow*)current_polyhedron)->condition+toadd, ((toprow*)current_polyhedron)->condition_order+1);