I am trying to implement a staggered thermal-elastic solver using the
matrix-free multigrid solver approach from step-37. I got an elasticity
solver and a thermal solver working and now I am trying to associate the
two. I have a shared triangulation, with two separate dof_handlers for the
thermal and stress problems. The question I have is when I do the assembly
loop, does the cell refer to the same cell for both the thermal and
elastic problems?
So from step-37, can I just call phi_temp.reinit(cell), and
phi_stress.reinit(cell), and have cell refer to the same physical cell in
the triangulation?
for (unsigned int cell = 0;
cell < system_matrix.get_matrix_free()->n_cell_batches();
++cell)
{
phi.reinit(cell);
for (const unsigned int q : phi.quadrature_point_indices())
phi.submit_value(make_vectorized_array
<https://www.dealii.org/current/doxygen/deal.II/vectorization_8h.html#ad7d7e08942faeecf438c75a254e06cbe><double>(1.0),
q);
phi.integrate(EvaluationFlags::values
<https://www.dealii.org/current/doxygen/deal.II/namespaceEvaluationFlags.html#a9b7c6d689cb76386839d0d13640f59aeaf9825c682f693a6a200094641a0d6a58>
);
phi.distribute_local_to_global(system_rhs);
}
Thanks!
Subramanya.
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion visit
https://groups.google.com/d/msgid/dealii/d8f59b24-2a30-4ade-83e4-75bc72112afen%40googlegroups.com.
