]> https://gitweb.dealii.org/ - dealii.git/commitdiff
fix tests 1836/head
authorTimo Heister <timo.heister@gmail.com>
Thu, 5 Nov 2015 14:02:50 +0000 (09:02 -0500)
committerTimo Heister <timo.heister@gmail.com>
Thu, 5 Nov 2015 14:03:38 +0000 (09:03 -0500)
- use std_cxx11::bind instead of std::bind
- require Trilinos if we use it
- fe/bdm_*: do not set matrix entries <1e-14 to zero (different on
different machines) and always print 0 entries
- use check_solver_within_range where needed

15 files changed:
tests/fe/bdm_8.cc
tests/fe/bdm_8.output
tests/fe/bdm_9.cc
tests/fe/bdm_9.output
tests/mpi/cell_weights_01_back_and_forth_01.cc
tests/mpi/cell_weights_01_back_and_forth_02.cc
tests/mpi/cell_weights_02.cc
tests/mpi/cell_weights_03.cc
tests/mpi/cell_weights_04.cc
tests/mpi/cell_weights_05.cc
tests/mpi/cell_weights_06.cc
tests/multigrid/transfer_04.with_trilinos=true.mpirun=1.output [moved from tests/multigrid/transfer_04.mpirun=1.output with 100% similarity]
tests/multigrid/transfer_04.with_trilinos=true.mpirun=2.output [moved from tests/multigrid/transfer_04.mpirun=2.output with 100% similarity]
tests/multigrid/transfer_04a.with_trilinos=true.mpirun=1.debug.output [moved from tests/multigrid/transfer_04a.mpirun=1.debug.output with 100% similarity]
tests/multigrid/transfer_04a.with_trilinos=true.mpirun=2.debug.output [moved from tests/multigrid/transfer_04a.mpirun=2.debug.output with 100% similarity]

index de8f98129dc5909551ed8f88799b30303567b0b4..7ba8dfabcf6112c044ac86c72e04b1a9f3c8c917 100644 (file)
@@ -73,26 +73,23 @@ test (const unsigned int degree)
           += (fe.shape_value_component(i,q_point,d) *
               fe.shape_value_component(j,q_point,d) *
               fe.JxW(q_point));
-  for (unsigned int i=0; i<dofs_per_cell; ++i)
-    for (unsigned int j=0; j<dofs_per_cell; ++j)
-      if (std::fabs(mass_matrix(i,j)) < 1e-14)
-        mass_matrix(i,j) = 0;
-  mass_matrix.print_formatted (logfile, 3, false, 0, " ", 1);
+
+  mass_matrix.print_formatted (logfile, 3, false, 0, "0", 1);
 
   SolverControl           solver_control (2*dofs_per_cell,
                                           1e-8);
   PrimitiveVectorMemory<> vector_memory;
-  SolverCG<>              cg (solver_control, vector_memory);
+  SolverCG<>              solver (solver_control, vector_memory);
 
   Vector<double> tmp1(dofs_per_cell), tmp2(dofs_per_cell);
   for (unsigned int i=0; i<dofs_per_cell; ++i)
     tmp1(i) = 1.*Testing::rand()/RAND_MAX;
-  cg.solve (mass_matrix, tmp2, tmp1, PreconditionIdentity());
 
-  deallog << "Degree=" << degree
-          << ": " << solver_control.last_step()
-          << " iterations to obtain convergence."
-          << std::endl;
+  deallog << "solving degree = " << degree << std::endl;
+  check_solver_within_range(
+    solver.solve (mass_matrix, tmp2, tmp1,
+                  PreconditionIdentity()),
+    solver_control.last_step(), 3, 45);
 }
 
 
index 0409547b66d6080d6cf259e472db05bed47f678d..42a90f21ab6b2ca0843c106ffd817074349a550f 100644 (file)
 
-0.208 -0.036 0.058 0.020 0.068             0.068 
--0.036 0.208 0.020 0.058       -0.068 -0.068       
-0.058 0.020 0.208 -0.036       -0.068 -0.068       
-0.020 0.058 -0.036 0.208 0.068             0.068 
-0.068             0.068 0.208 -0.036 0.058 0.020 
-      -0.068 -0.068       -0.036 0.208 0.020 0.058 
-      -0.068 -0.068       0.058 0.020 0.208 -0.036 
-0.068             0.068 0.020 0.058 -0.036 0.208 
-DEAL:cg::Starting value 1.91960
-DEAL:cg::Convergence step 4 value 0
-DEAL::Degree=1: 4 iterations to obtain convergence.
+0.208 -0.036 0.058 0.020 0.068 0.000 0.000 0.068 
+-0.036 0.208 0.020 0.058 0.000 -0.068 -0.068 0.000 
+0.058 0.020 0.208 -0.036 0.000 -0.068 -0.068 0.000 
+0.020 0.058 -0.036 0.208 0.068 0.000 0.000 0.068 
+0.068 0.000 0.000 0.068 0.208 -0.036 0.058 0.020 
+0.000 -0.068 -0.068 0.000 -0.036 0.208 0.020 0.058 
+0.000 -0.068 -0.068 0.000 0.058 0.020 0.208 -0.036 
+0.068 0.000 0.000 0.068 0.020 0.058 -0.036 0.208 
+DEAL::solving degree = 1
+DEAL::Solver stopped within 3 - 45 iterations
 3.262 -2.918 -2.918 0.613 -0.457 -0.457 0.768 -0.598 -0.469 -0.299 0.598 0.469 0.768 -0.469 -0.598 -0.299 0.469 0.598 
--2.918 5.836       -0.457 0.914       -0.469             -0.469             -0.598       1.195 0.598       -1.195 
--2.918       5.836 -0.457       0.914 -0.598 1.195       0.598 -1.195       -0.469             -0.469             
+-2.918 5.836 0.000 -0.457 0.914 0.000 -0.469 0.000 0.000 -0.469 0.000 0.000 -0.598 0.000 1.195 0.598 0.000 -1.195 
+-2.918 0.000 5.836 -0.457 0.000 0.914 -0.598 1.195 0.000 0.598 -1.195 0.000 -0.469 0.000 0.000 -0.469 0.000 0.000 
 0.613 -0.457 -0.457 3.262 -2.918 -2.918 -0.299 0.598 -0.469 -0.170 -0.598 0.469 -0.299 -0.469 0.598 -0.170 0.469 -0.598 
--0.457 0.914       -2.918 5.836       0.469             0.469             0.598       -1.195 -0.598       1.195 
--0.457       0.914 -2.918       5.836 0.598 -1.195       -0.598 1.195       0.469             0.469             
+-0.457 0.914 0.000 -2.918 5.836 0.000 0.469 0.000 0.000 0.469 0.000 0.000 0.598 0.000 -1.195 -0.598 0.000 1.195 
+-0.457 0.000 0.914 -2.918 0.000 5.836 0.598 -1.195 0.000 -0.598 1.195 0.000 0.469 0.000 0.000 0.469 0.000 0.000 
 0.768 -0.469 -0.598 -0.299 0.469 0.598 3.262 -2.918 -2.918 0.613 -0.457 -0.457 0.768 -0.598 -0.469 -0.299 0.598 0.469 
--0.598       1.195 0.598       -1.195 -2.918 5.836       -0.457 0.914       -0.469             -0.469             
--0.469             -0.469             -2.918       5.836 -0.457       0.914 -0.598 1.195       0.598 -1.195       
+-0.598 0.000 1.195 0.598 0.000 -1.195 -2.918 5.836 0.000 -0.457 0.914 0.000 -0.469 0.000 0.000 -0.469 0     0.000 
+-0.469 0.000 0.000 -0.469 0.000 0.000 -2.918 0.000 5.836 -0.457 0.000 0.914 -0.598 1.195 0.000 0.598 -1.195 0.000 
 -0.299 -0.469 0.598 -0.170 0.469 -0.598 0.613 -0.457 -0.457 3.262 -2.918 -2.918 -0.299 0.598 -0.469 -0.170 -0.598 0.469 
