BAIJ matrix type does not work with multi-process runs

[Slaedr]

When any case is run with implicit time-stepping using more than one MPI rank, using the BAIJ matrix type for the Jacobian matrices results in a runtime error by PETSc. It says "PETSc has generated inconsistent data" and "invalid column map". However, using the regular AIJ format works fine.