Question about BDF2

Dear PyFR developers,
recently I read some papers about researching incompressible flow with PyFR, almost all researchers consider BDF2 as the physical time discrete format, which is no longer available in the latest version(PyFRv1.15.0), is there any defect about this format ? could I use sdirk33/sdirk43 as the physical time scheme when I test the cases in the papers?
Best regards, Mao

Yep, we got rid of BDF support a few versions ago. The problem with BDF is the initialisation of the time history. There are several ways to overcome this, but they all affect the solution for several physical time steps. Alternatively, DIRK schemes don’t have this issue, which makes them much more attractive as a method.