QMCPACK
ForceCeperley Struct Reference
+ Inheritance diagram for ForceCeperley:
+ Collaboration diagram for ForceCeperley:

Public Member Functions

 ForceCeperley (ParticleSet &ions, ParticleSet &elns)
 
std::string getClassName () const override
 return class name More...
 
Return_t evaluate (ParticleSet &P) override
 Evaluate the local energy contribution of this component. More...
 
void InitMatrix ()
 
void registerObservables (std::vector< ObservableHelper > &h5list, hdf_archive &file) const override
 add to observable descriptor for hdf5 The default implementation is to register a scalar for this->value_ More...
 
void addObservables (PropertySetType &plist, BufferType &collectables) override
 named values to the property list Default implementaton uses addValue(plist_) More...
 
void setObservables (PropertySetType &plist) override
 Set the values evaluated by this object to plist Default implementation is to assign Value which is updated by evaluate function using my_index_. More...
 
void resetTargetParticleSet (ParticleSet &P) override
 Reset the data with the target ParticleSet. More...
 
void evaluate_IonIon (ParticleSet::ParticlePos &forces) const
 
void setParticlePropertyList (PropertySetType &plist, int offset) override
 
std::unique_ptr< OperatorBasemakeClone (ParticleSet &qp, TrialWaveFunction &psi) final
 
bool put (xmlNodePtr cur) override
 Read the input parameter. More...
 
bool get (std::ostream &os) const override
 write about the class More...
 
- Public Member Functions inherited from OperatorBase
 OperatorBase ()
 Construct a new Operator Base object Default and unique empty constructor. More...
 
virtual ~OperatorBase ()=default
 
virtual bool dependsOnWaveFunction () const
 return true if this operator depends on a wavefunction More...
 
std::bitset< 8 > & getUpdateMode () noexcept
 get update_mode_ reference More...
 
Return_t getValue () const noexcept
 get a copy of value_ More...
 
std::string getName () const noexcept
 getter a copy of my_name_, rvalue small string optimization More...
 
void setName (const std::string name) noexcept
 Set my_name member, uses small string optimization (pass by value) More...
 
TraceRequestgetRequest () noexcept
 Get request_ member. More...
 
virtual void registerCollectables (std::vector< ObservableHelper > &h5desc, hdf_archive &file) const
 
virtual void setHistories (Walker_t &ThisWalker)
 
virtual Return_t evaluateDeterministic (ParticleSet &P)
 Evaluate the local energy contribution of this component, deterministically based on current state. More...
 
virtual void mw_evaluate (const RefVectorWithLeader< OperatorBase > &o_list, const RefVectorWithLeader< TrialWaveFunction > &wf_list, const RefVectorWithLeader< ParticleSet > &p_list) const
 Evaluate the contribution of this component of multiple walkers. More...
 
virtual void mw_evaluatePerParticle (const RefVectorWithLeader< OperatorBase > &o_list, const RefVectorWithLeader< TrialWaveFunction > &wf_list, const RefVectorWithLeader< ParticleSet > &p_list, const std::vector< ListenerVector< RealType >> &listeners, const std::vector< ListenerVector< RealType >> &listeners_ions) const
 Evaluate the contribution of this component of multiple walkers per particle and report to registerd listeners from objects in Estimators. More...
 
virtual void mw_evaluateWithParameterDerivatives (const RefVectorWithLeader< OperatorBase > &o_list, const RefVectorWithLeader< ParticleSet > &p_list, const opt_variables_type &optvars, const RecordArray< ValueType > &dlogpsi, RecordArray< ValueType > &dhpsioverpsi) const
 TODO: add docs. More...
 
virtual Return_t rejectedMove (ParticleSet &P)
 TODO: add docs. More...
 
virtual Return_t evaluateWithToperator (ParticleSet &P)
 Evaluate the local energy contribution of this component with Toperators updated if requested. More...
 
virtual void mw_evaluateWithToperator (const RefVectorWithLeader< OperatorBase > &o_list, const RefVectorWithLeader< TrialWaveFunction > &wf_list, const RefVectorWithLeader< ParticleSet > &p_list) const
 Evaluate the contribution of this component of multiple walkers. More...
 
