Given a collection of MTS, the function returns the pairwise distance matrix,
where the distance between two MTS \(\boldsymbol X_T\) and \(\boldsymbol Y_T\) is defined
as \(d_{Eros}(\boldsymbol X_T, \boldsymbol Y_T)=\sqrt{2-2Eros(\boldsymbol X_T, \boldsymbol Y_T)}\),
where $$Eros(\boldsymbol X_T, \boldsymbol Y_T)=\sum_{i=1}^{d}w_i|<\boldsymbol x_i,\boldsymbol y_i>|=
\sum_{i=1}^{d}w_i|\cos \theta_i|,$$
where \(\{\boldsymbol x_1, \ldots, \boldsymbol x_d\}\), \(\{\boldsymbol y_1, \ldots, \boldsymbol y_d\}\)
are sets of eigenvectors concerning the covariance or correlation matrix of series \(\boldsymbol X_T\) and
\(\boldsymbol Y_T\), respectively, \(<\boldsymbol x_i,\boldsymbol y_i>\) is the inner product of
\(\boldsymbol x_i\) and \(\boldsymbol y_i\), \(\boldsymbol w=(w_1, \ldots, w_d)\)
is a vector of weights which is based on the eigenvalues of the MTS dataset with \(\sum_{i=1}^{d}w_i=1\)
and \(\theta_i\) is the angle between \(\boldsymbol x_i\) and \(\boldsymbol y_i\).