I have a code that takes a gmesh file and attempts to assemble subsystems from partitioned meshes. I can assemble the right side of the PDE of interest, but when I try to assemble the stiffness matrix, I get the following error:
*** Error: Unable to successfully call PETSc function 'MatAssemblyEnd'.
*** Reason: PETSc error code is: 63.
*** Where: This error was encountered inside /build/dolfin-k_QrtL/dolfin-1.6.0/dolfin/la/PETScMatrix.cpp.
*** Process: 0
*** DOLFIN version: 1.6.0
*** Git changeset: unknown
Where the PETSc error code 63 is:
#define PETSC_ERR_ARG_OUTOFRANGE 63 /* input argument, out of range *
I can't imagine why there would be a difference between the calls to assemble since either uses the test and trial functions from the same mesh!
