#ifndef WR_MATRIXTRAITS_HH_INCLUDED
#define WR_MATRIXTRAITS_HH_INCLUDED

#include <dune/common/stdstreams.hh>
#include <dune/fem/operator/matrix/spmatrix.hh>

template <class RowSpaceImp, class ColSpaceImp = RowSpaceImp>
struct MatrixTraits : public Dune::SparseRowMatrixTraits<RowSpaceImp,ColSpaceImp> {
    struct StencilType {
        int nonZerosEstimate( ColSpaceImp& rangeSpace_ ) {
            return 0;
        }
    };
};

/* usage:
    typedef Dune::SparseRowMatrixObject<  DiscreteFunctionSpaceTypeA, DiscreteFunctionSpaceTypeB,
                                            MatrixTraits<DiscreteFunctionSpaceTypeA, DiscreteFunctionSpaceTypeB >
        MassMatrixType;
*/
#endif