I'm about to start writing multistep time integrator based on FEniCS, probably some variation of BDF. And I'm wondering if anyone else is/was working on such subject?
Gryphon project did tremendous job for one-step Runge-Kutta methods, but I can't find anything about multistep. Maybe I'm looking wrong, or maybe it's just weird idea to use BDF?
My problem is very stiff, non-linear, modified diffusion equation modelling bioheat transfer.
Any comments will be very appreciated :)