(defun riemann-left-fast (a b step)
  (declare (type double-float a b)
           (optimize (speed 3)
                     (safety 0)))
  (* step (loop for n double-float from a to (- b step)
        by step
        summing (sin n))))