Re: [petsc-users] multi-threaded assembly

2014-06-23 Thread Jed Brown
Wen Jiang writes: > Dear all, > > I am trying to change my MPI finite element code to OPENMP one. I am not > familiar with the usage of OPENMP in PETSc and could anyone give me some > suggestions? > > To assemble the matrix in parallel using OpenMP pragmas, can I directly > call MATSETVALUES(ADD_

[petsc-users] multi-threaded assembly

2014-06-23 Thread Wen Jiang
Dear all, I am trying to change my MPI finite element code to OPENMP one. I am not familiar with the usage of OPENMP in PETSc and could anyone give me some suggestions? To assemble the matrix in parallel using OpenMP pragmas, can I directly call MATSETVALUES(ADD_VALUES) or do I need to add some l

Re: [petsc-users] pardiso interface (Barry Smith)

2014-06-23 Thread Jed Brown
Barry Smith writes: >BTW: when someone puts up a publicly available branch like this, >ANYONE can git it and try it out: it is just the regular PETSc plus >their additions. You don’t need to wait until it is merged into >PETSc (though it may be a bit more stable once it has been me

Re: [petsc-users] pardiso interface (Barry Smith)

2014-06-23 Thread Barry Smith
; > > https://bitbucket.org/petsc/petsc/pull-request/105/added-support-for-mkl-pardiso-solver/commits > > Since this uses the Pardiso in MKL, we should be able to test it without > the licensing problems of upstream Pardiso. > -- next part -- > A non-

Re: [petsc-users] pardiso interface (Barry Smith)

2014-06-23 Thread Wen Jiang
;Hmm, I don?t remember seeing it. Could you please let us know where > we can download it or send it again. > > The PR is on bitbucket: > > > https://bitbucket.org/petsc/petsc/pull-request/105/added-support-for-mkl-pardiso-solver/commits > > Since this uses the Pardiso

Re: [petsc-users] PETSC ERROR: PetscInitialize()

2014-06-23 Thread Barry Smith
On Jun 23, 2014, at 1:32 PM, Oo wrote: > > Hi Barry, > There is no other error in the following. > I get a result. Did an xterm open? Did the program just run to completion? Is this your machine? org.macosforge.xquartzAre you running on some remote machine? Barry > >

Re: [petsc-users] PETSC ERROR: PetscInitialize()

2014-06-23 Thread Oo
Hi Barry, There is no other error in the following. I get a result. Do you mean that this error must don't impact on the result? Thanks for your reply. :) Meng -- Original -- From: "Barry Smith";; Date: Jun 24, 2014 To: "Oo "; Cc: "petsc-users";

Re: [petsc-users] pardiso interface

2014-06-23 Thread Jed Brown
Barry Smith writes: >Hmm, I don’t remember seeing it. Could you please let us know where we can > download it or send it again. The PR is on bitbucket: https://bitbucket.org/petsc/petsc/pull-request/105/added-support-for-mkl-pardiso-solver/commits Since this uses the Pardiso in MKL, we

Re: [petsc-users] pardiso interface

2014-06-23 Thread Barry Smith
On Jun 23, 2014, at 11:48 AM, Michael Povolotskyi wrote: > On 6/23/2014 11:33 AM, Wen Jiang wrote: >> Hi, >> >> Could anyone tell me whether the current version of petsc supports pardiso? >> If so, are there any tricks to configure and use it? Thanks. >> >> Regards, >> Wen > We have built the

Re: [petsc-users] pardiso interface

2014-06-23 Thread Michael Povolotskyi
On 6/23/2014 11:33 AM, Wen Jiang wrote: Hi, Could anyone tell me whether the current version of petsc supports pardiso? If so, are there any tricks to configure and use it? Thanks. Regards, Wen We have built the interface and sent it to PETSc team. It works nicely for us. Dear Petsc develope

Re: [petsc-users] PETSC ERROR: PetscInitialize()

2014-06-23 Thread Satish Balay
On Mon, 23 Jun 2014, Barry Smith wrote: >If you are running with only a single process you can do > -start_in_debugger noxterm and it will open the debugger in the current > window instead of starting a new xterm. > or just: gdb executable Satish

Re: [petsc-users] PETSC ERROR: PetscInitialize()

2014-06-23 Thread Barry Smith
What happens at this point? Did the xterm open with the debugger or not? What is printed below seems like only a warning? If you are running with only a single process you can do -start_in_debugger noxterm and it will open the debugger in the current window instead of starting a new x

[petsc-users] PETSC ERROR: PetscInitialize()

2014-06-23 Thread Oo
Dear all, Recently, I try to write a simple PETSC code for solving linear system problem. At the beginning of main function, I use PetscInitialize(). However, when I run this program with "-start_in_debugger", there is a PETSC ERROR caused by this line. The PETSC ERROR is: PETSC: Attaching

[petsc-users] pardiso interface

2014-06-23 Thread Wen Jiang
Hi, Could anyone tell me whether the current version of petsc supports pardiso? If so, are there any tricks to configure and use it? Thanks. Regards, Wen

Re: [petsc-users] (no subject)

2014-06-23 Thread Dave May
And use the same sequential communicator when you create the KSP on each rank On Monday, 23 June 2014, Dave May wrote: > Yes, just assemble the same sequential matrices on each rank. To do this, > create the matrix using the communicator PETSC_COMM_SELF > > Cheers > Dave > > On Monday, 23 June

Re: [petsc-users] (no subject)

2014-06-23 Thread Dave May
Yes, just assemble the same sequential matrices on each rank. To do this, create the matrix using the communicator PETSC_COMM_SELF Cheers Dave On Monday, 23 June 2014, Bogdan Dita wrote: > > > Hello, > > I wanted to see how well Umfpack performs in PETSc but i encountered a > problem regardin

[petsc-users] (no subject)

2014-06-23 Thread Bogdan Dita
Hello, I wanted to see how well Umfpack performs in PETSc but i encountered a problem regarding the matrix type used and the way is distributed. I'm trying to solve for multiple frequencies in parallel, where system matrix A = omega*E-C, and omega is a frequency dependent scalar. Can i send matr