Inheritance diagram for FreeOrbital:

Collaboration diagram for FreeOrbital:

Public Member Functions
	FreeOrbital (const std::string &my_name, const std::vector< PosType > &kpts_cart)

	~FreeOrbital ()

std::string	getClassName () const override
	return class name More...

void	evaluate_notranspose (const ParticleSet &P, int first, int last, ValueMatrix &phi, GradMatrix &dphi, ValueMatrix &d2phi) override
	evaluate the values, gradients and laplacians of this single-particle orbital for [first,last) particles More...

void	evaluateVGL (const ParticleSet &P, int i, ValueVector &pvec, GradVector &dpvec, ValueVector &d2pvec) override
	evaluate the values, gradients and laplacians of this single-particle orbital set More...

void	evaluateValue (const ParticleSet &P, int iat, ValueVector &pvec) override
	evaluate the values of this single-particle orbital set More...

void	evaluate_notranspose (const ParticleSet &P, int first, int last, ValueMatrix &phi, GradMatrix &dphi, HessMatrix &d2phi_mat) override
	evaluate the values, gradients and hessians of this single-particle orbital for [first,last) particles More...

void	evaluate_notranspose (const ParticleSet &P, int first, int last, ValueMatrix &phi, GradMatrix &dphi, HessMatrix &d2phi_mat, GGGMatrix &d3phi_mat) override
	evaluate the values, gradients, hessians and third derivatives of this single-particle orbital for [first,last) particles More...

void	report (const std::string &pad) const override
	print SPOSet information More...

std::unique_ptr< SPOSet >	makeClone () const override
	make a clone of itself every derived class must implement this to have threading working correctly. More...

void	setOrbitalSetSize (int norbs) override
	set the OrbitalSetSize More...

Public Member Functions inherited from SPOSet
	SPOSet (const std::string &my_name)
	constructor More...

virtual	~SPOSet ()=default
	destructor More...

int	size () const
	return the size of the orbital set Ye: this needs to be replaced by getOrbitalSetSize(); More...

void	basic_report (const std::string &pad="") const
	print basic SPOSet information More...

int	getOrbitalSetSize () const
	return the size of the orbitals More...

virtual bool	isOptimizable () const
	Query if this SPOSet is optimizable. More...

virtual void	extractOptimizableObjectRefs (UniqueOptObjRefs &opt_obj_refs)
	extract underlying OptimizableObject references More...

virtual void	checkOutVariables (const opt_variables_type &active)
	check out variational optimizable variables More...

virtual bool	isOMPoffload () const
	Query if this SPOSet uses OpenMP offload. More...

virtual bool	hasIonDerivs () const
	Query if this SPOSet has an explicit ion dependence. More...

virtual void	checkObject () const
	check a few key parameters before putting the SPO into a determinant More...

virtual bool	isRotationSupported () const
	return true if this SPOSet can be wrappered by RotatedSPO More...

virtual void	storeParamsBeforeRotation ()
	store parameters before getting destroyed by rotation. More...

virtual void	applyRotation (const ValueMatrix &rot_mat, bool use_stored_copy=false)
	apply rotation to all the orbitals More...

virtual void	evaluateDerivatives (ParticleSet &P, const opt_variables_type &optvars, Vector< ValueType > &dlogpsi, Vector< ValueType > &dhpsioverpsi, const int &FirstIndex, const int &LastIndex)
	Parameter derivatives of the wavefunction and the Laplacian of the wavefunction. More...

virtual void	evaluateDerivativesWF (ParticleSet &P, const opt_variables_type &optvars, Vector< ValueType > &dlogpsi, int FirstIndex, int LastIndex)
	Parameter derivatives of the wavefunction. More...

