In case the call to Cholesky fails, we now switch to Gram-Schmidt orthogonalization (for this call only). This avoid crashes of the code and makes it more stable.
Making the Cholesky orthogonalization more stable.
- I have checked that my code follows the Octopus coding standards
- I have added tests for all the new features added in this request.