solution = 0;
- try
- {
cg.solve (system_matrix, solution, system_rhs,
preconditioner);
- }
- catch (...)
- {
- output_results(42);
- constraints.distribute (solution);
- output_results(43);
- MPI_Barrier(MPI_COMM_WORLD);
- exit(0);
- }
-
constraints.distribute (solution);
}
GridRefinement::refine_and_coarsen_fixed_number (triangulation,
estimated_error_per_cell,
0.3, 0.03);
-
triangulation.execute_coarsening_and_refinement ();
}
TrilinosWrappers::MPI::Vector res_ghosted = temp_solution;
res_ghosted = temp;
-
data_out.attach_dof_handler (mg_dof_handler);
data_out.add_data_vector (temp_solution, "solution");
data_out.add_data_vector (res_ghosted, "res");
using namespace dealii;
using namespace Step50;
- LaplaceProblem<2> laplace_problem(3);
+ LaplaceProblem<2> laplace_problem(1);
laplace_problem.run ();
}
catch (std::exception &exc)