virtual void mw_evaluatePerParticleWithToperator (const RefVectorWithLeader< OperatorBase > &o_list, const RefVectorWithLeader< TrialWaveFunction > &wf_list, const RefVectorWithLeader< ParticleSet > &p_list, const std::vector< ListenerVector< RealType >> &listeners, const std::vector< ListenerVector< RealType >> &listeners_ions) const
 Evaluate the contribution of this component of multiple walkers per particle and report to registerd listeners from objects in Estimators. More...
 
virtual Return_t evaluateValueAndDerivatives (ParticleSet &P, const opt_variables_type &optvars, const Vector< ValueType > &dlogpsi, Vector< ValueType > &dhpsioverpsi)
 Evaluate value and derivatives wrt the optimizables. More...
 
virtual Return_t evaluateWithIonDerivs (ParticleSet &P, ParticleSet &ions, TrialWaveFunction &psi, ParticleSet::ParticlePos &hf_term, ParticleSet::ParticlePos &pulay_term)
 Evaluate contribution to local energy and derivatives w.r.t ionic coordinates from OperatorBase. More...
 
virtual Return_t evaluateWithIonDerivsDeterministic (ParticleSet &P, ParticleSet &ions, TrialWaveFunction &psi, ParticleSet::ParticlePos &hf_term, ParticleSet::ParticlePos &pulay_term)
 Evaluate contribution to local energy and derivatives w.r.t ionic coordinates from OperatorBase. More...
 
virtual void evaluateOneBodyOpMatrix (ParticleSet &P, const TWFFastDerivWrapper &psi, std::vector< ValueMatrix > &B)
 Evaluate "B" matrix for observable. More...
 
virtual void evaluateOneBodyOpMatrixForceDeriv (ParticleSet &P, ParticleSet &source, const TWFFastDerivWrapper &psi, const int iat, std::vector< std::vector< ValueMatrix >> &Bforce)
 Evaluate "dB/dR" matrices for observable. More...
 
virtual void updateSource (ParticleSet &s)
 Update data associated with a particleset. More...
 
virtual Return_t getEnsembleAverage ()
 Return an average value by collective operation. More...
 
virtual void createResource (ResourceCollection &collection) const
 Initialize a shared resource and hand it to a collection. More...
 
virtual void acquireResource (ResourceCollection &collection, const RefVectorWithLeader< OperatorBase > &o_list) const
 Acquire a shared resource from a collection. More...
 
virtual void releaseResource (ResourceCollection &collection, const RefVectorWithLeader< OperatorBase > &o_list) const
 Return a shared resource to a collection. More...
 
virtual void setRandomGenerator (RandomBase< FullPrecRealType > *rng)
 Set the Random Generator object TODO: add docs. More...
 
virtual void add2Hamiltonian (ParticleSet &qp, TrialWaveFunction &psi, QMCHamiltonian &targetH)
 TODO: add docs. More...
 
virtual void getRequiredTraces (TraceManager &tm)
 TODO: add docs. More...
 
virtual void informOfPerParticleListener ()
 
bool isClassical () const noexcept
 
bool isQuantum () const noexcept
 
bool isClassicalClassical () const noexcept
 
bool isQuantumClassical () const noexcept
 
bool isQuantumQuantum () const noexcept
 
bool getMode (const int i) const noexcept
 Return the mode i. More...
 
bool isNonLocal () const noexcept
 TODO: add docs. More...
 
bool hasListener () const noexcept
 
void contributeTraceQuantities ()
 Make trace quantities available. More...
 
void checkoutTraceQuantities (TraceManager &tm)
 Checkout trace arrays Derived classes must guard individual checkouts using request info. More...
 
void collectScalarTraces ()
 Collect scalar trace data. More...
 
void deleteTraceQuantities ()
 delete trace arrays More...
 
- Public Member Functions inherited from ForceBase
Real g (Real r)
 
void initVarReduction (Real rcut, int m, int numFuncs)
 
 ForceBase (ParticleSet &ions, ParticleSet &elns)
 
virtual ~ForceBase ()
 
void registerObservablesF (std::vector< ObservableHelper > &h5list, hdf_archive &file) const
 
void addObservablesF (QMCTraits::PropertySetType &plist)
 
void addObservablesStress (QMCTraits::PropertySetType &plist)
 
void setObservablesF (QMCTraits::PropertySetType &plist)
 
void setObservablesStress (QMCTraits::PropertySetType &plist)
 
void setParticleSetF (QMCTraits::PropertySetType &plist, int offset)
 
