Re: [petsc-users] hypre-ILU vs hypre Euclid

[Jed Brown]

It looks like Victor is working on hypre-ILU so it is active. PETSc used to 
have PILUT support, but it was so buggy/leaky that we removed the interface.

Alexander Lindsay  writes:

> Haha no I am not sure. There are a few other preconditioning options I will 
> explore before knocking on this door some more. 
>
>  On Jun 22, 2023, at 6:49 PM, Matthew Knepley  wrote:
>
>  On Thu, Jun 22, 2023 at 8:37 PM Alexander Lindsay 
>  wrote:
>
>  I know that PETSc has hooks for Euclid but I discovered today that it does 
> not support 64 bit indices, which many MOOSE
>  applications need. This would probably be more appropriate for a hypre 
> support forum (does anyone know if such a forum
>  exists other than opening GitHub issues?), but does anyone here know what 
> the difference between hypre-ILU and
>  hypre-Euclid are? From the docs it seems they are both supposed to be 
> parallel ILU solvers.
>
>  If hypre-ILU worked with 64 bit indices (I can probably check this sifting 
> through the sources), then I would probably add
>  hooks for it in PETSc (AFAICT those don't exist at present).
>
>  My understanding was that two different people were working on them. I do 
> not know if either is still actively supported. We
>  would of course like a binding to whatever is supported.
>
>  Are you sure you want to run ILU?
>
>THanks,
>
>   Matt
>  -- 
>  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/

Re: [petsc-users] hypre-ILU vs hypre Euclid

[Alexander Lindsay]

Haha no I am not sure. There are a few other preconditioning options I will explore before knocking on this door some more. On Jun 22, 2023, at 6:49 PM, Matthew Knepley  wrote:On Thu, Jun 22, 2023 at 8:37 PM Alexander Lindsay  wrote:I know that PETSc has hooks for Euclid but I discovered today that it does not support 64 bit indices, which many MOOSE applications need. This would probably be more appropriate for a hypre support forum (does anyone know if such a forum exists other than opening GitHub issues?), but does anyone here know what the difference between hypre-ILU and hypre-Euclid are? From the docs it seems they are both supposed to be parallel ILU solvers.If hypre-ILU worked with 64 bit indices (I can probably check this sifting through the sources), then I would probably add hooks for it in PETSc (AFAICT those don't exist at present).
My understanding was that two different people were working on them. I do not know if either is still actively supported. We would of course like a binding to whatever is supported.Are you sure you want to run ILU?  THanks,     Matt-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.-- Norbert Wienerhttps://www.cse.buffalo.edu/~knepley/

Re: [petsc-users] hypre-ILU vs hypre Euclid

[Matthew Knepley]

On Thu, Jun 22, 2023 at 8:37 PM Alexander Lindsay 
wrote:

> I know that PETSc has hooks for Euclid but I discovered today that it does
> not support 64 bit indices, which many MOOSE applications need. This would
> probably be more appropriate for a hypre support forum (does anyone know if
> such a forum exists other than opening GitHub issues?), but does anyone
> here know what the difference between hypre-ILU and hypre-Euclid are? From
> the docs it seems they are both supposed to be parallel ILU solvers.
>
> If hypre-ILU worked with 64 bit indices (I can probably check this sifting
> through the sources), then I would probably add hooks for it in PETSc
> (AFAICT those don't exist at present).
>

My understanding was that two different people were working on them. I do
not know if either is still actively supported. We would of course like a
binding to whatever is supported.

Are you sure you want to run ILU?

  THanks,

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

[petsc-users] hypre-ILU vs hypre Euclid

[Alexander Lindsay]

I know that PETSc has hooks for Euclid but I discovered today that it does
not support 64 bit indices, which many MOOSE applications need. This would
probably be more appropriate for a hypre support forum (does anyone know if
such a forum exists other than opening GitHub issues?), but does anyone
here know what the difference between hypre-ILU and hypre-Euclid are? From
the docs it seems they are both supposed to be parallel ILU solvers.

If hypre-ILU worked with 64 bit indices (I can probably check this sifting
through the sources), then I would probably add hooks for it in PETSc
(AFAICT those don't exist at present).

Re: [petsc-users] snes_type aspin without DA example

[Barry Smith]

  A test would be great.

> On Jun 22, 2023, at 3:20 PM, Karin  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  > 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 >> > 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
>> 

Re: [petsc-users] snes_type aspin without DA example

[Karin]

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

Re: [petsc-users] snes_type aspin without DA example

[Barry Smith]

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

[petsc-users] snes_type aspin without DA example

[Karin]

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

Re: [petsc-users] petsc configure problem

[Gael]

Thanks for your help! Using 'mpiicc, mpiicpc, mpiifort' works well! And yes, 
intel compiler is used intentionally for my installation.


Jiahong CHEN




--Original--
From:   
 "petsc-users"