in default alg branch on ```julia if SciMLBase.has_ldiv!(A) # do ldiv!(u, A, b) via solve_function elseif SciMLBase.has_ldiv(A) # do u .= A \ b via solve_function end ```