PetscFEIntegrateResidual_Basic seems to have a redundant argument
list. It takes a single PetscFE, an array of PetscFE's, and a field
index into the array. The single PetscFE argument is ignored. I
would imagine that either the field index or the single PetscFE should
be eliminated.
Also, PetscFEIntegrateResidual_Basic allocates the f0 array with size
quad.numPoints * (spatial dimension of fe[0])
but accesses it as if it had size
quad.numPoints * (number of components of fe[field])
Am I reading this right?
Geoffrey
