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