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

Public Member Functions

 HarmonicExternalPotential (ParticleSet &P)
 
 ~HarmonicExternalPotential () override
 
std::string getClassName () const override
 return class name More...
 
void resetTargetParticleSet (ParticleSet &P) override
 Reset the data with the target ParticleSet. More...
 
bool put (xmlNodePtr cur) override
 Read the input parameter. More...
 
bool get (std::ostream &os) const override
 write about the class More...
 
std::unique_ptr< OperatorBasemakeClone (ParticleSet &P, TrialWaveFunction &psi) final
 
Return_t evaluate (ParticleSet &P) override
 Evaluate the local energy contribution of this component. More...
 
Return_t evaluate (ParticleSet &P, std::vector< NonLocalData > &Txy)
 
void contributeParticleQuantities () override
 
void checkoutParticleQuantities (TraceManager &tm) override
 
void deleteParticleQuantities () override
 
Return_t evaluate_sp (ParticleSet &P)
 
- 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 addObservables (PropertySetType &plist, BufferType &collectables)
 named values to the property list Default implementaton uses addValue(plist_) More...
 
virtual void registerObservables (std::vector< ObservableHelper > &h5desc, hdf_archive &file) const
 add to observable descriptor for hdf5 The default implementation is to register a scalar for this->value_ More...
 
virtual void registerCollectables (std::vector< ObservableHelper > &h5desc, hdf_archive &file) const
 
virtual void setObservables (PropertySetType &plist)
 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...
 
virtual void setParticlePropertyList (PropertySetType &plist, int offset)
 
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 Attributes

RealType mass
 
RealType energy
 
RealType length
 
PosType center
 
const ParticleSetPs
 
Array< TraceReal, 1 > * V_sample
 single particle trace sample array More...
 

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 > >
 
- Protected Member Functions inherited from OperatorBase
virtual void contributeScalarQuantities ()
 
virtual void checkoutScalarQuantities (TraceManager &tm)
 
virtual void collectScalarQuantities ()
 
virtual void deleteScalarQuantities ()
 
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_
 

Detailed Description

Definition at line 22 of file HarmonicExternalPotential.h.

Constructor & Destructor Documentation

◆ HarmonicExternalPotential()

Definition at line 37 of file HarmonicExternalPotential.h.

References OperatorBase::oneBodyQuantumDomain(), OperatorBase::POTENTIAL, and OperatorBase::setEnergyDomain().

37  : Ps(P)
38  {
41  }
void oneBodyQuantumDomain(const ParticleSet &P)
set quantum domain for one-body operator
void setEnergyDomain(EnergyDomains edomain)
Set the Energy Domain.

◆ ~HarmonicExternalPotential()

~HarmonicExternalPotential ( )
inlineoverride

Definition at line 43 of file HarmonicExternalPotential.h.

43 {}

Member Function Documentation

◆ checkoutParticleQuantities()

void checkoutParticleQuantities ( TraceManager tm)
inlineoverridevirtual

Reimplemented from OperatorBase.

Definition at line 62 of file HarmonicExternalPotential.h.

References TraceManager::checkout_real(), OperatorBase::name_, HarmonicExternalPotential::Ps, OperatorBase::request_, TraceRequest::streaming_array(), OperatorBase::streaming_particles_, and HarmonicExternalPotential::V_sample.

63  {
66  V_sample = tm.checkout_real<1>(name_, Ps);
67  }
bool streaming_array(const std::string &name)
Definition: TraceManager.h:249
TraceRequest request_
whether traces are being collected
Definition: OperatorBase.h:531
Array< TraceReal, 1 > * V_sample
single particle trace sample array
std::string name_
name of this object
Definition: OperatorBase.h:527

◆ contributeParticleQuantities()

void contributeParticleQuantities ( )
inlineoverridevirtual

Reimplemented from OperatorBase.

Definition at line 60 of file HarmonicExternalPotential.h.

References TraceRequest::contribute_array(), OperatorBase::name_, and OperatorBase::request_.

TraceRequest request_
whether traces are being collected
Definition: OperatorBase.h:531
std::string name_
name of this object
Definition: OperatorBase.h:527
void contribute_array(const std::string &name, bool default_quantity=false)
Definition: TraceManager.h:198

◆ deleteParticleQuantities()

void deleteParticleQuantities ( )
inlineoverridevirtual

Reimplemented from OperatorBase.

Definition at line 69 of file HarmonicExternalPotential.h.

References OperatorBase::streaming_particles_, and HarmonicExternalPotential::V_sample.

70  {
72  delete V_sample;
73  }
Array< TraceReal, 1 > * V_sample
single particle trace sample array

◆ evaluate() [1/2]

HarmonicExternalPotential::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 63 of file HarmonicExternalPotential.cpp.

References HarmonicExternalPotential::center, qmcplusplus::dot(), HarmonicExternalPotential::energy, HarmonicExternalPotential::evaluate_sp(), ParticleSet::getTotalNum(), HarmonicExternalPotential::length, ParticleSet::R, OperatorBase::streaming_particles_, and OperatorBase::value_.

64 {
65 #if !defined(REMOVE_TRACEMANAGER)
67  value_ = evaluate_sp(P);
68  else
69  {
70 #endif
71  value_ = 0.0;
72  RealType prefactor = .5 * energy / (length * length);
73  for (int i = 0; i < P.getTotalNum(); ++i)
74  {
75  PosType r = P.R[i] - center;
76  value_ += prefactor * dot(r, r);
77  }
78 #if !defined(REMOVE_TRACEMANAGER)
79  }
80 #endif
81  return value_;
82 }
QMCTraits::PosType PosType
QMCTraits::RealType RealType
Return_t value_
current value
Definition: OperatorBase.h:524
Tensor< typename BinaryReturn< T1, T2, OpMultiply >::Type_t, D > dot(const AntiSymTensor< T1, D > &lhs, const AntiSymTensor< T2, D > &rhs)

◆ evaluate() [2/2]

Return_t evaluate ( ParticleSet P,
std::vector< NonLocalData > &  Txy 
)
inline

Definition at line 56 of file HarmonicExternalPotential.h.

References HarmonicExternalPotential::evaluate().

Referenced by HarmonicExternalPotential::evaluate().

56 { return evaluate(P); }
Return_t evaluate(ParticleSet &P) override
Evaluate the local energy contribution of this component.

◆ evaluate_sp()

HarmonicExternalPotential::Return_t evaluate_sp ( ParticleSet P)
inline

Definition at line 86 of file HarmonicExternalPotential.cpp.

References HarmonicExternalPotential::center, qmcplusplus::dot(), HarmonicExternalPotential::energy, ParticleSet::getTotalNum(), HarmonicExternalPotential::length, ParticleSet::R, HarmonicExternalPotential::V_sample, and OperatorBase::value_.

Referenced by HarmonicExternalPotential::evaluate().

87 {
88  Array<TraceReal, 1>& V_samp = *V_sample;
89  value_ = 0.0;
90  RealType prefactor = .5 * energy / (length * length);
91  for (int i = 0; i < P.getTotalNum(); ++i)
92  {
93  PosType r = P.R[i] - center;
94  RealType v1 = prefactor * dot(r, r);
95  V_samp(i) = v1;
96  value_ += v1;
97  }
98  return value_;
99 }
QMCTraits::PosType PosType
Array< TraceReal, 1 > * V_sample
single particle trace sample array
QMCTraits::RealType RealType
Return_t value_
current value
Definition: OperatorBase.h:524
Tensor< typename BinaryReturn< T1, T2, OpMultiply >::Type_t, D > dot(const AntiSymTensor< T1, D > &lhs, const AntiSymTensor< T2, D > &rhs)

◆ get()

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

write about the class

Implements OperatorBase.

Definition at line 50 of file HarmonicExternalPotential.cpp.

51 {
52  os << "External harmonic potential" << std::endl;
53  return true;
54 }

◆ getClassName()

std::string getClassName ( ) const
inlineoverridevirtual

return class name

Implements OperatorBase.

Definition at line 45 of file HarmonicExternalPotential.h.

45 { return "HarmonicExternalPotential"; }

◆ makeClone()

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

Implements OperatorBase.

Definition at line 57 of file HarmonicExternalPotential.cpp.

58 {
59  return std::make_unique<HarmonicExternalPotential>(*this);
60 }

◆ put()

bool put ( xmlNodePtr  cur)
overridevirtual

Read the input parameter.

Parameters
curxml node for a OperatorBase object

Implements OperatorBase.

Definition at line 20 of file HarmonicExternalPotential.cpp.

References OhmmsAttributeSet::add(), HarmonicExternalPotential::center, HarmonicExternalPotential::energy, HarmonicExternalPotential::length, HarmonicExternalPotential::mass, OhmmsAttributeSet::put(), and qmcplusplus::sqrt().

21 {
22  using std::sqrt;
23 
24  mass = -1.0;
25  energy = -1.0;
26  length = -1.0;
27  center = 0.0;
28 
29  OhmmsAttributeSet attrib;
30  attrib.add(mass, "mass");
31  attrib.add(energy, "frequency");
32  attrib.add(energy, "energy");
33  attrib.add(length, "length");
34  attrib.add(center, "center");
35  attrib.put(cur);
36 
37  if (energy < 0.0)
38  energy = 1.0;
39  if (mass < 0.0 && length < 0.0)
40  length = 1.0;
41  if (mass < 0.0)
42  mass = 1.0 / (energy * length * length);
43  else if (length < 0.0)
44  length = 1.0 / sqrt(mass * energy);
45 
46  return true;
47 }
bool put(xmlNodePtr cur)
assign attributes to the set
Definition: AttributeSet.h:55
class to handle a set of attributes of an xmlNode
Definition: AttributeSet.h:24
MakeReturn< UnaryNode< FnSqrt, typename CreateLeaf< Vector< T1, C1 > >::Leaf_t > >::Expression_t sqrt(const Vector< T1, C1 > &l)
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

◆ resetTargetParticleSet()

void resetTargetParticleSet ( ParticleSet P)
inlineoverridevirtual

Reset the data with the target ParticleSet.

Parameters
Pnew target ParticleSet

Implements OperatorBase.

Definition at line 47 of file HarmonicExternalPotential.h.

47 {}

Member Data Documentation

◆ center

◆ energy

◆ length

◆ mass

RealType mass

Definition at line 25 of file HarmonicExternalPotential.h.

Referenced by HarmonicExternalPotential::put().

◆ Ps

◆ V_sample


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