[petsc-users] Valgrind Errors with KSPSolve

2018-10-26 Thread Phil Tooley
insights are appreciated. Many Thanks Phil Tooley Typical errors are as follows: ==7368== Conditional jump or move depends on uninitialised value(s) ==7368==    at 0x58894EE: PetscAbsScalar(double) (petscmath.h:389) ==7368==    by 0x5891702: MatLUFactorNumeric_SeqAIJ (aijfact.c:553) ==7368

Re: [petsc-users] Valgrind Errors with KSPSolve

2018-10-26 Thread Phil Tooley
help. Phil On 26/10/2018 15:01, Matthew Knepley wrote: > On Fri, Oct 26, 2018 at 9:54 AM Phil Tooley > mailto:phil.too...@sheffield.ac.uk>> wrote: > > Hi All, > > Running valgrind over my code reveals a huge number (it stops > collecting > after 1e6)

Re: [petsc-users] DMDA with dimension of size 1

2018-10-01 Thread Phil Tooley
Message > -- > [0]PETSC ERROR: Overflow in integer operation: > [0]PETSC ERROR: Mesh of 10 by 10 by 1 (dof) is too large for > 32 bit indices > > Is there a workaround other than "

[petsc-users] DMDA with dimension of size 1

2018-10-01 Thread Phil Tooley
;? Thanks Phil -- Phil Tooley Research Software Engineering University of Sheffield

[petsc-users] Better understanding VecScatter

2018-10-04 Thread Phil Tooley
please? Many Thanks Phil -- Phil Tooley Research Software Engineering University of Sheffield

[petsc-users] DMCreateXXXVector vs VecDuplicate

2018-09-24 Thread Phil Tooley
class all with the same distributed layout. It makes life much tidier if I can just pass the DMDA rather than also having to pass vectors to be copied. Cheers Phil -- Phil Tooley Research Software Engineering University of Sheffield

[petsc-users] Checking if a vector is a localvector of a given DMDA

2018-09-25 Thread Phil Tooley
that do what I need, but I am likely missing something obvious. Is there a way to achieve the above? Thanks Phil -- Phil Tooley Research Software Engineering University of Sheffield

Re: [petsc-users] Checking if a vector is a localvector of a given DMDA

2018-09-25 Thread Phil Tooley
with garbage data. Phil On 25/09/18 11:55, Matthew Knepley wrote: > On Tue, Sep 25, 2018 at 6:49 AM Phil Tooley > mailto:phil.too...@sheffield.ac.uk>> wrote: > > Hi all, > > Given a vector I know I can get an associated DM (if there is one) by > calling VecGetDM,

Re: [petsc-users] Checking if a vector is a localvector of a given DMDA

2018-09-25 Thread Phil Tooley
n Tue, 25 Sep 2018 at 13:20, Matthew Knepley <mailto:knep...@gmail.com>> wrote: > > On Tue, Sep 25, 2018 at 7:03 AM Dave May <mailto:dave.mayhe...@gmail.com>> wrote: > > On Tue, 25 Sep 2018 at 11:49, Phil Tooley > &

Re: [petsc-users] Indexing/using a 3D DMDA as a 1D vector

2018-09-12 Thread Phil Tooley
/09/18 14:27, Phil Tooley wrote: > I will preface this by saying I am new to PETSc and am still trying to > get my head around all of the layout mapping that is done.  That means I > may well have fundamentally misunderstood something, but hopefully > someone will be able to to

[petsc-users] Indexing/using a 3D DMDA as a 1D vector

2018-09-12 Thread Phil Tooley
point me in the correct direction please? Many Thanks -- Phil Tooley Research Software Engineering University of Sheffield