Re: [petsc-users] TS and petscFE

[Matthew Knepley]

On Tue, Aug 2, 2016 at 8:22 AM, Maximilian Hartig <imilian.har...@gmail.com>
wrote:

> Hello all,
>
> I would like to run a transient problem with PetscFE. Example ex11.c seems
> relevant since it uses the PestcFV context to create boundary conditions
> and RHS Functions for the TS.
> Is there an easy way to do transient analysis with TS and petscFE or do I
> have to code my own time-stepping routine?
>

You can use

    ierr = DMTSSetBoundaryLocal(adaptedDM,  DMPlexTSComputeBoundary,
user);CHKERRQ(ierr);
    ierr = DMTSSetIFunctionLocal(adaptedDM, DMPlexTSComputeIFunctionFEM,
user);CHKERRQ(ierr);
    ierr = DMTSSetIJacobianLocal(adaptedDM, DMPlexTSComputeIJacobianFEM,
user);CHKERRQ(ierr);

I have been meaning to write a heat equation example, but I have not
finished yet,

  Thanks,

     Matt


> Thanks,
> Max




-- 
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