Sorry for the long delay in getting back to you. I think there is not
directly a problem with the code; rather, with these Q0 elements, Galerkin is
not a good approach.
I went back to src/ksp/ksp/tutorials/ex32.c and ran a little deeper and noticed
that when running without Galerkin it gives
Ben,
This is interesting; it seems like it may be useful inside the PETSc
library code. Please take a quick look at
https://urldefense.us/v3/__https://bitbucket.org/saws/saws/src/master/__;!!G_uCfscf7eWS!eBNC69y99tie5v_KSbJ3PosiD5GwewD_jS-U2-1sd7l6XOaXwSy97N9job28xHc_2UBnSz-cuOUpWS-u6j5
Thank you Satish and Junchao,
Actually following Satish's instructions I found that I am indeed using
the same compiler but for some reason my command does not include the
mpi include path.
I am going to check now my cmake files because it was initially found so
I should be doing something w