Balancing may reduce the norm of the matrix or the condition number of some
eigenvalues. It applies to the ST operator, (A-sigma*B)^{-1}*B in case of
shift-and-invert. The case that sigma is close to an eigenvalue is usually not
a problem, provided that you use a robust direct solver (MUMPS). In
Good afternoon Mr Roman,
Thank you very much for your detailed and quick answer. I'll make the
use of eps_view and log_view and see if I can optimize the preallocation
of my matrix.
Good to know that it is possible to play with the "mpd" parameter, I
thought it was only for when large values
> El 19 feb 2018, a las 19:15, Thibaut Appel
> escribió:
>
> Good afternoon,
>
> I am solving generalized eigenvalue problems {Ax = omegaBx} in complex
> arithmetic, where A is non-hermitian and B is singular. I think the only way
> to get round the singularity is to employ a shift-and-inve
Good afternoon,
I am solving generalized eigenvalue problems {Ax = omegaBx} in complex
arithmetic, where A is non-hermitian and B is singular. I think the only
way to get round the singularity is to employ a shift-and-invert method,
where I am using MUMPS to invert the shifted matrix.
I am u