I have just noticed that the function linsolve_by_lu is not documented. Can somebody fix this please? Thanks, Stefano