<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 style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">

Hello,<br>

<br>

I have a general non hermitian eigenvalue problem arising from the 3D helmholtz equation.<br>

The form of the helmholtz equaton is:<br>

<br>

(S - k^2M)v = lambda k^2 M v<br>

<br>

Where S is the stiffness/curl-curl matrix and M is the mass matrix associated with edge elements used to discretize the problem.<br>

The helmholtz equation creates eigenvalues of -1.0, which I believe are eigenvectors that are part of the null space of the curl-curl operator S.

<br>

<br>

For my application, I would like to compute eigenvalues > -1.0, and avoid computation of eigenvalues of -1.0.<br>

I am currently using shift invert ST with mumps LU direct solver. By increasing the shift away from lambda=-1.0. I get faster computation of eigenvectors, and the lambda=-1.0 eigenvectors appear to slow down the computation by about a factor of two.<br>

Is there a way to avoid these lambda = -1.0 eigenpairs with a GNHEP problem type?<br>

<br>

Regards,

<div>Lucas</div>

<br>

<br>

</div>

</body>

</html>