WaveMultiField don't support inhomogenuous Dirichlet BC
Following !4 (merged) , WMF support homogenuous Dirichlet BCs. However, when calling WaveAuxiliaryField::assignValues()
, the constrained DOFs get value 0, irrelevant of the value of the constraint. However, the Auxiliary Field produces the good result, so it's just that results on the boundary will be off but correct elsewhere
Potential fixes:
- When calling
WaveMultiField::operator[]()
, use some FEM call to re-apply the BC afterwards. - In the
_dofValues
of the WMF, store at initialization the boundary values ?