comp_with_fs.h

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

#ifndef comp_with_fsH
#define comp_with_fsH

#include "../davidson/matrixbase.h"

/*
This class make possible comparison this implementation of Davidson method
with Andreas Stathopoulos and Charlotte F. Fischer implementation.

For more informations see:
	A DAVIDSON PROGRAM FOR FINDING A FEW SELECTED EXTREME
	EIGENPAIRS OF A LARGE, SPARSE, REAL, SYMMETRIC MATRIX.
          		A. STATHOPOULOS, C.F. FISCHER.
          REF. IN COMP. PHYS. COMMUN. 79 (1994) 268

*/

class CompareWithFS: public MatrixBase
{
    private:

		ntyp Nmax,
        	 IBAND,
             NZERmax;

        ntyp N;

        ftyp *A, *DIAG;
        ntyp *IndCol, *IndRow;

        bool lupper;

        void dinit(ntyp , ftyp , ftyp *, ntyp ) const;
		void gather(ntyp , ftyp *, ftyp *, ntyp *) const;
        void scatter(ntyp , ftyp *, ntyp *, ftyp *) const;
    	ftyp ddot(ntyp , ftyp *, ntyp , ftyp *, ntyp) const;
        void daxpy(ntyp , ftyp , ftyp *, ntyp , ftyp *, ntyp) const;

        void Init();

    public:
        CompareWithFS();

        ftyp *Diagonal() const;

        void dssbmv(ntyp , ftyp *, ftyp *, ntyp ) const;
        ntyp GetN() const;

        void print(Davidson *) const;

        ~CompareWithFS();

};

#endif