[petsc-users] Meaning of error message (gamg & fieldsplit related)

Hello, I am running some code that employs gamg preconditioning within a fieldsplit, and for sufficiently large/refined meshes, I am getting the following error:

Re: [petsc-users] Nested Fieldsplit for custom index sets

Lawrence, Thank you for your help, the program works perfectly now. Artur From: Lawrence Mitchell Sent: Thursday, July 28, 2016 3:35 AM To: Safin, Artur; Barry Smith Cc: petsc-users@mcs.anl.gov Subject: Re:

Re: [petsc-users] block matrix without MatCreateNest

On Mon, Aug 1, 2016 at 9:36 AM, Klaij, Christiaan wrote: > Matt, > > > 1) great! > > > 2) ??? that's precisely why I paste the output of "cat mattry.F90" in > the emails, so you have a small example that produces the errors I mention. > Now I'm also attaching it to this email.

[petsc-users] Segmentation faults: Derived types

Hello there, I'm having problems defining some variables into derived types in Fortran. Before, I had a similar problems with an allocatable array "PetsInt" but I solved it just doing a non-collective Petsc Vec. Today I'm having troubles with "PetscBool" or "Logical": In a module which define

Re: [petsc-users] vec norm for local portion of a vector

On Fri, Jul 29, 2016 at 10:41 AM, Barry Smith wrote: > > > On Jul 27, 2016, at 4:42 PM, Xiangdong wrote: > > > > Hello everyone, > > > > I have a global dmda vector vg. On each processor, if I want to know the > norm of local portion of vg, which function

Re: [petsc-users] Petsc mesh scalability issue with iterative solver and direct solver

On Mon, Aug 1, 2016 at 12:52 PM, Jinlei Shen wrote: > Hi Barry, > > Thanks for your reply. > > Firstly, as you suggested, I checked my program under valgrind. The > results for both sequential and parallel cases showed there are no memory > errors detected. > > Second, I coded a

Re: [petsc-users] Petsc mesh scalability issue with iterative solver and direct solver

Hi Barry, Thanks for your reply. Firstly, as you suggested, I checked my program under valgrind. The results for both sequential and parallel cases showed there are no memory errors detected. Second, I coded a sequential program without using PETSC to generate the global matrix of small mesh

Re: [petsc-users] block matrix without MatCreateNest

Matt, 1) great! 2) ??? that's precisely why I paste the output of "cat mattry.F90" in the emails, so you have a small example that produces the errors I mention. Now I'm also attaching it to this email. Thanks, Chris dr. ir. Christiaan Klaij | CFD Researcher | Research & Development

Re: [petsc-users] block matrix without MatCreateNest

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

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

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] Ignore command line arguments with fortran code using PETSc

I made a pull request. Cheers, Volker 2016-07-29 16:49 GMT+02:00 Barry Smith : > >> On Jul 28, 2016, at 2:35 AM, Tim Steinhoff wrote: >> >> 2016-07-27 21:42 GMT+02:00 Barry Smith : >>> >>> Actually there is currently no way to

Re: [petsc-users] block matrix without MatCreateNest

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] Ignore command line arguments with fortran code using PETSc

Okay, I'll let you know when I make the pull request. Thanks. 2016-07-29 16:49 GMT+02:00 Barry Smith : > >> On Jul 28, 2016, at 2:35 AM, Tim Steinhoff wrote: >> >> 2016-07-27 21:42 GMT+02:00 Barry Smith : >>> >>> Actually there is