Re: [petsc-users] GAMG and zero pivots follow up

2015-11-10 Thread David Knezevic
On Tue, Nov 10, 2015 at 10:24 PM, Matthew Knepley wrote: > On Tue, Nov 10, 2015 at 9:21 PM, David Knezevic < > david.kneze...@akselos.com> wrote: > >> On Tue, Nov 10, 2015 at 10:00 PM, Matthew Knepley >> wrote: >> >>> On Tue, Nov 10, 2015 at 8:39 PM, David Knezevic < >>> david.kneze...@akselos.c

Re: [petsc-users] GAMG and zero pivots follow up

2015-11-10 Thread Matthew Knepley
On Tue, Nov 10, 2015 at 9:21 PM, David Knezevic wrote: > On Tue, Nov 10, 2015 at 10:00 PM, Matthew Knepley > wrote: > >> On Tue, Nov 10, 2015 at 8:39 PM, David Knezevic < >> david.kneze...@akselos.com> wrote: >> >>> I'm looking into using GAMG, so I wanted to start with a simple 3D >>> elasticit

Re: [petsc-users] GAMG and zero pivots follow up

2015-11-10 Thread David Knezevic
On Tue, Nov 10, 2015 at 10:00 PM, Matthew Knepley wrote: > On Tue, Nov 10, 2015 at 8:39 PM, David Knezevic < > david.kneze...@akselos.com> wrote: > >> I'm looking into using GAMG, so I wanted to start with a simple 3D >> elasticity problem. When I first tried this, I got the following "zero >> pi

Re: [petsc-users] GAMG and zero pivots follow up

2015-11-10 Thread Matthew Knepley
On Tue, Nov 10, 2015 at 8:39 PM, David Knezevic wrote: > I'm looking into using GAMG, so I wanted to start with a simple 3D > elasticity problem. When I first tried this, I got the following "zero > pivot" error: > > --- > > [0]P

[petsc-users] GAMG and zero pivots follow up

2015-11-10 Thread David Knezevic
I'm looking into using GAMG, so I wanted to start with a simple 3D elasticity problem. When I first tried this, I got the following "zero pivot" error: --- [0]PETSC ERROR: Zero pivot in LU factorization: http://www.mcs.anl.gov/pe

Re: [petsc-users] CG+GAMG convergence issues in GHEP Krylov-Schur for some MPI runs

2015-11-10 Thread Barry Smith
Please send me the full output. This is nuts and should be reported once we understand it better to NERSc as something to be fixed. When I pay $60 million in taxes to a computing center I expect something that works fine for free on my laptop to work also there. Barry > On Nov 10, 2015, a

Re: [petsc-users] CG+GAMG convergence issues in GHEP Krylov-Schur for some MPI runs

2015-11-10 Thread Mark Adams
valgrind on Edison does seem to give a lot of false positives. The line number are accurate (not always the case). "assert" triggers it, as does SETERRQ. On Tue, Nov 10, 2015 at 8:51 AM, Mark Adams wrote: > I ran an 8 processor job on Edison of a small code for a short run (just a > linear sol

Re: [petsc-users] petsc-users Digest, Vol 83, Issue 34

2015-11-10 Thread Åsmund Ervik
> From: TAY wee-beng > To: PETSc users list > Subject: [petsc-users] Memory usage with DMDACreate3d and > DMDAGetCorners > Message-ID: <56419dce.4090...@gmail.com> > Content-Type: text/plain; charset=utf-8; format=flowed > > Hi, > > I need a subroutine in Fortran to partition a subset of

Re: [petsc-users] Memory usage with DMDACreate3d and DMDAGetCorners

2015-11-10 Thread Matthew Knepley
On Tue, Nov 10, 2015 at 6:34 AM, TAY wee-beng wrote: > > On 10/11/2015 8:25 PM, Matthew Knepley wrote: > > On Tue, Nov 10, 2015 at 1:33 AM, TAY wee-beng wrote: > >> Hi, >> >> I need a subroutine in Fortran to partition a subset of my grid in the 3 >> x,y,z directions for MPI. I thought of using

Re: [petsc-users] Use of DMDAVecGetArrayF90 and DMDAVecRestoreArrayF90

2015-11-10 Thread Matthew Knepley
On Tue, Nov 10, 2015 at 6:30 AM, TAY wee-beng wrote: > > On 10/11/2015 8:27 PM, Matthew Knepley wrote: > > On Tue, Nov 10, 2015 at 3:27 AM, TAY wee-beng wrote: > >> Hi, >> >> Inside my subroutine, I need to access the DA variable cu_types_array >> frequently. >> >> So I need to call DMDAVecGetAr

Re: [petsc-users] Memory usage with DMDACreate3d and DMDAGetCorners

2015-11-10 Thread TAY wee-beng
On 10/11/2015 8:25 PM, Matthew Knepley wrote: On Tue, Nov 10, 2015 at 1:33 AM, TAY wee-beng > wrote: Hi, I need a subroutine in Fortran to partition a subset of my grid in the 3 x,y,z directions for MPI. I thought of using DMDACreate3d and DMDAGetCorner

Re: [petsc-users] Use of DMDAVecGetArrayF90 and DMDAVecRestoreArrayF90

2015-11-10 Thread TAY wee-beng
On 10/11/2015 8:27 PM, Matthew Knepley wrote: On Tue, Nov 10, 2015 at 3:27 AM, TAY wee-beng > wrote: Hi, Inside my subroutine, I need to access the DA variable cu_types_array frequently. So I need to call DMDAVecGetArrayF90 and DMDAVecRestoreArrayF90

Re: [petsc-users] Use of DMDAVecGetArrayF90 and DMDAVecRestoreArrayF90

2015-11-10 Thread Matthew Knepley
On Tue, Nov 10, 2015 at 3:27 AM, TAY wee-beng wrote: > Hi, > > Inside my subroutine, I need to access the DA variable cu_types_array > frequently. > > So I need to call DMDAVecGetArrayF90 and DMDAVecRestoreArrayF90 before and > after frequently. > > Is this necessary? Can I call DMDAVecGetArrayF9

Re: [petsc-users] Memory usage with DMDACreate3d and DMDAGetCorners

2015-11-10 Thread Matthew Knepley
On Tue, Nov 10, 2015 at 1:33 AM, TAY wee-beng wrote: > Hi, > > I need a subroutine in Fortran to partition a subset of my grid in the 3 > x,y,z directions for MPI. I thought of using DMDACreate3d and > DMDAGetCorners to get the starting and width of the partitioned grid. > > Because I need to par

[petsc-users] Use of DMDAVecGetArrayF90 and DMDAVecRestoreArrayF90

2015-11-10 Thread TAY wee-beng
Hi, Inside my subroutine, I need to access the DA variable cu_types_array frequently. So I need to call DMDAVecGetArrayF90 and DMDAVecRestoreArrayF90 before and after frequently. Is this necessary? Can I call DMDAVecGetArrayF90 at the start and only call DMDAVecRestoreArrayF90 towards the