I'm using an affine transformation matrix to transform 2D coordinates from screen (magnitude 10e3) to small parts of fractal sets (magnitude as little as 10e-15).
I also need to map the other way round, so I simply invert the matrix. It works when the magnitudes are not too far apart, but it fails when the determinant gets too close to zero to be represented with standard double numbers.
The matrices I use are of the form:
a c e b d f 0 0 1
And the inversion algorithm I have is:
var dt = a * d - b * c; return new Matrix(d/dt, -b/dt, -c/dt, a/dt, (c * f - d * e) / dt, -(a * f - b * e) / dt);
Is there an alternate way of inverting the matrix? I'm not very good at mathematics, so I would need a simple solution, and one I could find the code, or an algorithm, for.
An example of such a matrix:
3.26091378894248e-9 -1.882689453850103e-9 -0.7172216437740687 -1.882689453814925e-9 -3.2609137888815494e-9 -0.23371832131832268 0 0 1