pygaia/astrometry/vectorastrometry.py
def elementaryRotationMatrix(axis, rotationAngle):
  """
  Construct an elementary rotation matrix describing a rotation around the x, y, or z-axis.

  Parameters
  ----------

  axis          - Axis around which to rotate ("x", "y", or "z")
  rotationAngle - the rotation angle in radians

  Returns
  -------

  The rotation matrix

  Example usage
  -------------

  rotmat = elementaryRotationMatrix("y", pi/6.0)
  """
  if (axis=="x" or axis=="X"):
    return array([[1.0, 0.0, 0.0], [0.0, cos(rotationAngle), sin(rotationAngle)], [0.0,
      -sin(rotationAngle), cos(rotationAngle)]])
  elif (axis=="y" or axis=="Y"):
    return array([[cos(rotationAngle), 0.0, -sin(rotationAngle)], [0.0, 1.0, 0.0], [sin(rotationAngle),
      0.0, cos(rotationAngle)]])
  elif (axis=="z" or axis=="Z"):
    return array([[cos(rotationAngle), sin(rotationAngle), 0.0], [-sin(rotationAngle),
      cos(rotationAngle), 0.0], [0.0, 0.0, 1.0]])
  else:
    raise Exception("Unknown rotation axis "+axis+"!")
def elementaryRotationMatrix(axis, rotationAngle):
  if (axis=="x" or axis=="X"):
    return array([[1.0, 0.0, 0.0], [0.0, cos(rotationAngle), sin(rotationAngle)], [0.0,
      -sin(rotationAngle), cos(rotationAngle)]])
  elif (axis=="y" or axis=="Y"):
    return array([[cos(rotationAngle), 0.0, -sin(rotationAngle)], [0.0, 1.0, 0.0], [sin(rotationAngle),
      0.0, cos(rotationAngle)]])
  elif (axis=="z" or axis=="Z"):
    return array([[cos(rotationAngle), sin(rotationAngle), 0.0], [-sin(rotationAngle),
      cos(rotationAngle), 0.0], [0.0, 0.0, 1.0]])
  else:
    raise Exception("Unknown rotation axis "+axis+"!")
def elementaryRotationMatrix ( axis , rotationAngle ) : if ( axis == "x" or axis == "X" ) : return array ( [ [ 1.0 , 0.0 , 0.0 ] , [ 0.0 , cos ( rotationAngle ) , sin ( rotationAngle ) ] , [ 0.0 , - sin ( rotationAngle ) , cos ( rotationAngle ) ] ] ) elif ( axis == "y" or axis == "Y" ) : return array ( [ [ cos ( rotationAngle ) , 0.0 , - sin ( rotationAngle ) ] , [ 0.0 , 1.0 , 0.0 ] , [ sin ( rotationAngle ) , 0.0 , cos ( rotationAngle ) ] ] ) elif ( axis == "z" or axis == "Z" ) : return array ( [ [ cos ( rotationAngle ) , sin ( rotationAngle ) , 0.0 ] , [ - sin ( rotationAngle ) , cos ( rotationAngle ) , 0.0 ] , [ 0.0 , 0.0 , 1.0 ] ] ) else : raise Exception ( "Unknown rotation axis " + axis + "!" )
Construct an elementary rotation matrix describing a rotation around the x, y, or z-axis.

  Parameters
  ----------

  axis          - Axis around which to rotate ("x", "y", or "z")
  rotationAngle - the rotation angle in radians

  Returns
  -------

  The rotation matrix

  Example usage
  -------------

  rotmat = elementaryRotationMatrix("y", pi/6.0)
Construct an elementary rotation matrix describing a rotation around the x y or z - axis .