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!
> >
> >
: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:
> >
-users@mcs.anl.gov>
> Subject: Re: [petsc-users] block matrix without MatCreateNest
> Message-ID:
>
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
.@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:
>
>
--
> *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
___
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:
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
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
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
>
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
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
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
13 matches
Mail list logo