Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
Reindent everything (possibly breaking all branches, but for the last time)
[simgrid.git] / testsuite / surf / lmm_usage.c
index 3b534d9..2339813 100644 (file)
@@ -17,8 +17,7 @@
 #include "xbt/module.h"
 #include <math.h>
 
-XBT_LOG_NEW_DEFAULT_CATEGORY(surf_test,
-                            "Messages specific for surf example");
+XBT_LOG_NEW_DEFAULT_CATEGORY(surf_test, "Messages specific for surf example");
 
 #define PRINT_VAR(var) DEBUG1(#var " = %g",lmm_variable_getvalue(var));
 #define SHOW_EXPR(expr) DEBUG1(#expr " = %g",expr);
@@ -37,7 +36,7 @@ typedef enum {
 } method_t;
 
 static double dichotomy(double func(double), double min, double max,
-                       double min_error)
+                        double min_error)
 {
   double middle;
   double min_func, max_func, middle_func;
@@ -59,7 +58,7 @@ static double dichotomy(double func(double), double min, double max,
   while (overall_error > min_error) {
     SHOW_EXPR(overall_error);
     if ((min_func > 0 && max_func > 0) ||
-       (min_func < 0 && max_func < 0) || (min_func > 0 && max_func < 0)) {
+        (min_func < 0 && max_func < 0) || (min_func > 0 && max_func < 0)) {
       abort();
     }
 
@@ -83,7 +82,7 @@ static double dichotomy(double func(double), double min, double max,
     } else if (middle_func > 0) {
       max = middle;
       max_func = middle_func;
-      overall_error = middle_func-min_func;
+      overall_error = middle_func - min_func;
     } else {
       overall_error = 0;
     }
@@ -96,10 +95,10 @@ double b_test_1 = 0;
 static double diff_lagrange_test_1(double x)
 {
   return -(3 / (1 + 3 * x * x / 2) -
-          3 / (2 * (3 * (a_test_1 - x) * (a_test_1 - x) / 2 + 1)) +
-          3 / (2 *
-               (3 * (b_test_1 - a_test_1 + x) *
-                (b_test_1 - a_test_1 + x) / 2 + 1)));
+           3 / (2 * (3 * (a_test_1 - x) * (a_test_1 - x) / 2 + 1)) +
+           3 / (2 *
+                (3 * (b_test_1 - a_test_1 + x) *
+                 (b_test_1 - a_test_1 + x) / 2 + 1)));
 }
 
 void test1(method_t method);
@@ -119,10 +118,10 @@ void test1(method_t method)
 
   if (method == LAGRANGE_VEGAS)
     lmm_set_default_protocol_function(func_vegas_f, func_vegas_fp,
-                                     func_vegas_fpi);
+                                      func_vegas_fpi);
   else if (method == LAGRANGE_RENO)
     lmm_set_default_protocol_function(func_reno_f, func_reno_fpi,
-                                     func_reno_fpi);
+                                      func_reno_fpi);
 
   Sys = lmm_system_new();
   L1 = lmm_constraint_new(Sys, (void *) "L1", a);
