Computing the numerical range of Krein space operators
Consider the Hilbert space (H,〈• , •〉) equipped with the indefinite inner product[u,v]=v*J u,u,v∈ H, where J is an indefinite self-adjoint involution acting on H. The Krein space numerical range WJ(T) of an operator T acting on H is the set of all the values attained by the quadratic form [Tu,u], with u ∈H satisfying [u,u]=± 1. We develop, implement and test an alternative algorithm to compute WJ(T) in the finite dimensional case, constructing 2 by 2 matrix compressions of T and their easily determined...