Re: [petsc-users] How to parallelize a linear system solver A*x=q in a scalar code

2014-10-13 Thread Jed Brown
Valerio Grazioso writes: > Dear All, > I have to parallelize an old fortran77-fortran90 cfd scalar By "scalar", do you mean "serial"? > code. Before starting from scratch, I decided to parallelize “just” > the pressure solver (for instance a linear system solver A*x=q). > Because of the subrout

Re: [petsc-users] -log_summary output in PETSc 3.5.2

2014-10-13 Thread Constantine Khroulev
Barry, Thanks! Would you mind adding this to the manual? A minor comment: even though -log_view :filename:format *is* mentioned, supported formats are not documented in the PetscLogView manual page. -- Constantine On Mon, Oct 13, 2014 at 3:26 PM, Barry Smith wrote: > > -log_view :filename:A

Re: [petsc-users] -log_summary output in PETSc 3.5.2

2014-10-13 Thread Barry Smith
-log_view :filename:ASCII_INFO_DETAIL It saves information from each process in one file in Python so it is easy to process the summaries you want. Barry On Oct 13, 2014, at 6:04 PM, Constantine Khroulev wrote: > Hello PETSc developers, > > Is there a way to get profiling stats simil

[petsc-users] How to parallelize a linear system solver A*x=q in a scalar code

2014-10-13 Thread Valerio Grazioso
Dear All, I have to parallelize an old fortran77-fortran90 cfd scalar code. Before starting from scratch, I decided to parallelize “just” the pressure solver (for instance a linear system solver A*x=q). Because of the subroutines call tree, and in particular for the fact a lot of subroutines are

[petsc-users] -log_summary output in PETSc 3.5.2

2014-10-13 Thread Constantine Khroulev
Hello PETSc developers, Is there a way to get profiling stats similar to -log_summary, but in a format that can be processed programmatically? As far as I can tell, current -log_summary output is tricky to parse: whether output is space-separated depends on values that are printed. In addition to

Re: [petsc-users] PetscFEIntegrateResidual_OpenCL w/ viennacl

2014-10-13 Thread Matthew Knepley
On Mon, Oct 13, 2014 at 2:16 PM, Karl Rupp wrote: > Hi David, > > > when using the snes/ex12 example > >> https://bitbucket.org/petsc/petsc/src/f5476a2b443e0eedb22d31cb1b2d85 >> b9b918f73f/src/snes/examples/tutorials/ex12.c?at=master >> with >> >> -variable_coefficient field -petscfe_type opencl

Re: [petsc-users] PetscFEIntegrateResidual_OpenCL w/ viennacl

2014-10-13 Thread Karl Rupp
Hi David, > when using the snes/ex12 example https://bitbucket.org/petsc/petsc/src/f5476a2b443e0eedb22d31cb1b2d85b9b918f73f/src/snes/examples/tutorials/ex12.c?at=master with -variable_coefficient field -petscfe_type opencl -mat_petscfe_type opencl and -dm_vec_type viennacl -dm_mat_type seqa

Re: [petsc-users] PetscFEIntegrateResidual_OpenCL w/ viennacl

2014-10-13 Thread Matthew Knepley
On Mon, Oct 13, 2014 at 12:59 PM, David Fuentes wrote: > Hi, > > when using the snes/ex12 example > https://bitbucket.org/petsc/petsc/src/f5476a2b443e0eedb22d31cb1b2d85b9b918f73f/src/snes/examples/tutorials/ex12.c?at=master > with > > -variable_coefficient field -petscfe_type opencl -mat_petscfe

[petsc-users] PetscFEIntegrateResidual_OpenCL w/ viennacl

2014-10-13 Thread David Fuentes
Hi, when using the snes/ex12 example https://bitbucket.org/petsc/petsc/src/f5476a2b443e0eedb22d31cb1b2d85b9b918f73f/src/snes/examples/tutorials/ex12.c?at=master with -variable_coefficient field -petscfe_type opencl -mat_petscfe_type opencl and -dm_vec_type viennacl -dm_mat_type seqaijviennacl

