1
#include <Rcpp.h>
2

3
using namespace Rcpp;
4

5
// [[Rcpp::export]]
6 2
NumericMatrix dependCurFlow(NumericMatrix Tmat,IntegerMatrix el, int m, int n) {
7 2
  NumericMatrix betmat(n,n);
8 2
  for(int e=0; e<m; ++e){
9 2
    int i = el(e,0);
10 2
    int j = el(e,1);
11 2
    for(int s=0;s<n;++s){
12 2
      for(int t=0;t<n;++t){
13 2
        if((i!=s) & (i!=t) & (s!=t)){
14 2
          betmat(i,s) += 0.5*std::abs(Tmat(i,s)-Tmat(i,t)-Tmat(j,s)+Tmat(j,t));
15
        }
16 2
        if((j!=s) & (j!=t) & (s!=t)){
17 2
          betmat(j,s) += 0.5*std::abs(Tmat(j,s)-Tmat(j,t)-Tmat(i,s)+Tmat(i,t));
18
        }
19
      }
20
    }
21
  }
22 2
  return betmat;
23
}
24
  

Read our documentation on viewing source code .

Loading