INNER CODE UNIT · C++
diff_mat
GreenleafLab/chromVAR · src/pwm_similarity.cpp:7
arma::mat diff_mat = arma::square(mat1 - mat2);
int width = mat1.n_cols;
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;