Introduce internal::MatrixFreeFunctions::ConstraintInfo #13519
Conversation
There was a problem hiding this comment.
@kronbichler As discussed, I have copied together the most relevant functions to be able to use them in MGTwoLevelTransfer. This is not a perfect solution and should be regarded as intermediate step, but should give us a feeling what to generalize in the other classes.
| } | ||
|
|
||
| void | ||
| read_dof_indices( |
There was a problem hiding this comment.
This function and the next one are adopted from DoFInfo and MatrixFree.
|
|
||
| template <typename T, typename VectorType> | ||
| void | ||
| read_write_operation(const T & operation, |
There was a problem hiding this comment.
This and the next one are adopted from FEEvaluationBase.
| } | ||
|
|
||
| void | ||
| finalize() |
There was a problem hiding this comment.
A difference to the setup in the original code is that I am first collecting everything in vectors of vectors and after that convert them to CRS data structures. I do this to put the data to the right position during the cell loop in the setup of the global-coarsening transfer. But what could do is to loop over all cells category by category or alternatively to sort the data in a second step. Both needs some more involved refactoring in MGTwoLevelTransfer so that I would like to postpone this to a follow-up PR!
There was a problem hiding this comment.
Yes, that makes sense. In any case, this looks like a good intermediate step what could later be unified with how FEEvaluation/DoFInfo treats this case with one set of code.
|
/rebuild |
peterrum
force-pushed
the
constraint_info
branch
from
c4a13d0
to
eb907a7
Compare
|
@kronbichler I have split up the declaration and interpolation and have added a test. |
peterrum
force-pushed
the
constraint_info
branch
2 times, most recently
from
def708e
to
ccb8c1d
Compare
peterrum
force-pushed
the
constraint_info
branch
from
ccb8c1d
to
653832d
Compare
Add dummy test
peterrum
force-pushed
the
constraint_info
branch
from
653832d
to
d3ea260
Compare
A class that allows to efficiently evaluate constraints in matrix-free loops in user codes. The first application is the global-coarsening transfer operator in #13520.