-0.598       -1.195 -0.598       1.195 -0.457 0.914       -2.918 5.836       0.469             0.469             
-0.469             0.469             -0.457       0.914 -2.918       5.836 0.598 -1.195       -0.598 1.195       
+0.598 0.000 -1.195 -0.598 0.000 1.195 -0.457 0.914 0.000 -2.918 5.836 0     0.469 0.000 0.000 0.469 0     0.000 
+0.469 0.000 0.000 0.469 0.000 0.000 -0.457 0.000 0.914 -2.918 0     5.836 0.598 -1.195 0     -0.598 1.195 0     
 0.768 -0.598 -0.469 -0.299 0.598 0.469 0.768 -0.469 -0.598 -0.299 0.469 0.598 3.262 -2.918 -2.918 0.613 -0.457 -0.457 
--0.469             -0.469             -0.598       1.195 0.598       -1.195 -2.918 5.836       -0.457 0.914       
--0.598 1.195       0.598 -1.195       -0.469             -0.469             -2.918       5.836 -0.457       0.914 
+-0.469 0.000 0.000 -0.469 0.000 0.000 -0.598 0.000 1.195 0.598 0.000 -1.195 -2.918 5.836 0.000 -0.457 0.914 0.000 
+-0.598 1.195 0.000 0.598 -1.195 0.000 -0.469 0.000 0.000 -0.469 0.000 0     -2.918 0.000 5.836 -0.457 0.000 0.914 
 -0.299 0.598 -0.469 -0.170 -0.598 0.469 -0.299 -0.469 0.598 -0.170 0.469 -0.598 0.613 -0.457 -0.457 3.262 -2.918 -2.918 
-0.469             0.469             0.598       -1.195 -0.598       1.195 -0.457 0.914       -2.918 5.836       
-0.598 -1.195       -0.598 1.195       0.469             0.469             -0.457       0.914 -2.918       5.836 
-DEAL:cg::Starting value 2.31677
-DEAL:cg::Convergence step 11 value 0
-DEAL::Degree=1: 11 iterations to obtain convergence.
-0.105 -0.050 -0.003 0.036 -0.030 -0.013 0.004 -0.002 0.004 -0.004 0.002 -0.004 -0.020       
--0.050 0.169 -0.050 -0.030 0.058 -0.030 -0.002 -0.010 -0.002 0.002 0.010 0.002 -0.090       
--0.003 -0.050 0.105 -0.013 -0.030 0.036 0.004 -0.002 0.004 -0.004 0.002 -0.004 -0.020       
-0.036 -0.030 -0.013 0.105 -0.050 -0.003 -0.004 0.002 -0.004 0.004 -0.002 0.004 -0.020       
--0.030 0.058 -0.030 -0.050 0.169 -0.050 0.002 0.010 0.002 -0.002 -0.010 -0.002 -0.090       
--0.013 -0.030 0.036 -0.003 -0.050 0.105 -0.004 0.002 -0.004 0.004 -0.002 0.004 -0.020       
-0.004 -0.002 0.004 -0.004 0.002 -0.004 0.105 -0.050 -0.003 0.036 -0.030 -0.013       -0.020 
--0.002 -0.010 -0.002 0.002 0.010 0.002 -0.050 0.169 -0.050 -0.030 0.058 -0.030       -0.090 
-0.004 -0.002 0.004 -0.004 0.002 -0.004 -0.003 -0.050 0.105 -0.013 -0.030 0.036       -0.020 
--0.004 0.002 -0.004 0.004 -0.002 0.004 0.036 -0.030 -0.013 0.105 -0.050 -0.003       -0.020 
-0.002 0.010 0.002 -0.002 -0.010 -0.002 -0.030 0.058 -0.030 -0.050 0.169 -0.050       -0.090 
--0.004 0.002 -0.004 0.004 -0.002 0.004 -0.013 -0.030 0.036 -0.003 -0.050 0.105       -0.020 
--0.020 -0.090 -0.020 -0.020 -0.090 -0.020                                     1.195       
-                                    -0.020 -0.090 -0.020 -0.020 -0.090 -0.020       1.195 
-DEAL:cg::Starting value 2.24859
-DEAL:cg::Convergence step 9 value 0.00000
-DEAL::Degree=2: 9 iterations to obtain convergence.
-13.347 -29.750 -29.750 19.696 17.582 17.582 4.028 -9.390 -9.390 3.085 6.793 6.793 1.091 -2.351 -0.681 0.000 0.769 0.681 0.491 -0.813 0.681 0.000 -0.769 -0.681 1.091 -2.263 -2.351 3.164 0.681 0.769 -1.091 2.263 2.351 -3.164 -0.681 -0.769 0.164 0.000       
--29.750 127.389 20.905 -39.393 -103.052 -1.208 -9.390 46.283 3.634 -6.170 -41.089 -0.549 -2.263 3.823 1.318 0.000 -0.659 -1.318 -0.901 2.505 -1.318 0.000 0.659 1.318 -2.351 3.823 7.119 -6.328 -0.659 -3.955 2.351 -3.823 -7.119 6.328 0.659 3.955 -0.879       0.000 
--29.750 20.905 127.389 -39.393 -1.208 -103.052 -9.390 3.634 46.283 -6.170 -0.549 -41.089 -2.351 7.119 0.659 0.000 -3.955 -0.659 -0.813 -0.791 -0.659 0.000 3.955 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -0.879 0.000       
-19.696 -39.393 -39.393 78.785 0.000 0.000 3.085 -6.170 -6.170 12.340 0.000 0.000 3.164 -6.328 0.000 0.000 0.000 0.000 3.164 -6.328 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000       0.000       
+0.469 0.000 0.000 0.469 0.000 0.000 0.598 0     -1.195 -0.598 0     1.195 -0.457 0.914 0.000 -2.918 5.836 0     
+0.598 -1.195 0.000 -0.598 1.195 0.000 0.469 0.000 0.000 0.469 0.000 0     -0.457 0.000 0.914 -2.918 0     5.836 
+DEAL::solving degree = 1
+DEAL::Solver stopped within 3 - 45 iterations
+0.105 -0.050 -0.003 0.036 -0.030 -0.013 0.004 -0.002 0.004 -0.004 0.002 -0.004 -0.020 0.000 
+-0.050 0.169 -0.050 -0.030 0.058 -0.030 -0.002 -0.010 -0.002 0.002 0.010 0.002 -0.090 0.000 
+-0.003 -0.050 0.105 -0.013 -0.030 0.036 0.004 -0.002 0.004 -0.004 0.002 -0.004 -0.020 0.000 
+0.036 -0.030 -0.013 0.105 -0.050 -0.003 -0.004 0.002 -0.004 0.004 -0.002 0.004 -0.020 0.000 
+-0.030 0.058 -0.030 -0.050 0.169 -0.050 0.002 0.010 0.002 -0.002 -0.010 -0.002 -0.090 0.000 
+-0.013 -0.030 0.036 -0.003 -0.050 0.105 -0.004 0.002 -0.004 0.004 -0.002 0.004 -0.020 0.000 
+0.004 -0.002 0.004 -0.004 0.002 -0.004 0.105 -0.050 -0.003 0.036 -0.030 -0.013 0.000 -0.020 
+-0.002 -0.010 -0.002 0.002 0.010 0.002 -0.050 0.169 -0.050 -0.030 0.058 -0.030 0.000 -0.090 
+0.004 -0.002 0.004 -0.004 0.002 -0.004 -0.003 -0.050 0.105 -0.013 -0.030 0.036 0.000 -0.020 
+-0.004 0.002 -0.004 0.004 -0.002 0.004 0.036 -0.030 -0.013 0.105 -0.050 -0.003 0.000 -0.020 
+0.002 0.010 0.002 -0.002 -0.010 -0.002 -0.030 0.058 -0.030 -0.050 0.169 -0.050 0.000 -0.090 
+-0.004 0.002 -0.004 0.004 -0.002 0.004 -0.013 -0.030 0.036 -0.003 -0.050 0.105 0.000 -0.020 
+-0.020 -0.090 -0.020 -0.020 -0.090 -0.020 0.000 0.000 0.000 0.000 0.000 0.000 1.195 0     
+0.000 0.000 0.000 0.000 0.000 0.000 -0.020 -0.090 -0.020 -0.020 -0.090 -0.020 0     1.195 
+DEAL::solving degree = 2
+DEAL::Solver stopped within 3 - 45 iterations
+13.347 -29.750 -29.750 19.696 17.582 17.582 4.028 -9.390 -9.390 3.085 6.793 6.793 1.091 -2.351 -0.681 0.000 0.769 0.681 0.491 -0.813 0.681 0.000 -0.769 -0.681 1.091 -2.263 -2.351 3.164 0.681 0.769 -1.091 2.263 2.351 -3.164 -0.681 -0.769 0.164 0.000 0.000 
+-29.750 127.389 20.905 -39.393 -103.052 -1.208 -9.390 46.283 3.634 -6.170 -41.089 -0.549 -2.263 3.823 1.318 0.000 -0.659 -1.318 -0.901 2.505 -1.318 0.000 0.659 1.318 -2.351 3.823 7.119 -6.328 -0.659 -3.955 2.351 -3.823 -7.119 6.328 0.659 3.955 -0.879 0.000 0.000 
+-29.750 20.905 127.389 -39.393 -1.208 -103.052 -9.390 3.634 46.283 -6.170 -0.549 -41.089 -2.351 7.119 0.659 0.000 -3.955 -0.659 -0.813 -0.791 -0.659 0.000 3.955 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -0.879 0.000 0.000 
+19.696 -39.393 -39.393 78.785 0.000 0.000 3.085 -6.170 -6.170 12.340 0.000 0.000 3.164 -6.328 0.000 0.000 0.000 0.000 3.164 -6.328 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
 17.582 -103.052 -1.208 0.000 103.052 1.208 6.793 -41.089 -0.549 0.000 41.089 0.549 0.681 -0.659 -1.318 0.000 0.659 1.318 -0.681 0.659 1.318 0.000 -0.659 -1.318 0.769 -0.659 -3.955 0.000 0.659 3.955 -0.769 0.659 3.955 0.000 -0.659 -3.955 0.879 0.000 0.000 
-17.582 -1.208 -103.052 0.000 1.208 103.052 6.793 -0.549 -41.089 0.000 0.549 41.089 0.769 -3.955 -0.659 0.000 3.955 0.659 -0.769 3.955 0.659 0.000 -3.955 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 0.879 0.000       
-4.028 -9.390 -9.390 3.085 6.793 6.793 13.347 -29.750 -29.750 19.696 17.582 17.582 -1.091 2.351 0.681 0.000 -0.769 -0.681 -0.491 0.813 -0.681 0.000 0.769 0.681 0.491 -0.901 -0.813 3.164 -0.681 -0.769 -0.491 0.901 0.813 -3.164 0.681 0.769 0.164 0.000       
--9.390 46.283 3.634 -6.170 -41.089 -0.549 -29.750 127.389 20.905 -39.393 -103.052 -1.208 2.263 -3.823 -1.318 0.000 0.659 1.318 0.901 -2.505 1.318 0.000 -0.659 -1.318 -0.813 2.505 -0.791 -6.328 0.659 3.955 0.813 -2.505 0.791 6.328 -0.659 -3.955 -0.879       0.000 
--9.390 3.634 46.283 -6.170 -0.549 -41.089 -29.750 20.905 127.389 -39.393 -1.208 -103.052 2.351 -7.119 -0.659 0.000 3.955 0.659 0.813 0.791 0.659 0.000 -3.955 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 -0.879 0.000       
-3.085 -6.170 -6.170 12.340 0.000 0.000 19.696 -39.393 -39.393 78.785 0.000 0.000 -3.164 6.328 0.000 0.000 0.000 0.000 -3.164 6.328 0.000 0.000 0.000       0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000       
+17.582 -1.208 -103.052 0.000 1.208 103.052 6.793 -0.549 -41.089 0.000 0.549 41.089 0.769 -3.955 -0.659 0.000 3.955 0.659 -0.769 3.955 0.659 0.000 -3.955 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 0.879 0.000 0.000 
+4.028 -9.390 -9.390 3.085 6.793 6.793 13.347 -29.750 -29.750 19.696 17.582 17.582 -1.091 2.351 0.681 0.000 -0.769 -0.681 -0.491 0.813 -0.681 0.000 0.769 0.681 0.491 -0.901 -0.813 3.164 -0.681 -0.769 -0.491 0.901 0.813 -3.164 0.681 0.769 0.164 0.000 0.000 
+-9.390 46.283 3.634 -6.170 -41.089 -0.549 -29.750 127.389 20.905 -39.393 -103.052 -1.208 2.263 -3.823 -1.318 0.000 0.659 1.318 0.901 -2.505 1.318 0.000 -0.659 -1.318 -0.813 2.505 -0.791 -6.328 0.659 3.955 0.813 -2.505 0.791 6.328 -0.659 -3.955 -0.879 0.000 0.000 
+-9.390 3.634 46.283 -6.170 -0.549 -41.089 -29.750 20.905 127.389 -39.393 -1.208 -103.052 2.351 -7.119 -0.659 0.000 3.955 0.659 0.813 0.791 0.659 0.000 -3.955 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 -0.879 0.000 0.000 
+3.085 -6.170 -6.170 12.340 0.000 0.000 19.696 -39.393 -39.393 78.785 0.000 0.000 -3.164 6.328 0.000 0.000 0.000 0.000 -3.164 6.328 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
 6.793 -41.089 -0.549 0.000 41.089 0.549 17.582 -103.052 -1.208 0.000 103.052 1.208 -0.681 0.659 1.318 0.000 -0.659 -1.318 0.681 -0.659 -1.318 0.000 0.659 1.318 -0.769 0.659 3.955 0.000 -0.659 -3.955 0.769 -0.659 -3.955 0.000 0.659 3.955 0.879 0.000 0.000 
-6.793 -0.549 -41.089 0.000 0.549 41.089 17.582 -1.208 -103.052 0.000 1.208 103.052 -0.769 3.955 0.659 0.000 -3.955 -0.659 0.769 -3.955 -0.659 0.000 3.955 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 0.879 0.000       
-1.091 -2.263 -2.351 3.164 0.681 0.769 -1.091 2.263 2.351 -3.164 -0.681 -0.769 13.347 -29.750 -29.750 19.696 17.582 17.582 4.028 -9.390 -9.390 3.085 6.793 6.793 1.091 -2.351 -0.681 0.000 0.769 0.681 0.491 -0.813 0.681 0.000 -0.769 -0.681       0.164       
--2.351 3.823 7.119 -6.328 -0.659 -3.955 2.351 -3.823 -7.119 6.328 0.659 3.955 -29.750 127.389 20.905 -39.393 -103.052 -1.208 -9.390 46.283 3.634 -6.170 -41.089 -0.549 -2.263 3.823 1.318 0.000 -0.659 -1.318 -0.901 2.505 -1.318 0.000 0.659 1.318       -0.879 0.000 
--0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -29.750 20.905 127.389 -39.393 -1.208 -103.052 -9.390 3.634 46.283 -6.170 -0.549 -41.089 -2.351 7.119 0.659 0.000 -3.955 -0.659 -0.813 -0.791 -0.659 0.000 3.955 0.659       -0.879 0.000 
-0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 19.696 -39.393 -39.393 78.785 0.000 0.000 3.085 -6.170 -6.170 12.340 0.000 0.000 3.164 -6.328 0.000 0.000 0.000 0.000 3.164 -6.328 0.000 0.000       0.000                   
-0.769 -0.659 -3.955 0.000 0.659 3.955 -0.769 0.659 3.955 0.000 -0.659 -3.955 17.582 -103.052 -1.208 0.000 103.052 1.208 6.793 -41.089 -0.549 0.000 41.089 0.549 0.681 -0.659 -1.318 0.000 0.659 1.318 -0.681 0.659 1.318 0.000 -0.659 -1.318       0.879       
-0.681 -1.318 -0.659 0.000 1.318 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 17.582 -1.208 -103.052 0.000 1.208 103.052 6.793 -0.549 -41.089 0.000 0.549 41.089 0.769 -3.955 -0.659 0.000 3.955 0.659 -0.769 3.955 0.659 0.000 -3.955 -0.659       0.879 0.000 
-0.491 -0.901 -0.813 3.164 -0.681 -0.769 -0.491 0.901 0.813 -3.164 0.681 0.769 4.028 -9.390 -9.390 3.085 6.793 6.793 13.347 -29.750 -29.750 19.696 17.582 17.582 -1.091 2.351 0.681 0.000 -0.769 -0.681 -0.491 0.813 -0.681 0.000 0.769 0.681       0.164       
--0.813 2.505 -0.791 -6.328 0.659 3.955 0.813 -2.505 0.791 6.328 -0.659 -3.955 -9.390 46.283 3.634 -6.170 -41.089 -0.549 -29.750 127.389 20.905 -39.393 -103.052 -1.208 2.263 -3.823 -1.318 0.000 0.659 1.318 0.901 -2.505 1.318 0.000 -0.659 -1.318       -0.879       
-0.681 -1.318 -0.659 0.000 1.318 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 -9.390 3.634 46.283 -6.170 -0.549 -41.089 -29.750 20.905 127.389 -39.393 -1.208 -103.052 2.351 -7.119 -0.659 0.000 3.955 0.659 0.813 0.791 0.659 0.000 -3.955 -0.659       -0.879 0.000 
-0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 3.085 -6.170 -6.170 12.340 0.000 0.000 19.696 -39.393 -39.393 78.785 0.000 0.000 -3.164 6.328 0.000 0.000 0.000 0.000 -3.164 6.328 0.000 0.000 0.000 0.000                   
--0.769 0.659 3.955 0.000 -0.659 -3.955 0.769 -0.659 -3.955 0.000 0.659 3.955 6.793 -41.089 -0.549 0.000 41.089 0.549 17.582 -103.052 -1.208 0.000 103.052 1.208 -0.681 0.659 1.318 0.000 -0.659 -1.318 0.681 -0.659 -1.318 0.000 0.659 1.318       0.879       
--0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659       1.318 0.659 6.793 -0.549 -41.089 0.000 0.549 41.089 17.582 -1.208 -103.052 0.000 1.208 103.052 -0.769 3.955 0.659 0.000 -3.955 -0.659 0.769 -3.955 -0.659 0.000 3.955 0.659       0.879 0.000 
-1.091 -2.351 -0.681 0.000 0.769 0.681 0.491 -0.813 0.681 0.000 -0.769 -0.681 1.091 -2.263 -2.351 3.164 0.681 0.769 -1.091 2.263 2.351 -3.164 -0.681 -0.769 13.347 -29.750 -29.750 19.696 17.582 17.582 4.028 -9.390 -9.390 3.085 6.793 6.793             0.164 
--2.263 3.823 1.318 0.000 -0.659 -1.318 -0.901 2.505 -1.318 0.000 0.659 1.318 -2.351 3.823 7.119 -6.328 -0.659 -3.955 2.351 -3.823 -7.119 6.328 0.659 3.955 -29.750 127.389 20.905 -39.393 -103.052 -1.208 -9.390 46.283 3.634 -6.170 -41.089 -0.549 0.000       -0.879 
--2.351 7.119 0.659 0.000 -3.955 -0.659 -0.813 -0.791 -0.659 0.000 3.955 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -29.750 20.905 127.389 -39.393 -1.208 -103.052 -9.390 3.634 46.283 -6.170 -0.549 -41.089             -0.879 
-3.164 -6.328 0.000 0.000 0.000 0.000 3.164 -6.328 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 19.696 -39.393 -39.393 78.785 0.000 0.000 3.085 -6.170 -6.170 12.340 0.000 0.000                   
-0.681 -0.659 -1.318 0.000 0.659 1.318 -0.681 0.659 1.318 0.000 -0.659 -1.318 0.769 -0.659 -3.955 0.000 0.659 3.955 -0.769 0.659 3.955 0.000 -0.659 -3.955 17.582 -103.052 -1.208 0.000 103.052 1.208 6.793 -41.089 -0.549 0.000 41.089 0.549             0.879 
-0.769 -3.955 -0.659 0.000 3.955 0.659 -0.769 3.955 0.659 0.000 -3.955 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 17.582 -1.208 -103.052 0.000 1.208 103.052 6.793 -0.549 -41.089 0.000 0.549 41.089             0.879 
--1.091 2.351 0.681 0.000 -0.769 -0.681 -0.491 0.813 -0.681 0.000 0.769 0.681 0.491 -0.901 -0.813 3.164 -0.681 -0.769 -0.491 0.901 0.813 -3.164 0.681 0.769 4.028 -9.390 -9.390 3.085 6.793 6.793 13.347 -29.750 -29.750 19.696 17.582 17.582             0.164 
-2.263 -3.823 -1.318 0.000 0.659 1.318 0.901 -2.505 1.318 0.000 -0.659 -1.318 -0.813 2.505 -0.791 -6.328 0.659 3.955 0.813 -2.505 0.791 6.328 -0.659 -3.955 -9.390 46.283 3.634 -6.170 -41.089 -0.549 -29.750 127.389 20.905 -39.393 -103.052 -1.208             -0.879 
-2.351 -7.119 -0.659 0.000 3.955 0.659 0.813 0.791 0.659 0.000 -3.955 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 -9.390 3.634 46.283 -6.170 -0.549 -41.089 -29.750 20.905 127.389 -39.393 -1.208 -103.052             -0.879 
--3.164 6.328 0.000 0.000 0.000 0.000 -3.164 6.328 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 3.085 -6.170 -6.170 12.340 0.000 0.000 19.696 -39.393 -39.393 78.785 0.000 0.000                   
--0.681 0.659 1.318 0.000 -0.659 -1.318 0.681 -0.659 -1.318 0.000 0.659 1.318 -0.769 0.659 3.955       -0.659 -3.955 0.769 -0.659 -3.955 0.000 0.659 3.955 6.793 -41.089 -0.549 0.000 41.089 0.549 17.582 -103.052 -1.208 0.000 103.052 1.208             0.879 
--0.769 3.955 0.659 0.000 -3.955 -0.659 0.769 -3.955 -0.659 0.000 3.955 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 6.793 -0.549 -41.089 0.000 0.549 41.089 17.582 -1.208 -103.052 0.000 1.208 103.052             0.879 
-0.164 -0.879 -0.879       0.879 0.879 0.164 -0.879 -0.879 0.000 0.879 0.879                                                                               0.000                                                             1.195             
-0.000       0.000 0.000 0.000 0.000 0.000       0.000 0.000 0.000 0.000 0.164 -0.879 -0.879       0.879 0.879 0.164 -0.879 -0.879       0.879 0.879                                                                               1.195       
-      0.000             0.000             0.000             0.000             0.000 0.000             0.000             0.000             0.000 0.164 -0.879 -0.879       0.879 0.879 0.164 -0.879 -0.879       0.879 0.879             1.195 
-DEAL:cg::Starting value 3.79719
-DEAL:cg::Convergence step 40 value 0.00000
-DEAL::Degree=2: 40 iterations to obtain convergence.
-0.072 -0.058 0.015 0.002 0.016 -0.031 0.001 0.011 0.000 0.000 0.001             0.001 0.000 0.000 -0.040 0.084 -0.016 -0.006             
--0.058 0.151 -0.056 0.015 -0.031 0.067 -0.039 0.001 0.000 -0.002       -0.001 -0.001       -0.002 0.000 -0.404 0.224 0.468 0.016             
-0.015 -0.056 0.151 -0.058 0.001 -0.039 0.067 -0.031 0.001       0.002 0.000 0.000 0.002       0.001 0.064 0.224 -0.468 -0.016             
-0.002 0.015 -0.058 0.072 0.011 0.001 -0.031 0.016       -0.001 0.000 0.000 0.000 0.000 -0.001       -0.057 0.084 0.016 0.006             
-0.016 -0.031 0.001 0.011 0.072 -0.058 0.015 0.002       -0.001 0.000 0.000 0.000 0.000 -0.001       0.043 -0.084 -0.016 0.006             
--0.031 0.067 -0.039 0.001 -0.058 0.151 -0.056 0.015 0.001       0.002 0.000 0.000 0.002       0.001 -0.180 -0.224 0.468 -0.016             
-0.001 -0.039 0.067 -0.031 0.015 -0.056 0.151 -0.058 0.000 -0.002       -0.001 -0.001       -0.002 0.000 0.288 -0.224 -0.468 0.016             
-0.011 0.001 -0.031 0.016 0.002 0.015 -0.058 0.072 0.000 0.000 0.001             0.001 0.000 0.000 0.027 -0.084 0.016 -0.006             
-0.000 0.000 0.001             0.001 0.000 0.000 0.072 -0.058 0.015 0.002 0.016 -0.031 0.001 0.011 -0.006             -0.040 -0.016 0.084 
-0.000 -0.002       -0.001 -0.001       -0.002 0.000 -0.058 0.151 -0.056 0.015 -0.031 0.067 -0.039 0.001 0.016             -0.404 0.468 0.224 
-0.001       0.002 0.000 0.000 0.002       0.001 0.015 -0.056 0.151 -0.058 0.001 -0.039 0.067 -0.031 -0.016             0.064 -0.468 0.224 
-      -0.001 0.000 0.000 0.000 0.000 -0.001       0.002 0.015 -0.058 0.072 0.011 0.001 -0.031 0.016 0.006             -0.057 0.016 0.084 
-      -0.001 0.000 0.000 0.000 0.000 -0.001       0.016 -0.031 0.001 0.011 0.072 -0.058 0.015 0.002 0.006             0.043 -0.016 -0.084 
-0.001       0.002 0.000 0.000 0.002       0.001 -0.031 0.067 -0.039 0.001 -0.058 0.151 -0.056 0.015 -0.016             -0.180 0.468 -0.224 
-0.000 -0.002       -0.001 -0.001       -0.002 0.000 0.001 -0.039 0.067 -0.031 0.015 -0.056 0.151 -0.058 0.016             0.288 -0.468 -0.224 
-0.000 0.000 0.001             0.001 0.000 0.000 0.011 0.001 -0.031 0.016 0.002 0.015 -0.058 0.072 -0.006             0.027 0.016 -0.084 
--0.040 -0.404 0.064 -0.057 0.043 -0.180 0.288 0.027 -0.006 0.016 -0.016 0.006 0.006 -0.016 0.016 -0.006 9.185 -7.910 -8.068                   
-0.084 0.224 0.224 0.084 -0.084 -0.224 -0.224 -0.084                                                 -7.910 15.820                         
--0.016 0.468 -0.468 0.016 -0.016 0.468 -0.468 0.016                                                 -8.068       16.137                   
--0.006 0.016 -0.016 0.006 0.006 -0.016 0.016 -0.006 -0.040 -0.404 0.064 -0.057 0.043 -0.180 0.288 0.027                   9.185 -8.068 -7.910 
-                                                -0.016 0.468 -0.468 0.016 -0.016 0.468 -0.468 0.016                   -8.068 16.137       
-                                                0.084 0.224 0.224 0.084 -0.084 -0.224 -0.224 -0.084                   -7.910       15.820 
-DEAL:cg::Starting value 3.10915
-DEAL:cg::Convergence step 14 value 0.00000
-DEAL::Degree=3: 14 iterations to obtain convergence.
+6.793 -0.549 -41.089 0.000 0.549 41.089 17.582 -1.208 -103.052 0.000 1.208 103.052 -0.769 3.955 0.659 0.000 -3.955 -0.659 0.769 -3.955 -0.659 0.000 3.955 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 0.879 0.000 0.000 
+1.091 -2.263 -2.351 3.164 0.681 0.769 -1.091 2.263 2.351 -3.164 -0.681 -0.769 13.347 -29.750 -29.750 19.696 17.582 17.582 4.028 -9.390 -9.390 3.085 6.793 6.793 1.091 -2.351 -0.681 0.000 0.769 0.681 0.491 -0.813 0.681 0.000 -0.769 -0.681 0.000 0.164 0.000 
+-2.351 3.823 7.119 -6.328 -0.659 -3.955 2.351 -3.823 -7.119 6.328 0.659 3.955 -29.750 127.389 20.905 -39.393 -103.052 -1.208 -9.390 46.283 3.634 -6.170 -41.089 -0.549 -2.263 3.823 1.318 0.000 -0.659 -1.318 -0.901 2.505 -1.318 0.000 0.659 1.318 0.000 -0.879 0.000 
+-0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -29.750 20.905 127.389 -39.393 -1.208 -103.052 -9.390 3.634 46.283 -6.170 -0.549 -41.089 -2.351 7.119 0.659 0.000 -3.955 -0.659 -0.813 -0.791 -0.659 0.000 3.955 0.659 0.000 -0.879 0.000 
+0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 19.696 -39.393 -39.393 78.785 0.000 0.000 3.085 -6.170 -6.170 12.340 0.000 0.000 3.164 -6.328 0.000 0.000 0.000 0.000 3.164 -6.328 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
+0.769 -0.659 -3.955 0.000 0.659 3.955 -0.769 0.659 3.955 0.000 -0.659 -3.955 17.582 -103.052 -1.208 0.000 103.052 1.208 6.793 -41.089 -0.549 0.000 41.089 0.549 0.681 -0.659 -1.318 0.000 0.659 1.318 -0.681 0.659 1.318 0.000 -0.659 -1.318 0.000 0.879 0.000 
+0.681 -1.318 -0.659 0.000 1.318 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 17.582 -1.208 -103.052 0.000 1.208 103.052 6.793 -0.549 -41.089 0.000 0.549 41.089 0.769 -3.955 -0.659 0.000 3.955 0.659 -0.769 3.955 0.659 0.000 -3.955 -0.659 0.000 0.879 0.000 
+0.491 -0.901 -0.813 3.164 -0.681 -0.769 -0.491 0.901 0.813 -3.164 0.681 0.769 4.028 -9.390 -9.390 3.085 6.793 6.793 13.347 -29.750 -29.750 19.696 17.582 17.582 -1.091 2.351 0.681 0.000 -0.769 -0.681 -0.491 0.813 -0.681 0.000 0.769 0.681 0.000 0.164 0.000 
+-0.813 2.505 -0.791 -6.328 0.659 3.955 0.813 -2.505 0.791 6.328 -0.659 -3.955 -9.390 46.283 3.634 -6.170 -41.089 -0.549 -29.750 127.389 20.905 -39.393 -103.052 -1.208 2.263 -3.823 -1.318 0.000 0.659 1.318 0.901 -2.505 1.318 0.000 -0.659 -1.318 0.000 -0.879 0.000 
+0.681 -1.318 -0.659 0.000 1.318 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 -9.390 3.634 46.283 -6.170 -0.549 -41.089 -29.750 20.905 127.389 -39.393 -1.208 -103.052 2.351 -7.119 -0.659 0.000 3.955 0.659 0.813 0.791 0.659 0.000 -3.955 -0.659 0.000 -0.879 0.000 
+0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 3.085 -6.170 -6.170 12.340 0.000 0.000 19.696 -39.393 -39.393 78.785 0.000 0.000 -3.164 6.328 0.000 0.000 0.000 0.000 -3.164 6.328 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
+-0.769 0.659 3.955 0.000 -0.659 -3.955 0.769 -0.659 -3.955 0.000 0.659 3.955 6.793 -41.089 -0.549 0.000 41.089 0.549 17.582 -103.052 -1.208 0.000 103.052 1.208 -0.681 0.659 1.318 0.000 -0.659 -1.318 0.681 -0.659 -1.318 0.000 0.659 1.318 0.000 0.879 0.000 
+-0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 6.793 -0.549 -41.089 0.000 0.549 41.089 17.582 -1.208 -103.052 0.000 1.208 103.052 -0.769 3.955 0.659 0.000 -3.955 -0.659 0.769 -3.955 -0.659 0.000 3.955 0.659 0.000 0.879 0.000 
+1.091 -2.351 -0.681 0.000 0.769 0.681 0.491 -0.813 0.681 0.000 -0.769 -0.681 1.091 -2.263 -2.351 3.164 0.681 0.769 -1.091 2.263 2.351 -3.164 -0.681 -0.769 13.347 -29.750 -29.750 19.696 17.582 17.582 4.028 -9.390 -9.390 3.085 6.793 6.793 0.000 0.000 0.164 
+-2.263 3.823 1.318 0.000 -0.659 -1.318 -0.901 2.505 -1.318 0.000 0.659 1.318 -2.351 3.823 7.119 -6.328 -0.659 -3.955 2.351 -3.823 -7.119 6.328 0.659 3.955 -29.750 127.389 20.905 -39.393 -103.052 -1.208 -9.390 46.283 3.634 -6.170 -41.089 -0.549 0.000 0.000 -0.879 
+-2.351 7.119 0.659 0.000 -3.955 -0.659 -0.813 -0.791 -0.659 0.000 3.955 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -29.750 20.905 127.389 -39.393 -1.208 -103.052 -9.390 3.634 46.283 -6.170 -0.549 -41.089 0.000 0.000 -0.879 
+3.164 -6.328 0.000 0.000 0.000 0.000 3.164 -6.328 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 19.696 -39.393 -39.393 78.785 0.000 0.000 3.085 -6.170 -6.170 12.340 0.000 0.000 0.000 0.000 0.000 
+0.681 -0.659 -1.318 0.000 0.659 1.318 -0.681 0.659 1.318 0.000 -0.659 -1.318 0.769 -0.659 -3.955 0.000 0.659 3.955 -0.769 0.659 3.955 0.000 -0.659 -3.955 17.582 -103.052 -1.208 0.000 103.052 1.208 6.793 -41.089 -0.549 0.000 41.089 0.549 0.000 0.000 0.879 
+0.769 -3.955 -0.659 0.000 3.955 0.659 -0.769 3.955 0.659 0.000 -3.955 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 17.582 -1.208 -103.052 0.000 1.208 103.052 6.793 -0.549 -41.089 0.000 0.549 41.089 0.000 0.000 0.879 
+-1.091 2.351 0.681 0.000 -0.769 -0.681 -0.491 0.813 -0.681 0.000 0.769 0.681 0.491 -0.901 -0.813 3.164 -0.681 -0.769 -0.491 0.901 0.813 -3.164 0.681 0.769 4.028 -9.390 -9.390 3.085 6.793 6.793 13.347 -29.750 -29.750 19.696 17.582 17.582 0.000 0.000 0.164 
+2.263 -3.823 -1.318 0.000 0.659 1.318 0.901 -2.505 1.318 0.000 -0.659 -1.318 -0.813 2.505 -0.791 -6.328 0.659 3.955 0.813 -2.505 0.791 6.328 -0.659 -3.955 -9.390 46.283 3.634 -6.170 -41.089 -0.549 -29.750 127.389 20.905 -39.393 -103.052 -1.208 0.000 0.000 -0.879 
+2.351 -7.119 -0.659 0.000 3.955 0.659 0.813 0.791 0.659 0.000 -3.955 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 -9.390 3.634 46.283 -6.170 -0.549 -41.089 -29.750 20.905 127.389 -39.393 -1.208 -103.052 0.000 0.000 -0.879 
+-3.164 6.328 0.000 0.000 0.000 0.000 -3.164 6.328 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 3.085 -6.170 -6.170 12.340 0.000 0.000 19.696 -39.393 -39.393 78.785 0.000 0.000 0.000 0.000 0.000 
+-0.681 0.659 1.318 0.000 -0.659 -1.318 0.681 -0.659 -1.318 0.000 0.659 1.318 -0.769 0.659 3.955 0.000 -0.659 -3.955 0.769 -0.659 -3.955 0.000 0.659 3.955 6.793 -41.089 -0.549 0.000 41.089 0.549 17.582 -103.052 -1.208 0.000 103.052 1.208 0.000 0.000 0.879 
+-0.769 3.955 0.659 0.000 -3.955 -0.659 0.769 -3.955 -0.659 0.000 3.955 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 6.793 -0.549 -41.089 0.000 0.549 41.089 17.582 -1.208 -103.052 0.000 1.208 103.052 0.000 0.000 0.879 
+0.164 -0.879 -0.879 0.000 0.879 0.879 0.164 -0.879 -0.879 0.000 0.879 0.879 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.195 0.000 0.000 
+0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.164 -0.879 -0.879 0.000 0.879 0.879 0.164 -0.879 -0.879 0.000 0.879 0.879 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.195 0.000 
+0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.164 -0.879 -0.879 0.000 0.879 0.879 0.164 -0.879 -0.879 0.000 0.879 0.879 0.000 0.000 1.195 
+DEAL::solving degree = 2
+DEAL::Solver stopped within 3 - 45 iterations
+0.072 -0.058 0.015 0.002 0.016 -0.031 0.001 0.011 0.000 0.000 0.001 0.000 0.000 0.001 0.000 0.000 -0.040 0.084 -0.016 -0.006 0.000 0.000 
+-0.058 0.151 -0.056 0.015 -0.031 0.067 -0.039 0.001 0.000 -0.002 0.000 -0.001 -0.001 0.000 -0.002 0.000 -0.404 0.224 0.468 0.016 0.000 0.000 
+0.015 -0.056 0.151 -0.058 0.001 -0.039 0.067 -0.031 0.001 0.000 0.002 0.000 0.000 0.002 0.000 0.001 0.064 0.224 -0.468 -0.016 0.000 0.000 
+0.002 0.015 -0.058 0.072 0.011 0.001 -0.031 0.016 0.000 -0.001 0.000 0.000 0.000 0.000 -0.001 0.000 -0.057 0.084 0.016 0.006 0.000 0.000 
+0.016 -0.031 0.001 0.011 0.072 -0.058 0.015 0.002 0.000 -0.001 0.000 0.000 0.000 0.000 -0.001 0.000 0.043 -0.084 -0.016 0.006 0.000 0.000 
+-0.031 0.067 -0.039 0.001 -0.058 0.151 -0.056 0.015 0.001 0.000 0.002 0.000 0.000 0.002 0.000 0.001 -0.180 -0.224 0.468 -0.016 0.000 0.000 
+0.001 -0.039 0.067 -0.031 0.015 -0.056 0.151 -0.058 0.000 -0.002 0.000 -0.001 -0.001 0.000 -0.002 0.000 0.288 -0.224 -0.468 0.016 0.000 0.000 
+0.011 0.001 -0.031 0.016 0.002 0.015 -0.058 0.072 0.000 0.000 0.001 0.000 0.000 0.001 0.000 0.000 0.027 -0.084 0.016 -0.006 0.000 0.000 
+0.000 0.000 0.001 0.000 0.000 0.001 0.000 0.000 0.072 -0.058 0.015 0.002 0.016 -0.031 0.001 0.011 -0.006 0.000 0.000 -0.040 -0.016 0.084 
+0.000 -0.002 0.000 -0.001 -0.001 0.000 -0.002 0.000 -0.058 0.151 -0.056 0.015 -0.031 0.067 -0.039 0.001 0.016 0.000 0.000 -0.404 0.468 0.224 
+0.001 0.000 0.002 0.000 0.000 0.002 0.000 0.001 0.015 -0.056 0.151 -0.058 0.001 -0.039 0.067 -0.031 -0.016 0.000 0.000 0.064 -0.468 0.224 
+0.000 -0.001 0.000 0.000 0.000 0.000 -0.001 0.000 0.002 0.015 -0.058 0.072 0.011 0.001 -0.031 0.016 0.006 0.000 0.000 -0.057 0.016 0.084 
+0.000 -0.001 0.000 0.000 0.000 0.000 -0.001 0.000 0.016 -0.031 0.001 0.011 0.072 -0.058 0.015 0.002 0.006 0.000 0.000 0.043 -0.016 -0.084 
+0.001 0.000 0.002 0.000 0.000 0.002 0.000 0.001 -0.031 0.067 -0.039 0.001 -0.058 0.151 -0.056 0.015 -0.016 0.000 0.000 -0.180 0.468 -0.224 
+0.000 -0.002 0.000 -0.001 -0.001 0.000 -0.002 0.000 0.001 -0.039 0.067 -0.031 0.015 -0.056 0.151 -0.058 0.016 0.000 0.000 0.288 -0.468 -0.224 
+0.000 0.000 0.001 0.000 0.000 0.001 0.000 0.000 0.011 0.001 -0.031 0.016 0.002 0.015 -0.058 0.072 -0.006 0.000 0.000 0.027 0.016 -0.084 
+-0.040 -0.404 0.064 -0.057 0.043 -0.180 0.288 0.027 -0.006 0.016 -0.016 0.006 0.006 -0.016 0.016 -0.006 9.185 -7.910 -8.068 0.000 0.000 0.000 
+0.084 0.224 0.224 0.084 -0.084 -0.224 -0.224 -0.084 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -7.910 15.820 0.000 0.000 0.000 0.000 
+-0.016 0.468 -0.468 0.016 -0.016 0.468 -0.468 0.016 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -8.068 0.000 16.137 0.000 0     0.000 
+-0.006 0.016 -0.016 0.006 0.006 -0.016 0.016 -0.006 -0.040 -0.404 0.064 -0.057 0.043 -0.180 0.288 0.027 0.000 0.000 0.000 9.185 -8.068 -7.910 
+0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.016 0.468 -0.468 0.016 -0.016 0.468 -0.468 0.016 0.000 0.000 0     -8.068 16.137 0.000 
+0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.084 0.224 0.224 0.084 -0.084 -0.224 -0.224 -0.084 0.000 0.000 0.000 -7.910 0.000 15.820 
+DEAL::solving degree = 3
+DEAL::Solver stopped within 3 - 45 iterations
index 8a606c6f6a7569a6b86235a848782ca6dfef6c39..3d38b31a0651e63914ea97d6432659fa026edfe7 100644 (file)
@@ -68,22 +68,22 @@ test (const unsigned int degree)
 
   MatrixTools::create_mass_matrix (dof, q, mass_matrix);
 
