code: def vecdiff(htilde, hinterp, sample_points, psd=None):
      
      vecdiffs = numpy.zeros(sample_points.size-1, dtype=float)
      for kk,thisf in enumerate(sample_points[:-1]):
          nextf = sample_points[kk+1]
          vecdiffs[kk] = abs(_vecdiff(htilde, hinterp, thisf, nextf, psd=psd))
      return vecdiffs
docs: Computes a statistic indicating between which sample points a waveform
      and the interpolated waveform differ the most.