void setParticleSetStress (QMCTraits::PropertySetType &plist, int offset)
 
bool getAddIonIon () const noexcept
 
void setAddIonIon (bool val) noexcept
 
const ParticleSet::ParticlePosgetForces () const noexcept
 
void setForces (const ParticleSet::ParticlePos &forces)
 
void setForces (Real val)
 
const ParticleSet::ParticlePosgetForcesIonIon () const noexcept
 
void setForcesIonIon (const ParticleSet::ParticlePos &forces_ion_ion)
 
const SymTensor< Real, OHMMS_DIM > & getStressIonIon () const noexcept
 
const SymTensor< Real, OHMMS_DIM > & getStressEE () const noexcept
 
const SymTensor< Real, OHMMS_DIM > & getStressEI () const noexcept
 
const SymTensor< Real, OHMMS_DIM > & getStressKin () const noexcept
 
const SymTensor< Real, OHMMS_DIM > & getStress () const noexcept
 

Public Attributes

double Rcut
 
int m_exp
 
int N_basis
 
Matrix< FullPrecRealTypeSinv
 
Vector< FullPrecRealTypeh
 
Vector< FullPrecRealTypec
 

Private Attributes

const int d_aa_ID
 
const int d_ei_ID
 

Additional Inherited Members

- Public Types inherited from OperatorBase
enum  EnergyDomains { KINETIC = 0, POTENTIAL, NO_ENERGY_DOMAIN }
 enum to denote energy domain of operators More...
 
enum  QuantumDomains {
  NO_QUANTUM_DOMAIN = 0, CLASSICAL, QUANTUM, CLASSICAL_CLASSICAL,
  QUANTUM_CLASSICAL, QUANTUM_QUANTUM
}
 
enum  {
  PRIMARY = 0, OPTIMIZABLE = 1, RATIOUPDATE = 2, PHYSICAL = 3,
  COLLECTABLE = 4, NONLOCAL = 5
}
 enum for update_mode More...
 
using Return_t = FullPrecRealType
 type of return value of evaluate More...
 
using ValueMatrix = SPOSet::ValueMatrix
 For fast derivative evaluation. More...
 
using GradMatrix = SPOSet::GradMatrix
 
using BufferType = ParticleSet::Buffer_t
 typedef for the serialized buffer More...
 
using Walker_t = ParticleSet::Walker_t
 typedef for the walker More...
 
using ParticleScalar = ParticleSet::Scalar_t
 typedef for the ParticleScalar More...
 
using SPOMap = SPOSet::SPOMap
 typedef for SPOMap More...
 
- Public Types inherited from QMCTraits
enum  { DIM = OHMMS_DIM, DIM_VGL = OHMMS_DIM + 2 }
 
using QTBase = QMCTypes< OHMMS_PRECISION, DIM >
 
using QTFull = QMCTypes< OHMMS_PRECISION_FULL, DIM >
 
using RealType = QTBase::RealType
 
using ComplexType = QTBase::ComplexType
 
using ValueType = QTBase::ValueType
 
using PosType = QTBase::PosType
 
using GradType = QTBase::GradType
 
using TensorType = QTBase::TensorType
 
using IndexType = OHMMS_INDEXTYPE
 define other types More...
 
using FullPrecRealType = QTFull::RealType
 
using FullPrecValueType = QTFull::ValueType
 
using PropertySetType = RecordNamedProperty< FullPrecRealType >
 define PropertyList_t More...
 
using PtclGrpIndexes = std::vector< std::pair< int, int > >
 
- Public Types inherited from ForceBase
using Real = QMCTraits::RealType
 cheat, need to use virtual inheriance to clean up More...
 
- Protected Member Functions inherited from OperatorBase
virtual void contributeScalarQuantities ()
 
virtual void checkoutScalarQuantities (TraceManager &tm)
 
virtual void collectScalarQuantities ()
 
virtual void deleteScalarQuantities ()
 
virtual void contributeParticleQuantities ()
 
virtual void checkoutParticleQuantities (TraceManager &tm)
 
virtual void deleteParticleQuantities ()
 
virtual void setComputeForces (bool compute)
 
void setEnergyDomain (EnergyDomains edomain)
 Set the Energy Domain. More...
 
void setQuantumDomain (QuantumDomains qdomain)
 set quantum domain More...
 
