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

Reply via email to