Re: [petsc-users] PETSC ERROR: Petsc has generated inconsistent data with PCGAMG and KSPGRMES

2019-07-11 Thread 我 via petsc-users
Thank you very much! Daye At 2019-07-10 02:36:55, "Smith, Barry F." wrote: > > Most preconditioners don't care if the matrix has non-symmetric non-zero > structure. GAMG requires it have a symmetric nonzero structure. It lists some > options you can try to make it symmetric. > > >>

Re: [petsc-users] problem with CISS in SLEPC

2019-07-11 Thread Jose E. Roman via petsc-users
Internally the solver builds a subspace of dimension L*M, which defaults to 128. See the details here http://slepc.upv.es/documentation/reports/str11.pdf I guess the method is not prepared for smallish matrices. Try reducing the values of L and M. Jose > El 11 jul 2019, a las 19:43, Matthew Kne

[petsc-users] problem with CISS in SLEPC

2019-07-11 Thread Povolotskyi, Mykhailo via petsc-users
Hello, I want to use CISS for my eigenvalue problem. To test it I tried a simple (2x2) matrix case, but the code failed with the following message: [0]PETSC ERROR: - Error Message -- [0]PETSC ERROR: Error in externa

Re: [petsc-users] What is the best way to do domain decomposition with petsc?

2019-07-11 Thread Jed Brown via petsc-users
Matthew Knepley via petsc-users writes: >> I am just wondering which way is better, or do you have any other >> suggestion? >> > If you plan on doing a lot of mesh manipulation by hand and want to control > everything, the first option might be better. > On the other hand, if you use Plex, you co

[petsc-users] What is the best way to do domain decomposition with petsc?

2019-07-11 Thread Dongyu Liu - CITG via petsc-users
Hi, We want to incorporate the PETSc into our in-house FEM package. We found two ways to do the domain decomposition. The first one is to read the mesh with partitioning, and the partitioning is done with gmsh. For this one, we need to do the index mapping (renumbering). We found that the mes