bare-ham.h

Go to the documentation of this file.

/* Copyright (C) 2014  Ward Poelmans

This file is part of Hubbard-GPU.

Hubbard-GPU is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.

Hubbard-GPU is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.

You should have received a copy of the GNU General Public License
along with Hubbard-GPU.  If not, see <http://www.gnu.org/licenses/>.
*/

#ifndef BARE_HAM_H
#define BARE_HAM_H

#include <vector>
#include <string>

typedef unsigned int myint;

extern "C" {
    void dsyevd_( char* jobz, char* uplo, int* n, double* a, int* lda, double* w, double* work, int* lwork, int* iwork, int* liwork, int* info);
    double ddot_(int *n,double *x,int *incx,double *y,int *incy);
    void dscal_(int *n,double *alpha,double *x,int *incx);
    void dsymv_(char *uplo, const int *n, const double *alpha, const double *a, const int *lda, const double *x, const int *incx, const double *beta, double *y, const int *incy);
    void daxpy_(int *n,const double *alpha,double *x,int *incx,double *y,int *incy);
    void dstev_( const char* jobz, const int* n, double* d, double* e, double* z, const int* ldz, double* work, int* info );
    void dgesvd_( char* jobu, char* jobvt, int* m, int* n, double* a, int* lda, double* s, double* u, int* ldu, double* vt, int* ldvt, double* work, int* lwork, int* info );
    void dgemm_(char *transA,char *transB,const int *m,const int *n,const int *k,double *alpha,double *A,const int *lda,double *B,const int *ldb,double *beta,double *C,const int *ldc);

    void dsaupd_(int *ido, char *bmat, int *n, char *which,
            int *nev, double *tol, double *resid, int *ncv,
            double *v, int *ldv, int *iparam, int *ipntr,
            double *workd, double *workl, int *lworkl, int *info);

    void dseupd_(int *rvec, char *All, int *select, double *d,
            double *z, int *ldz, double *sigma,
            char *bmat, int *n, char *which, int *nev,
            double *tol, double *resid, int *ncv, double *v,
            int *ldv, int *iparam, int *ipntr, double *workd,
            double *workl, int *lworkl, int *info);
}

class BareHamiltonian
{
    public:
        BareHamiltonian(int L, int Nu, int Nd, double J, double U);
        virtual ~BareHamiltonian();

        static int CalcDim(int L, int N);

        static int CountBits(myint bits);

        std::string print_bin(myint num) const;

        static std::string print_bin(myint num,int bitcount);

        virtual void BuildBase();

        virtual void BuildFullHam() = 0;

        virtual void BuildHam();

        int getL() const;
        int getNu() const;
        int getNd() const;
        int getDim() const;

        double getJ() const;
        double getU() const;

        void setU(double);

        myint getBaseUp(unsigned int i) const;
        myint getBaseDown(unsigned int i) const;

        virtual std::vector<double> ExactDiagonalizeFull(bool calc_eigenvectors=false);
        virtual double LanczosDiagonalize(int m=0);
        virtual double arpackDiagonalize();

        void Print(bool list=false) const;

        void PrintBase() const;

        void PrintGroundstateVector() const;

        virtual void mvprod(double *x, double *y, double alpha) const = 0;

        double MemoryNeededFull() const;

        double MemoryNeededLanczos() const;

        double MemoryNeededArpack() const;

        virtual void SaveToFile(const std::string filename) const;

    protected:
        int CalcSign(int i,int j,myint a) const;

        static void Diagonalize(int dim, double *mat, double *eigs, bool calc_eigenvectors);

        void SaveToFile(const std::string filename, double *data, int dim) const;

        int L;
        int Nu;
        int Nd;
        double J;
        double U;
        double *ham;
        int dim;

        myint Hb;
        int Hbc;

        std::vector<myint> baseUp;
        std::vector<myint> baseDown;
};

#endif /* BARE_HAM_H */

/* vim: set ts=8 sw=4 tw=0 expandtab :*/