@@ -158,7 +157,7 @@ void test1(method_t method)
 #endif
   } else if (method == LAGRANGE_VEGAS) {
     double x = 3 * a / 4 - 3 * b / 8 +
-       sqrt(9 * b * b + 4 * a * a - 4 * a * b) / 8;
+      sqrt(9 * b * b + 4 * a * a - 4 * a * b) / 8;
     /* Computed with mupad and D_f=1.0 */
     double max_deviation = 0.0;
     if (x > a) {
@@ -170,26 +169,24 @@ void test1(method_t method)
 
     lagrange_solve(Sys);
 
+    max_deviation = MAX(max_deviation, fabs(lmm_variable_getvalue(R_1) - x));
+    max_deviation = MAX(max_deviation, fabs(lmm_variable_getvalue(R_3) - x));
     max_deviation =
-       MAX(max_deviation, fabs(lmm_variable_getvalue(R_1) - x));
+      MAX(max_deviation, fabs(lmm_variable_getvalue(R_2) - (b - a + x)));
     max_deviation =
-       MAX(max_deviation, fabs(lmm_variable_getvalue(R_3) - x));
-    max_deviation =
-       MAX(max_deviation, fabs(lmm_variable_getvalue(R_2) - (b - a + x)));
-    max_deviation =
-       MAX(max_deviation, fabs(lmm_variable_getvalue(R_1_2_3) - (a - x)));
+      MAX(max_deviation, fabs(lmm_variable_getvalue(R_1_2_3) - (a - x)));
 
     if (max_deviation > MAXMIN_PRECISION) {
       WARN1("Max Deviation from optimal solution : %g", max_deviation);
       WARN1("Found x = %1.20f", x);
       WARN2("Deviation from optimal solution (R_1 = %g): %1.20f", x,
-           lmm_variable_getvalue(R_1) - x);
+            lmm_variable_getvalue(R_1) - x);
       WARN2("Deviation from optimal solution (R_2 = %g): %1.20f",
-           b - a + x, lmm_variable_getvalue(R_2) - (b - a + x));
+            b - a + x, lmm_variable_getvalue(R_2) - (b - a + x));
       WARN2("Deviation from optimal solution (R_3 = %g): %1.20f", x,
-           lmm_variable_getvalue(R_3) - x);
+            lmm_variable_getvalue(R_3) - x);
       WARN2("Deviation from optimal solution (R_1_2_3 = %g): %1.20f",
-           a - x, lmm_variable_getvalue(R_1_2_3) - (a - x));
+            a - x, lmm_variable_getvalue(R_1_2_3) - (a - x));
     }
   } else if (method == LAGRANGE_RENO) {
     double x;
@@ -205,26 +202,24 @@ void test1(method_t method)
       x = a;
     lagrange_solve(Sys);
 
+    max_deviation = MAX(max_deviation, fabs(lmm_variable_getvalue(R_1) - x));
+    max_deviation = MAX(max_deviation, fabs(lmm_variable_getvalue(R_3) - x));
     max_deviation =
-       MAX(max_deviation, fabs(lmm_variable_getvalue(R_1) - x));
-    max_deviation =
-       MAX(max_deviation, fabs(lmm_variable_getvalue(R_3) - x));
-    max_deviation =
-       MAX(max_deviation, fabs(lmm_variable_getvalue(R_2) - (b - a + x)));
+      MAX(max_deviation, fabs(lmm_variable_getvalue(R_2) - (b - a + x)));
     max_deviation =
-       MAX(max_deviation, fabs(lmm_variable_getvalue(R_1_2_3) - (a - x)));
+      MAX(max_deviation, fabs(lmm_variable_getvalue(R_1_2_3) - (a - x)));
 
     if (max_deviation > MAXMIN_PRECISION) {
       WARN1("Max Deviation from optimal solution : %g", max_deviation);
       WARN1("Found x = %1.20f", x);
       WARN2("Deviation from optimal solution (R_1 = %g): %1.20f", x,
-           lmm_variable_getvalue(R_1) - x);
+            lmm_variable_getvalue(R_1) - x);
       WARN2("Deviation from optimal solution (R_2 = %g): %1.20f",
-           b - a + x, lmm_variable_getvalue(R_2) - (b - a + x));
+            b - a + x, lmm_variable_getvalue(R_2) - (b - a + x));
       WARN2("Deviation from optimal solution (R_3 = %g): %1.20f", x,
-           lmm_variable_getvalue(R_3) - x);
+            lmm_variable_getvalue(R_3) - x);
       WARN2("Deviation from optimal solution (R_1_2_3 = %g): %1.20f",
-           a - x, lmm_variable_getvalue(R_1_2_3) - (a - x));
+            a - x, lmm_variable_getvalue(R_1_2_3) - (a - x));
     }
   } else {
     xbt_assert0(0, "Invalid method");
@@ -251,10 +246,10 @@ void test2(method_t method)
 
   if (method == LAGRANGE_VEGAS)
     lmm_set_default_protocol_function(func_vegas_f, func_vegas_fp,
-                                     func_vegas_fpi);
+                                      func_vegas_fpi);
   else if (method == LAGRANGE_RENO)
     lmm_set_default_protocol_function(func_reno_f, func_reno_fp,
-                                     func_reno_fpi);
+                                      func_reno_fpi);
 
   Sys = lmm_system_new();
   CPU1 = lmm_constraint_new(Sys, (void *) "CPU1", 200.0);
