Dear Sambit,
If the result of the two cases is different and
ConstraintMatrix::distribute() was called in both cases, I expect there
to be some confusion regarding the indices of ghost entries. In debug
mode, there should be a check that the parallel partitioner of the
vector inside FEEvaluation::read_dof_values* does match with the
expected index numbering. Did you run in debug mode? To localize the
issue, can you check whether you called
"matrix_free_data.initialize_dof_vector(VECTOR_NAME, 1);" in the second
case? Note the optional argument "1" that must match with the "1" passed
to FEEvaluation. If the issue still appears, it must be because
ConstraintMatrix::distribute() does not do all updates. In that case, I
would appreciate if you can give us a workable example.
Best,
Martin
--
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.
