H_Ewald_pw.h

Go to the documentation of this file.

#ifndef H_EWALD_PW_H
#define H_EWALD_PW_H
 
#include "source_base/global_function.h"
#include "source_base/global_variable.h"
#include "source_cell/unitcell.h"
#include "source_basis/module_pw/pw_basis.h"
#include "source_pw/module_pwdft/forces.h"
#include "source_pw/module_pwdft/stress_func.h"
 
class H_Ewald_pw 
{
  public:
    H_Ewald_pw();
    ~H_Ewald_pw();
 
    // compute the Ewald energy
    static double compute_ewald(const UnitCell& cell,
                                const ModulePW::PW_Basis* rho_basis,
                                const ModuleBase::ComplexMatrix& strucFac);
 
  public:
    static void rgen(
        const ModuleBase::Vector3<double> &dtau,
        const double &rmax,
        int *irr,
        const ModuleBase::Matrix3 &at,
        const ModuleBase::Matrix3 &bg,
        ModuleBase::Vector3<double> *r,
        double *r2,
        int  &nrm
    );
 
    // the coefficient of ewald method
    static double alpha;
    static int mxr;
 
};
 
#endif //ewald energy