R/isAnc_fast_LE.R

In ericstrobl/CCI: Cyclic Causal Inference

isAnc_fast_LE <- function(graph,a,b,visited=rep(FALSE,nrow(graph)))
{

  if (a %in% b){
    return(TRUE)
  }

  visited[a] = TRUE;
    
  adj = which(graph[a,] & !visited);
  
  out=FALSE;
  for (j in adj){
    out=isAnc_fast_LE(graph, j, b, visited);
    if(out==TRUE){
      break;
    }
  }
  
  return(out)

}