36 namespace sixDoFSolvers
77 v() = tConstraints() & (v0() + aDamp()*0.5*deltaT0*
a0());
78 pi() = rConstraints() & (pi0() + aDamp()*0.5*deltaT0*tau0());
80 centreOfRotation() = centreOfRotation0() + deltaT*v();
84 pi() = rConstraints() & Qpi.
second();
87 updateAcceleration(fGlobal, tauGlobal);
95 v() += tConstraints() & aDamp()*0.5*deltaT*a();
96 pi() += rConstraints() & aDamp()*0.5*deltaT*tau();
Six degree of freedom motion for a rigid body.
virtual ~symplectic()
Destructor.
A list of keyword definitions, which are a keyword followed by a number of values (eg...
A 2-tuple for storing two objects of dissimilar types. The container is similar in purpose to std::pa...
defineTypeNameAndDebug(CrankNicolson, 0)
virtual void solve(bool firstIter, const vector &fGlobal, const vector &tauGlobal, scalar deltaT, scalar deltaT0)
Drag coefficient.
Macros for easy insertion into run-time selection tables.
addToRunTimeSelectionTable(sixDoFSolver, CrankNicolson, dictionary)
constexpr scalar pi(M_PI)
symplectic(const dictionary &dict, sixDoFRigidBodyMotion &body)
Construct from a dictionary and the body.
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...
const T2 & second() const noexcept
Access the second element.
const dimensionedScalar a0
Bohr radius: default SI units: [m].
const T1 & first() const noexcept
Access the first element.