<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Thank you for the fast reply. <br>
</p>
<p><br>
</p>
<p>I've tried running my program (using the defaul Krylov-Schur method for sparse MPI matrices) with the additional options:</p>
<p><br>
</p>
<p><span>-st_type filter -st_filter_degree 2 -st_filter_interval 2.,2.7 -st_filter_range -2000,2000</span></p>
<p><span><br>
</span></p>
<p><span>and I get the following error message:</span></p>
<p><span><br>
</span></p>
<p><span></p>
<div>[0]PETSC ERROR: STFILTER cannot get the filter specified; please adjust your filter parameters (e.g. increasing the polynomial degree)<br>
</div>
....</span>
<p></p>
<p><span></p>
<div>[0]PETSC ERROR: #1 FILTLAN_GetIntervals() line 451 in /home/applications/sources/libraries/slepc-3.9.2/src/sys/classes/st/impls/filter/filtlan.c<br>
[0]PETSC ERROR: #2 STFilter_FILTLAN_setFilter() line 1016 in /home/applications/sources/libraries/slepc-3.9.2/src/sys/classes/st/impls/filter/filtlan.c<br>
[0]PETSC ERROR: #3 STSetUp_Filter() line 42 in /home/applications/sources/libraries/slepc-3.9.2/src/sys/classes/st/impls/filter/filter.c<br>
[0]PETSC ERROR: #4 STSetUp() line 271 in /home/applications/sources/libraries/slepc-3.9.2/src/sys/classes/st/interface/stsolve.c<br>
[0]PETSC ERROR: #5 EPSSetUp() line 263 in /home/applications/sources/libraries/slepc-3.9.2/src/eps/interface/epssetup.c<br>
[0]PETSC ERROR: #6 EPSSolve() line 135 in /home/applications/sources/libraries/slepc-3.9.2/src/eps/interface/epssolve.c<br>
<br>
</div>
<br>
</span>
<p></p>
<p><span><br>
</span></p>
<p><span>Do you have a clue what I've missed?</span></p>
<p><span><br>
</span></p>
<p><span><br>
</span></p>
<p><span>Moritz<br>
</span></p>
<br>
<br>
<div style="color:rgb(0,0,0)">
<div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> Jose E. Roman <jroman@dsic.upv.es><br>
<b>Sent:</b> Wednesday, October 10, 2018 2:30 PM<br>
<b>To:</b> Moritz Cygorek<br>
<b>Cc:</b> petsc-users@mcs.anl.gov<br>
<b>Subject:</b> Re: [petsc-users] STFILTER in slepc</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt">
<div class="PlainText"><br>
> El 10 oct 2018, a las 19:54, Moritz Cygorek <mcygorek@uottawa.ca> escribió:<br>
> <br>
> Hi,<br>
> <br>
> in the list of changes to SLEPc version 3.8, it is stated that there is a preliminary implementation of polynomial filtering using STFILTER.
<br>
> <br>
> 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, but I could not find any documentation yet.<br>
> <br>
> Does anybody have an example of how to use STFILTER or any documentation about it?<br>
> <br>
> Thanks in advance,<br>
> Moritz<br>
<br>
There are no examples. You just set the type to STFILTER and set some parameters such as the interval of interest or the polynomial degree. See functions starting with STFilter here:
<a href="http://slepc.upv.es/documentation/current/docs/manualpages/ST/index.html" id="LPlnk589364" previewremoved="true">
http://slepc.upv.es/documentation/current/docs/manualpages/ST/index.html</a><br>
<br>
In some problems it works well, but don't expect too much. It is still in our to-do list to make it more usable. It will be good to have your feedback. If you want, send results to slepc-maint, maybe we can help tuning the parameters.<br>
<br>
Jose<br>
<br>
</div>
</span></font></div>
</div>
</body>
</html>