virtual void	evaluateDerivatives (ParticleSet &P, const opt_variables_type &optvars, Vector< ValueType > &dlogpsi, Vector< ValueType > &dhpsioverpsi, const ValueType &psiCurrent, const std::vector< ValueType > &Coeff, const std::vector< size_t > &C2node_up, const std::vector< size_t > &C2node_dn, const ValueVector &detValues_up, const ValueVector &detValues_dn, const GradMatrix &grads_up, const GradMatrix &grads_dn, const ValueMatrix &lapls_up, const ValueMatrix &lapls_dn, const ValueMatrix &M_up, const ValueMatrix &M_dn, const ValueMatrix &Minv_up, const ValueMatrix &Minv_dn, const GradMatrix &B_grad, const ValueMatrix &B_lapl, const std::vector< int > &detData_up, const size_t N1, const size_t N2, const size_t NP1, const size_t NP2, const std::vector< std::vector< int >> &lookup_tbl)
	Evaluate the derivative of the optimized orbitals with respect to the parameters this is used only for MSD, to be refined for better serving both single and multi SD. More...

virtual void	evaluateDerivativesWF (ParticleSet &P, const opt_variables_type &optvars, Vector< ValueType > &dlogpsi, const QTFull::ValueType &psiCurrent, const std::vector< ValueType > &Coeff, const std::vector< size_t > &C2node_up, const std::vector< size_t > &C2node_dn, const ValueVector &detValues_up, const ValueVector &detValues_dn, const ValueMatrix &M_up, const ValueMatrix &M_dn, const ValueMatrix &Minv_up, const ValueMatrix &Minv_dn, const std::vector< int > &detData_up, const std::vector< std::vector< int >> &lookup_tbl)
	Evaluate the derivative of the optimized orbitals with respect to the parameters this is used only for MSD, to be refined for better serving both single and multi SD. More...

virtual void	evaluateDetRatios (const VirtualParticleSet &VP, ValueVector &psi, const ValueVector &psiinv, std::vector< ValueType > &ratios)
	evaluate determinant ratios for virtual moves, e.g., sphere move for nonlocalPP More...

virtual void	evaluateDetSpinorRatios (const VirtualParticleSet &VP, ValueVector &psi, const std::pair< ValueVector, ValueVector > &spinor_multiplier, const ValueVector &invrow, std::vector< ValueType > &ratios)
	evaluate determinant ratios for virtual moves, specifically for Spinor SPOSets More...

virtual void	evaluateDerivRatios (const VirtualParticleSet &VP, const opt_variables_type &optvars, ValueVector &psi, const ValueVector &psiinv, std::vector< ValueType > &ratios, Matrix< ValueType > &dratios, int FirstIndex, int LastIndex)
	Determinant ratios and parameter derivatives of the wavefunction for virtual moves. More...

virtual void	mw_evaluateDetRatios (const RefVectorWithLeader< SPOSet > &spo_list, const RefVectorWithLeader< const VirtualParticleSet > &vp_list, const RefVector< ValueVector > &psi_list, const std::vector< const ValueType *> &invRow_ptr_list, std::vector< std::vector< ValueType >> &ratios_list) const
	evaluate determinant ratios for virtual moves, e.g., sphere move for nonlocalPP, of multiple walkers More...

virtual void	evaluateVGL_spin (const ParticleSet &P, int iat, ValueVector &psi, GradVector &dpsi, ValueVector &d2psi, ValueVector &dspin)
	evaluate the values, gradients and laplacians and spin gradient of this single-particle orbital set More...

virtual void	mw_evaluateValue (const RefVectorWithLeader< SPOSet > &spo_list, const RefVectorWithLeader< ParticleSet > &P_list, int iat, const RefVector< ValueVector > &psi_v_list) const
	evaluate the values this single-particle orbital sets of multiple walkers More...

virtual void	mw_evaluateVGL (const RefVectorWithLeader< SPOSet > &spo_list, const RefVectorWithLeader< ParticleSet > &P_list, int iat, const RefVector< ValueVector > &psi_v_list, const RefVector< GradVector > &dpsi_v_list, const RefVector< ValueVector > &d2psi_v_list) const
	evaluate the values, gradients and laplacians of this single-particle orbital sets of multiple walkers More...

virtual void	mw_evaluateVGLWithSpin (const RefVectorWithLeader< SPOSet > &spo_list, const RefVectorWithLeader< ParticleSet > &P_list, int iat, const RefVector< ValueVector > &psi_v_list, const RefVector< GradVector > &dpsi_v_list, const RefVector< ValueVector > &d2psi_v_list, OffloadMatrix< ComplexType > &mw_dspin) const
	evaluate the values, gradients and laplacians and spin gradient of this single-particle orbital sets of multiple walkers More...

