Changeset 1262 for applications/robust/robustlib.h

Timestamp:

12/13/10 19:25:45 (14 years ago)

Author:

sindj

Message:

Dodelavani korekcnich faktoru pri stjene hodnote funkce ve dvou bodech polyhedronu. Rozdelano. JS

Files:

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

applications/robust/robustlib.h

@@ -12 +12 @@
 #include <vector>
 #include <list>
+#include <map>
 #include <set>
 #include <algorithm>
+#include <string>
         
 using namespace bdm;
@@ -413 +415 @@
 class emlig // : eEF
 {
+        vector<set<ivec>> correction_factors;
 
         /// A statistic in a form of a Hasse diagram representing a complex of convex polyhedrons obtained as a result
@@ -425 +428 @@
 
         double normalization_factor;
+
+        
 
         void alter_toprow_conditions(vec condition, bool should_be_added)
@@ -915 +920 @@
                         sizevector.push_back(statistic.row_size(s));
                 }*/
-        }
+
+                
+        }
+
+        vector<list<ivec>> get_correction_factors(int order)
+                {
+                        for(int remaining_order = correction_factors.size();remaining_order>order;remaining_order++)
+                        {
+                                set<ivec> factor_templates;
+                                set<ivec> final_factors;
+
+                                for(int i = 1;i==number_of_parameters-order+1;i++)
+                                {
+                                        for(int j = 1;j==remaining_order;j++)
+                                        {
+                                                factor_templates.insert(zeros(number_of_parameters-order+2));
+
+                                                for(set<ivec>::iterator fac_ref = factor_templates.begin();fac_ref!=factor_templates.end();fac_ref++)
+                                                {
+                                                        ivec current_template = (*fac_ref);
+
+                                                        current_template[0]+=1;
+                                                        current_template[i]+=1;
+
+                                                        if(current_template[0]==remaining_order)
+                                                        {
+                                                                final_factors.insert(current_template.right(current_template.size()-1);
+                                                        }
+                                                        else
+                                                        {
+                                                                factor_templates.insert(current_template);
+                                                        }
+                                                }                                       
+                                        }
+                                }       
+
+                                correction_factors.push_back(final_factors);
+
+                        }
+                }
 
 protected:
@@ -1234 +1278 @@
 
 
+
 #endif //TRAGE_H