Added
**Dynamic time warping (`curvesimilarities.dtw`)**
- Dynamic time warping distance `dtw()`, its accumulated cost matrix `dtw_acm()` and its optimal warping path `dtw_owp()`.
Changed
**Fréchet distance (`curvesimilarities.frechet`)**
- Parametric search tolerance of `fd()` is now the machine epsilon of float64.