virtual void	mw_evaluateVGLandDetRatioGrads (const RefVectorWithLeader< SPOSet > &spo_list, const RefVectorWithLeader< ParticleSet > &P_list, int iat, const std::vector< const ValueType *> &invRow_ptr_list, OffloadMWVGLArray &phi_vgl_v, std::vector< ValueType > &ratios, std::vector< GradType > &grads) const
	evaluate the values, gradients and laplacians of this single-particle orbital sets and determinant ratio and grads of multiple walkers. More...

virtual void	mw_evaluateVGLandDetRatioGradsWithSpin (const RefVectorWithLeader< SPOSet > &spo_list, const RefVectorWithLeader< ParticleSet > &P_list, int iat, const std::vector< const ValueType *> &invRow_ptr_list, OffloadMWVGLArray &phi_vgl_v, std::vector< ValueType > &ratios, std::vector< GradType > &grads, std::vector< ValueType > &spingrads) const
	evaluate the values, gradients and laplacians of this single-particle orbital sets and determinant ratio and grads of multiple walkers. More...

virtual void	evaluateVGH (const ParticleSet &P, int iat, ValueVector &psi, GradVector &dpsi, HessVector &grad_grad_psi)
	evaluate the values, gradients and hessians of this single-particle orbital set More...

virtual void	evaluateVGHGH (const ParticleSet &P, int iat, ValueVector &psi, GradVector &dpsi, HessVector &grad_grad_psi, GGGVector &grad_grad_grad_psi)
	evaluate the values, gradients, hessians, and grad hessians of this single-particle orbital set More...

virtual void	evaluate_spin (const ParticleSet &P, int iat, ValueVector &psi, ValueVector &dpsi)
	evaluate the values of this single-particle orbital set More...

virtual void	evaluateThirdDeriv (const ParticleSet &P, int first, int last, GGGMatrix &grad_grad_grad_logdet)
	evaluate the third derivatives of this single-particle orbital set More...

virtual void	evaluate_notranspose_spin (const ParticleSet &P, int first, int last, ValueMatrix &logdet, GradMatrix &dlogdet, ValueMatrix &d2logdet, ValueMatrix &dspinlogdet)
	evaluate the values, gradients and laplacians of this single-particle orbital for [first,last) particles, including the spin gradient More...

virtual void	mw_evaluate_notranspose (const RefVectorWithLeader< SPOSet > &spo_list, const RefVectorWithLeader< ParticleSet > &P_list, int first, int last, const RefVector< ValueMatrix > &logdet_list, const RefVector< GradMatrix > &dlogdet_list, const RefVector< ValueMatrix > &d2logdet_list) const

virtual void	evaluateGradSource (const ParticleSet &P, int first, int last, const ParticleSet &source, int iat_src, GradMatrix &gradphi)
	evaluate the gradients of this single-particle orbital for [first,last) target particles with respect to the given source particle More...

virtual void	evaluateGradSource (const ParticleSet &P, int first, int last, const ParticleSet &source, int iat_src, GradMatrix &grad_phi, HessMatrix &grad_grad_phi, GradMatrix &grad_lapl_phi)
	evaluate the gradients of values, gradients, laplacians of this single-particle orbital for [first,last) target particles with respect to the given source particle More...

virtual void	evaluateGradSourceRow (const ParticleSet &P, int iel, const ParticleSet &source, int iat_src, GradVector &gradphi)
	Returns a row of d/dR_iat phi_j(r) evaluated at position r. More...

virtual PosType	get_k (int orb)
	access the k point related to the given orbital More...

virtual void	createResource (ResourceCollection &collection) const
	initialize a shared resource and hand it to collection More...

virtual void	acquireResource (ResourceCollection &collection, const RefVectorWithLeader< SPOSet > &spo_list) const
	acquire a shared resource from collection More...

virtual void	releaseResource (ResourceCollection &collection, const RefVectorWithLeader< SPOSet > &spo_list) const
	return a shared resource to collection More...

