+
+ i = pair[1] ; j = pair[2]
+ if (verbose && j==i+1 && !parll)
+ cat(paste(" Distances (",i,",",j,"), (",i,",",j+1,") ...\n", sep=""))
+
+ # Compute CWT of columns i and j in synchrones
+ L = nrow(synchrones)
+ cwt_i <- getSynchroneCWT(i, L)
+ cwt_j <- getSynchroneCWT(j, L)
+
+ # Compute the ratio of integrals formula 5.6 for WER^2
+ # in https://arxiv.org/abs/1101.4744v2 §5.3
+ num <- filterMA(Mod(cwt_i * Conj(cwt_j)))
+ WX <- filterMA(Mod(cwt_i * Conj(cwt_i)))
+ WY <- filterMA(Mod(cwt_j * Conj(cwt_j)))
+ wer2 <- sum(colSums(num)^2) / sum(colSums(WX) * colSums(WY))
+
+ Xwer_dist[i,j] <- sqrt(L * ncol(cwt_i) * (1 - wer2))
+ Xwer_dist[j,i] <- Xwer_dist[i,j]
+ Xwer_dist[i,i] <- 0.