Re: [petsc-users] Code break when PETSc not debugging

2018-08-14 Thread Manuel Valera
Thanks, i was able to find the bug and correct it, it is working now, I was calling a wrong DM for some DA, Regards, On Tue, Aug 14, 2018 at 2:13 PM, Jed Brown wrote: > Manuel Valera writes: > > > Thanks Jed, > > > > I got the attached, it looks is coming from one of my routines > > CorrectU4

Re: [petsc-users] Code break when PETSc not debugging

2018-08-14 Thread Jed Brown
Manuel Valera writes: > Thanks Jed, > > I got the attached, it looks is coming from one of my routines > CorrectU4Pressure.F90, what other information can i get from this log? If you compile with debug info (COPTFLAGS=-g FOPTFLAGS=-g) you'll be able to see the source line that's causing this err

Re: [petsc-users] Code break when PETSc not debugging

2018-08-14 Thread Manuel Valera
Thanks Jed, I got the attached, it looks is coming from one of my routines CorrectU4Pressure.F90, what other information can i get from this log? Thanks, On Tue, Aug 14, 2018 at 12:57 PM, Jed Brown wrote: > Manuel Valera writes: > > > Hello everyone, > > > > I am working on running part of my

Re: [petsc-users] Code break when PETSc not debugging

2018-08-14 Thread Jed Brown
Manuel Valera writes: > Hello everyone, > > I am working on running part of my code in a GPU, recently i was able to > run the whole model using one P100 GPU and one processor with good timing > but using --with-debugging=1 as configure argument, > > With this in mind i compiled PETSc in a separ

[petsc-users] Code break when PETSc not debugging

2018-08-14 Thread Manuel Valera
Hello everyone, I am working on running part of my code in a GPU, recently i was able to run the whole model using one P100 GPU and one processor with good timing but using --with-debugging=1 as configure argument, With this in mind i compiled PETSc in a separate folder with the same exact flags

[petsc-users] an error message related to nest matrix

2018-08-14 Thread Mike Wick
Hi: I have a simple code myself to learn the basic usage of nest matrix. It assembles a 3x3 nest matrix and solve with that nest matrix. I can run it without a preconditioner. Running the code with -pc_type fieldsplit just leads to an error message. I have no clue of what goes wrong. The message i

Re: [petsc-users] Max color value exeeds 65535 limit

2018-08-14 Thread Smith, Barry F.
Your "sparse" matrix is too dense to realistically use coloring to compute the Jacobian. Use instead SNESComputeJacobianDefault Barry > On Aug 14, 2018, at 9:06 AM, Florian Kauer wrote: > > Hi, > I was just testing out if my calculations also work for several 1 > variables a

Re: [petsc-users] Max color value exeeds 65535 limit

2018-08-14 Thread Matthew Knepley
On Tue, Aug 14, 2018 at 10:12 AM Florian Kauer wrote: > Hi, > I was just testing out if my calculations also work for several 1 > variables and got the following error. Is this a fundamental limit with > very large (or complex) systems of equations or does it hint to a bug/wrong > configurati

[petsc-users] Max color value exeeds 65535 limit

2018-08-14 Thread Florian Kauer
Hi, I was just testing out if my calculations also work for several 1 variables and got the following error. Is this a fundamental limit with very large (or complex) systems of equations or does it hint to a bug/wrong configuration in my program? Everything works fine for smaller systems. [