Improve `energy_preserving_order`

[ranocha]

• It should share more functionality with is_energy_preserving
• It should be implemented for a general B-series