This covers both, the serial and parallel cases: seqaij and mpiaij.
The latter matrix type will not work for serial cases---#3203.