本文整理汇总了C++中ParticleSet::G方法的典型用法代码示例。如果您正苦于以下问题：C++ ParticleSet::G方法的具体用法？C++ ParticleSet::G怎么用？C++ ParticleSet::G使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类ParticleSet的用法示例。
在下文中一共展示了ParticleSet::G方法的2个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: alternateRatio
RNDiracDeterminantBase::ValueType RNDiracDeterminantBase::alternateRatio(ParticleSet& P)
{
//returns psi_T/psi_G
for (int i=0, iat=FirstIndex; i<NumPtcls; i++, iat++)
{
P.G(iat) += myG_alternate(iat) - myG(iat);
P.L(iat) += myL_alternate(iat) - myL(iat);
}
RealType sgn = std::cos(alternatePhaseValue);
return sgn*std::exp(alternateLogValue-LogValue);
}示例2: grad
void
ThreeBodyGeminal::evaluateLogAndStore(ParticleSet& P)
{
GeminalBasis->evaluateForWalkerMove(P);
MatrixOperators::product(GeminalBasis->Y, Lambda, V);
Y=GeminalBasis->Y;
dY=GeminalBasis->dY;
d2Y=GeminalBasis->d2Y;
Uk=0.0;
LogValue=RealType();
for(int i=0; i< NumPtcls-1; i++)
{
const RealType* restrict yptr=GeminalBasis->Y[i];
for(int j=i+1; j<NumPtcls; j++)
{
RealType x= simd::dot(V[j],yptr,BasisSize);
LogValue += x;
Uk[i]+= x;
Uk[j]+= x;
}
}
for(int i=0; i<NumPtcls; i++)
{
const PosType* restrict dptr=GeminalBasis->dY[i];
const RealType* restrict d2ptr=GeminalBasis->d2Y[i];
const RealType* restrict vptr=V[0];
BasisSetType::GradType grad(0.0);
BasisSetType::ValueType lap(0.0);
for(int j=0; j<NumPtcls; j++, vptr+=BasisSize)
{
if(j==i)
{
dUk(i,j) = 0.0;
d2Uk(i,j)= 0.0;
}
else
{
grad+= (dUk(i,j) = simd::dot(vptr,dptr,BasisSize));
lap += (d2Uk(i,j)= simd::dot(vptr,d2ptr,BasisSize));
}
}
P.G(i)+=grad;
P.L(i)+=lap;
}
}