[petsc-dev] Problem mixing sieve and TS

Jed Brown jedbrown at mcs.anl.gov
Thu Oct 27 10:15:13 CDT 2011


On Thu, Oct 27, 2011 at 08:54, Paul Sicsic <paul.sicsic at polytechnique.edu>wrote:

>    SubRoutine IJacobianPoisson(dummyTS, t, U, Udot, a, Jac, PreJac, AppCtx,
> iErr)
>       Type(TS)                                     :: dummyTS
>       PetscReal                                    :: t, a
>       Type(Vec)                                    :: U, Udot
>       Type(Mat)                                    :: Jac, PreJac
>        Type(AppCtx_Type)                            :: AppCtx
>       PetscInt                                     :: iErr
>

You are missing the MatStructure argument here. Add that and your AppCtx
should come through correctly.

(This produces a warning in sane languages like C because the function type
is checked.)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111027/a64a556f/attachment.html>


More information about the petsc-dev mailing list