Computing the numerical range of Krein space operators