matrixbase.h

An object-oriented C++ implementation of Davidson method for finding a few selected extreme eigenpai

#ifndef __MATRIX_BASE
#define __MATRIX_BASE

#include "../davidson/defaulttd.h"

class MatrixBase
{
	public:
        // the method returns pointer for the vector of the diagonal elements
    	virtual ftyp *Diagonal() const=0;

        /* computes the product C of matrix A (of size n) and a block m of
           vectors B (of size n), that is C=A B*/
        virtual void dssbmv(ntyp , ftyp *, ftyp *, ntyp ) const=0;

        /* returns the dimension of the table containing the diagonal elements
           (equal to dimension of matrix A) */
        virtual ntyp GetN() const=0;

        //this operator returns  the value of the i-th diagonal element
        ftyp operator[](const ntyp &i) const {return *(Diagonal()+i);};
};

#endif