Characteristic-based numerical algorithms for stiff hyperbolic relaxation systems