Better reverse mode propagation through linear solver, #889