psi_init_nao.h

Go to the documentation of this file.

#ifndef PSI_INIT_NAO_H
#define PSI_INIT_NAO_H
#include "source_base/cubic_spline.h"
#include "source_base/realarray.h"
#include "source_base/spherical_bessel_transformer.h"
#include "psi_initializer.h"
 
#include <memory>
/*
Psi (planewave based wavefunction) initializer: numerical atomic orbital method
*/
template <typename T>
class psi_init_nao : public psi_initializer<T>
{
  private:
    using Real = typename GetTypeReal<T>::type;
 
  public:
    psi_init_nao()
    {
        this->method_ = "nao";
    };
    ~psi_init_nao(){};
 
    virtual void init_psig(T* psig, const int& ik) override;
 
    virtual void initialize(const Structure_Factor*,             //< structure factor
                            const ModulePW::PW_Basis_K*,         //< planewave basis
                            const UnitCell*,                     //< unit cell
                            const K_Vectors*,                    //< kpoints
                            const int& = 1,                      //< random seed
                            const pseudopot_cell_vnl* = nullptr, //< nonlocal pseudopotential
                            const int& = 0) override;            //< MPI rank
 
    void read_external_orbs(const std::string* orbital_files, const int& rank);
    virtual void tabulate() override;
    std::vector<std::string> external_orbs() const
    {
        return orbital_files_;
    }
    std::vector<std::vector<int>> nr() const
    {
        return nr_;
    }
    std::vector<int> nr(const int& itype) const
    {
        return nr_[itype];
    }
    int nr(const int& itype, const int& ichi) const
    {
        return nr_[itype][ichi];
    }
    std::vector<std::vector<std::vector<double>>> chi() const
    {
        return chi_;
    }
    std::vector<std::vector<double>> chi(const int& itype) const
    {
        return chi_[itype];
    }
    std::vector<double> chi(const int& itype, const int& ichi) const
    {
        return chi_[itype][ichi];
    }
    double chi(const int& itype, const int& ichi, const int& ir) const
    {
        return chi_[itype][ichi][ir];
    }
    std::vector<std::vector<std::vector<double>>> rgrid() const
    {
        return rgrid_;
    }
    std::vector<std::vector<double>> rgrid(const int& itype) const
    {
        return rgrid_[itype];
    }
    std::vector<double> rgrid(const int& itype, const int& ichi) const
    {
        return rgrid_[itype][ichi];
    }
    double rgrid(const int& itype, const int& ichi, const int& ir) const
    {
        return rgrid_[itype][ichi][ir];
    }
 
  protected:
    void allocate_ao_table();
    std::vector<std::string> orbital_files_;
    std::unique_ptr<ModuleBase::CubicSpline> cubspl_;
    ModuleBase::realArray projmap_;
    std::vector<std::vector<int>> nr_;
    std::vector<std::vector<std::vector<double>>> chi_;
    std::vector<std::vector<std::vector<double>>> rgrid_;
    ModuleBase::SphericalBesselTransformer sbt;
};
#endif