[petsc-users] Modification to ex11.c

[Mark Adams]

Line 1528 of dmplexsnes.c ?

On Fri, Apr 17, 2020 at 10:41 AM Matthew Knepley <knepley at gmail.com> wrote:

> Tried it again and works for me.
>
>    Matt
>
> On Fri, Apr 17, 2020 at 10:39 AM Mark Adams <mfadams at lbl.gov> wrote:
>
>> Matt, that kink seems messed up.
>>
>> On Fri, Apr 17, 2020 at 9:18 AM Matthew Knepley <knepley at 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 <mfadams at 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 <mukkundsunjii at 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/>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200417/9e21f1af/attachment.html>