From approximating subdivision schemes for exponential splines to high-performance interpolating algorithms