RotationMatrix Namespace Reference

Utility functions to return rotations matrics. More...

Functions
template<bool is_ad = false>
GenericRealTensorValue< is_ad >	rotVecToZ (GenericRealVectorValue< is_ad > vec)
	provides a rotation matrix that will rotate the vector vec to the z axis (the "2" direction) More...
template<bool is_ad = false>
GenericRealTensorValue< is_ad >	rotVec1ToVec2 (GenericRealVectorValue< is_ad > vec1, GenericRealVectorValue< is_ad > vec2)
	provides a rotation matrix that will rotate the vector vec1 to vec2 More...
template<bool is_ad = false>
GenericRealTensorValue< is_ad >	rotVec2DToX (const GenericRealVectorValue< is_ad > &vec)
	provides a rotation matrix that will rotate the vector vec1 to the [1,0,0], assuming vec1[2]==0 More...

Detailed Description

Utility functions to return rotations matrics.

Function Documentation

rotVec1ToVec2()

template<bool is_ad = false>

GenericRealTensorValue<is_ad> RotationMatrix::rotVec1ToVec2	(	GenericRealVectorValue< is_ad >	vec1,
		GenericRealVectorValue< is_ad >	vec2
	)

provides a rotation matrix that will rotate the vector vec1 to vec2

Definition at line 66 of file RotationMatrix.h.

{
  GenericRealTensorValue<is_ad> rot1_to_z = rotVecToZ<is_ad>(vec1);
  GenericRealTensorValue<is_ad> rot2_to_z = rotVecToZ<is_ad>(vec2);
  return rot2_to_z.transpose() * rot1_to_z;
}

rotVec2DToX()

template<bool is_ad = false>

GenericRealTensorValue<is_ad> RotationMatrix::rotVec2DToX

(

const GenericRealVectorValue< is_ad > &

vec

)

provides a rotation matrix that will rotate the vector vec1 to the [1,0,0], assuming vec1[2]==0

Definition at line 77 of file RotationMatrix.h.

{
  using std::atan2, std::sin, std::cos;
  const GenericReal<is_ad> theta = atan2(vec(1), vec(0));
  const GenericReal<is_ad> st = sin(theta);
  const GenericReal<is_ad> ct = cos(theta);
  return GenericRealTensorValue<is_ad>(ct, st, 0., -st, ct, 0., 0., 0., 1.);
}

rotVecToZ()

template<bool is_ad = false>

GenericRealTensorValue<is_ad> RotationMatrix::rotVecToZ

(

GenericRealVectorValue< is_ad >

vec

)

provides a rotation matrix that will rotate the vector vec to the z axis (the "2" direction)

Definition at line 24 of file RotationMatrix.h.

Referenced by SolutionUserObjectBase::SolutionUserObjectBase().

{
  using std::sqrt, std::abs;
  // ensure that vec is normalised
  vec /= sqrt(vec * vec);

  // construct v0 and v1 to be orthonormal to vec
  // and form a RH basis, that is, so v1 x vec = v0

  // Use Gram-Schmidt method to find v1.
  GenericRealVectorValue<is_ad> v1;
  // Need a prototype for v1 first, and this is done by looking at the smallest component of vec
  GenericRealVectorValue<is_ad> w(abs(vec(0)), abs(vec(1)), abs(vec(2)));
  if ((w(2) >= w(1) && w(1) >= w(0)) || (w(1) >= w(2) && w(2) >= w(0)))
    // vec(0) is the smallest component
    v1(0) = 1;
  else if ((w(2) >= w(0) && w(0) >= w(1)) || (w(0) >= w(2) && w(2) >= w(1)))
    // vec(1) is the smallest component
    v1(1) = 1;
  else
    // vec(2) is the smallest component
    v1(2) = 1;
  // now Gram-Schmidt
  v1 -= (v1 * vec) * vec;
  v1 /= sqrt(v1 * v1);

  // now use v0 = v1 x vec
  GenericRealVectorValue<is_ad> v0;
  v0(0) = v1(1) * vec(2) - v1(2) * vec(1);
  v0(1) = v1(2) * vec(0) - v1(0) * vec(2);
  v0(2) = v1(0) * vec(1) - v1(1) * vec(0);

  // the desired rotation matrix is just
  GenericRealTensorValue<is_ad> rot(
      v0(0), v0(1), v0(2), v1(0), v1(1), v1(2), vec(0), vec(1), vec(2));
  return rot;
}