Slow performance while solving 1-d poisson equation using GeneralMIRK4() from BoundaryValueDiffEq

[ChrisRackauckas]

https://github.com/luchr/ODEInterface.jl has some wrappers of BVP solvers that might do a little bit better until we improve in this domain.