virtual bool	transformSPOSet ()
	Used only by cusp correction in AOS LCAO. More...

virtual void	finalizeConstruction ()
	finalize the construction of SPOSet More...

const std::string &	getName () const
	return object name More...

Private Attributes
const std::vector< PosType >	kvecs

const int	mink

const int	maxk

std::vector< RealType >	k2neg

Additional Inherited Members
Public Types inherited from SPOSet
using	ValueVector = OrbitalSetTraits< ValueType >::ValueVector

using	ValueMatrix = OrbitalSetTraits< ValueType >::ValueMatrix

using	GradVector = OrbitalSetTraits< ValueType >::GradVector

using	GradMatrix = OrbitalSetTraits< ValueType >::GradMatrix

using	HessVector = OrbitalSetTraits< ValueType >::HessVector

using	HessMatrix = OrbitalSetTraits< ValueType >::HessMatrix

using	GGGVector = OrbitalSetTraits< ValueType >::GradHessVector

using	GGGMatrix = OrbitalSetTraits< ValueType >::GradHessMatrix

using	SPOMap = std::map< std::string, const std::unique_ptr< const SPOSet > >

using	OffloadMWVGLArray = Array< ValueType, 3, OffloadPinnedAllocator< ValueType > >

using	OffloadMWVArray = Array< ValueType, 2, OffloadPinnedAllocator< ValueType > >

template<typename DT >
using	OffloadMatrix = Matrix< DT, OffloadPinnedAllocator< DT > >

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 Attributes inherited from SPOSet
const std::string	my_name_
	name of the object, unique identifier More...

IndexType	OrbitalSetSize
	number of Single-particle orbitals More...

opt_variables_type	myVars
	Optimizable variables. More...

Detailed Description

Definition at line 24 of file FreeOrbital.h.

Constructor & Destructor Documentation

◆ FreeOrbital()

FreeOrbital	(	const std::string &	my_name,
		const std::vector< PosType > &	kpts_cart
	)

Definition at line 6 of file FreeOrbital.cpp.

References qmcplusplus::dot(), FreeOrbital::k2neg, FreeOrbital::kvecs, FreeOrbital::maxk, and SPOSet::OrbitalSetSize.

     : SPOSet(my_name),
       kvecs(kpts_cart),
 #ifdef QMC_COMPLEX
       mink(0), // first k at twist may not be 0
 #else
       mink(1),                   // treat k=0 as special case
 #endif
       maxk(kpts_cart.size())
 {
 #ifdef QMC_COMPLEX
   OrbitalSetSize = maxk;
 #else
   OrbitalSetSize = 2 * maxk - 1; // k=0 has no (cos, sin) split
 #endif
   k2neg.resize(maxk);
   for (int ik = 0; ik < maxk; ik++)
     k2neg[ik] = -dot(kvecs[ik], kvecs[ik]);
 }

◆ ~FreeOrbital()

~FreeOrbital ( )

Definition at line 26 of file FreeOrbital.cpp.

26 {}

Member Function Documentation

◆ evaluate_notranspose() [1/3]

void evaluate_notranspose	(	const ParticleSet &	P,
		int	first,
		int	last,
		ValueMatrix &	logdet,
		GradMatrix &	dlogdet,
		ValueMatrix &	d2logdet
	)

overridevirtual

evaluate the values, gradients and laplacians of this single-particle orbital for [first,last) particles

Parameters

[in]	P	current ParticleSet
[in]	first	starting index of the particles
[in]	last	ending index of the particles
[out]	logdet	determinant matrix to be inverted
[out]	dlogdet	gradients
[out]	d2logdet	laplacians

Implements SPOSet.

Definition at line 78 of file FreeOrbital.cpp.

References FreeOrbital::evaluateVGL(), and SPOSet::OrbitalSetSize.

 {
   for (int iat = first, i = 0; iat < last; iat++, i++)
   {
     ValueVector p(phi[i], OrbitalSetSize);
     GradVector dp(dphi[i], OrbitalSetSize);
     ValueVector d2p(d2phi[i], OrbitalSetSize);
     evaluateVGL(P, iat, p, dp, d2p);
   }
 }

◆ evaluate_notranspose() [2/3]

