Hello Javier,
what does happen, if you replace the "1000" in the function
AdvectionProblem<dim>::solve with "system_rhs.size ()"?
Best Regards,
Markus
Am 12.01.2012 16:16, schrieb Javier Muñoz:
Hi
I am practicing the tutorials examples of Deal.ii. In step-9, which
deals with the advection equation and shows the use of refinement
meshing, I am trying to use a different mesh, generated with Gmsh (I
attach the file).
if (cycle == 0)
{
// GridGenerator::hyper_cube (triangulation, -1, 1);
// triangulation.refine_global (4);
GridIn<dim> grid_in;
grid_in.attach_triangulation (triangulation);
std::ifstream input_file ("square.msh");
grid_in.read_msh (input_file);
}
else
{
refine_grid ();
};
The domain is the same and no other modifications were made. After 3
cycles the program gives me an error about convergence:
Cycle 0:
Number of active cells: 416
Number of degrees of freedom: 453
Cycle 1:
Number of active cells: 1043
Number of degrees of freedom: 1204
Cycle 2:
Number of active cells: 2657
Number of degrees of freedom: 3077
----------------------------------------------------
Exception on processing:
Iterative method reported convergence failure in step 1000 with residual
3.12503e-07
Aborting!
----------------------------------------------------
I have no idea why I'm receiving this error. Can you give me any
suggestions?
Best Regards
Javier Muñoz
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
