diff --git a/dune/istl/matrix.hh b/dune/istl/matrix.hh index 1dc4fb9f62b028ffd6c98b1beb52e7d5c66067ae..5b32bc37e1fd7800ba0477df71d1f25901fe3470 100644 --- a/dune/istl/matrix.hh +++ b/dune/istl/matrix.hh @@ -320,7 +320,7 @@ namespace Dune { /// Generic matrix multiplication. friend Matrix<T> operator*(const Matrix<T>& m1, const Matrix<T>& m2) { Matrix<T> out(m1.N(), m2.M()); - out.clear(); + out = 0; for (size_type i=0; i<out.N(); i++ ) { for ( size_type j=0; j<out.M(); j++ )