[petsc-users] PETSC MPI error on PetscGatherMessageLengths2

2015-07-14 Thread Sun, Hui
I'm running some code written by myself, using PETSC with MPI. It runs fine with less than or equal to 12 cores. However, if I ran it with 16 cores, it gives me an error. By looking at the error message, it seems that the error happens at a line where I have something like this: ierr =

Re: [petsc-users] Double free in KSPDestroy

2015-07-14 Thread Florian Lindner
Hello, Am Montag, 13. Juli 2015, 12:26:21 schrieb Barry Smith: Run under valgrind first, see if it gives any more details about the memory issue http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind I tried running it like that: valgrind --tool=memcheck ./pmpi A -malloc off

[petsc-users] A combined petsc4py, FEniCS, and Firedrake issue

2015-07-14 Thread Justin Chang
Hi everyone, As the title says, I have an issue regarding all three. Not sure which mailing to post this to but i figured somebody here may help. I was attempting to install Firedrake on my Ubuntu PC, so far it has not worked yet. However, after my failed attempt, i tried to run some of my old

Re: [petsc-users] A combined petsc4py, FEniCS, and Firedrake issue

2015-07-14 Thread Justin Chang
Hi Lawrence, I just resolved this issue, you're right I did have two petsc4py's (manually removing the one in /usr/local/... did the trick) and I undid every installation step i took for Firedrake. Seems the issue was me installing FEniCS with sudo apt-get install fenics and installing Firedrake

Re: [petsc-users] A combined petsc4py, FEniCS, and Firedrake issue

2015-07-14 Thread Lawrence Mitchell
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi Justin, On 14/07/15 09:24, Justin Chang wrote: Hi everyone, As the title says, I have an issue regarding all three. Not sure which mailing to post this to but i figured somebody here may help. I was attempting to install Firedrake on my

[petsc-users] Question about Nested matrices

2015-07-14 Thread Alejandro D Otero
Hi all, I am trying to nest to square matrices to form a rectangular one with dims 2n x n (with n the size of the square matrices) I want the two matrices to have interlaced rows so I passed the following 2 list of indices to MatCreateNest: [[ 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34

Re: [petsc-users] PETSC MPI error on PetscGatherMessageLengths2

2015-07-14 Thread Hong
Hui : It is impossible to know what causes the problem from these error messages. I notice that you use Petsc-3.5.2. I would suggest 1. update to the latest petsc-3.6; 2. if problem remains, send us your matrices C and B in petsc binary format, we can investigate. With matrix files, one can test

[petsc-users] How to disable KSP message

2015-07-14 Thread Jin, Shuangshuang
Hi, does anyone know what is the command line option to disable the following printout message about KSP convergence? 0: PETSc KSP converged after 1 iterations, reason: 3 0: PETSc KSP converged after 0 iterations, reason: 3 0: PETSc KSP converged after 1 iterations, reason: 3 0: PETSc KSP

Re: [petsc-users] Fieldsplit PC null pointer error on getksp

2015-07-14 Thread Sander Land
Hi Matt, Doing things to the PC like this has worked for me for a long time. It just happened to be formerly KSPPREONLY/PCLU. Setting things like GMRES/PCNONE works as well with no errors, though it converges far too slowly. The actual preconditioner I am trying to form is one from recent

Re: [petsc-users] Fieldsplit PC null pointer error on getksp

2015-07-14 Thread Matthew Knepley
On Tue, Jul 14, 2015 at 11:37 AM, Sander Land sander.l...@gmail.com wrote: Hi Matt, Doing things to the PC like this has worked for me for a long time. It just happened to be formerly KSPPREONLY/PCLU. Setting things like GMRES/PCNONE works as well with no errors, though it converges far too

Re: [petsc-users] How to disable KSP message

2015-07-14 Thread Matthew Knepley
I believe this comes from your code rather than PETSc. Matt On Tue, Jul 14, 2015 at 10:55 AM, Jin, Shuangshuang shuangshuang@pnnl.gov wrote: Hi, does anyone know what is the command line option to disable the following printout message about KSP convergence? 0: PETSc KSP

Re: [petsc-users] How to disable KSP message

2015-07-14 Thread Abhyankar, Shrirang G.
This seems to be a result of some print statement in GridPACK, not PETSc. Shri From: Jin, Shuangshuang shuangshuang@pnnl.govmailto:shuangshuang@pnnl.gov Date: Tuesday, July 14, 2015 at 10:55 AM To: petsc-users@mcs.anl.govmailto:petsc-users@mcs.anl.gov

Re: [petsc-users] How to disable KSP message

2015-07-14 Thread Jin, Shuangshuang
Yes, found it in GridPACK's math module Petsc linear solver part. Thanks, Shuangshuang From: Abhyankar, Shrirang G. [mailto:abhy...@anl.gov] Sent: Tuesday, July 14, 2015 10:00 AM To: Jin, Shuangshuang; petsc-users@mcs.anl.gov Cc: Palmer, Bruce J Subject: Re: [petsc-users] How to disable KSP

Re: [petsc-users] Fwd: Fwd: same petsc problem

2015-07-14 Thread Barry Smith
Version should not matter. Do other command line options work? You can also put PetscOptionsSetValue(-ksp_view_mat,binary:filename); in the code just before the call to KSPSolve(). Or run in the debugger and put a break point in KSPSolve and then next/step through to see why the

Re: [petsc-users] PETSc / GAMG user counts

2015-07-14 Thread Jed Brown
Mark Adams mfad...@lbl.gov writes: Do we have any data with respect to the number of PETSc users? Better yet number of GAMG users? Mailing list statistics, downstream software packages, papers, and commits. We don't spy on users, so just measure products that are made public. signature.asc

Re: [petsc-users] Question about Nested matrices

2015-07-14 Thread Jed Brown
Alejandro D Otero aot...@fi.uba.ar writes: Hi all, I am trying to nest to square matrices to form a rectangular one with dims 2n x n (with n the size of the square matrices) I want the two matrices to have interlaced rows so I passed the following 2 list of indices to MatCreateNest: [[ 0

Re: [petsc-users] Building a rectangular MPI matrix relating to PetscSections

2015-07-14 Thread Jed Brown
Alejandro D Otero aot...@fi.uba.ar writes: Hi all, is it posible to build a rectangular matrix that relates two fields with (posible) different number of DoF per point defining the nonzero structure from 2 petscsections representing one field each? There isn't a convenience function, but of

Re: [petsc-users] Nodes coordinates in distributed dmplex

2015-07-14 Thread Jed Brown
Alejandro D Otero aot...@fi.uba.ar writes: But when I set the viewer to use VTK_VTU (in this case ASCII just for make it readable, but the same in native vtk format when read with paraview) I got: The ASCII viewer is unrelated I never use it. The VTU (XML/binary-appended) format writes

Re: [petsc-users] Fieldsplit PC null pointer error on getksp

2015-07-14 Thread Barry Smith
On Jul 14, 2015, at 11:37 AM, Sander Land sander.l...@gmail.com wrote: Hi Matt, Doing things to the PC like this has worked for me for a long time. It just happened to be formerly KSPPREONLY/PCLU. Setting things like GMRES/PCNONE works as well with no errors, though it converges far