pulay_force_stress.h

Go to the documentation of this file.

#pragma once
#include "source_basis/module_nao/two_center_bundle.h"
#include "source_estate/module_dm/density_matrix.h"
#include "source_lcao/module_gint/gint_gamma.h"
#include "source_lcao/module_gint/gint_k.h"
#include "source_estate/module_pot/potential_new.h"
#include "source_cell/unitcell.h"
#include "source_lcao/stress_tools.h"
#ifndef TGINT_H
#define TGINT_H
template <typename T>
struct TGint;
template <> struct TGint<double> { using type = Gint_Gamma; };
template <> struct TGint<std::complex<double>> { using type = Gint_k; };
#endif
 
namespace PulayForceStress
{
    template<typename TK, typename TR>
    void cal_pulay_fs(
        ModuleBase::matrix& f,  
        ModuleBase::matrix& s,  
        const elecstate::DensityMatrix<TK, TR>& dm,  
        const UnitCell& ucell,  
        const Parallel_Orbitals& pv,  
        const double* (&dHSx)[3],  
        const double* (&dHSxy)[6],  
        const bool& isforce,
        const bool& isstress,
        Record_adj* ra = nullptr,
        const double& factor_force = 1.0,
        const double& factor_stress = 1.0);
 
    template<typename TK, typename TR>
    void cal_pulay_fs(
        ModuleBase::matrix& f,  
        ModuleBase::matrix& s,  
        const elecstate::DensityMatrix<TK, TR>& dm,  
        const UnitCell& ucell,  
        const Parallel_Orbitals& pv,  
        const double* (&dHSx)[3],  
        const double* dtau,  
        const bool& isforce,
        const bool& isstress,
        Record_adj* ra = nullptr,
        const double& factor_force = 1.0,
        const double& factor_stress = 1.0);
 
    template<typename TK, typename TR>
    void cal_pulay_fs(
        ModuleBase::matrix& f,  
        ModuleBase::matrix& s,  
        const elecstate::DensityMatrix<TK, TR>& dm,  
        const UnitCell& ucell,  
        const elecstate::Potential* pot, 
        typename TGint<TK>::type& gint, 
        const bool& isforce,
        const bool& isstress,
        const bool& set_dmr_gint = true);
}
#include "pulay_force_stress_center2_template.hpp"
#include "pulay_force_stress_gint.hpp"