Reconstruction/MdcPatRec/MdcTrkRecon/include/MdcTrkRecon/mdcTwoVec.h

Go to the documentation of this file.

/* Multiply a symmetric 2-d matrix by a vector. */
 
inline void mdcTwoVec( const double matrix[3], const double invect[2], double outvect[2] ) {
 
  double buffvec[2]; /* Buffer the output in case outvect = invect. */
 
  /************************************************************************/
 
  buffvec[0] = matrix[0] * invect[0] + matrix[1] * invect[1];
  buffvec[1] = matrix[1] * invect[0] + matrix[2] * invect[1];
 
  outvect[0] = buffvec[0];
  outvect[1] = buffvec[1];
 
  return;
}