void oneBodyQuantumDomain (const ParticleSet &P)
 set quantum domain for one-body operator More...
 
void twoBodyQuantumDomain (const ParticleSet &P)
 set quantum domain for two-body operator More...
 
void twoBodyQuantumDomain (const ParticleSet &P1, const ParticleSet &P2)
 set quantum domain for two-body operator More...
 
void addValue (PropertySetType &plist)
 named values to the property list More...
 
- Protected Attributes inherited from OperatorBase
std::bitset< 8 > update_mode_
 set the current update mode More...
 
Return_t value_
 current value More...
 
std::string name_
 name of this object More...
 
TraceRequest request_
 whether traces are being collected More...
 
int my_index_
 starting index of this object More...
 
Return_t new_value_
 a new value for a proposed move More...
 
Walker_tt_walker_
 reference to the current walker More...
 
bool streaming_particles_
 
bool have_required_traces_
 
- Protected Attributes inherited from ForceBase
int first_force_index_
 
int n_nuc_
 
int n_el_
 
int tries_
 
bool first_time_
 
bool add_ion_ion_
 Determines if ion-ion force will be added to electron-ion force in derived force estimators. If false, forces_ion_ion_=0.0. More...
 
ParticleSetions_
 
ParticleSet::ParticlePos forces_
 
ParticleSet::ParticlePos forces_ion_ion_
 
SymTensor< Real, OHMMS_DIMstress_ion_ion_
 
SymTensor< Real, OHMMS_DIMstress_ee_
 
SymTensor< Real, OHMMS_DIMstress_ei_
 
SymTensor< Real, OHMMS_DIMstress_kin_
 
SymTensor< Real, OHMMS_DIMstress_
 
std::string prefix_
 
std::string pair_name_
 
Real rcut_
 
int m_
 
std::vector< Realck_
 

Detailed Description

Definition at line 28 of file ForceCeperley.h.

Constructor & Destructor Documentation

◆ ForceCeperley()

ForceCeperley ( ParticleSet ions,
ParticleSet elns 
)

Definition at line 28 of file ForceCeperley.cpp.

References ForceCeperley::evaluate_IonIon(), ForceBase::forces_, ForceBase::forces_ion_ion_, ForceCeperley::m_exp, ForceCeperley::N_basis, OperatorBase::name_, ForceBase::prefix_, ForceCeperley::Rcut, and ParticleSet::update().

29  : ForceBase(ions, elns), d_aa_ID(ions.addTable(ions)), d_ei_ID(elns.addTable(ions))
30 {
31  ReportEngine PRE("ForceCeperley", "ForceCeperley");
32  name_ = "Ceperley_Force_Base";
33  prefix_ = "HFCep";
34  // Defaults
35  Rcut = 0.4;
36  m_exp = 2;
37  N_basis = 4;
38  forces_ = 0.0;
39  ///////////////////////////////////////////////////////////////
40  ions.update();
42 }
ParticleSet::ParticlePos forces_
Definition: ForceBase.h:87
ForceBase(ParticleSet &ions, ParticleSet &elns)
Definition: ForceBase.cpp:28
std::string prefix_
Definition: ForceBase.h:96
std::string name_
name of this object
Definition: OperatorBase.h:527
void evaluate_IonIon(ParticleSet::ParticlePos &forces) const
ParticleSet::ParticlePos forces_ion_ion_
Definition: ForceBase.h:88

Member Function Documentation

◆ addObservables()

void addObservables ( PropertySetType plist,
BufferType collectables 
)
inlineoverridevirtual

named values to the property list Default implementaton uses addValue(plist_)

Parameters
plistRecordNameProperty
collectablesObservables that are accumulated by evaluate

Reimplemented from OperatorBase.

Definition at line 56 of file ForceCeperley.h.

References ForceBase::addObservablesF().

56 { addObservablesF(plist); }
void addObservablesF(QMCTraits::PropertySetType &plist)
Definition: ForceBase.cpp:47

◆ evaluate()

ForceCeperley::Return_t evaluate ( ParticleSet P)
overridevirtual

Evaluate the local energy contribution of this component.

Parameters
Pinput configuration containing N particles
Returns
the value of the Hamiltonian component

Implements OperatorBase.

Definition at line 82 of file ForceCeperley.cpp.