-  mass_matrix.print_formatted (logfile, 3, false, 0, " ", 1);
+  mass_matrix.print_formatted (logfile, 3, false, 0, "0", 1);
 
   SolverControl           solver_control (3*dofs_per_cell,
                                           1e-8);
   PrimitiveVectorMemory<> vector_memory;
-  SolverCG<>              cg (solver_control, vector_memory);
+  SolverCG<>              solver (solver_control, vector_memory);
 
   Vector<double> tmp1(dofs_per_cell), tmp2(dofs_per_cell);
   for (unsigned int i=0; i<dofs_per_cell; ++i)
     tmp1(i) = 1.*Testing::rand()/RAND_MAX;
-  cg.solve (mass_matrix, tmp2, tmp1, PreconditionIdentity());
 
-  deallog << "Degree=" << degree
-          << ": " << solver_control.last_step()
-          << " iterations to obtain convergence."
-          << std::endl;
+  deallog << "solving degree = " << degree << std::endl;
+  check_solver_within_range(
+    solver.solve (mass_matrix, tmp2, tmp1,
+                  PreconditionIdentity()),
+    solver_control.last_step(), 3, 45);
 }
 
 
index 3db37251abb9b44453b22c09419aa07d083baedc..c7d97f8d0d2ffd702c984c8498fad225782507b7 100644 (file)
@@ -7,9 +7,8 @@
 0.000 -0.068 -0.068 0.000 -0.036 0.208 0.020 0.058 
 0.000 -0.068 -0.068 0.000 0.058 0.020 0.208 -0.036 
 0.068 0.000 0.000 0.068 0.020 0.058 -0.036 0.208 
