Dear Martin,
Thank you for your quick reply. I understand my mistake now- the issue
stems from "matrix_free_data.initialize_dof_vector(VECTOR_NAME, 1);" in the
second
case. I have used matrix_free_data.initialize_dof_vector(VECTOR_NAME, 0)
in both cases which causes a mismatch with the
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
Hi All,
I have reduced my bug to the following minimal example: I create a
MatrixFree object for two different ConstraintMatrices (provided as vector
of ConstraintMatrices to the reinit(..)).
matrix_free_data.reinit(dofHandlerVector, d_constraintsVector,
quadratureVector,