Matrix preconditioning
The Doolittle LU matrix decomposition that is used to invert a matrix does not behave well in case of (very) badly conditioned matrices. Therefore it would be necessary to implement a conditioning step (partial pivoting). This is necessary in order to solve the regressions for the new despiking filter (issue 630).