Extracts the VJP implementing Christianson's two-phase algorithm from the definition of the function "two_phase_solver" so as to be able to call the VJP manually if desired.
This commit does not affect the way that two_phase_solver is called. The relevant tests still seem to hold.
Extracts the VJP implementing Christianson's two-phase algorithm from the definition of the function "two_phase_solver" so as to be able to call the VJP manually if desired.
This commit does not affect the way that two_phase_solver is called. The relevant tests still seem to hold.