References ForceBase::add_ion_ion_, ForceCeperley::c, ForceCeperley::d_ei_ID, ForceBase::forces_, ForceBase::forces_ion_ion_, ParticleSet::getDistTableAB(), ForceBase::ions_, ForceCeperley::m_exp, ForceCeperley::N_basis, ForceBase::n_el_, ForceBase::n_nuc_, qmcplusplus::pow(), ForceCeperley::Rcut, and ParticleSet::Z.

Referenced by qmcplusplus::TEST_CASE().

83 {
84  if (add_ion_ion_ == true)
86  else
87  forces_ = 0.0;
88  const auto& d_ab = P.getDistTableAB(d_ei_ID);
89  const ParticleScalar* restrict Zat = ions_.Z.first_address();
90  const ParticleScalar* restrict Qat = P.Z.first_address();
91  for (int jat = 0; jat < n_el_; jat++)
92  {
93  const auto& dist = d_ab.getDistRow(jat);
94  const auto& displ = d_ab.getDisplRow(jat);
95  for (int iat = 0; iat < n_nuc_; iat++)
96  {
97  // electron contribution (special treatment if distance is inside cutoff!)
98  RealType r = dist[iat];
99  RealType zoverr3 = Qat[jat] * Zat[iat] / (r * r * r);
100  if (r >= Rcut)
101  {
102  forces_[iat] += zoverr3 * displ[iat];
103  }
104  else
105  {
106  RealType g_q = 0.0;
107  for (int q = 0; q < N_basis; q++)
108  {
109  g_q += c[q] * std::pow(r, m_exp + q + 1);
110  }
111  g_q *= zoverr3;
112  // negative sign accounts for definition of target as electrons
113  forces_[iat] += g_q * displ[iat];
114  }
115  }
116  }
117  return 0.0;
118 }
ParticleSet::ParticlePos forces_
Definition: ForceBase.h:87
MakeReturn< BinaryNode< FnPow, typename CreateLeaf< Vector< T1, C1 > >::Leaf_t, typename CreateLeaf< Vector< T2, C2 > >::Leaf_t > >::Expression_t pow(const Vector< T1, C1 > &l, const Vector< T2, C2 > &r)
bool add_ion_ion_
Determines if ion-ion force will be added to electron-ion force in derived force estimators. If false, forces_ion_ion_=0.0.
Definition: ForceBase.h:84
ParticleScalar Z
charge of each particle
Definition: ParticleSet.h:89
ParticleSet::Scalar_t ParticleScalar
typedef for the ParticleScalar
Definition: OperatorBase.h:81
QMCTraits::RealType RealType
ParticleSet::ParticlePos forces_ion_ion_
Definition: ForceBase.h:88
Vector< FullPrecRealType > c
Definition: ForceCeperley.h:40
ParticleSet & ions_
Definition: ForceBase.h:86

◆ evaluate_IonIon()

void evaluate_IonIon ( ParticleSet::ParticlePos forces) const

Definition at line 44 of file ForceCeperley.cpp.

References ForceCeperley::d_aa_ID, ParticleSet::getDistTableAA(), ForceBase::ions_, ForceBase::n_nuc_, and ParticleSet::Z.

Referenced by ForceCeperley::ForceCeperley().

45 {
46  forces = 0.0;
47  const auto& d_aa(ions_.getDistTableAA(d_aa_ID));
48  const ParticleScalar* restrict Zat = ions_.Z.first_address();
49  for (size_t ipart = 1; ipart < n_nuc_; ipart++)
50  {
51  const auto& dist = d_aa.getDistRow(ipart);
52  const auto& displ = d_aa.getDisplRow(ipart);
53  for (size_t jpart = 0; jpart < ipart; ++jpart)
54  {
55  RealType rinv = 1.0 / dist[jpart];
56  RealType r3zz = Zat[jpart] * Zat[ipart] * rinv * rinv * rinv;
57  forces[jpart] += r3zz * displ[jpart];
58  forces[ipart] -= r3zz * displ[jpart];
59  }
60  }
61 }
const DistanceTableAA & getDistTableAA(int table_ID) const
get a distance table by table_ID and dyanmic_cast to DistanceTableAA
ParticleScalar Z
charge of each particle
Definition: ParticleSet.h:89
ParticleSet::Scalar_t ParticleScalar
typedef for the ParticleScalar
Definition: OperatorBase.h:81
QMCTraits::RealType RealType
ParticleSet & ions_
Definition: ForceBase.h:86

◆ get()

bool get ( std::ostream &  os) const
inlineoverridevirtual

write about the class

Implements OperatorBase.

Definition at line 70 of file ForceCeperley.h.

References ForceBase::pair_name_.

71  {
72  os << "Ceperley Force Estimator Hamiltonian: " << pair_name_;
73  return true;
74  }
std::string pair_name_
Definition: ForceBase.h:97

◆ getClassName()

std::string getClassName ( ) const
inlineoverridevirtual

return class name

Implements OperatorBase.

Definition at line 45 of file ForceCeperley.h.

45 { return "ForceCeperley"; }

◆ InitMatrix()

void InitMatrix ( )

Definition at line 63 of file ForceCeperley.cpp.

References ForceCeperley::c, ForceCeperley::h, qmcplusplus::invert_matrix(), ForceCeperley::m_exp, ForceCeperley::N_basis, qmcplusplus::pow(), qmcplusplus::MatrixOperators::product(), ForceCeperley::Rcut, Matrix< T, Alloc >::resize(), Vector< T, Alloc >::resize(), and ForceCeperley::Sinv.

Referenced by ForceCeperley::put(), and qmcplusplus::TEST_CASE().

64 {
66  h.resize(N_basis);
67  c.resize(N_basis);
68  for (int k = 0; k < N_basis; k++)
69  {
70  h[k] = std::pow(Rcut, (k + 2)) / static_cast<RealType>(k + 2);
71  for (int j = 0; j < N_basis; j++)
72  {
73  Sinv(k, j) = std::pow(Rcut, (m_exp + k + j + 3)) / static_cast<RealType>(m_exp + k + j + 3);
74  }
75  }
76  // in Numerics/DeterminantOperators.h
77  invert_matrix(Sinv, false);
78  // in Numerics/MatrixOperators.h
80 }
void resize(size_type n, Type_t val=Type_t())
Resize the container.
Definition: OhmmsVector.h:166
MatrixA::value_type invert_matrix(MatrixA &M, bool getdet=true)
invert a matrix
Matrix< FullPrecRealType > Sinv
Definition: ForceCeperley.h:38
void resize(size_type n, size_type m)
Resize the container.
Definition: OhmmsMatrix.h:99
MakeReturn< BinaryNode< FnPow, typename CreateLeaf< Vector< T1, C1 > >::Leaf_t, typename CreateLeaf< Vector< T2, C2 > >::Leaf_t > >::Expression_t pow(const Vector< T1, C1 > &l, const Vector< T2, C2 > &r)
void product(const Matrix< T > &A, const Matrix< T > &B, Matrix< T > &C)
static function to perform C=AB for real matrices
QMCTraits::RealType RealType
Vector< FullPrecRealType > h
Definition: ForceCeperley.h:39
Vector< FullPrecRealType > c
Definition: ForceCeperley.h:40

◆ makeClone()

std::unique_ptr< OperatorBase > makeClone ( ParticleSet qp,
TrialWaveFunction psi 
)
finalvirtual

Implements OperatorBase.

Definition at line 144 of file ForceCeperley.cpp.

145 {
146  return std::make_unique<ForceCeperley>(*this);
147 }

◆ put()

bool put ( xmlNodePtr  cur)
overridevirtual

Read the input parameter.

Parameters
curxml node for a OperatorBase object

Implements OperatorBase.

Definition at line 120 of file ForceCeperley.cpp.

References OhmmsAttributeSet::add(), ParameterSet::add(), ForceBase::add_ion_ion_, qmcplusplus::app_log(), ForceBase::first_time_, ForceCeperley::InitMatrix(), ForceCeperley::m_exp, ForceCeperley::N_basis, ForceBase::prefix_, ParameterSet::put(), OhmmsAttributeSet::put(), and ForceCeperley::Rcut.

121 {
122  std::string ionionforce("yes");
123  OhmmsAttributeSet attr;
124  attr.add(prefix_, "name");
125  attr.add(ionionforce, "add_ion_ion_");
126  attr.put(cur);
127  add_ion_ion_ = (ionionforce == "yes") || (ionionforce == "true");
128  app_log() << "ionionforce = " << ionionforce << std::endl;
129  app_log() << "add_ion_ion_=" << add_ion_ion_ << std::endl;
130  app_log() << "first_time_= " << first_time_ << std::endl;
131  ParameterSet fcep_param_set;
132  fcep_param_set.add(Rcut, "rcut");
133  fcep_param_set.add(N_basis, "nbasis");
134  fcep_param_set.add(m_exp, "weight_exp");
135  fcep_param_set.put(cur);
136  app_log() << " ForceCeperley Parameters" << std::endl;
137  app_log() << " ForceCeperley::Rcut=" << Rcut << std::endl;
138  app_log() << " ForceCeperley::N_basis=" << N_basis << std::endl;
139  app_log() << " ForceCeperley::m_exp=" << m_exp << std::endl;
140  InitMatrix();
141  return true;
142 }
std::ostream & app_log()
Definition: OutputManager.h:65
bool put(xmlNodePtr cur)
assign attributes to the set
Definition: AttributeSet.h:55
std::string prefix_
Definition: ForceBase.h:96
bool put(std::istream &is) override
read from std::istream
Definition: ParameterSet.h:42
class to handle a set of parameters
Definition: ParameterSet.h:27
bool add_ion_ion_
Determines if ion-ion force will be added to electron-ion force in derived force estimators. If false, forces_ion_ion_=0.0.
Definition: ForceBase.h:84
class to handle a set of attributes of an xmlNode
Definition: AttributeSet.h:24
void add(PDT &aparam, const std::string &aname_in, std::vector< PDT > candidate_values={}, TagStatus status=TagStatus::OPTIONAL)
add a new parameter corresponding to an xmlNode <parameter>
void add(PDT &aparam, const std::string &aname, std::vector< PDT > candidate_values={}, TagStatus status=TagStatus::OPTIONAL)
add a new attribute
Definition: AttributeSet.h:42

◆ registerObservables()

void registerObservables ( std::vector< ObservableHelper > &  h5desc,
hdf_archive file 
) const
inlineoverridevirtual

add to observable descriptor for hdf5 The default implementation is to register a scalar for this->value_

Parameters
h5desccontains a set of hdf5 descriptors for a scalar observable
gidhdf5 group to which the observables belong

Reimplemented from OperatorBase.

Definition at line 51 of file ForceCeperley.h.

References ForceBase::registerObservablesF().

52  {
53  registerObservablesF(h5list, file);
54  }
void registerObservablesF(std::vector< ObservableHelper > &h5list, hdf_archive &file) const
Definition: ForceBase.cpp:75

◆ resetTargetParticleSet()

void resetTargetParticleSet ( ParticleSet P)
inlineoverridevirtual

Reset the data with the target ParticleSet.

Parameters
Pnew target ParticleSet

Implements OperatorBase.

Definition at line 60 of file ForceCeperley.h.

60 {}

◆ setObservables()

void setObservables ( PropertySetType plist)
inlineoverridevirtual

Set the values evaluated by this object to plist Default implementation is to assign Value which is updated by evaluate function using my_index_.

Parameters
plistRecordNameProperty

Reimplemented from OperatorBase.

Definition at line 58 of file ForceCeperley.h.

References ForceBase::setObservablesF().

58 { setObservablesF(plist); }
void setObservablesF(QMCTraits::PropertySetType &plist)
Definition: ForceBase.cpp:86

◆ setParticlePropertyList()

void setParticlePropertyList ( PropertySetType plist,
int  offset 
)
inlineoverridevirtual

Reimplemented from OperatorBase.

Definition at line 65 of file ForceCeperley.h.

References ForceBase::setParticleSetF().

65 { setParticleSetF(plist, offset); }
void setParticleSetF(QMCTraits::PropertySetType &plist, int offset)
Definition: ForceBase.cpp:113

Member Data Documentation

◆ c

◆ d_aa_ID

const int d_aa_ID
private

Definition at line 31 of file ForceCeperley.h.

Referenced by ForceCeperley::evaluate_IonIon().

◆ d_ei_ID

const int d_ei_ID
private

Definition at line 32 of file ForceCeperley.h.

Referenced by ForceCeperley::evaluate().

◆ h

Definition at line 39 of file ForceCeperley.h.

Referenced by ForceCeperley::InitMatrix().

◆ m_exp

◆ N_basis

◆ Rcut

◆ Sinv

Definition at line 38 of file ForceCeperley.h.

Referenced by ForceCeperley::InitMatrix().


The documentation for this struct was generated from the following files: