[petsc-dev] PCREDUNDANT

2019-07-27 Thread Pierre Jolivet via petsc-dev
Hello, I’m facing multiple issues with PCREDUNDANT and MATMPISBAIJ: 1) https://www.mcs.anl.gov/petsc/petsc-current/src/mat/impls/sbaij/mpi/mpisbaij.c.html#line3354 shouldn’t this be sum != N? I’m runn

[petsc-dev] PCRedundant problem

2012-10-15 Thread Stefano Zampini
Hi, is it possible to use PCRedundant with PETSC_SUBCOMM_CONTIGUOUS as subcomm type? Regards, -- Stefano -- next part -- An HTML attachment was scrubbed... URL:

[petsc-dev] PCRedundant problem

2012-10-15 Thread Hong Zhang
Stefano : > > is it possible to use PCRedundant with PETSC_SUBCOMM_CONTIGUOUS as subcomm > type? > I think "PETSC_SUBCOMM_CONTIGUOUS" is the default for PCRedundant. Hong -- next part -- An HTML attachment was scrubbed... URL:

[petsc-dev] PCRedundant problem

2012-10-15 Thread Stefano Zampini
Hong, in petsc-dev, default is interlaced. I ran some tests, I got malloc errors when inserting mat values using subcomm of contiguous type. I think that PCRedundant with PETSC_SUBCOMM_CONTIGUOUS is buggy. 2012/10/15 Hong Zhang > Stefano : > >> >> is it possible to use PCRedundant with PETSC_SUB

[petsc-dev] PCRedundant problem

2012-10-15 Thread Stefano Zampini
Hong, INTERLACED is the default since http://petsc.cs.iit.edu/petsc/petsc-dev/rev/67df4fdba5bd See also http://petsc.cs.iit.edu/petsc/petsc-dev/rev/c04222d17a3d 2012/10/15 Stefano Zampini > Hong, in petsc-dev, default is interlaced. I ran some tests, I got malloc > errors when inserting mat va

[petsc-dev] PCRedundant problem

2012-10-15 Thread Hong Zhang
Stefano : > Hong, INTERLACED is the default since > http://petsc.cs.iit.edu/petsc/petsc-dev/rev/67df4fdba5bd > > See also http://petsc.cs.iit.edu/petsc/petsc-dev/rev/c04222d17a3d Hmm, you are correct :-( > Hong, in petsc-dev, default is interlaced. I ran some tests, I got malloc >> errors when

[petsc-dev] PCRedundant problem

2012-10-16 Thread Hong Zhang
Stefano, Checking PCRedundant, I see that current PCRedundant is hard-coded using PETSC_SUBCOMM_INTERLACED which 'avoid data scattering from subcomm back to original comm' See ~petsc/src/sys/objects/subcomm.c - PetscSubcommCreate_interlaced() Unless you write a new implementation, PETSC_SUBCOMM_C

Re: [petsc-dev] PCREDUNDANT

2019-07-28 Thread Mark Adams via petsc-dev
On Sun, Jul 28, 2019 at 2:54 AM Pierre Jolivet via petsc-dev < petsc-dev@mcs.anl.gov> wrote: > Hello, > I’m facing multiple issues with PCREDUNDANT and MATMPISBAIJ: > 1) > https://www.mcs.anl.gov/petsc/petsc-current/src/mat/impls/sbaij/mpi/mpisbaij.c.html#line3354 > shouldn’t > this be sum != N?

Re: [petsc-dev] PCREDUNDANT

2019-07-28 Thread Pierre Jolivet via petsc-dev
> On 28 Jul 2019, at 3:02 PM, Mark Adams wrote: > > > On Sun, Jul 28, 2019 at 2:54 AM Pierre Jolivet via petsc-dev > mailto:petsc-dev@mcs.anl.gov>> wrote: > Hello, > I’m facing multiple issues with PCREDUNDANT and MATMPISBAIJ: > 1) > https://www.mcs.anl.gov/petsc/petsc-current/src/mat/impls/