INNER CODE UNIT · C++

m

GreenleafLab/chromVAR · src/pwm_similarity.cpp:23

  arma::uword m = n1 + n2 - 1;

  arma::vec out(2);
  double tmp;
  //arma::uword min_overlap = 4;
  for (arma::uword i = min_overlap - 1; i < m-min_overlap + 1; i++){
    arma::mat mat1s,mat2s;
    if (i < n2){
      if (i < n1){
        mat1s = mat1.cols(0, i);
        mat2s = mat2.cols(n2 - 1 - i, n2 - 1);
      } else{
        mat1s = mat1.cols(0,  n1 - 1);
        mat2s = mat2.cols(n2 - 1 - i, m - i - 1);
      }

    } else{
      if (i < n1){

View source record →