Re: [petsc-users] How to know the cpu rank that hold the minimum vec values

2017-03-23 Thread Barry Smith
> On Mar 23, 2017, at 8:31 PM, 丁老师 wrote: > > Dear professor: > I got the location and value of the minimum in a vec with VecMin, But > how to know the rank of CPU which holds the minimum values. > This rank should be known in the whole communicator. You can call

[petsc-users] How to know the cpu rank that hold the minimum vec values

2017-03-23 Thread 丁老师
Dear professor: I got the location and value of the minimum in a vec with VecMin, But how to know the rank of CPU which holds the minimum values. This rank should be known in the whole communicator. Regards

Re: [petsc-users] coloring algorithms

2017-03-23 Thread Barry Smith
> On Mar 23, 2017, at 5:35 PM, Kong, Fande wrote: > > Thanks, Barry, > > On Thu, Mar 23, 2017 at 4:02 PM, Barry Smith wrote: > > Please send the matrix as a binary file. > >Are you computing a distance one coloring or distance 2. 2 is needed for

Re: [petsc-users] coloring algorithms

2017-03-23 Thread Kong, Fande
Thanks, Barry, On Thu, Mar 23, 2017 at 4:02 PM, Barry Smith wrote: > > Please send the matrix as a binary file. > >Are you computing a distance one coloring or distance 2. 2 is needed > for Jacobians. > The matrix does not come from PDE, and it is from a

Re: [petsc-users] coloring algorithms

2017-03-23 Thread Barry Smith
Please send the matrix as a binary file. Are you computing a distance one coloring or distance 2. 2 is needed for Jacobians. > On Mar 23, 2017, at 4:57 PM, Kong, Fande wrote: > > Thanks, Hong, > > I did some tests with a matrix (40x40): > > row 0: (0, 1.) (2, 1.)

Re: [petsc-users] coloring algorithms

2017-03-23 Thread Kong, Fande
Thanks, Hong, I did some tests with a matrix (40x40): *row 0: (0, 1.) (2, 1.) (3, 1.) (11, 1.) (14, 1.) (15, 1.) (19, 1.) (22, 1.) (23, 1.) (24, 1.) (27, 1.) (28, 1.) row 1: (1, 1.) (2, 1.) (3, 1.) (6, 1.) (16, 1.) (17, 1.) (18, 1.) (21, 1.)

Re: [petsc-users] Update to MUMPS 5.1.1

2017-03-23 Thread Satish Balay
On Thu, 23 Mar 2017, Tim Steinhoff wrote: > Hi all, > > there is new version of MUMPS with some nice improvements: > http://mumps.enseeiht.fr/index.php?page=dwnld#cl > Is it possible to update the package in the petsc repository? I have the changes in git branch 'balay/update-mumps-5.1.1'.

Re: [petsc-users] How to use f2py on a PETSc/SLEPc-based fortran code

2017-03-23 Thread Gaetan Kenway
We do this all the time. The trick is that you can't use f2py do actually do *any* of the compiling/linking for you. You just have to use it to get the module.c file and .f90 file. Then compile that yourself and use whatever compile flags, linking etc you would normally use to get a .so. We've had

Re: [petsc-users] How to use f2py on a PETSc/SLEPc-based fortran code

2017-03-23 Thread Austin Herrema
Sorry to self-bump this. Sounds like there are no quick and easy answers to get this solved, which I understand. I just wanted to say that if anybody has any suggestions about what I should be looking into specifically (C compiler? f2py? etc...) then even that would be helpful as I'm pretty stuck

[petsc-users] Update to MUMPS 5.1.1

2017-03-23 Thread Tim Steinhoff
Hi all, there is new version of MUMPS with some nice improvements: http://mumps.enseeiht.fr/index.php?page=dwnld#cl Is it possible to update the package in the petsc repository? Thanks and kind regards, Volker

Re: [petsc-users] coloring algorithms

2017-03-23 Thread Hong
> > Fande, > > I was wondering if the coloring approaches listed online are working? > Which ones are in parallel, and which ones are in sequential? > > http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/ > MatColoringType.html#MatColoringType > JP and Greedy are parallel. > > If

[petsc-users] coloring algorithms

2017-03-23 Thread Kong, Fande
Hi All, I was wondering if the coloring approaches listed online are working? Which ones are in parallel, and which ones are in sequential? http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatColoringType.html#MatColoringType If the coloring is in parallel, can it be used with

Re: [petsc-users] left and right preconditioning with a constant null space

2017-03-23 Thread Lawrence Mitchell
On 23/03/17 15:37, Klaij, Christiaan wrote: > Yes, that's clearer, thanks! I do have is0 and is1 so I can try > PetscObjectCompose and let you know. > > Note though that the viewer reports that both S and A11 have a > null space attached... My matrix is a matnest and I've attached a > null

Re: [petsc-users] left and right preconditioning with a constant null space

2017-03-23 Thread Klaij, Christiaan
Lawrence, Yes, that's clearer, thanks! I do have is0 and is1 so I can try PetscObjectCompose and let you know. Note though that the viewer reports that both S and A11 have a null space attached... My matrix is a matnest and I've attached a null space to A11, so the latter works as expected. But

Re: [petsc-users] About error

2017-03-23 Thread Matthew Knepley
On Thu, Mar 23, 2017 at 10:43 AM, Валентин Егоров wrote: > Hello! > My name is Valentin Egorov. I am from Russia. And I have a question for > you about PETSC. > > You recommended to add use petsc to .#nclude > "petsc/finclude/petscXXX.h", and for XXX

Re: [petsc-users] left and right preconditioning with a constant null space

2017-03-23 Thread Lawrence Mitchell
On 23/03/17 08:42, Klaij, Christiaan wrote: > Matt, Lawrence > > > The same problem happens when using gmres with rtol 1e-6 in the > schur complement (attachment "left_schur"). I'm not sure what > this tells us. If I understand Lawrence correctly, the null space > may be attached to the wrong

[petsc-users] About error

2017-03-23 Thread Валентин Егоров
Hello! My name is Valentin Egorov. I am from Russia. And I have a question for you about PETSC. You recommended to add use petsc to .#nclude "petsc/finclude/petscXXX.h", and for XXX variablename - type(tXXX) variablename for Fortran. However, I have "Error in opening the Library