Integrate robust linear solver for stack_solve_algo=0
Supersedes #984
To be done: add interface option for robust_lin_solve
to set options_.simul.robust_lin_solve=1
Supersedes #984
To be done: add interface option for robust_lin_solve
to set options_.simul.robust_lin_solve=1