[petsc-users] numerical quadrature

[Gideon Simpson]

Suppose I have a function f at sample points x, with x and f both stored as Vec distributed structures.  What I would like to do is compute an estimate of the anti derivative of f,

\int_a^x f(s)ds

for a<= x <=b.

One way I can see how to compute this efficiently is to do the numerical quadrature on each node, and then use standard MPI to send the successive cumulative quantity from processor 0 to 1 to 2, and so on.  I am wondering if there is a “PETSc” way to do this kind of calculation, as opposed to relying on MPI code.

-gideon

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150204/43330d7c/attachment.html>