INNER CODE UNIT · C++
raw_dist
GreenleafLab/chromVAR · src/pwm_similarity.cpp:9
double raw_dist = arma::sum(arma::sqrt((arma::rowvec)arma::sum(diff_mat,0)));
double norm_dist = raw_dist / sqrt(2) / width;
return norm_dist;
}
// [[Rcpp::export]]
arma::vec pwm_dist_single(arma::mat mat1, arma::mat mat2,
arma::uword min_overlap){
arma::uword n1 = mat1.n_cols;
arma::uword n2 = mat2.n_cols;
arma::uword m = n1 + n2 - 1;
arma::vec out(2);
double tmp;