Re: Legendre Transform

[Sanjay Govindjee <sanjay@xxxxxxxxxxxxxxx>]

Can you define what you mean by Legendre transform?  The usual 
definition in physics
is L{f}(p)  = max_x (p.x - f(x)) where f(x) is convex.    When you say 
parallel, then I presume that
x lies in a very high dimensional space.  In which case you are simply 
looking at a root
finding problem in high dimensions and you could certainly use PETSc to 
solve the intermediate
linear solves of a Newton scheme (or just invoke SNES).

-sg

Timothy Stitt wrote:
> Hi all,
>
> I am not sure if this query is directly related to the PETSc library 
> per se, but I was wondering if anyone knows of an efficient 
> implementation of the Legendre Transform. A parallel implementation 
> would be even more ideal.
> I don't know that much about the implementation of integral 
> transforms...could it be done in PETSc for instance, if no suitable 
> library exists? A web search doesn't throw up to many relevant hits.
>
> Thanks in advance for any guidance given,
>
> Best,
>
> Tim.