On Sun, Jan 19, 2014 at 5:03 PM, Jed Brown <[email protected]> wrote: > Geoffrey Irving <[email protected]> writes: >> It isn't, the only difference is that the function I'm imagining needs >> an extra quadrature rule argument. Nothing here is complicated, so if >> there are no particular opinions I'll write the function and send a >> patch. Can I write it hard coded to MPI_SUM first, and we can discuss >> how to generalize it after that? I am less clear on how to best write >> the generalized version. > > That's fine with me.
Should the result be an array of PetscScalar or PetscReal? My application is for objectives, so PetscReal, but PetscScalar is more generic. Geoffrey