-DEAL:cg::Starting value 1.92
-DEAL:cg::Convergence step 4 value 0
-DEAL::Degree=1: 4 iterations to obtain convergence.
+DEAL::solving degree = 1
+DEAL::Solver stopped within 3 - 45 iterations
 3.262 -2.918 -2.918 0.613 -0.457 -0.457 0.768 -0.598 -0.469 -0.299 0.598 0.469 0.768 -0.469 -0.598 -0.299 0.469 0.598 
 -2.918 5.836 0.000 -0.457 0.914 0.000 -0.469 0.000 0.000 -0.469 0.000 0.000 -0.598 0.000 1.195 0.598 0.000 -1.195 
 -2.918 0.000 5.836 -0.457 0.000 0.914 -0.598 1.195 0.000 0.598 -1.195 0.000 -0.469 0.000 0.000 -0.469 0.000 0.000 
@@ -28,9 +27,8 @@ DEAL::Degree=1: 4 iterations to obtain convergence.
 -0.299 0.598 -0.469 -0.170 -0.598 0.469 -0.299 -0.469 0.598 -0.170 0.469 -0.598 0.613 -0.457 -0.457 3.262 -2.918 -2.918 
 0.469 0.000 0.000 0.469 0.000 0.000 0.598 0.000 -1.195 -0.598 0.000 1.195 -0.457 0.914 0.000 -2.918 5.836 0.000 
 0.598 -1.195 0.000 -0.598 1.195 0.000 0.469 0.000 0.000 0.469 0.000 0.000 -0.457 0.000 0.914 -2.918 0.000 5.836 
