  The effective numerical method is developed performing the test of the
hyperbolicity of chaotic dynamics. The method employs ideas of algorithms for
covariant Lyapunov vectors but avoids their explicit computation. The outcome
is a distribution of a characteristic value which is bounded within the unit
interval and whose zero indicate the presence of tangency between expanding and
contracting subspaces. To perform the test one needs to solve several copies of
equations for infinitesimal perturbations whose amount is equal to the sum of
numbers of positive and zero Lyapunov exponents. Since for high-dimensional
system this amount is normally much less then the full phase space dimension,
this method provide the fast and memory saving way for numerical hyperbolicity
test of such systems.
