Re: [petsc-users] mesh refinement PETSc3.5.4 with HYPRE

Please upgrade to the latest 3.8 version of PETSc. and then run with -fp_trap This message is usually an indication that a NaN or Inf got into the numerical computation. Using the latest PETSc will make it much easier to track down. Barry > On Sep 28, 2017, at 8:57 PM, Hao Zhang

[petsc-users] mesh refinement PETSc3.5.4 with HYPRE

I'm experiencing error when doing mesh refinement(mesh x2 for X, Y and Z direction) for 3D poisson problem. PETSc produces no errors when no mesh refinement. I've pasted some debugging info here. Please advise. [0]PETSC ERROR: - Error Message

Re: [petsc-users] Incorrect Eigenvalues when Setting KSP and PC types

Greg: > Thanks so much for the detailed response. I am glad to know the reason > behind it--hopefully we eventually figure out why the solvers have this > behavior! Hong, I really appreciate you working on a patch to throw an > error in this case. It definitely bit me and some people using my

Re: [petsc-users] Using higher order ILU preconditioner

> On Sep 28, 2017, at 2:27 PM, Rachit Prasad wrote: > > I guess I missed out an important piece of information, due to legacy code > issues, I am working with Petsc version 2.3.2. Unfortunately, it doesn't have > PCFactorSetMatOrderingType(). Thanks for the advice though. I

Re: [petsc-users] Using higher order ILU preconditioner

I guess I missed out an important piece of information, due to legacy code issues, I am working with Petsc version 2.3.2. Unfortunately, it doesn't have PCFactorSetMatOrderingType(). Thanks for the advice though. I still don't understand, why would it ask for so much memory, it is almost 16

Re: [petsc-users] Using higher order ILU preconditioner

Calling PCFactorSetMatOrderingType() (or command line option: -pc_factor_mat_ordering_type) to use RCM or 1WD usually helps me a lot. My application is based on highly-unstructured meshes. Fande, On Thu, Sep 28, 2017 at 2:53 PM, Rachit Prasad wrote: > Hi, > > I am trying to

Re: [petsc-users] How to copy vectors to/from petsc?

On Thu, 28 Sep 2017, Karl Rupp wrote: > Hi Klaus, > > there are several options, e.g.: > * > http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/VecCreateSeqWithArray.html > (sequential) For parallel

[petsc-users] Using higher order ILU preconditioner

Hi, I am trying to solve a highly ill-conditioned matrix and am unable to get convergence when using ILU(0) or Jacobi and SOR preconditioners. I tried to implement a higher order ILU by increasing the level of fill-in to 1 by calling call PCFactorSetLevels(pc,1,ierr) However, when I do that I

Re: [petsc-users] How to copy vectors to/from petsc?

Hi Klaus, there are several options, e.g.: * http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/VecCreateSeqWithArray.html (sequential) * http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/VecSetValues.html (use with multiple MPI ranks if you cannot easily