-DEAL:cg::Starting value 2.32
-DEAL:cg::Convergence step 11 value 0
-DEAL::Degree=1: 11 iterations to obtain convergence.
+DEAL::solving degree = 1
+DEAL::Solver stopped within 3 - 45 iterations
 0.105 -0.050 -0.003 0.036 -0.030 -0.013 0.004 -0.002 0.004 -0.004 0.002 -0.004 -0.020 0.000 
 -0.050 0.169 -0.050 -0.030 0.058 -0.030 -0.002 -0.010 -0.002 0.002 0.010 0.002 -0.090 0.000 
 -0.003 -0.050 0.105 -0.013 -0.030 0.036 0.004 -0.002 0.004 -0.004 0.002 -0.004 -0.020 0.000 
@@ -45,9 +43,8 @@ DEAL::Degree=1: 11 iterations to obtain convergence.
 -0.004 0.002 -0.004 0.004 -0.002 0.004 -0.013 -0.030 0.036 -0.003 -0.050 0.105 0.000 -0.020 
 -0.020 -0.090 -0.020 -0.020 -0.090 -0.020 0.000 0.000 0.000 0.000 0.000 0.000 1.195 0.000 
 0.000 0.000 0.000 0.000 0.000 0.000 -0.020 -0.090 -0.020 -0.020 -0.090 -0.020 0.000 1.195 
-DEAL:cg::Starting value 2.25
-DEAL:cg::Convergence step 10 value 0
-DEAL::Degree=2: 10 iterations to obtain convergence.
+DEAL::solving degree = 2
+DEAL::Solver stopped within 3 - 45 iterations
 13.347 -29.750 -29.750 19.696 17.582 17.582 4.028 -9.390 -9.390 3.085 6.793 6.793 1.091 -2.351 -0.681 0.000 0.769 0.681 0.491 -0.813 0.681 0.000 -0.769 -0.681 1.091 -2.263 -2.351 3.164 0.681 0.769 -1.091 2.263 2.351 -3.164 -0.681 -0.769 0.164 0.000 0.000 
 -29.750 127.389 20.905 -39.393 -103.052 -1.208 -9.390 46.283 3.634 -6.170 -41.089 -0.549 -2.263 3.823 1.318 0.000 -0.659 -1.318 -0.901 2.505 -1.318 0.000 0.659 1.318 -2.351 3.823 7.119 -6.328 -0.659 -3.955 2.351 -3.823 -7.119 6.328 0.659 3.955 -0.879 0.000 0.000 
 -29.750 20.905 127.389 -39.393 -1.208 -103.052 -9.390 3.634 46.283 -6.170 -0.549 -41.089 -2.351 7.119 0.659 0.000 -3.955 -0.659 -0.813 -0.791 -0.659 0.000 3.955 0.659 -0.681 1.318 0.659 0.000 -1.318 -0.659 0.681 -1.318 -0.659 0.000 1.318 0.659 -0.879 0.000 0.000 
