本文整理汇总了C++中shared_ptr::AddForceVector方法的典型用法代码示例。如果您正苦于以下问题:C++ shared_ptr::AddForceVector方法的具体用法?C++ shared_ptr::AddForceVector怎么用?C++ shared_ptr::AddForceVector使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类shared_ptr
的用法示例。
在下文中一共展示了shared_ptr::AddForceVector方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1:
Eigen::Vector3d Atom::CalculateElectricalForce(const Eigen::Vector3d& chargePosition, Real charge, shared_ptr<ForceVectors> forceVectorStorage) {
Eigen::Vector3d forceDirection = chargePosition.array()-m_position.array();
Real rSquared = forceDirection.squaredNorm();
Real forceMagnitude = K_VACUUM*charge*m_effectiveNuclearCharge*Square(ELECTRIC_CHARGE)/rSquared;
// The force is in the direction of the charge.
forceDirection.normalize();
Eigen::Vector3d forceVector = forceDirection.array()*forceMagnitude;
if(forceVectorStorage) {
forceVectorStorage->AddForceVector(GetName(), forceVector);
}
// Add in forces due to sp orbitals.
for(int i=0; i<m_spOrbitals.size(); i++) {
shared_ptr<SpOrbital> orbital = m_spOrbitals[i];
Eigen::Vector3d spForce = orbital->CalculateElectricalForce(chargePosition, charge, forceVectorStorage);
forceVector.array() += spForce.array();
}
return forceVector;
}