void evaluate_notranspose	(	const ParticleSet &	P,
		int	first,
		int	last,
		ValueMatrix &	logdet,
		GradMatrix &	dlogdet,
		HessMatrix &	grad_grad_logdet
	)

overridevirtual

evaluate the values, gradients and hessians of this single-particle orbital for [first,last) particles

Parameters

P	current ParticleSet
first	starting index of the particles
last	ending index of the particles
logdet	determinant matrix to be inverted
dlogdet	gradients
grad_grad_logdet	hessians

Reimplemented from SPOSet.

Definition at line 94 of file FreeOrbital.cpp.

References ParticleSet::activeR(), qmcplusplus::dot(), FreeOrbital::kvecs, FreeOrbital::maxk, FreeOrbital::mink, OHMMS_DIM, SPOSet::OrbitalSetSize, and qmcplusplus::sincos().

 {
   RealType sinkr, coskr;
   ValueType phi_of_r;
   for (int iat = first, i = 0; iat < last; iat++, i++)
   {
     ValueVector p(phi[i], OrbitalSetSize);
     GradVector dp(dphi[i], OrbitalSetSize);
     HessVector hess(d2phi_mat[i], OrbitalSetSize);
 
     const PosType& r = P.activeR(iat);
     for (int ik = mink; ik < maxk; ik++)
     {
       sincos(dot(kvecs[ik], r), &sinkr, &coskr);
 #ifdef QMC_COMPLEX
       // phi(r) = cos(kr)+i*sin(kr)
       phi_of_r = ValueType(coskr, sinkr);
       p[ik]    = phi_of_r;
       // i*phi(r) = -sin(kr) + i*cos(kr)
       dp[ik] = ValueType(-sinkr, coskr) * kvecs[ik];
       for (int la = 0; la < OHMMS_DIM; la++)
       {
         hess[ik](la, la) = -phi_of_r * (kvecs[ik])[la] * (kvecs[ik])[la];
         for (int lb = la + 1; lb < OHMMS_DIM; lb++)
         {
           hess[ik](la, lb) = -phi_of_r * (kvecs[ik])[la] * (kvecs[ik])[lb];
           hess[ik](lb, la) = hess[ik](la, lb);
         }
       }
 #else
       const int j2 = 2 * ik;
       const int j1 = j2 - 1;
       p[j1]        = coskr;
       p[j2]        = sinkr;
       dp[j1]       = -sinkr * kvecs[ik];
       dp[j2]       = coskr * kvecs[ik];
       for (int la = 0; la < OHMMS_DIM; la++)
       {
         hess[j1](la, la) = -coskr * (kvecs[ik])[la] * (kvecs[ik])[la];
         hess[j2](la, la) = -sinkr * (kvecs[ik])[la] * (kvecs[ik])[la];
         for (int lb = la + 1; lb < OHMMS_DIM; lb++)
         {
           hess[j1](la, lb) = -coskr * (kvecs[ik])[la] * (kvecs[ik])[lb];
           hess[j2](la, lb) = -sinkr * (kvecs[ik])[la] * (kvecs[ik])[lb];
           hess[j1](lb, la) = hess[j1](la, lb);
           hess[j2](lb, la) = hess[j2](la, lb);
         }
       }
 #endif
     }
 #ifndef QMC_COMPLEX
     p[0]    = 1.0;
     dp[0]   = 0.0;
     hess[0] = 0.0;
 #endif
   }
 }

◆ evaluate_notranspose() [3/3]

void evaluate_notranspose	(	const ParticleSet &	P,
		int	first,
		int	last,
		ValueMatrix &	logdet,
		GradMatrix &	dlogdet,
		HessMatrix &	grad_grad_logdet,
		GGGMatrix &	grad_grad_grad_logdet
	)

overridevirtual

evaluate the values, gradients, hessians and third derivatives of this single-particle orbital for [first,last) particles

Parameters

P	current ParticleSet
first	starting index of the particles
last	ending index of the particles
logdet	determinant matrix to be inverted
dlogdet	gradients
grad_grad_logdet	hessians
grad_grad_grad_logdet	third derivatives

Reimplemented from SPOSet.

