Tests for solving systems of PDEs with various solvers added. In particular solving systems with a direct solver (MKL) was not done in a best way as a matrix pattern was created twice. This problem should be fixed now. A message is printed to inform about a pattern unrolling. I leave the message in there for now to check that unrolling called only if appropriate.
There is also a function now which turns face tags into node tags.
