Re: [petsc-users] DIVERGED_NANORING with PC GAMG

2018-11-05 Thread Mark Adams via petsc-users
On Mon, Nov 5, 2018 at 4:11 PM Appel, Thibaut wrote: > "Local" as in serial? > Block Jacobi with ILU as the solver on each block. Each block corresponds to an MPI process by default. So it is completely parallel it is just not a true ILU. I the limit of one equation per processor it is just

Re: [petsc-users] DIVERGED_NANORING with PC GAMG

2018-11-05 Thread Appel, Thibaut via petsc-users
"Local" as in serial? Thibaut On 5 Nov 2018, at 20:12, Mark Adams mailto:mfad...@lbl.gov>> wrote: On Mon, Nov 5, 2018 at 12:50 PM Thibaut Appel mailto:t.appe...@imperial.ac.uk>> wrote: Hi Mark, Yes it doesn't seem to be usable. Unfortunately we're aiming to do 3D so direct solvers are

Re: [petsc-users] TAOIPM for AC optimal power flow

2018-11-05 Thread Matthew Knepley via petsc-users
On Mon, Nov 5, 2018 at 3:23 PM Justin Chang via petsc-users < petsc-users@mcs.anl.gov> wrote: > Hi everyone, > > I am working on a generic AC optimal power flow solver, and I hope to use > DMNetwork's data structure and TAO's optimization solvers for this purpose. > Last time I inquired about IPM

Re: [petsc-users] DIVERGED_NANORING with PC GAMG

2018-11-05 Thread Mark Adams via petsc-users
On Mon, Nov 5, 2018 at 12:50 PM Thibaut Appel wrote: > Hi Mark, > > Yes it doesn't seem to be usable. Unfortunately we're aiming to do 3D so > direct solvers are not a viable solution and PETSc' ILU is not parallel and > we can't use HYPRE (complex arithmetic) > I think SuperLU has a parallel

Re: [petsc-users] DIVERGED_NANORING with PC GAMG

2018-11-05 Thread Thibaut Appel via petsc-users
Hi Mark, Yes it doesn't seem to be usable. Unfortunately we're aiming to do 3D so direct solvers are not a viable solution and PETSc' ILU is not parallel and we can't use HYPRE (complex arithmetic) Thibaut On 01/11/2018 20:42, Mark Adams wrote: On Wed, Oct 31, 2018 at 8:11 PM Smith,

Re: [petsc-users] PetscViewerFileSetName() wrong error code in Fortran

2018-11-05 Thread Smith, Barry F. via petsc-users
Fixed in the branch barry/fix-fortran-petscfileviewersetname soon to be fixed in maint and then the next patch release of PETSc. Thanks for the report Barry > On Nov 5, 2018, at 7:29 AM, Tim Steinhoff via petsc-users > wrote: > > Dear PETSc Team, > > I am having the issue that

Re: [petsc-users] Vec, Mat and binaryfiles.

2018-11-05 Thread Jed Brown via petsc-users
Sal Am via petsc-users writes: > Hi, > > I am trying to solve a Ax=b complex system. the vector b and "matrix" A are > both binary and NOT created by PETSc. So I keep getting error messages that > they are not correct format when I read the files with PetscViewBinaryOpen, > after some digging it

Re: [petsc-users] Force SNES diverge

2018-11-05 Thread Matthew Knepley via petsc-users
On Mon, Nov 5, 2018 at 8:41 AM Karol Lewandowski via petsc-users < petsc-users@mcs.anl.gov> wrote: > Hi, > > I am solving a highly nonlinear problem using SNES solver. Under certain > conditions during the iterations I already know that the step will diverge > (in the next few iterations). Is

Re: [petsc-users] Problems about Assemble DMComposite Precondition Matrix

2018-11-05 Thread Jed Brown via petsc-users
Yingjie Wu via petsc-users writes: > Thank you very much for your reply. > My equation is a neutron diffusion equation with eigenvalues, which is why > I use DMConposite because there is a single non-physical field variable, > eigenvalue. I am not very familiar with FieldSplit. I can understand

Re: [petsc-users] Problems about Assemble DMComposite Precondition Matrix

2018-11-05 Thread Mark Adams via petsc-users
On Mon, Nov 5, 2018 at 10:37 AM Yingjie Wu wrote: > Thank you very much for your reply. > My equation is a neutron diffusion equation with eigenvalues, which is why > I use DMConposite because there is a single non-physical field variable, > eigenvalue. > OK, DMComposite might be your best

Re: [petsc-users] Problems about Assemble DMComposite Precondition Matrix

2018-11-05 Thread Yingjie Wu via petsc-users
Thank you very much for your reply. My equation is a neutron diffusion equation with eigenvalues, which is why I use DMConposite because there is a single non-physical field variable, eigenvalue. I am not very familiar with FieldSplit. I can understand it first. It seems not the problem of

Re: [petsc-users] Problems about Assemble DMComposite Precondition Matrix

2018-11-05 Thread Mark Adams via petsc-users
DMComposite is not very mature, the last time I checked and I don't of anyone having worked on it recently, and it is probably not what you want anyway. FieldSplit is most likely what you want. What are your equations and discretization? eg, Stokes with cell centered pressure? There are probably