Re: [petsc-users] block matrix without MatCreateNest

2016-09-14 Thread Klaij, Christiaan
c-users@mcs.anl.gov; Jed Brown > > Subject: Re: [petsc-users] block matrix without MatCreateNest > > > > On Mon, Aug 1, 2016 at 9:36 AM, Klaij, Christiaan <c.kl...@marin.nl> > wrote: > > > > Matt, > > > > > > 1) great! > > > >

Re: [petsc-users] block matrix without MatCreateNest

2016-08-04 Thread Matthew Knepley
:36 -0500 > > From: Matthew Knepley <knep...@gmail.com> > > To: "Klaij, Christiaan" <c.kl...@marin.nl> > > Cc: "petsc-users@mcs.anl.gov" <petsc-users@mcs.anl.gov> > > Subject: Re: [petsc-users] block matrix without MatCreateNest > > Message-ID: > >

Re: [petsc-users] block matrix without MatCreateNest

2016-08-04 Thread Klaij, Christiaan
-users@mcs.anl.gov> > Subject: Re: [petsc-users] block matrix without MatCreateNest > Message-ID: >

Re: [petsc-users] block matrix without MatCreateNest

2016-08-02 Thread Matthew Knepley
I know who needed it had explicitly coded it in rather than trying to use it from options. It should not take long to get this fixed. Thanks, Matt > Chris > > > > From: Matthew Knepley <knep...@gmail.com> > > Sent: Tuesday, August 02, 2016 12:28 AM > > To: Kla

Re: [petsc-users] block matrix without MatCreateNest

2016-08-02 Thread Klaij, Christiaan
.@gmail.com> > Sent: Tuesday, August 02, 2016 12:28 AM > To: Klaij, Christiaan > Cc: petsc-users@mcs.anl.gov; Jed Brown > Subject: Re: [petsc-users] block matrix without MatCreateNest > > On Mon, Aug 1, 2016 at 9:36 AM, Klaij, Christiaan <c.kl...@marin.nl> wrote: > >

Re: [petsc-users] block matrix without MatCreateNest

2016-08-01 Thread Matthew Knepley
-- > *From:* Matthew Knepley <knep...@gmail.com> > *Sent:* Monday, August 01, 2016 4:15 PM > *To:* Klaij, Christiaan > *Cc:* petsc-users@mcs.anl.gov > *Subject:* Re: [petsc-users] block matrix without MatCreateNest > > On Mon, Aug 1, 2016 at 8:59 AM, Kla

Re: [petsc-users] block matrix without MatCreateNest

2016-08-01 Thread Klaij, Christiaan
___ From: Matthew Knepley <knep...@gmail.com> Sent: Monday, August 01, 2016 4:15 PM To: Klaij, Christiaan Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] block matrix without MatCreateNest On Mon, Aug 1, 2016 at 8:59 AM, Klaij, Christiaan <c.kl...@marin.nl<mailto:

Re: [petsc-users] block matrix without MatCreateNest

2016-08-01 Thread Matthew Knepley
On Mon, Aug 1, 2016 at 8:59 AM, Klaij, Christiaan wrote: > Matt, > > Thanks for your suggestions. Here's the outcome: > > 1) without the "if type=nest" protection, I get a "Cannot > locate function MatNestSetSubMats_C" error when using > type mpiaij, see below. > That is a

Re: [petsc-users] block matrix without MatCreateNest

2016-08-01 Thread Klaij, Christiaan
Matt, Thanks for your suggestions. Here's the outcome: 1) without the "if type=nest" protection, I get a "Cannot locate function MatNestSetSubMats_C" error when using type mpiaij, see below. 2) with the isg in a proper array, I get the same "Invalid Pointer to Object" error, see below. Chris

Re: [petsc-users] block matrix without MatCreateNest

2016-08-01 Thread Matthew Knepley
On Mon, Aug 1, 2016 at 3:00 AM, Klaij, Christiaan wrote: > Matt, Barry > > Thanks for your replies! I've added a call to MatNestSetSubMats() > but something's still wrong, see below. > > Chris > > > $ cat mattry.F90 > program mattry > > use petscksp > implicit none >

Re: [petsc-users] block matrix without MatCreateNest

2016-08-01 Thread Klaij, Christiaan
Matt, Barry Thanks for your replies! I've added a call to MatNestSetSubMats() but something's still wrong, see below. Chris $ cat mattry.F90 program mattry use petscksp implicit none #include PetscInt :: n=4 ! setting 4 cells per process PetscErrorCode :: ierr PetscInt

Re: [petsc-users] block matrix without MatCreateNest

2016-07-30 Thread Klaij, Christiaan
Anyone? (my guess is an if-statement, something like "if type nest then setup nest"...) > Date: Thu, 28 Jul 2016 08:38:54 + > From: "Klaij, Christiaan" <c.kl...@marin.nl> > To: "petsc-users@mcs.anl.gov" <petsc-users@mcs.anl.gov> > Sub

[petsc-users] block matrix without MatCreateNest

2016-07-28 Thread Klaij, Christiaan
I'm trying to understand how to assemble a block matrix in a format-independent manner, so that I can switch between types mpiaij and matnest. The manual states that the key to format-independent assembly is to use MatGetLocalSubMatrix. So, in the code below, I'm using this to assemble a 3-by-3