void DiagonalMatrixTemplate<T>::mulVector(const VectorT& a, VectorT& b) const { b.componentMul(a,*this); }