#1565 Implement a ParameterOperator to compute the global coordinates of dofs.
This is useful for the computation of the volume within a closed surface when we are dealing with non isoparametric formulations (for instance a P2 unknown with P1 geometry). We need to compute:
V = - \frac{1}{3} (\underline{X}_0 + \underline{Y}) \cdot \int_{\partial \Omega} \underline{n}_t \cdot \underline{y}^* \, \textrm{d} S
Where \underline{X}_0
is the global position of the dofs associated to \underline{Y}
.
We have in fact 3 possible implementations:
- Use a free function taking in the relevant finite element space and the relevant unkown
- Implement a ParameterOperator which implies a rework of them as they are AtQuadraturePoint oriented
- Create a Linear VariationalOperator but we need to extend the API of the GlobalOperator so that we can choose the strategy used for the insertion of the values (
INSERT_VALUES
instead of the usualADD_VALUES
for instance).
Edited by DIAZ Jerome