The following numeric routines should be used in preference to others that do
the same function, as these are reasonably optimized. Generally these are
tailored to the target machine (e.g. on the Cray these are
based on Cray SCILIB routines, on the Convex, they are partially based
on VECLIB routines).