Line 1528 of dmplexsnes.c ? On Fri, Apr 17, 2020 at 10:41 AM Matthew Knepley <knep...@gmail.com> wrote:
> Tried it again and works for me. > > Matt > > On Fri, Apr 17, 2020 at 10:39 AM Mark Adams <mfad...@lbl.gov> wrote: > >> Matt, that kink seems messed up. >> >> On Fri, Apr 17, 2020 at 9:18 AM Matthew Knepley <knep...@gmail.com> >> wrote: >> >>> Mark is right, you could just add in a source term for each cell this >>> way. >>> >>> If you really wanted a lower level interface, it would go here I think: >>> >>> >>> https://gitlab.com/petsc/petsc/-/blob/master/src/snes/utils/dmplexsnes.c#L1528 >>> >>> What would you want the interface to look like? >>> >>> Thanks, >>> >>> Matt >>> >>> On Fri, Apr 17, 2020 at 7:55 AM Mark Adams <mfad...@lbl.gov> wrote: >>> >>>> You can write your own DMPlexTSComputeRHSFunctionFVM method like this >>>> and give that to the dm. >>>> >>>> PetscErrorCode foo(DM dm, PetscReal time, Vec locX, Vec F, void *user) >>>> { >>>> PetscErrorCode ierr; >>>> PetscFunctionBegin; >>>> ierr = DMPlexTSComputeRHSFunctionFVM(dm, time, locX, F, >>>> user);CHKERRQ(ierr); >>>> ... >>>> >>>> PetscFunctionReturn(0); >>>> } >>>> >>>> >>>> On Fri, Apr 17, 2020 at 5:28 AM MUKKUND SUNJII <mukkundsun...@gmail.com> >>>> wrote: >>>> >>>>> Greetings, >>>>> >>>>> I had been working on ts/tutorials/ex11.c as part of my master’s >>>>> thesis. >>>>> >>>>> In ex11.c, there is a 2D Shallow Water Model. I want to modify the >>>>> model by adding a source term to the model. This source term would factor >>>>> in the effects of bathymetry or bed elevation in addition to water height. >>>>> >>>>> I know how I will be able to achieve this if just use >>>>> TSSetRHSFunction(). However ex11.c makes use of the Riemann Solver for the >>>>> computation fluxes at the interfaces (therefore, it uses >>>>> DMTSSetRHSFunctionLocal(dm, DMPlexTSComputeRHSFunctionFVM, user) ). >>>>> >>>>> Is there a routine that allows me to *add* a source term *for every >>>>> cell* to the RHS function? >>>>> >>>>> Thank you in advance! >>>>> >>>>> Regards, >>>>> >>>>> Mukkund >>>>> >>>> >>> >>> -- >>> 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 >>> >>> https://www.cse.buffalo.edu/~knepley/ >>> <http://www.cse.buffalo.edu/~knepley/> >>> >> > > -- > 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 > > https://www.cse.buffalo.edu/~knepley/ > <http://www.cse.buffalo.edu/~knepley/> >
