int main()
{
ofstream logfile("mglocal.output");
-// logfile.setf(ios::fixed);
-// logfile.precision (3);
+ logfile.setf(ios::fixed);
+ logfile.precision (3);
deallog.attach(logfile);
- deallog.log_execution_time(true);
-// deallog.depth_console(0);
+// deallog.log_execution_time(true);
+ deallog.depth_console(0);
Helmholtz equation;
RHSFunction<2> rhs;
hanging_nodes.condense(A);
hanging_nodes.condense(f);
- if (true)
+ if (false)
{
ofstream out_file("MGf");
DataOut<2> out;
u.reinit(f);
PrimitiveVectorMemory<> mem;
- SolverControl control(20, 1.e-12, true);
+ SolverControl control(20, 1.e-12, false);
SolverCG<> solver(control, mem);
MGLevelObject<SparsityPattern> mgstruct(0, tr.n_levels()-1);
solver.solve(A, u, f, mgprecondition);
hanging_nodes.distribute(u);
- DataOut<2> out;
- char* name = new char[100];
+ if (false)
+ {
+ DataOut<2> out;
+ char* name = new char[100];
- sprintf(name, "MG-Q%d-%d", degree, step);
+ sprintf(name, "MG-Q%d-%d", degree, step);
- ofstream ofile(name);
- out.attach_dof_handler(dof);
- out.add_data_vector(u,"u");
- out.add_data_vector(f,"f");
- out.build_patches(5);
- out.write_gnuplot(ofile);
-
- delete[] name;
+ ofstream ofile(name);
+ out.attach_dof_handler(dof);
+ out.add_data_vector(u,"u");
+ out.add_data_vector(f,"f");
+ out.build_patches(5);
+ out.write_gnuplot(ofile);
+ delete[] name;
+ }
}
deallog.pop();
}
-DEAL:Q1::DoFs 137
-DEAL:Q1::Levels: 5
+DEAL:Q1::DoFs 41
+DEAL:Q1::Levels: 4
DEAL:Q1::smoothing-steps1
-DEAL:Q1:Richardson::Check 0 1.811
-DEAL:Q1:Richardson::Starting
-DEAL:Q1:Richardson::Check 1 0.317
-DEAL:Q1:Richardson::Check 2 0.058
-DEAL:Q1:Richardson::Check 3 0.009
-DEAL:Q1:Richardson::Check 4 0.001
-DEAL:Q1:Richardson::Check 5 0.000
-DEAL:Q1:Richardson::Check 6 0.000
-DEAL:Q1:Richardson::Check 7 0.000
-DEAL:Q1:Richardson::Check 8 0.000
-DEAL:Q1:Richardson::Check 9 0.000
-DEAL:Q1:Richardson::Check 10 0.000
-DEAL:Q1:Richardson::Check 11 0.000
-DEAL:Q1:Richardson::Check 12 0.000
-DEAL:Q1:Richardson::Check 13 0.000
-DEAL:Q1:Richardson::Check 14 0.000
-DEAL:Q1:Richardson::Check 15 0.000
-DEAL:Q1:Richardson::Check 16 0.000
-DEAL:Q1:Richardson::Convergence step 16
+DEAL:Q1:cg::Starting
+DEAL:Q1:cg::Convergence step 8
DEAL:Q1::smoothing-steps2
-DEAL:Q1:Richardson::Check 0 1.811
-DEAL:Q1:Richardson::Starting
-DEAL:Q1:Richardson::Check 1 0.345
-DEAL:Q1:Richardson::Check 2 0.052
-DEAL:Q1:Richardson::Check 3 0.005
-DEAL:Q1:Richardson::Check 4 0.000
-DEAL:Q1:Richardson::Check 5 0.000
-DEAL:Q1:Richardson::Check 6 0.000
-DEAL:Q1:Richardson::Check 7 0.000
-DEAL:Q1:Richardson::Check 8 0.000
-DEAL:Q1:Richardson::Check 9 0.000
-DEAL:Q1:Richardson::Check 10 0.000
-DEAL:Q1:Richardson::Check 11 0.000
-DEAL:Q1:Richardson::Check 12 0.000
-DEAL:Q1:Richardson::Check 13 0.000
-DEAL:Q1:Richardson::Check 14 0.000
-DEAL:Q1:Richardson::Check 15 0.000
-DEAL:Q1:Richardson::Convergence step 15
-DEAL:Q1::smoothing-steps3
-DEAL:Q1:Richardson::Check 0 1.811
-DEAL:Q1:Richardson::Starting
-DEAL:Q1:Richardson::Check 1 0.367
-DEAL:Q1:Richardson::Check 2 0.055
-DEAL:Q1:Richardson::Check 3 0.004
-DEAL:Q1:Richardson::Check 4 0.000
-DEAL:Q1:Richardson::Check 5 0.000
-DEAL:Q1:Richardson::Check 6 0.000
-DEAL:Q1:Richardson::Check 7 0.000
-DEAL:Q1:Richardson::Check 8 0.000
-DEAL:Q1:Richardson::Check 9 0.000
-DEAL:Q1:Richardson::Check 10 0.000
-DEAL:Q1:Richardson::Check 11 0.000
-DEAL:Q1:Richardson::Check 12 0.000
-DEAL:Q1:Richardson::Check 13 0.000
-DEAL:Q1:Richardson::Check 14 0.000
-DEAL:Q1:Richardson::Check 15 0.000
-DEAL:Q1:Richardson::Convergence step 15
-DEAL:Q1::smoothing-steps4
-DEAL:Q1:Richardson::Check 0 1.811
-DEAL:Q1:Richardson::Starting
-DEAL:Q1:Richardson::Check 1 0.387
-DEAL:Q1:Richardson::Check 2 0.060
-DEAL:Q1:Richardson::Check 3 0.004
-DEAL:Q1:Richardson::Check 4 0.000
-DEAL:Q1:Richardson::Check 5 0.000
-DEAL:Q1:Richardson::Check 6 0.000
-DEAL:Q1:Richardson::Check 7 0.000
-DEAL:Q1:Richardson::Check 8 0.000
-DEAL:Q1:Richardson::Check 9 0.000
-DEAL:Q1:Richardson::Check 10 0.000
-DEAL:Q1:Richardson::Check 11 0.000
-DEAL:Q1:Richardson::Check 12 0.000
-DEAL:Q1:Richardson::Check 13 0.000
-DEAL:Q1:Richardson::Check 14 0.000
-DEAL:Q1:Richardson::Check 15 0.000
-DEAL:Q1:Richardson::Convergence step 15
+DEAL:Q1:cg::Starting
+DEAL:Q1:cg::Convergence step 7
+DEAL:Q2::DoFs 141
+DEAL:Q2::Levels: 4
+DEAL:Q2::smoothing-steps1
+DEAL:Q2:cg::Starting
+DEAL:Q2:cg::Failure step 20
+DEAL:Q2::smoothing-steps2
+DEAL:Q2:cg::Starting
+DEAL:Q2:cg::Failure step 20
+DEAL:Q3::DoFs 297
+DEAL:Q3::Levels: 4
+DEAL:Q3::smoothing-steps1
+DEAL:Q3:cg::Starting
+DEAL:Q3:cg::Failure step 20
+DEAL:Q3::smoothing-steps2
+DEAL:Q3:cg::Starting
+DEAL:Q3:cg::Failure step 20