  The spectral problem associated with the linearization about solitary waves
of spinor systems or optical coupled mode equations supporting gap solitons is
formulated in terms of the Evans function, a complex analytic function whose
zeros correspond to eigenvalues. These problems may exhibit oscillatory
instabilities where eigenvalues detach from the edges of the continuous
spectrum, so called edge bifurcations. A numerical framework, based on a fast
robust shooting algorithm using exterior algebra is described. The complete
algorithm is robust in the sense that it does not produce spurious unstable
eigenvalues. The algorithm allows to locate exactly where the unstable discrete
eigenvalues detach from the continuous spectrum. Moreover, the algorithm allows
for stable shooting along multi-dimensional stable and unstable manifolds. The
method is illustrated by computing the stability and instability of gap
solitary waves of a coupled mode model.
