I wonder if it possible to re-use FEEvaluation within each cell to
calculate integrals of multiple right-hand-side like vectors.
Something along these lines:
for (unsigned int cell=0; cell<n_cells; ++cell)
{
fe_eval.reinit(cell);
for (unsigned int q=0; q<n_q_points; ++q)
// evaluate RHS functions and store results in
values[q][i]
for (unsigned int i = 0; i < n_rhs; ++i)
{
for (unsigned int q= 0; q < n_q_points; ++q)
fe_eval.submit_value(values[q][i], q);
fe_eval.integrate (true,false);
fe_eval.distribute_local_to_global (rhs_vectors[i]);
}
}
Essentially, the question is whether or not integrate()
and/or distribute_local_to_global() will re-set the FEEvaluation object to
the state right after fe_eval.reinit(cell).
If not, then i suppose my best bet is to move the reinit within the loop
over RHS vectors.
p.s. n_rhs is large, more than just a few, so I would prefer not to keep a
vector of shared pointers to FEEvaluation or alike.
Regards,
Denis.
--
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].
For more options, visit https://groups.google.com/d/optout.
