The computed residuals are not good. Something went wrong during the
eigensolve, probably within your shell matrix multiply. You should check all
operations, specially within your shell matrix ops, wrapping them with
PetscCall(). Otherwise, failures go unnoticed.
Regarding shift-and-invert, you
got it, thanks Pierre & Jose.
On Mon, Aug 14, 2023 at 12:50 PM Jose E. Roman wrote:
> See for instance ex3.c and ex9.c
> https://slepc.upv.es/documentation/current/src/eps/tutorials/index.html
>
> Jose
>
>
> > El 14 ago 2023, a las 10:45, Pierre Jolivet
> escribió:
> >
> >
> >
> >> On 14 Aug 20
See for instance ex3.c and ex9.c
https://slepc.upv.es/documentation/current/src/eps/tutorials/index.html
Jose
> El 14 ago 2023, a las 10:45, Pierre Jolivet escribió:
>
>
>
>> On 14 Aug 2023, at 10:39 AM, maitri ksh wrote:
>>
>>
>> Hi,
>> I need to solve an eigenvalue problem Ax=lmbda*x
> On 14 Aug 2023, at 10:39 AM, maitri ksh wrote:
>
>
> Hi,
> I need to solve an eigenvalue problem Ax=lmbda*x, where A=(B^-H)*Q*B^-1 is a
> hermitian matrix, 'B^-H' refers to the hermitian of the inverse of the matrix
> B. Theoretically it would take around 1.8TB to explicitly compute the
Hi,
I need to solve an eigenvalue problem *Ax=lmbda*x*, where A=(B^-H)*Q*B^-1
is a hermitian matrix, 'B^-H' refers to the hermitian of the inverse of the
matrix B. Theoretically it would take around 1.8TB to explicitly compute
the matrix B^-1 . A feasible way to solve this eigenvalue problem would
