[petsc-dev] resurrecting the finite element energy functions thread

[Geoffrey Irving]

On Tue, Jan 21, 2014 at 10:04 AM, Jed Brown <jed at jedbrown.org> wrote:
> Geoffrey Irving <irving at naml.us> writes:
>
>> On Sun, Jan 19, 2014 at 5:03 PM, Jed Brown <jed at jedbrown.org> wrote:
>>> Geoffrey Irving <irving at naml.us> 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.
>
> I would use PetscScalar.

I uploaded the branch:

    irving/fe-scalars

There are no unit tests in petsc yet, but I do have unit tests outside
of petsc which pass.

Geoffrey