Stephane,
> Is it possible to perform calculus on non-conforming mesh with DEAL II ?
I have no solid idea about that, but...
> When I restart my calculus from the mesh I saved, results are wrong,
> like if hanging nodes were note taken into account correctly.
If you got as far as creating a triangulation and writing the
local-to-global matrix elements, something is very right!
Did you forget to apply hanging node constraints to your solution vector?
Something like: hanging_node_constraints.distribute (solution_vector);
just *after* getting the solution from your solver is generally needed.
See if doing that helps things move along...
Of course, this will only work if your hanging nodes really do exist in
the triangulation after GridIn excecutes. You should definately check that
first... :-)
Best,
Toby
-----
Toby D. Young
Assistant Professor
Philosophy-Physics
Polish Academy of Sciences
Warszawa, Polska
www: http://www.ippt.gov.pl/~tyoung
skype: stenografia
On Tue, 4 May 2010, veys wrote:
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
