19 template<
class t_OutputMatrixType,
20 class t_InputMatrixType>
21 static void compute(t_OutputMatrixType & U_inverse,
22 const t_InputMatrixType & L)
28 U_inverse(i,i) = 1.0 / L(i, i);
31 double tmp = L.transpose().row(j).segment(j+1, i-j)*U_inverse.col(i).segment(j+1, i-j);
32 U_inverse(j, i) = - tmp / L(j,j);
static void compute(t_OutputMatrixType &U_inverse, const t_InputMatrixType &L)