harpaham.h

Go to the documentation of this file.

#ifndef HARPAHAM_H
#define HARPAHAM_H

#include "harpavec.h"


class HaRPAHam 
{
public:
        HaRPAHam();
        virtual ~HaRPAHam();

        bool SetOpMode(int imode); 
        bool SetEnergy(const double NewEne); 

        bool LoadGauCom(const HaQCMod & qc_mod) const; 


bool Apply(vector<HaRPAvec> & RPAvec); 


bool Apply_init(vector<HaRPAvec> & RPAvec) const; 
      // Apply zero order RPA Hamiltonian 


vector<HaRPAvec> 
operator *(const vector<HaRPAvec> & RPAvec); 

protected: 
        double Ene;
        enum OperMode {FULL=0,MO_DIAG=1,E_MIN_H=2,E_MIN_H0=3} opmode; 
};


class HaRPAResolv 
{
public:
        HaRPAResolv();
        virtual ~HaRPAResolv();


bool SetOpMode(const int imode);

bool SetEnergy(const double NewEne); 

bool SetImag(const bool new_imag);


bool Apply_Init(vector<HaRPAvec> & RPAvec) const; 


bool Apply(vector<HaRPAvec> & RPAvec) const; 


vector<HaRPAvec> solve(const vector<HaRPAvec> & CISvec) const; 


protected: 
        bool imag;    
        double Ene;   
        enum OperMode {FULL=0, MO_DIAG} opmode; 
};



#endif /* !HARPAHAM_H */

---