  In this paper we consider the problem of computing generalized eigenvectors
of a matrix pencil in real Schur form. In exact arithmetic, this problem can be
solved using substitution. In practice, substitution is vulnerable to
floating-point overflow. The robust solvers xTGEVC in LAPACK prevent overflow
by dynamically scaling the eigenvectors. These subroutines are sequential
scalar codes which compute the eigenvectors one by one. In this paper we
discuss how to derive robust blocked algorithms. The new StarNEig library
contains a robust task-parallel solver Zazamoukh which runs on top of StarPU.
Our numerical experiments show that Zazamoukh achieves a super-linear speedup
compared with DTGEVC for sufficiently large matrices.