Re: [petsc-users] Crow Instability

2014-10-13 Thread Jed Brown
"Krishnappa, Sanjay" writes: >I am looking for simulating crow instability using PETSc. I am >very new to PETSc and I would like to know if there is any >work/ Code that already exist relating to crow >instability. Thanks in advance. PETSc is an algebraic solvers p

Re: [petsc-users] Crow Instability

2014-10-13 Thread Krishnappa, Sanjay
Hi, gentle reminder. regards Sanjay Krishnappa Graduate Research & Teaching Assistant Aerospace Department Embry-Riddle Aeronautical University Daytona Beach, FL? From: Krishnappa, Sanjay Sent: Saturday, October 11, 2014 9:21 PM To: petsc-users@mcs.a

Re: [petsc-users] wrap MATMPIAIJ in MatShellSetOperation and use in SLEPc

2014-10-13 Thread Jose E. Roman
El 13/10/2014, a las 19:16, Denis Davydov escribió: > Dear all, > > I am solving a generalised eigenvalue problem > >S \cdot v = a M \cdot v (1) > > using Krylov-Schur eigensolver with the shift-and-invert transformation. > Now I would like to expand the problem to the case > >

[petsc-users] wrap MATMPIAIJ in MatShellSetOperation and use in SLEPc

2014-10-13 Thread Denis Davydov
Dear all, I am solving a generalised eigenvalue problem S \cdot v = a M \cdot v (1) using Krylov-Schur eigensolver with the shift-and-invert transformation. Now I would like to expand the problem to the case [S + c^\ast \otimes c ] \cdot v = a M \cdot v (2), where “c” is s

Re: [petsc-users] Convergence problem in nonlinear problem, with PCBDDC.

2014-10-13 Thread Barry Smith
On Oct 13, 2014, at 11:24 AM, De Groof, Vincent Frans Maria wrote: > > But if it is indefinite, shouldn't the factorization algorithm complain about > a zero pivot? It may not. Cholesky generally just chugs along “pushing” the zero pivot to the final entry at the lower right corner of th

Re: [petsc-users] Convergence problem in nonlinear problem, with PCBDDC.

2014-10-13 Thread Matthew Knepley
On Mon, Oct 13, 2014 at 11:24 AM, De Groof, Vincent Frans Maria < vincent.de-gr...@uibk.ac.at> wrote: > > But if it is indefinite, shouldn't the factorization algorithm complain > about a zero pivot? > Not necessarily. There are plenty of factorizable indefinite systems. Here is one: an upper tr

Re: [petsc-users] Convergence problem in nonlinear problem, with PCBDDC.

2014-10-13 Thread De Groof, Vincent Frans Maria
But if it is indefinite, shouldn't the factorization algorithm complain about a zero pivot? Vincent Von: Matthew Knepley [knep...@gmail.com] Gesendet: Montag, 13. Oktober 2014 17:23 An: De Groof, Vincent Frans Maria Cc: petsc-users@mcs.anl.gov Betreff: Re: [pet

Re: [petsc-users] Convergence problem in nonlinear problem, with PCBDDC.

2014-10-13 Thread Matthew Knepley
On Mon, Oct 13, 2014 at 10:18 AM, De Groof, Vincent Frans Maria < vincent.de-gr...@uibk.ac.at> wrote: > Hi, > > > I am having a convergence problem in a nonlinear FE problem (we only use > ksp from petsc, not snes). We don't have any problems using a direct solver > with a cholesky factorization

[petsc-users] Convergence problem in nonlinear problem, with PCBDDC.

2014-10-13 Thread De Groof, Vincent Frans Maria
Hi, I am having a convergence problem in a nonlinear FE problem (we only use ksp from petsc, not snes). We don't have any problems using a direct solver with a cholesky factorization (Petsc factor, cholmod, mumps). Now we want to solve it using the bddc preconditioner. No problems on some prob