Using PETSC for all computations. Now to test in parallel.
Parallel working.
Parallel mesh loading is functioning
Next step is to use a higher order ODE solver
Heun's method implemented. RK4 is not perfect yet.