Right now OperatorBase does it all by swtiching around if-statements. The proper way to do this is to derive the various operators from OperatorBase. Operator base would then have some properties from finite element methods such as quadratures and basis functions, but the derived classes would be responsible to implement how those operators are applied, similarly to the concept of LinearOperator from deal.II.
Right now OperatorBase does it all by swtiching around if-statements. The proper way to do this is to derive the various operators from OperatorBase. Operator base would then have some properties from finite element methods such as quadratures and basis functions, but the derived classes would be responsible to implement how those operators are applied, similarly to the concept of LinearOperator from deal.II.