In thinking about encapsulation and uses-cases #423, the idea arose of writing a back-end in scikit-fem for scipy.integrate.solve_bvp; that might make a nice demonstration of scikti-fem.
The problems are generally nonlinear, so this might build on ex10 or more recent ideas of using scipy.optimize.root #119 for Newton iteration or even Jacobian-free Newton–Krylov.
In thinking about encapsulation and uses-cases #423, the idea arose of writing a back-end in scikit-fem for
scipy.integrate.solve_bvp
; that might make a nice demonstration of scikti-fem.The problems are generally nonlinear, so this might build on ex10 or more recent ideas of using
scipy.optimize.root
#119 for Newton iteration or even Jacobian-free Newton–Krylov.