[petsc-dev] questions about PetscFEIntegrateResidual_Basic

[Geoffrey Irving]

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