Re: [petsc-users] Solve KSP in parallel.

Ok, last output was from simulated multicores, in an actual cluster the errors are of the kind: [valera@cinci CSRMatrix]$ petsc -n 2 ./solvelinearmgPETSc TrivSoln loaded, size:4 / 4 TrivSoln loaded, size:4 / 4 RHS loaded, size:4 /

Re: [petsc-users] Solve KSP in parallel.

Ok, i created a tiny testcase just for this, The output from n# calls are as follows: n1: Mat Object: 1 MPI processes type: mpiaij row 0: (0, 1.) (1, 2.) (2, 4.) (3, 3.) row 1: (0, 2.) (1, 1.) (2, 3.) (3, 4.) row 2: (0, 4.) (1, 3.) (2, 1.) (3, 2.) row 3: (0, 3.) (1, 4.) (2, 2.) (3,

Re: [petsc-users] Solve KSP in parallel.

> On Sep 26, 2016, at 5:07 PM, Manuel Valera wrote: > > Ok i was using a big matrix before, from a smaller testcase i got the output > and effectively, it looks like is not well read at all, results are attached > for DRAW viewer, output is too big to use STDOUT even in

Re: [petsc-users] Solve KSP in parallel.

Ok i was using a big matrix before, from a smaller testcase i got the output and effectively, it looks like is not well read at all, results are attached for DRAW viewer, output is too big to use STDOUT even in the small testcase. n# is the number of processors requested. is there a way to create

Re: [petsc-users] Solve KSP in parallel.

How large is the matrix? It will take a very long time if the matrix is large. Debug with a very small matrix. Barry > On Sep 26, 2016, at 4:34 PM, Manuel Valera wrote: > > Indeed there is something wrong with that call, it hangs out indefinitely > showing only:

Re: [petsc-users] Solve KSP in parallel.

Indeed there is something wrong with that call, it hangs out indefinitely showing only: Mat Object: 1 MPI processes type: mpiaij It draws my attention that this program works for 1 processor but not more, but it doesnt show anything for that viewer in either case. Thanks for the insight on

Re: [petsc-users] Solve KSP in parallel.

The call to MatCreateMPIAIJWithArrays() is likely interpreting the values you pass in different than you expect. Put a call to MatView(Ap,PETSC_VIEWER_STDOUT_WORLD,ierr) after the MatCreateMPIAIJWithArray() to see what PETSc thinks the matrix is. > On Sep 26, 2016, at 3:42 PM, Manuel

[petsc-users] Solve KSP in parallel.

Hello, I'm working on solve a linear system in parallel, following ex12 of the ksp tutorial i don't see major complication on doing so, so for a working linear system solver with PCJACOBI and KSPGCR i did only the following changes: call MatCreate(PETSC_COMM_WORLD,Ap,ierr) ! call

Re: [petsc-users] petsc4py with complex numbers?

The output is: I installed petsc via macport and wrongly assumed complex number support would have been default. I've just reinstalled petsc with it and it seems to be working. thanks Lisandro ! aurelien Le 25/09/16 à 23:00, Lisandro Dalcin a écrit : Are you sure you built petsc4py with

Re: [petsc-users] DMPlex problem

Morten, have a look at the local-to-global mapping that is created for the MATIS case in DMCreateMatrix_Plex in src/dm/impls/plex/plex.c. and how DMPlexMatSetClosure is used in src/snes/utils/dmplexsnes.c in the MATIS case. MATIS is a format for non-overlapping domain decomposition when you

Re: [petsc-users] DMPlex problem

On Mon, Sep 26, 2016 at 8:41 AM, Morten Nobel-Jørgensen wrote: > Hi Matt > > We are trying to do a simple FE using DMPlex, but when assemble the global > stiffness matrix we get in problems when running NP>1 - that is the global > matrix differs when we move to a distributed system,

Re: [petsc-users] DMPlex problem

Hi Matt We are trying to do a simple FE using DMPlex, but when assemble the global stiffness matrix we get in problems when running NP>1 - that is the global matrix differs when we move to a distributed system, where it should not. In pseudo-code our CreateGlobalStiffnessMatrix does the

Re: [petsc-users] DMPlex problem

On Mon, Sep 26, 2016 at 7:00 AM, Morten Nobel-Jørgensen wrote: > Hi Matthew > > It seems like the problem is not fully fixed. I have changed the code to > now run on with both 2,3 and 4 cells. When I run the code using NP = 1..3 I > get different result both for NP=1 to NP=2/3 when

Re: [petsc-users] DMPlex problem

Hi Matthew It seems like the problem is not fully fixed. I have changed the code to now run on with both 2,3 and 4 cells. When I run the code using NP = 1..3 I get different result both for NP=1 to NP=2/3 when cell count is larger than 2. Kind regards, Morten mpiexec -np 1 ./ex18k cells 2

[petsc-users] Position for Junior Professorship in High Performance Computing, Ruhr University Bochum, Germany

Dear colleagues There is the position opening for Junior Professorship in HPC at my department. Please take a look if you are interested in such the position or forward this information to your colleagues. Thank you very much. --- MSc. Hoang-Giang Bui Ruhr University Bochum Institute for

Re: [petsc-users] petsc4py with complex numbers?

Are you sure you built petsc4py with a PETSc build with complex scalars? Whats the output of "print(PETSc.ScalarType)" ? On 26 September 2016 at 07:47, Aurelien Ponte wrote: > Hi, > > I am trying to solve a linear problem whose operator has complex > coefficients