Definition at line 157 of file FreeOrbital.cpp.

References ParticleSet::activeR(), qmcplusplus::dot(), FreeOrbital::kvecs, FreeOrbital::maxk, FreeOrbital::mink, OHMMS_DIM, SPOSet::OrbitalSetSize, and qmcplusplus::sincos().

 {
   RealType sinkr, coskr;
   ValueType phi_of_r;
   for (int iat = first, i = 0; iat < last; iat++, i++)
   {
     ValueVector p(phi[i], OrbitalSetSize);
     GradVector dp(dphi[i], OrbitalSetSize);
     HessVector hess(d2phi_mat[i], OrbitalSetSize);
     GGGVector ggg(d3phi_mat[i], OrbitalSetSize);
 
     const PosType& r = P.activeR(iat);
     for (int ik = mink; ik < maxk; ik++)
     {
       sincos(dot(kvecs[ik], r), &sinkr, &coskr);
 #ifdef QMC_COMPLEX
       const ValueType compi(0, 1);
       phi_of_r = ValueType(coskr, sinkr);
       p[ik]    = phi_of_r;
       dp[ik]   = compi * phi_of_r * kvecs[ik];
       for (int la = 0; la < OHMMS_DIM; la++)
       {
         hess[ik](la, la) = -phi_of_r * (kvecs[ik])[la] * (kvecs[ik])[la];
         for (int lb = la + 1; lb < OHMMS_DIM; lb++)
         {
           hess[ik](la, lb) = -phi_of_r * (kvecs[ik])[la] * (kvecs[ik])[lb];
           hess[ik](lb, la) = hess[ik](la, lb);
         }
       }
       for (int la = 0; la < OHMMS_DIM; la++)
       {
         ggg[ik][la] = compi * (kvecs[ik])[la] * hess[ik];
       }
 #else
       const int j2 = 2 * ik;
       const int j1 = j2 - 1;
       p[j1]        = coskr;
       p[j2]        = sinkr;
       dp[j1]       = -sinkr * kvecs[ik];
       dp[j2]       = coskr * kvecs[ik];
       for (int la = 0; la < OHMMS_DIM; la++)
       {
         hess[j1](la, la)    = -coskr * (kvecs[ik])[la] * (kvecs[ik])[la];
         hess[j2](la, la)    = -sinkr * (kvecs[ik])[la] * (kvecs[ik])[la];
         ggg[j1][la](la, la) = sinkr * (kvecs[ik])[la] * (kvecs[ik])[la] * (kvecs[ik])[la];
         ggg[j2][la](la, la) = -coskr * (kvecs[ik])[la] * (kvecs[ik])[la] * (kvecs[ik])[la];
         for (int lb = la + 1; lb < OHMMS_DIM; lb++)
         {
           hess[j1](la, lb)    = -coskr * (kvecs[ik])[la] * (kvecs[ik])[lb];
           hess[j2](la, lb)    = -sinkr * (kvecs[ik])[la] * (kvecs[ik])[lb];
           hess[j1](lb, la)    = hess[j1](la, lb);
           hess[j2](lb, la)    = hess[j2](la, lb);
           ggg[j1][la](lb, la) = sinkr * (kvecs[ik])[la] * (kvecs[ik])[lb] * (kvecs[ik])[la];
           ggg[j2][la](lb, la) = -coskr * (kvecs[ik])[la] * (kvecs[ik])[lb] * (kvecs[ik])[la];
           ggg[j1][la](la, lb) = ggg[j1][la](lb, la);
           ggg[j2][la](la, lb) = ggg[j2][la](lb, la);
           ggg[j1][lb](la, la) = ggg[j1][la](lb, la);
           ggg[j2][lb](la, la) = ggg[j2][la](lb, la);
           ggg[j1][la](lb, lb) = sinkr * (kvecs[ik])[la] * (kvecs[ik])[lb] * (kvecs[ik])[lb];
           ggg[j2][la](lb, lb) = -coskr * (kvecs[ik])[la] * (kvecs[ik])[lb] * (kvecs[ik])[lb];
           ggg[j1][lb](la, lb) = ggg[j1][la](lb, lb);
           ggg[j2][lb](la, lb) = ggg[j2][la](lb, lb);
           ggg[j1][lb](lb, la) = ggg[j1][la](lb, lb);
           ggg[j2][lb](lb, la) = ggg[j2][la](lb, lb);
           for (int lc = lb + 1; lc < OHMMS_DIM; lc++)
           {
             ggg[j1][la](lb, lc) = sinkr * (kvecs[ik])[la] * (kvecs[ik])[lb] * (kvecs[ik])[lc];
             ggg[j2][la](lb, lc) = -coskr * (kvecs[ik])[la] * (kvecs[ik])[lb] * (kvecs[ik])[lc];
             ggg[j1][la](lc, lb) = ggg[j1][la](lb, lc);
             ggg[j2][la](lc, lb) = ggg[j2][la](lb, lc);
             ggg[j1][lb](la, lc) = ggg[j1][la](lb, lc);
             ggg[j2][lb](la, lc) = ggg[j2][la](lb, lc);
             ggg[j1][lb](lc, la) = ggg[j1][la](lb, lc);
             ggg[j2][lb](lc, la) = ggg[j2][la](lb, lc);
             ggg[j1][lc](la, lb) = ggg[j1][la](lb, lc);
             ggg[j2][lc](la, lb) = ggg[j2][la](lb, lc);
             ggg[j1][lc](lb, la) = ggg[j1][la](lb, lc);
             ggg[j2][lc](lb, la) = ggg[j2][la](lb, lc);
           }
         }
       }
 #endif
     }
 #ifndef QMC_COMPLEX
     p[0]    = 1.0;
     dp[0]   = 0.0;
     hess[0] = 0.0;
     ggg[0]  = 0.0;
 #endif
   }
 }

