PETSc MPI gather matrix
I generated a matrix from a bilinear form with PETSc as the linear-algebra backend. To run a solver that supports complex numbers, I want to collect this distributed matrix into a single matrix at node 0 and then convert it to scipy. Is there a routine to do so (I was checking PETSc Matrix routines but didnot find such)? Another option would be to use PETSc with support for complex numbers, can I do so with fenics?
Thanks in advance! Christopher