Legendre Transform
Sanjay Govindjee
sanjay at ce.berkeley.edu
Thu Feb 7 09:22:01 CST 2008
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.
>
More information about the petsc-users
mailing list