1inline int mdcTwoInv(
double matrix[3],
double invmat[3] ) {
3 double det = matrix[0] * matrix[2] - matrix[1] * matrix[1];
6 invmat[0] = invmat[2] = invmat[1] = 0.00000;
11 double detinv = 1. / det;
12 invmat[0] = matrix[2] * detinv;
13 invmat[2] = matrix[0] * detinv;
14 invmat[1] = -matrix[1] * detinv;