Parallel constraint calculation on SerialMesh too

This will hopefully be a speedup for everybody, but more importantly
it should fix a bug where FEMFunctionBase boundary conditions could
not be used on a SerialMesh on more than one processor.