Lines Matching refs:LU
208 inverse_transform(const double *LU, const unsigned *p, float *v) in inverse_transform() argument
220 b[1] = (double)v[p[1]] - b[0] * LU[1 + 0 * 4]; in inverse_transform()
221 b[2] = (double)v[p[2]] - b[0] * LU[2 + 0 * 4]; in inverse_transform()
222 b[3] = (double)v[p[3]] - b[0] * LU[3 + 0 * 4]; in inverse_transform()
223 b[2] -= b[1] * LU[2 + 1 * 4]; in inverse_transform()
224 b[3] -= b[1] * LU[3 + 1 * 4]; in inverse_transform()
225 b[3] -= b[2] * LU[3 + 2 * 4]; in inverse_transform()
230 b[3] /= LU[3 + 3 * 4]; in inverse_transform()
231 b[0] -= b[3] * LU[0 + 3 * 4]; in inverse_transform()
232 b[1] -= b[3] * LU[1 + 3 * 4]; in inverse_transform()
233 b[2] -= b[3] * LU[2 + 3 * 4]; in inverse_transform()
235 b[2] /= LU[2 + 2 * 4]; in inverse_transform()
236 b[0] -= b[2] * LU[0 + 2 * 4]; in inverse_transform()
237 b[1] -= b[2] * LU[1 + 2 * 4]; in inverse_transform()
239 b[1] /= LU[1 + 1 * 4]; in inverse_transform()
240 b[0] -= b[1] * LU[0 + 1 * 4]; in inverse_transform()
242 b[0] /= LU[0 + 0 * 4]; in inverse_transform()
246 b[j] /= LU[j + j * 4]; in inverse_transform()
248 b[k] -= b[j] * LU[k + j * 4]; in inverse_transform()
251 b[0] /= LU[0 + 0 * 4]; in inverse_transform()
263 double LU[16]; /* column-major */ in weston_matrix_invert() local
267 if (matrix_invert(LU, perm, matrix) < 0) in weston_matrix_invert()
272 inverse_transform(LU, perm, &inverse->d[c * 4]); in weston_matrix_invert()