I have 4x4 matrix (GL_MODELVIEW of the current object) …
How to find inverse matrix?
(i want to get object space position of light)
This worked well for me, and is reasonably simple :
fmatrix4 &fmatrix4::invert(void)
{
fmatrix4 D;
int i, j, k;D.identity();
for( i = 0; i < 4; i++ ) {
float tmp = 1.0f / data[i*4+i]; for( j = 3; j >= 0; j-- ) { D.data[i*4+j] *= tmp; data[i*4+j] *= tmp; } for( j = 0; j < 4; j++ ) if( j != i ) { tmp = data[j*4+i]; for( k = 3; k >= 0; k-- ) { D.data[j*4+k] -= D.data[i*4+k] * tmp; data[j*4+k] -= data[i*4+k] * tmp; } }
}
*this = D;
return( *this );
}
Oh yes: and if your matrix is guaranteed to only encode a rotation in the upper 3x3 part, you could of course use a simple transpose. Optional translation can also be taken into account, without needing a full inversion. If you have scaling, shearing or even a projection (rather unlikely for the modelview matrix), then you need a full invert.
[This message has been edited by AlexH (edited 01-04-2004).]
Thanks! It work very well …