diff --git a/dune/solvers/transferoperators/genericmultigridtransfer.hh b/dune/solvers/transferoperators/genericmultigridtransfer.hh index 92f0720904fdb09915be1ac4d263ac76f2ecee70..d50319f583573c7a060ca0bc2ab4acfe3dac07b5 100644 --- a/dune/solvers/transferoperators/genericmultigridtransfer.hh +++ b/dune/solvers/transferoperators/genericmultigridtransfer.hh @@ -13,6 +13,7 @@ #include <dune/localfunctions/lagrange/pqkfactory.hh> #include "dune/solvers/common/staticmatrixtools.hh" +#include <dune/solvers/common/arithmetic.hh> /** \brief Restriction and prolongation operator for standard multigrid @@ -825,7 +826,7 @@ public: // Compute cm = im^T * m * jm if(TransferOperatorType::block_type::rows==1) - StaticMatrix::axpy(cm, (*im)[0][0] * (*jm)[0][0], *m); + Arithmetic::addProduct(cm, (*im)[0][0] * (*jm)[0][0], *m); else StaticMatrix::addTransformedMatrix(cm, *im, *m, *jm); } @@ -892,7 +893,7 @@ public: // Compute cm = im^T * m * jm if(TransferOperatorType::block_type::rows==1) - StaticMatrix::axpy(cm, (*im)[0][0] * (*jm)[0][0], *m); + Arithmetic::addProduct(cm, (*im)[0][0] * (*jm)[0][0], *m); else StaticMatrix::addTransformedMatrix(cm, *im, *m, *jm); }