Hello DealII's users,
I am trying to explain what I would like to get:
I have two triangulations that are read by "read_ucd()" function.
I distribute two different DoFHandlers, each one for each triangulation:
std::ifstream in ("file1.dat");
GridIn<DIMENSION> grid_leida1;
grid_leida1.attach_triangulation(triangulation1);
grid_leida1.read_ucd(in1);
dof_handler_1.distribute_dofs(fe_1);
DoFRenumbering::component_wise(dof_handler_1);
[
.
. // I solve a equation with this data and I obtain vector_solution_1
.
]
std::ifstream in ("file2.dat");
GridIn<DIMENSION> grid_leida2;
grid_leida2.attach_triangulation(triangulation2);
grid_leida2.read_ucd(in2);
dof_handler_2.distribute_dofs(fe_2);
DoFRenumbering::component_wise(dof_handler_2);
[
.
. // I solve a equation with this data and I obtain vector_solution_2
.
]
For each domain, with its triangulation, its dof_handler, its
FESystem..., I obtain its vector_solution_1 and vector_solution_2,
respectively.
Once I have these two vectors, totally independent, I would like to
transfer them to a common (DIMENSION-1) mesh, I mean a third
triangulation3 with its dof_handler_3, its FEsystem... I don't have
created this new triangulation3, but it should introduce the less
distortion as possible. And if, at the beginning, I am working with a 3D
triangulation, this third one would be a 2D triangulation. So, I would
like to obtain two new vectors: "vector_solution_from_1_to_3" and
"vector_solution_from_2_to_3", both belonging to triangulation3,
dofhandler3...
I know I cannot use "SolutionTransfer" method because the original grid
isn't the same for the two of them.
I suppose I would have to use some kinds of interpolation function...
Could you give me some ideas in order to create the new triangulation
and solutions in this one?
Thanks in advance!
Best regards!
Isa
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
