Re: [petsc-users] pc gamg did not converge in sinv where it used to

Mark, Look at the PetscFree() macro; it does set the pointer to null after freeing. Denis, Valgrind the code next > On Feb 7, 2017, at 4:57 PM, Mark Adams wrote: > > This is failing on a test for a data cache being NULL. It is not. It is in > the reset routine,

[petsc-users] pc gamg did not converge in sinv where it used to

Dear all, I rerun calculations (unit tests) which used to work with slightly older versions of PETSc/SLEPc (a year ago, or so) and see the "KSPSolve has not converged” error for shift and invert transformation with gamg preconditioner (below). Shifted matrix is SPD, could have bad condition

Re: [petsc-users] Replacement for Euclid ?

Hi Michel, I am not entirely sure if there is a PETSc version, but there is a preconditioner called DS factorization (details in High-Performance Scientific Computing: Algorithms and Applications). It shares few properties of Euclid (like being quite scalable). If you have to use a factorization

Re: [petsc-users] Replacement for Euclid ?

> On Feb 7, 2017, at 2:49 AM, Michel Kern wrote: > > I have a code that uses euclid for ILU preconditioning (as part of a CPR-AMG > preconditioner for a mulitphase flow solver). As noted in the release notes, > euclid is no longer supported starting from 3.6. What is the

Re: [petsc-users] missing types in petscsysdef.h

PETSC_USE_REAL_DOUBLE should be defined in petscconf.h You did mention earlier that you found it in configure.log. Please verify if this exists in petscconf.h for your install. And how are you compiling your code? Not using petsc makefiles? Perhaps your makefiles are picking up the wrong

Re: [petsc-users] Where to restrict MPI Communicator

On Tue, Feb 7, 2017 at 8:06 AM, Florian Lindner wrote: > Hello, > > so far we have happily set PETSC_COMM_WORLD anywhere in our code and we > were not having trouble with that. > > However http://www.mcs.anl.gov/petsc/petsc-current/docs/ >

[petsc-users] Where to restrict MPI Communicator

Hello, so far we have happily set PETSC_COMM_WORLD anywhere in our code and we were not having trouble with that. However http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Sys/PETSC_COMM_WORLD.html states that: "Notes: By default PETSC_COMM_WORLD and MPI_COMM_WORLD are identical

Re: [petsc-users] Replacement for Euclid ?

On Tue, Feb 7, 2017 at 7:53 AM, Michel Kern wrote: > Thanks Mat, > > according tot he hypre documentation, pilut is "no longer supported by the > hypre team", as less efficient than Euclid, so I wouldn't consider an > option. I wsa looking for a "natural" upgrade path. > I

Re: [petsc-users] Replacement for Euclid ?

Thanks Mat, according tot he hypre documentation, pilut is "no longer supported by the hypre team", as less efficient than Euclid, so I wouldn't consider an option. I wsa looking for a "natural" upgrade path. The problems we've looked at for the moment can be up to a few millions cells

Re: [petsc-users] missing types in petscsysdef.h

On Tue, Feb 7, 2017 at 6:09 AM, Lukas van de Wiel < lukas.drinkt.t...@gmail.com> wrote: > Hi Matt, > > thank you for your fast reply. > > I have tried this with the dev version. > The structure of the include files is a lot more transparent, now! Thanks! > > However, PetscReal is still not

Re: [petsc-users] missing types in petscsysdef.h

Hi Matt, thank you for your fast reply. I have tried this with the dev version. The structure of the include files is a lot more transparent, now! Thanks! However, PetscReal is still not defined. I can use my old method to change the third to last line from petscsys.h block #if defined

Re: [petsc-users] missing types in petscsysdef.h

On Tue, Feb 7, 2017 at 5:12 AM, Lukas van de Wiel < lukas.drinkt.t...@gmail.com> wrote: > Dear fellow PETSc users, > > I have been using PETSc for several years now, but I notice that after > every upgrade I have to make changes to include/finclude/petscsysdef.h > to get my Fortran code to work

Re: [petsc-users] Replacement for Euclid ?

On Tue, Feb 7, 2017 at 2:49 AM, Michel Kern wrote: > I have a code that uses euclid for ILU preconditioning (as part of a > CPR-AMG preconditioner for a mulitphase flow solver). As noted in the > release notes, euclid is no longer supported starting from 3.6. What is the >

[petsc-users] missing types in petscsysdef.h

Dear fellow PETSc users, I have been using PETSc for several years now, but I notice that after every upgrade I have to make changes to include/finclude/petscsysdef.h to get my Fortran code to work with it (or, in 3.7.5 include/petsc/finclude/petscsysdef.h) I have been running 3.4.2 for a few