53 #define M4(i, j)	((i << 2) + j)  macro
64             res[M4(i,j)] = m1[M4(i,0)] * m2[M4(0,j)] +  in gl_mult_matrixf()
65                         m1[M4(i,1)] * m2[M4(1,j)] +  in gl_mult_matrixf()
66                         m1[M4(i,2)] * m2[M4(2,j)] +  in gl_mult_matrixf()
67                         m1[M4(i,3)] * m2[M4(3,j)];  in gl_mult_matrixf()
202     float determinant = +m1[M4(0,0)]*(m1[M4(1,1)]*m1[M4(2,2)]-m1[M4(2,1)]*m1[M4(1,2)])  in get_inv_transpose_3x3()
203                         -m1[M4(0,1)]*(m1[M4(1,0)]*m1[M4(2,2)]-m1[M4(1,2)]*m1[M4(2,0)])  in get_inv_transpose_3x3()
204                         +m1[M4(0,2)]*(m1[M4(1,0)]*m1[M4(2,1)]-m1[M4(1,1)]*m1[M4(2,0)]);  in get_inv_transpose_3x3()
208     inv_transpose_result[M3(0,0)] =  (m1[M4(1,1)]*m1[M4(2,2)]-m1[M4(2,1)]*m1[M4(1,2)])*invdet;  in get_inv_transpose_3x3()
209     inv_transpose_result[M3(1,0)] = -(m1[M4(0,1)]*m1[M4(2,2)]-m1[M4(0,2)]*m1[M4(2,1)])*invdet;  in get_inv_transpose_3x3()
210     inv_transpose_result[M3(2,0)] =  (m1[M4(0,1)]*m1[M4(1,2)]-m1[M4(0,2)]*m1[M4(1,1)])*invdet;  in get_inv_transpose_3x3()
211     inv_transpose_result[M3(0,1)] = -(m1[M4(1,0)]*m1[M4(2,2)]-m1[M4(1,2)]*m1[M4(2,0)])*invdet;  in get_inv_transpose_3x3()
212     inv_transpose_result[M3(1,1)] =  (m1[M4(0,0)]*m1[M4(2,2)]-m1[M4(0,2)]*m1[M4(2,0)])*invdet;  in get_inv_transpose_3x3()
213     inv_transpose_result[M3(2,1)] = -(m1[M4(0,0)]*m1[M4(1,2)]-m1[M4(1,0)]*m1[M4(0,2)])*invdet;  in get_inv_transpose_3x3()
214     inv_transpose_result[M3(0,2)] =  (m1[M4(1,0)]*m1[M4(2,1)]-m1[M4(2,0)]*m1[M4(1,1)])*invdet;  in get_inv_transpose_3x3()
215     inv_transpose_result[M3(1,2)] = -(m1[M4(0,0)]*m1[M4(2,1)]-m1[M4(2,0)]*m1[M4(0,1)])*invdet;  in get_inv_transpose_3x3()
216     inv_transpose_result[M3(2,2)] =  (m1[M4(0,0)]*m1[M4(1,1)]-m1[M4(1,0)]*m1[M4(0,1)])*invdet;  in get_inv_transpose_3x3()