R/zzz_reachability.R

In markovchain: Easy Handling Discrete Time Markov Chains

# Normalize the matrix orientation before calling the C++ reachability kernel.
# The kernel interprets the matrix as row-stochastic, while a markovchain
# object may represent transitions by columns when byrow = FALSE.
setMethod("is.accessible", c("markovchain", "missing", "missing"),
  function(object, from, to) {
    if (!object@byrow) {
      object@transitionMatrix <- t(object@transitionMatrix)
      object@byrow <- TRUE
    }

    .reachabilityMatrixRcpp(object)
  }
)