@@ -87,9 +84,8 @@ DEAL::Degree=2: 10 iterations to obtain convergence.
 0.164 -0.879 -0.879 0.000 0.879 0.879 0.164 -0.879 -0.879 0.000 0.879 0.879 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.195 0.000 0.000 
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.164 -0.879 -0.879 0.000 0.879 0.879 0.164 -0.879 -0.879 0.000 0.879 0.879 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.195 0.000 
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.164 -0.879 -0.879 0.000 0.879 0.879 0.164 -0.879 -0.879 0.000 0.879 0.879 0.000 0.000 1.195 
-DEAL:cg::Starting value 3.80
-DEAL:cg::Convergence step 40 value 0.00
-DEAL::Degree=2: 40 iterations to obtain convergence.
+DEAL::solving degree = 2
+DEAL::Solver stopped within 3 - 45 iterations
 0.072 -0.058 0.015 0.002 0.016 -0.031 0.001 0.011 0.000 0.000 0.001 0.000 0.000 0.001 0.000 0.000 -0.040 0.084 -0.016 -0.006 0.000 0.000 
 -0.058 0.151 -0.056 0.015 -0.031 0.067 -0.039 0.001 0.000 -0.002 0.000 -0.001 -0.001 0.000 -0.002 0.000 -0.404 0.224 0.468 0.016 0.000 0.000 
 0.015 -0.056 0.151 -0.058 0.001 -0.039 0.067 -0.031 0.001 0.000 0.002 0.000 0.000 0.002 0.000 0.001 0.064 0.224 -0.468 -0.016 0.000 0.000 
@@ -112,6 +108,5 @@ DEAL::Degree=2: 40 iterations to obtain convergence.
 -0.006 0.016 -0.016 0.006 0.006 -0.016 0.016 -0.006 -0.040 -0.404 0.064 -0.057 0.043 -0.180 0.288 0.027 0.000 0.000 0.000 9.185 -8.068 -7.910 
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 -0.016 0.468 -0.468 0.016 -0.016 0.468 -0.468 0.016 0.000 0.000 0.000 -8.068 16.137 0.000 
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.084 0.224 0.224 0.084 -0.084 -0.224 -0.224 -0.084 0.000 0.000 0.000 -7.910 0.000 15.820 
-DEAL:cg::Starting value 3.11
-DEAL:cg::Convergence step 14 value 0.00
-DEAL::Degree=3: 14 iterations to obtain convergence.
+DEAL::solving degree = 3
+DEAL::Solver stopped within 3 - 45 iterations
index c3befbb39b1ac1cdfdc04a70cb7da73ee0e2b92a..b2ab5acbc798d00b0f68fc50db9982d64c4b2879 100644 (file)
@@ -71,16 +71,16 @@ void test()
 
   // repartition the mesh as described above, first in some arbitrary
   // way, and then with all equal weights
-  tr.signals.cell_weight.connect(std::bind(&cell_weight_1<dim>,
-                                           std_cxx11::_1,
-                                           std_cxx11::_2));
+  tr.signals.cell_weight.connect(std_cxx11::bind(&cell_weight_1<dim>,
+                                                 std_cxx11::_1,
+                                                 std_cxx11::_2));
   tr.repartition();
 
   tr.signals.cell_weight.disconnect_all_slots();
 
-  tr.signals.cell_weight.connect(std::bind(&cell_weight_2<dim>,
-                                           std_cxx11::_1,
-                                           std_cxx11::_2));
+  tr.signals.cell_weight.connect(std_cxx11::bind(&cell_weight_2<dim>,
+                                                 std_cxx11::_1,
+                                                 std_cxx11::_2));
   tr.repartition();
 
 
index 6486bb8dd2b18df38046b9947004089f4b9c5566..85046d085042d39f8fb27607eef20ca389a3979e 100644 (file)
@@ -61,9 +61,9 @@ void test()
 
   // repartition the mesh as described above, first in some arbitrary
   // way, and then with no weights
-  tr.signals.cell_weight.connect(std::bind(&cell_weight_1<dim>,
-                                           std_cxx11::_1,
-                                           std_cxx11::_2));
+  tr.signals.cell_weight.connect(std_cxx11::bind(&cell_weight_1<dim>,
+                                                 std_cxx11::_1,
+                                                 std_cxx11::_2));
   tr.repartition();
 
   tr.signals.cell_weight.disconnect_all_slots();
index 9be95619af42eefc75bc26843b678649f3e77edc..d1e54904f438e57aca20fc47743992c9e00a7b31 100644 (file)
@@ -53,9 +53,9 @@ void test()
 
   GridGenerator::subdivided_hyper_cube(tr, 16);
 
-  tr.signals.cell_weight.connect(std::bind(&cell_weight<dim>,
-                                           std_cxx11::_1,
-                                           std_cxx11::_2));
+  tr.signals.cell_weight.connect(std_cxx11::bind(&cell_weight<dim>,
+                                                 std_cxx11::_1,
+                                                 std_cxx11::_2));
   tr.refine_global(1);
 
   if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)
index 8b124c50531b7d4e99dc982e0f6ae4c78e1b6a3f..2806b13b6c83e222ca993354c000f918e355c776 100644 (file)
@@ -64,9 +64,9 @@ void test()
 
   tr.refine_global(1);
 
-  tr.signals.cell_weight.connect(std::bind(&cell_weight<dim>,
-                                           std_cxx11::_1,
-                                           std_cxx11::_2));
+  tr.signals.cell_weight.connect(std_cxx11::bind(&cell_weight<dim>,
+                                                 std_cxx11::_1,
+                                                 std_cxx11::_2));
 
   tr.repartition();
 
index 64de6a11fa08db51a1caf03242e59e5a5b5173d3..0be7901f233fa4b32aa396d7234cba348fa8af76 100644 (file)
@@ -59,9 +59,9 @@ void test()
 
   GridGenerator::subdivided_hyper_cube(tr, 16);
 
-  tr.signals.cell_weight.connect(std::bind(&cell_weight<dim>,
-                                           std_cxx11::_1,
-                                           std_cxx11::_2));
+  tr.signals.cell_weight.connect(std_cxx11::bind(&cell_weight<dim>,
+                                                 std_cxx11::_1,
+                                                 std_cxx11::_2));
 
   tr.refine_global(1);
 
index 7d4865799fc560a133ad21445918edc82b1c5918..547bc9de95a17419f0047d3f20cf225a928359e1 100644 (file)
@@ -79,9 +79,9 @@ void test()
 
 
   // repartition the mesh; attach different weights to all cells
-  tr.signals.cell_weight.connect(std::bind(&cell_weight<dim>,
-                                           std_cxx11::_1,
-                                           std_cxx11::_2));
+  tr.signals.cell_weight.connect(std_cxx11::bind(&cell_weight<dim>,
+                                                 std_cxx11::_1,
+                                                 std_cxx11::_2));
 
   tr.repartition ();
 
index 55386a5dfb250bee787729aec2474ed6fd8ea6b4..9d5c095f8e42bf53d60dda0549f2b4a4cf01ef1b 100644 (file)
@@ -79,9 +79,9 @@ void test()
 
   // repartition the mesh; attach different weights to all cells
   n_global_active_cells = tr.n_global_active_cells();
-  tr.signals.cell_weight.connect(std::bind(&cell_weight<dim>,
-                                           std_cxx11::_1,
-                                           std_cxx11::_2));
+  tr.signals.cell_weight.connect(std_cxx11::bind(&cell_weight<dim>,
+                                                 std_cxx11::_1,
+                                                 std_cxx11::_2));
   tr.repartition ();
 
   if (Utilities::MPI::this_mpi_process (MPI_COMM_WORLD) == 0)

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.