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).