Corrected mat array index (im to i)
Added first parallel test case, tweaked writes

Some write statements had to be pre-pended to ensure that only processor
0 is producing them.  The GNUPlot output was set to only be produced by
processor 0.
Corrected xs array index
Added execution time output
Fixed fidas call so p>0 gets input

Also set up nnull to redirect p>0 output to (as needed).
Actually implement /dev/null output
Added test case for 4 processors
Updated test matrix readme
Added mpi_startup.f90 to the repo
Added spatial decomp and corresponding output

Calculation is performed to distribute spatial mesh to all processors
and perform load balancing.
Merge branch 'parallel' of github.com:jkulesza/EECS587_Proj into parallel