derivative and adjoint_derivative in solve_and_derivative_internal do not actually pass their kwargs (particularly tolerances) to lsqr or lsmr solvers. lsmr appears to use different, hardcoded, tolerances for forward and reverse modes.
|
dz = _diffcp.lsqr(M, rhs).solution |