IMP: faster analyse step by avoiding matrix inversion
In the analysis step: To avoid the inversion of :math:\left(C^{l}_{DD}+\alpha_{l+1} C_{D}\right)
, the product :math:\left(C^{l}_{DD}+\alpha_{l+1} C_{D}\right) ^{-1}\left(d^{l}_{uc,j} - d^{l}_{j} \right)
is solved linearly as :math:A^{-1}b = x
which is equivalent to solve :math:Ax = b
.
See Issue #11 (closed)
Edited by Antoine COLLET