◆ evaluateValue()

void evaluateValue	(	const ParticleSet &	P,
		int	iat,
		ValueVector &	psi
	)

overridevirtual

evaluate the values of this single-particle orbital set

Parameters

P	current ParticleSet
iat	active particle
psi	values of the SPO

Implements SPOSet.

Definition at line 57 of file FreeOrbital.cpp.

References ParticleSet::activeR(), qmcplusplus::dot(), FreeOrbital::kvecs, FreeOrbital::maxk, FreeOrbital::mink, and qmcplusplus::sincos().

 {
   const PosType& r = P.activeR(iat);
   RealType sinkr, coskr;
   for (int ik = mink; ik < maxk; ik++)
   {
     sincos(dot(kvecs[ik], r), &sinkr, &coskr);
 #ifdef QMC_COMPLEX
     pvec[ik] = ValueType(coskr, sinkr);
 #else
     const int j2 = 2 * ik;
     const int j1 = j2 - 1;
     pvec[j1]     = coskr;
     pvec[j2]     = sinkr;
 #endif
   }
 #ifndef QMC_COMPLEX
   pvec[0] = 1.0;
 #endif
 }

◆ evaluateVGL()

void evaluateVGL	(	const ParticleSet &	P,
		int	iat,
		ValueVector &	psi,
		GradVector &	dpsi,
		ValueVector &	d2psi
	)

overridevirtual

evaluate the values, gradients and laplacians of this single-particle orbital set

Parameters

P	current ParticleSet
iat	active particle
psi	values of the SPO
dpsi	gradients of the SPO
d2psi	laplacians of the SPO

Implements SPOSet.

Definition at line 28 of file FreeOrbital.cpp.

References ParticleSet::activeR(), qmcplusplus::dot(), FreeOrbital::k2neg, FreeOrbital::kvecs, FreeOrbital::maxk, FreeOrbital::mink, and qmcplusplus::sincos().

