"calc.graph.pars.fast" <-
function (sp.mat, rank, sub, maxWgt = 1, verbose = FALSE) 
{
  N <- length(rank)
  DE <- rank <= sub
  sp.mat.A <- sp.mat[DE, DE]
  sp.mat.B <- sp.mat[!DE, !DE]
  deg <- rowSumsSp(sp.mat)
  Eall <- mean(deg)
  EA <- mean(deg[DE])
  EB <- mean(deg[!DE])
  beta <- EA/EB
  pA <- sum(sp.mat.A@ra)/2/(choose(sub,2) * maxWgt)
  pB <- sum(sp.mat.B@ra)/2/(choose(N-sub,2) * maxWgt)
  pAB <- (EA/maxWgt - (sub - 1) * pA)/(N - sub)
  relcon <- pA/pAB
  if (verbose) {
    cat(paste("E_deg(all) = ", Eall, "\n", sep = ""))
    cat(paste("E_deg(A) = ", EA, "\n", sep = ""))
    cat(paste("E_deg(B) = ", EB, "\n", sep = ""))
    cat(paste("beta = ", beta, "\n", sep = ""))
    cat(paste("p_A = ", pA, "\n", sep = ""))
    cat(paste("p_B = ", pB, "\n", sep = ""))
    cat(paste("p_AB = ", pAB, "\n", sep = ""))
    cat(paste("Rel. conn. = ", relcon, "\n", sep = ""))
  }
  result <- list(Eall = Eall, EA = EA, EB = EB, beta = beta,
      pA = pA, pB = pB, pAB = pAB, rconn = relcon)
  return(result)
}