anonymous_namespace{ibm_triel.cpp} Namespace Reference

Functions
std::tuple< Utils::Vector3d, Utils::Vector3d, Utils::Vector3d >	RotateForces (Utils::Vector2d const &f1_rot, Utils::Vector2d const &f2_rot, Utils::Vector3d const &v12, Utils::Vector3d const &v13)
	Rotate calculated trielastic forces in the 2d plane back to the 3d plane.

Function Documentation

RotateForces()

std::tuple< Utils::Vector3d, Utils::Vector3d, Utils::Vector3d > anonymous_namespace{ibm_triel.cpp}::RotateForces	(	Utils::Vector2d const &	f1_rot,
		Utils::Vector2d const &	f2_rot,
		Utils::Vector3d const &	v12,
		Utils::Vector3d const &	v13
	)

Rotate calculated trielastic forces in the 2d plane back to the 3d plane.

Use knowledge that the x-axis in rotated system is parallel to r(p1->p2) in original system; To find the corresponding unit vector to y in the rotated system, construct vector perpendicular to r(p1->p2); note that f3 is not calculated here but is implicitly calculated by f3 = -(f1+f2) which is consistent with the literature

Returns

forces on particles 1, 2, 3

Definition at line 42 of file ibm_triel.cpp.