diff --git a/dune/solvers/common/staticmatrixtools.hh b/dune/solvers/common/staticmatrixtools.hh index af26041b3952a87d89c2f5517f1bc638d4eff883..041baf14203af30ff95c742490c86a72c233d429 100644 --- a/dune/solvers/common/staticmatrixtools.hh +++ b/dune/solvers/common/staticmatrixtools.hh @@ -88,6 +88,11 @@ class StaticMatrix x.diagonal()[i] += a * y[i][i]; } + template <typename FieldType, int n, typename Scalar> + static void axpy(Dune::FieldVector<FieldType,n>& x, const Scalar a, const Dune::FieldVector<FieldType,n>& y) + { + x.axpy(a, y); + } // add transformed matrix X = A^T*Y*B ******************************************************