Set up fieldsplit ISes correctly for parallel running

Numbering of fields on process N is offset by the sum of the rows of the
big matrix on processes 0, ..., N-1 (inclusive).
Test mixed poisson problem in parallel with fieldsplit pc