[petsc-dev] explicit FEM

Matthew Knepley knepley at gmail.com
Tue May 16 19:41:12 CDT 2017


On Tue, May 16, 2017 at 7:22 PM, Mark Adams <mfadams at lbl.gov> wrote:

> I am having problems with explicit + FEM. I set the DM with something like:
>
>     ierr = DMTSSetRHSFunctionLocal(dm, DMPlexTSComputeRHSFunctionFVM,
> &ctx);CHKERRQ(ierr);
>
> But DMPlexTSComputeRHSFunctionFVM calls DMPlexComputeResidual_Internal
> without a "locX_t" so there is no time derivative and so the FEM point
> functions do not get a time derivative (u_t).
>
> I'm not sure if I am missing something or if PETSc is missing this
> functionality.
>

Its missing. We will have to put in a  DMPlexTSComputeRHSFunctionFEM().

   Matt

Mark
>
-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20170516/d31b93cf/attachment.html>


More information about the petsc-dev mailing list