only output exelem files for first timestep and species (big IO / disk space saver since we don't currently have problems that require time-dependent exelem files. This is set with a flag so it can be easily changed if we have future equations sets or problem setups that require this)
BDF solver
Check whether array has already been allocated in PETSc context set
In RHS evaluate use known stateData for contiguous state data (the temporary states used in the non-contiguous case is intended for use with a non-contiguous offset)
Test now passing: Python_Bindings_CellMLModelIntegration_n98_BDF
Reaction diffusion
BDF solver