Re: [petsc-users] Trouble with ISEmbed

2022-10-21 Thread Pierre Jolivet
> On 21 Oct 2022, at 10:46 AM, TARDIEU Nicolas wrote: > > Dear Pierre, > > To complete my last post, in fact, the initial code (playing with the LGMap) > was correct. It was my test case that was wrong. > Once fixed according to your suggestion, everything turns out to be OK. > I am neverthel

Re: [petsc-users] Does flagging a matrix as symmetric improving performances?

2022-10-21 Thread Barry Smith
Set the flag, it cannot harm. I cannot see that hypre boomerAMG has options to indicate matrix symmetry to take advantage of it. Stefano or Mark would know better than I. No, do not change MPIAIJ it offers the most preconditioners. Barry > On Oct 21, 2022, at 11:36 AM, Edoardo ali

Re: [petsc-users] Does flagging a matrix as symmetric improving performances?

2022-10-21 Thread Edoardo alinovi
Hi Berry, Thank you for the hints. Actually, I am using MPIAIJ, does this mean I need to change the matrix format? Most of the time I am using CG + boomerAMG, do you think that using a sym matrix can be any good for the performance? Cheers

Re: [petsc-users] Does flagging a matrix as symmetric improving performances?

2022-10-21 Thread Barry Smith
For most solvers, just setting this flag will not, by itself, will not improve the setup time. An exception to this is PCGAMG in the latest release where it will improve the solution time. For some preconditioner situations you can use MATSBAIJ to store the matrix, this can help noticeably

[petsc-users] Does flagging a matrix as symmetric improving performances?

2022-10-21 Thread Edoardo alinovi
Hi PETSc friends, As per object, do you think that flagging a matrix as symmetric might improve setup times of the preconditioner? Thank you as always.

Re: [petsc-users] Periodic boundary condition

2022-10-21 Thread Matthew Knepley
On Thu, Oct 20, 2022 at 10:27 PM Barry Smith wrote: > > Some of the valgrind information does not appear to make sense > > PetscMemcpy() is not calling SNESSolve() so I suspect there must be some > serious corruption of something to this impossible stack trace > I ran Valgrind on it, and it lo

Re: [petsc-users] Trouble with ISEmbed

2022-10-21 Thread TARDIEU Nicolas via petsc-users
Dear Pierre, To complete my last post, in fact, the initial code (playing with the LGMap) was correct. It was my test case that was wrong. Once fixed according to your suggestion, everything turns out to be OK. I am nevertheless wondering if this IS embeding according to the global numbering sho