Identifying the convective/absolute instability nature of a local base flow requires an analysis of its linear impulse response. One must find the appropriate singularity in the eigenvalue problem with complex frequencies and wavenumbers and prove causality. One way to do so is to show that the appropriate integration contour of this response, a steepest decent path through the relevant singularity, exists. Due to the inherent difficulties of such a proof, one often verifies instead whether this singularity satisfies the collision criterion. In other words, one must show that the branches involved in the formation of this singularity come from distinct halves of the complex wavenumber plane. However, this graphical search is computationally intensive in a single plane and essentially prohibitive in two planes. A significant computational cost reduction can be achieved when root finding procedures are applied instead of graphical ones to search for singularities. They focus on locating these points, with causality being verified graphically a posteriori for a small parametric sample size. The use of root-finding procedures require auxiliary equations, often derived by applying the zero group velocity conditions to the dispersion relation. This relation, in turn, is derived by applying matrix forming to the differential eigenvalue problem and taking the determinant of the resulting system of algebraic equations. Taking the derivative of the dispersion relation with respect to the wavenumbers generates the auxiliary equations. If the algebraic system is decoupled, this derivation is straightforward. However, its computational cost is often prohibitive when the algebraic system is coupled. Other methods exist, but often they can also be too costly and/or not reliable for two wavenumber plane searches. This paper describes an alternative methodology based on sensitivity analysis and adjoints that allow the zero group velocity conditions to be applied directly to the differential eigenvalue problem. In doing so, the direct and auxiliary differential eigenvalue problems can be solved simultaneously using standard shooting methods to directly locate singularities. Auxiliary dispersion relations no longer have to be derived, although it is shown that they are the algebraic equivalent of the auxiliary differential eigenvalue problems obtained by this alternative methodology. Using the latter dramatically reduces computational costs. The search for arbitrary singularities is then not only accelerated in single wavenumber planes but it also becomes viable in two wavenumber planes. Finally, the new method also allows group velocity calculations, greatly facilitating the verification of causality. Several test cases are presented to illustrate the capabilities of this new method.