[petsc-users] Generalized eigenvalue problem using quad precision

2018-03-04 Thread Santiago Andres Triana
Dear all, A rather general question, is there any possibility of solving a complex-valued generalized eigenvalue problem using quad (or extended) precision when the 'B' matrix is singular? So far I have been using MUMPS with double precision with good results but I require eventually extended prec

Re: [petsc-users] how to check if cell is local owned in DMPlex

2018-03-04 Thread Matthew Knepley
On Fri, Mar 2, 2018 at 3:22 PM, Danyang Su wrote: > Hi Matt, > I use the latest Fortran style in PETSc 3.8. Enclosed are the PETSc > configuration, code compiling log and the function that causes compiling > error. The compiling error happens after I include petscsf.h in the > following section.

Re: [petsc-users] how to check if cell is local owned in DMPlex

2018-03-04 Thread Danyang Su
On 18-03-04 08:08 AM, Matthew Knepley wrote: On Fri, Mar 2, 2018 at 3:22 PM, Danyang Su > wrote: Hi Matt, I use the latest Fortran style in PETSc 3.8. Enclosed are the PETSc configuration, code compiling log and the function that causes compiling e

Re: [petsc-users] [petsc-maint] how to check if cell is local owned in DMPlex

2018-03-04 Thread Smith, Barry F.
See src/vec/is/sf/examples/tutorials/ex1f.F90 in the master branch of the PETSc git repository BTW: git grep -i petscsfgetgraph will show every use of the function in the source code. Very useful tool Barry > On Mar 4, 2018, at 1:05 PM, Danyang Su wrote: > > > > On 18-03-04

Re: [petsc-users] Generalized eigenvalue problem using quad precision

2018-03-04 Thread Jose E. Roman
Why do you want to move to quad precision? Double precision is usually enough. The fact that B is singular should not be a problem, provided that you do shift-and-invert with a nonzero target value. Can you send the output of -eps_view so that I can get a better idea what you are doing? Jose >

Re: [petsc-users] Scaling problem when cores > 600

2018-03-04 Thread TAY wee-beng
On 1/3/2018 12:14 PM, Smith, Barry F. wrote: On Feb 28, 2018, at 8:01 PM, TAY wee-beng wrote: On 1/3/2018 12:10 AM, Matthew Knepley wrote: On Wed, Feb 28, 2018 at 10:45 AM, TAY wee-beng wrote: Hi, I have a CFD code which uses PETSc and HYPRE. I found that for a certain case with grid si

Re: [petsc-users] Scaling problem when cores > 600

2018-03-04 Thread Smith, Barry F.
360 process KSPSolve 99 1.0 2.6403e+02 1.0 6.67e+10 1.1 2.7e+05 9.9e+05 5.1e+02 15100 17 42 19 15100 17 42 19 87401 1920 processes KSPSolve 99 1.0 2.3184e+01 1.0 1.32e+10 1.2 1.5e+06 4.3e+05 5.1e+02 4100 17 42 19 4100 17 42 19 967717 Ratio of number of processes