Reconstruction/MdcPatRec/MdcTrkRecon/include/MdcTrkRecon/mdcTwoInv.h

Go to the documentation of this file.

inline int mdcTwoInv( double matrix[3], double invmat[3] ) {
 
  double det = matrix[0] * matrix[2] - matrix[1] * matrix[1];
  if ( det == 0.0 )
  {
    invmat[0] = invmat[2] = invmat[1] = 0.00000;
    return 1;
  }
  else
  {
    double detinv = 1. / det;
    invmat[0]     = matrix[2] * detinv;
    invmat[2]     = matrix[0] * detinv;
    invmat[1]     = -matrix[1] * detinv;
    return 0;
  }
}