Hello!
I would like to interpolate the solution, "vector1", obtained in a
"triangulation1" with its "dof_handler1" into a solution, "vector2" over
a "triangulation2" with its "dof_handler2".
Both meshes make reference to the same domain, so I have decided, as you
advised me, using the "Functions::FEFieldFunction class", and I do:
******
Functions::FEFieldFunction<DIMENSION> fe_function1
(dof_handler_total1,vector1);
VectorTools::interpolate(dof_handler_total2,fe_function1,vector2);
******
where "vector1" has a size of "dof_handler_total1.n_dofs()"and it
contains the solution over the "triangulation1" and
"vector2" has "dof_handler_total2.n_dofs()" and it is the vector over
which I want to interpolate the value.
The problem is when I execute the program, because "fe_function1" is
created without problems but, later, when running
"VectorTools::interpolate(...)", the program breaks down showing me the
following message:
******
p0_29592: p4_error: interrupt SIGSEGV: 11
******
I would be really grateful if you could tell me what I am doing wrong.
Thanks in advance.
Best
Isa
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii