[petsc-users] snes_type aspin without DA example

Karin&NiKo niko.karin at gmail.com
Fri Jun 23 02:40:10 CDT 2023


Dear Barry,
I have started looking at the code but I miss an example using
SNESNASMSetSubdomains. In fact I do not even find a single use of the
function in PETSc.
Could someone provide me with an example ?
Thanks,
Nicolas

Le ven. 23 juin 2023 à 02:17, Barry Smith <bsmith at petsc.dev> a écrit :

>
>   A test would be great.
>
> On Jun 22, 2023, at 3:20 PM, Karin&NiKo <niko.karin at gmail.com> wrote:
>
> Thank you Barry. I will try this.
> Should I provide a test in src/binding/petsc4py/test/test_snes.py ?
>
> Le jeu. 22 juin 2023 à 20:41, Barry Smith <bsmith at petsc.dev> a écrit :
>
>>
>>   You are not missing anything. The petsc4py stub for
>> SNESNASMSetSubdomains() has not been written. You could add it by adding to
>> src/petsc4py/PETSc/SNES.pyx and src/petsc4py/PETSc/petscsnes.pxi and
>> then make a merge request
>> https://petsc.org/release/developers/contributing/ to get it into PETSc.
>>
>> On Jun 22, 2023, at 1:54 PM, Karin&NiKo <niko.karin at gmail.com> wrote:
>>
>> Dear PETSc team,
>>
>> I would like to play with aspin-type nonlinear solvers. I have found
>> several tests like snes/tutorials/ex19.c but they all use DA, which I don't
>> want to use since I need to stick at the algebraic level.
>> Then, I started looking at petsc4py/demo/ode/heat.py and tried to set up
>> things.
>> Unfortunately, I get the error "DM has no default decomposition defined.
>> Set subsolves manually with SNESNASMSetSubdomains()" which, I think, I do
>> understand.
>> But I do not find any implementation of the SNESNASMSetSubdomains in
>> petsc4py.
>> Am I missing something ?
>> Thanks,
>> Nicolas
>>
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230623/d77e49d6/attachment.html>


More information about the petsc-users mailing list