GEM can now calculate with multiple grid (matrix size 8000 rows +), better matrix reading and memory operations implemented. Text almost done.

Went back in commits for old GPU code. Now MPI turns on MPI code and without MPI turns on old GPU code.

and subtraction of diagonals

ODE solvers

See merge request !125