LocalFiniteElementShapeFunctionSet is not thread safe.
The evaluateEach (https://gitlab.dune-project.org/dune-fem/dune-fem/-/blob/master/dune/fem/space/localfiniteelement/shapefunctionset.hh#L81), jacobianEach, hessianEach are not thread safe. It does not seem enough to make the local values_
variable thread safe. There is also an issue with localBasis. It is very difficult to find out which code in dune-localfunctions actually is called here.
This becomes a problem when threading
is used with ElementQuadrature
.
Edited by Robert K