@@ -316,9 +311,9 @@ void test3(method_t method)
 
     1, 1, 1, 1, 1
   };
-       
-  /*A = xbt_new0(double*, links + 5);*/
-  A = xbt_new0(double*, links + 5);
+
+  /*A = xbt_new0(double*, links + 5); */
+  A = xbt_new0(double *, links + 5);
 
   for (i = 0; i < links + 5; i++) {
     A[i] = xbt_new0(double, flows + 5);
@@ -326,7 +321,7 @@ void test3(method_t method)
       A[i][j] = 0.0;
 
       if (i >= links || j >= flows) {
-       A[i][j] = 0.0;
+        A[i][j] = 0.0;
       }
     }
   }
@@ -403,21 +398,21 @@ void test3(method_t method)
 
   if (method == LAGRANGE_VEGAS)
     lmm_set_default_protocol_function(func_vegas_f, func_vegas_fp,
-                                     func_vegas_fpi);
+                                      func_vegas_fpi);
   else if (method == LAGRANGE_RENO)
     lmm_set_default_protocol_function(func_reno_f, func_reno_fp,
-                                     func_reno_fpi);
+                                      func_reno_fpi);
 
   Sys = lmm_system_new();
 
 
 
-  tmp_name = xbt_new0(char *,31);
+  tmp_name = xbt_new0(char *, 31);
 
   /*
    * Creates the constraints
    */
-  tmp_cnst = xbt_new0(lmm_constraint_t,15);
+  tmp_cnst = xbt_new0(lmm_constraint_t, 15);
   for (i = 0; i < 15; i++) {
     tmp_name[i] = bprintf("C_%03d", i);
     tmp_cnst[i] = lmm_constraint_new(Sys, (void *) tmp_name[i], B[i]);
@@ -427,11 +422,11 @@ void test3(method_t method)
   /*
    * Creates the variables
    */
-  tmp_var = xbt_new0(lmm_variable_t,16);
+  tmp_var = xbt_new0(lmm_variable_t, 16);
   for (j = 0; j < 16; j++) {
     tmp_name[i + j] = bprintf("X_%03d", j);
     tmp_var[j] =
-       lmm_variable_new(Sys, (void *) tmp_name[i + j], 1.0, -1.0, 15);
+      lmm_variable_new(Sys, (void *) tmp_name[i + j], 1.0, -1.0, 15);
     lmm_update_variable_weight(Sys, tmp_var[j], 1.0);
   }
 
@@ -441,7 +436,7 @@ void test3(method_t method)
   for (i = 0; i < 15; i++) {
     for (j = 0; j < 16; j++) {
       if (A[i][j]) {
-       lmm_expand(Sys, tmp_cnst[i], tmp_var[j], 1.0);
+        lmm_expand(Sys, tmp_cnst[i], tmp_var[j], 1.0);
       }
     }
   }
@@ -491,8 +486,8 @@ int main(int argc, char **argv)
   INFO0("***** Test 1 (Lagrange - Reno)");
   test1(LAGRANGE_RENO);
 
-   
-   
+
+
   INFO0("***** Test 2 (Max-Min)");
   test2(MAXMIN);
 #ifdef HAVE_SDP