I have found that Pan has a question about system matrix and solution as follows:
at this point I would like to ask, how many system
matrix and solution will be generated.
I mean if I run parallel code in a cluster, does it
mean there is one system matrix and system solution?
And each node manages one part of the matrix and
solution vector?


I want to know whether I may generate multiple system matrixes and obtain multiple solutions? In addition, if the problem to be solved includes several PDEs (they are coupled), whether is libmesh avalible? If it ok, how to deal with them?