Initial implementation of LU decomposition for matrix4f

this version uses pivoting in order to provide high stability for
inversion algorithms. still not one of the fastest approaches out there
- O(NM).
but since we are dealing with at mot 4x4 matrices (3d homogeneous
transformations, this should work fine for us).