Peter Munch
Peter Munch
Question: why do you have a look at the distributed implementation? The distributed trianulations do not support anisotropy. I guess the function that needs to be fixed is: https://github.com/dealii/dealii/blob/91bff6e62a6d9fc1fe71e61dd5eeab497ad618ef/source/grid/tria.cc#L13926-L13935
My best guess is that https://github.com/dealii/dealii/blob/b0114df040e3873837fa5ac149dafb05dc4dfabe/include/deal.II/matrix_free/operators.h#L2073 needs to call the other constructor which also takes the range. This is needed for hp.
@kronbichler I don't think we will be able to implement a fix till the release. Postponing.
@nabw I'll give the PR another round in the next days.
> Could we put these in MatrixTools instead? We don't use that namespace much but I would prefer not to introduce a new one that is similar. I know the...
> Regarding the namespace, I have no strong opinion. SparseMatrixTools sounds OK in my opinion. One question I'm wondering, without having looked at it: Did you compare this to the...
I have updated PR and removed the dependency to #14124. Ready from my side.
@kronbichler I have fixed the typos and updated the documentation.
@marcfehling You could also call `AffineConstraints::make_consistent_in_parallel()`; after that you should be on the safe side.
Postponing and setting this (and related topics) onto the agenda of the workshop.