Here’s a version of batty’s answer, but this computes the correct inverse. batty’s version computes the transpose of the inverse. // computes the inverse of a matrix m double det = m(0, 0) * (m(1, 1) * m(2, 2) – m(2, 1) * m(1, 2)) – m(0, 1) * (m(1, 0) * m(2, 2) – … Read more
You should have a look at numpy if you do matrix manipulation. This is a module mainly written in C, which will be much faster than programming in pure python. Here is an example of how to invert a matrix, and do other matrix manipulation. from numpy import matrix from numpy import linalg A = … Read more
here: bool gluInvertMatrix(const double m[16], double invOut[16]) { double inv[16], det; int i; inv[0] = m[5] * m[10] * m[15] – m[5] * m[11] * m[14] – m[9] * m[6] * m[15] + m[9] * m[7] * m[14] + m[13] * m[6] * m[11] – m[13] * m[7] * m[10]; inv[4] = -m[4] * m[10] … Read more
solve(c) does give the correct inverse. The issue with your code is that you are using the wrong operator for matrix multiplication. You should use solve(c) %*% c to invoke matrix multiplication in R. R performs element by element multiplication when you invoke solve(c) * c.