Re: [petsc-users] STFILTER in slepc

2018-10-10 Thread Matthew Knepley
On Wed, Oct 10, 2018 at 8:41 PM Moritz Cygorek wrote: > Thank you very much. Apparently, I've misunderstood what the filter > actually does. I thought about the much simpler process, where you > diagonalize > > > -(A- tau*I)^2 +offset*I > > > where tau is my target an offset is large enough so

Re: [petsc-users] STFILTER in slepc

2018-10-10 Thread Moritz Cygorek
Thank you very much. Apparently, I've misunderstood what the filter actually does. I thought about the much simpler process, where you diagonalize -(A- tau*I)^2 +offset*I where tau is my target an offset is large enough so that the global maximum is reached for eigenvalues around tau. Then

Re: [petsc-users] STFILTER in slepc

2018-10-10 Thread Jose E. Roman
This type of method requires a very high degree polynomial; suggest using degree=100 at least (this is the default value), but larger values may be necessary. Also, for this particular filter the "range" must be approximately equal to the numerical range; if you have no clue where your first

Re: [petsc-users] STFILTER in slepc

2018-10-10 Thread Moritz Cygorek
Thank you for the fast reply. I've tried running my program (using the defaul Krylov-Schur method for sparse MPI matrices) with the additional options: -st_type filter -st_filter_degree 2 -st_filter_interval 2.,2.7 -st_filter_range -2000,2000 and I get the following error message:

Re: [petsc-users] STFILTER in slepc

2018-10-10 Thread Jose E. Roman
> El 10 oct 2018, a las 19:54, Moritz Cygorek escribió: > > Hi, > > in the list of changes to SLEPc version 3.8, it is stated that there is a > preliminary implementation of polynomial filtering using STFILTER. > > Because I am struggling to obtain interior eigenvalues and harmonic >

[petsc-users] STFILTER in slepc

2018-10-10 Thread Moritz Cygorek
Hi, in the list of changes to SLEPc version 3.8, it is stated that there is a preliminary implementation of polynomial filtering using STFILTER. Because I am struggling to obtain interior eigenvalues and harmonic extraction seems not to be stable enough in my case, I wanted to give it a try,

Re: [petsc-users] error: ‘DMSetSection’ was not declared in this scope

2018-10-10 Thread Ling Zou
BTW, if I move the source file to: petsc-3.10.2/src/dm/impls/plex/examples/tutorials Using the makefile there, it compiles just fine. Ling From: petsc-users on behalf of Ling Zou Sent: Wednesday, October 10, 2018 10:35:37 AM To: PETSc users list Subject:

[petsc-users] error: ‘DMSetSection’ was not declared in this scope

2018-10-10 Thread Ling Zou
Hi All, I am starting playing with DMPlex, but get this error: error: ‘DMSetSection’ was not declared in this scope However, other DMPlex basic functions work just fine, e.g., DMPlexCreateFluentFromFile DMPlexCheckSymmetry DMPlexCheckSkeleton DMPlexCheckFaces DMPlexGetChart

Re: [petsc-users] Problems about Block Jocabi and Matrix-Free Method in SNES

2018-10-10 Thread Matthew Knepley
On Wed, Oct 10, 2018 at 8:50 AM Yingjie Wu wrote: > Thank you very much, Matt. > I followed your advice and learned some MATMFFD related functions. Since > I really wanted to use the Matrix-Free method, I learned the few examples > of this method. I would like to enhance my understanding of the

Re: [petsc-users] Problems about Block Jocabi and Matrix-Free Method in SNES

2018-10-10 Thread Yingjie Wu
Thank you very much, Matt. I followed your advice and learned some MATMFFD related functions. Since I really wanted to use the Matrix-Free method, I learned the few examples of this method. I would like to enhance my understanding of the Matrix-Free method in PETSc according to my questions about

Re: [petsc-users] PETSc/SLEPc: Memory consumption, particularly during solver initialization/solve

2018-10-10 Thread Matthew Knepley
On Wed, Oct 10, 2018 at 3:38 AM Ale Foggia wrote: > Jed, Jose and Matthew, > I've finally managed to make massif (it gives pretty detailed information, > I like it) work in the correct way in the cluster and I'm able to track > down the memory consumption, and what's more important (for me), I

Re: [petsc-users] PETSc/SLEPc: Memory consumption, particularly during solver initialization/solve

2018-10-10 Thread Ale Foggia
Jed, Jose and Matthew, I've finally managed to make massif (it gives pretty detailed information, I like it) work in the correct way in the cluster and I'm able to track down the memory consumption, and what's more important (for me), I think now I'm able to make a more accurate prediction of the