Re: [petsc-users] Understanding MatCreate bsize parameter

2015-03-27 Thread Florian Lindner
Am Donnerstag, 26. März 2015, 07:34:27 schrieb Jed Brown: Florian Lindner mailingli...@xgm.de writes: Hello, I'm using petsc with petsc4py. A matrix is created like that MPIrank = MPI.COMM_WORLD.Get_rank() MPIsize = MPI.COMM_WORLD.Get_size() print(MPI Rank = ,

Re: [petsc-users] Column #j is wrong in parallel from message Inserting a new nonzero (i, j) into matrix

2015-03-27 Thread Eric Chamberland
On 03/26/2015 09:52 PM, Barry Smith wrote: Eric, I have now updated all the standard MPI matrix types AIJ, BAIJ, SBAIJ to print the correct global indices in the error messages when a new nonzero location is generated thus making debugging this issue easier. In the branches

Re: [petsc-users] Understanding MatCreate bsize parameter

2015-03-27 Thread Matthew Knepley
On Fri, Mar 27, 2015 at 4:28 AM, Florian Lindner mailingli...@xgm.de wrote: Am Donnerstag, 26. März 2015, 07:34:27 schrieb Jed Brown: Florian Lindner mailingli...@xgm.de writes: Hello, I'm using petsc with petsc4py. A matrix is created like that MPIrank =

Re: [petsc-users] Understanding MatCreate bsize parameter

2015-03-27 Thread Florian Lindner
Am Freitag, 27. März 2015, 07:26:11 schrieb Matthew Knepley: On Fri, Mar 27, 2015 at 4:28 AM, Florian Lindner mailingli...@xgm.de wrote: Am Donnerstag, 26. März 2015, 07:34:27 schrieb Jed Brown: Florian Lindner mailingli...@xgm.de writes: Hello, I'm using petsc with

Re: [petsc-users] Understanding MatCreate bsize parameter

2015-03-27 Thread Matthew Knepley
On Fri, Mar 27, 2015 at 7:31 AM, Florian Lindner mailingli...@xgm.de wrote: Am Freitag, 27. März 2015, 07:26:11 schrieb Matthew Knepley: On Fri, Mar 27, 2015 at 4:28 AM, Florian Lindner mailingli...@xgm.de wrote: Am Donnerstag, 26. März 2015, 07:34:27 schrieb Jed Brown: Florian

[petsc-users] installation question

2015-03-27 Thread Gideon Simpson
My admin did an install of petsc on ubuntu 14.0 using the package manager. Things seem to be ok, expect I noticed the following when I compile using the petsc rules: mpicxx -o ex19.o -c -Wall -Wwrite-strings -Wno-strict-aliasing -Wno-unknown-pragmas -O -fPIC

Re: [petsc-users] Mixed matrices MatMatMult

2015-03-27 Thread Matthew Knepley
On Thu, Mar 26, 2015 at 6:08 PM, Luc Berger-Vergiat lb2...@columbia.edu wrote: The difference between IS(1), IS(2) and IS(1)',IS(2)' is that they operate on two different matrices. The former operate on A whereas the operates on S and these matrices have different sizes so it's not obvious

[petsc-users] Problem with MPI, MatAXPY and SAME_NONZERO_PATTERN

2015-03-27 Thread Klaus Kaiser
Hallo, I have a strange behavior in my code concerning the function MatAXPY. I create 3 different Matrices ierr = MatCreateBAIJ(PETSC_COMM_WORLD, block_size, local_size, local_size, system_size, system_size, 0, d_nnz, 0, o_nnz,A); ierr =

Re: [petsc-users] Problem with MPI, MatAXPY and SAME_NONZERO_PATTERN

2015-03-27 Thread Matthew Knepley
On Fri, Mar 27, 2015 at 9:48 AM, Klaus Kaiser kai...@igpm.rwth-aachen.de wrote: Hallo, I have a strange behavior in my code concerning the function MatAXPY. I create 3 different Matrices ierr = MatCreateBAIJ(PETSC_COMM_WORLD, block_size, local_size, local_size, system_size,

Re: [petsc-users] Problem with MPI, MatAXPY and SAME_NONZERO_PATTERN

2015-03-27 Thread Barry Smith
On Mar 27, 2015, at 11:14 AM, Klaus Kaiser kai...@igpm.rwth-aachen.de wrote: Hallo Barry, I think you're right, but this is the setup I used for more than a year, so I first want to check if there is an error in my code and then switch with an correct code to a newer version when I do

Re: [petsc-users] Problem with MPI, MatAXPY and SAME_NONZERO_PATTERN

2015-03-27 Thread Barry Smith
Klaus, You would really benefit by upgrading to PETSc 3.5.3, we added much faster symbolic MatAXPY() for the DIFFERENT_NONZERO_STRUCTURE flag. Plus it is much easier for us to support the newest version. Barry On Mar 27, 2015, at 9:48 AM, Klaus Kaiser kai...@igpm.rwth-aachen.de

Re: [petsc-users] installation question

2015-03-27 Thread Satish Balay
On Fri, 27 Mar 2015, Gideon Simpson wrote: My admin did an install of petsc on ubuntu 14.0 using the package manager. Things seem to be ok, expect I noticed the following when I compile using the petsc rules: mpicxx -o ex19.o -c -Wall -Wwrite-strings -Wno-strict-aliasing

Re: [petsc-users] Problem with MPI, MatAXPY and SAME_NONZERO_PATTERN

2015-03-27 Thread Matthew Knepley
On Fri, Mar 27, 2015 at 10:16 AM, Klaus Kaiser kai...@igpm.rwth-aachen.de wrote: Hallo Matthew, thanks for your fast response. With no beside the nonzero structure I meant, that I do not create a different non-zero structure while I'm adding values to my matrices. I also tried

Re: [petsc-users] running Petsc programs without a connection.

2015-03-27 Thread Satish Balay
On Thu, 26 Mar 2015, Matthew Knepley wrote: On Thu, Mar 26, 2015 at 2:15 AM, Sanjay Kharche sanjay.khar...@manchester.ac.uk wrote: Dear All I have a fedora pc which some times cannot have an internet connection. I found that I can run non-Petsc MPI programs when it is not connected

Re: [petsc-users] Problem with MPI, MatAXPY and SAME_NONZERO_PATTERN

2015-03-27 Thread Klaus Kaiser
Hallo Matthew, thanks for your fast response. With no beside the nonzero structure I meant, that I do not create a different non-zero structure while I'm adding values to my matrices. I also tried MAT_COPY_VALUES and the result was the same. I do not get a error message but the resulting

Re: [petsc-users] Problem with MPI, MatAXPY and SAME_NONZERO_PATTERN

2015-03-27 Thread Klaus Kaiser
Hallo Matthew, here is a short example I used for 8 cores: int d_nnz[3]; int o_nnz[3]; d_nnz[0] = 3; d_nnz[1] = 3; d_nnz[2] = 3; o_nnz[0] = 6; o_nnz[1] = 6; o_nnz[2] = 6; ierr = MatCreateBAIJ(PETSC_COMM_WORLD, 3, 9, 9, 72, 72, 0, d_nnz, 0, o_nnz,A); ierr =

Re: [petsc-users] Problem with MPI, MatAXPY and SAME_NONZERO_PATTERN

2015-03-27 Thread Klaus Kaiser
Hallo Barry, I think you're right, but this is the setup I used for more than a year, so I first want to check if there is an error in my code and then switch with an correct code to a newer version when I do have the time for the installation. Is the faster symbolic MatAXPY for the