Referenced by FreeOrbital::evaluate_notranspose().

 {
   const PosType& r = P.activeR(iat);
   RealType sinkr, coskr;
   for (int ik = mink; ik < maxk; ik++)
   {
     sincos(dot(kvecs[ik], r), &sinkr, &coskr);
 #ifdef QMC_COMPLEX
     pvec[ik]   = ValueType(coskr, sinkr);
     dpvec[ik]  = ValueType(-sinkr, coskr) * kvecs[ik];
     d2pvec[ik] = ValueType(k2neg[ik] * coskr, k2neg[ik] * sinkr);
 #else
     const int j2 = 2 * ik;
     const int j1 = j2 - 1;
     pvec[j1]     = coskr;
     pvec[j2]     = sinkr;
     dpvec[j1]    = -sinkr * kvecs[ik];
     dpvec[j2]    = coskr * kvecs[ik];
     d2pvec[j1]   = k2neg[ik] * coskr;
     d2pvec[j2]   = k2neg[ik] * sinkr;
 #endif
   }
 #ifndef QMC_COMPLEX
   pvec[0]   = 1.0;
   dpvec[0]  = 0.0;
   d2pvec[0] = 0.0;
 #endif
 }

◆ getClassName()

std::string getClassName ( ) const

inlineoverridevirtual

return class name

Implements SPOSet.

Definition at line 30 of file FreeOrbital.h.

30 { return "FreeOrbital"; }

◆ makeClone()

std::unique_ptr<SPOSet> makeClone ( ) const

inlineoverridevirtual

make a clone of itself every derived class must implement this to have threading working correctly.

Reimplemented from SPOSet.

Definition at line 64 of file FreeOrbital.h.

64 { return std::make_unique<FreeOrbital>(*this); }

◆ report()

void report ( const std::string & pad ) const

overridevirtual

print SPOSet information

Reimplemented from SPOSet.

Definition at line 255 of file FreeOrbital.cpp.

References qmcplusplus::app_log(), and FreeOrbital::kvecs.

 {
   app_log() << pad << "FreeOrbital report" << std::endl;
   for (int ik = 0; ik < kvecs.size(); ik++)
   {
     app_log() << pad << ik << " " << kvecs[ik] << std::endl;
   }
   app_log() << pad << "end FreeOrbital report" << std::endl;
   app_log().flush();
 }

◆ setOrbitalSetSize()

void setOrbitalSetSize ( int norbs )

inlineoverridevirtual

set the OrbitalSetSize

Parameters

norbs number of single-particle orbitals Ye: I prefer to remove this interface in the future. SPOSet builders need to handle the size correctly. It doesn't make sense allowing to set the value at any place in the code.

Implements SPOSet.

Definition at line 65 of file FreeOrbital.h.

65 { throw std::runtime_error("not implemented"); }

Member Data Documentation

◆ k2neg

std::vector<RealType> k2neg

private

Definition at line 71 of file FreeOrbital.h.

Referenced by FreeOrbital::evaluateVGL(), and FreeOrbital::FreeOrbital().

◆ kvecs

const std::vector<PosType> kvecs

private

Definition at line 68 of file FreeOrbital.h.

Referenced by FreeOrbital::evaluate_notranspose(), FreeOrbital::evaluateValue(), FreeOrbital::evaluateVGL(), FreeOrbital::FreeOrbital(), and FreeOrbital::report().

◆ maxk

const int maxk

private

Definition at line 70 of file FreeOrbital.h.

Referenced by FreeOrbital::evaluate_notranspose(), FreeOrbital::evaluateValue(), FreeOrbital::evaluateVGL(), and FreeOrbital::FreeOrbital().

◆ mink

const int mink

private

Definition at line 69 of file FreeOrbital.h.

Referenced by FreeOrbital::evaluate_notranspose(), FreeOrbital::evaluateValue(), and FreeOrbital::evaluateVGL().

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

/home/pk7/projects/qmc/for_cron_doxygen/qmcpack/src/QMCWaveFunctions/ElectronGas/FreeOrbital.h
/home/pk7/projects/qmc/for_cron_doxygen/qmcpack/src/QMCWaveFunctions/ElectronGas/FreeOrbital.cpp

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ FreeOrbital()

◆ ~FreeOrbital()

Member Function Documentation

◆ evaluate_notranspose() [1/3]

◆ evaluate_notranspose() [2/3]

◆ evaluate_notranspose() [3/3]

◆ evaluateValue()

◆ evaluateVGL()

◆ getClassName()

◆ makeClone()

◆ report()

◆ setOrbitalSetSize()

Member Data Documentation

◆ k2neg

◆ kvecs

◆ maxk

◆ mink