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;
  }
}