<div dir="ltr">Hi, thanks for the quick answer.<br><br>I just did it, and it does not work. My problem is GNHEP and I use the default solver (Krylov-Schur). Moreover I run the code with the options: -st_ksp_type preonly -st_pc_type lu -st_pc_factor_mat_solver_type mumps<br><div><br></div><div>Any other suggestions?</div>Kind regards.</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El lun., 17 de feb. de 2020 a la(s) 12:39, Jeremy Theler (<a href="mailto:jeremy@seamplex.com" target="_blank">jeremy@seamplex.com</a>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">The usual trick is to set ones in one matrix and zeros in the other<br>
one.<br>
<br>
<br>
On Mon, 2020-02-17 at 12:35 -0600, Emmanuel Ayala wrote:<br>
> Hi everyone,<br>
> <br>
> I have an eigenvalue problem where I need to apply BCs to the<br>
> stiffness and mass matrix. <br>
> <br>
> Usually, for KSP solver, it is enough to set to zero the rows and<br>
> columns related to the boundary conditions. I used to apply it with<br>
> MatZeroRowsColumns, with a 1s on the diagonal. Then the solver works<br>
> well.<br>
> <br>
> There is something similar to KSP for EPS solver ? <br>
> <br>
> I already used MatZeroRowsColumns (for EPS solver), with a 1s on the<br>
> diagonal, and I got wrong result.<br>
> <br>
> Kind regards.<br>
> <br>
> <br>
> <br>
> <br>
<br>
</blockquote></div>