Warhammer 40k, R - rextester

Language: Layout:

# Set number of dice and hit and wound roll needed
hitdice     <- 100  ### how many dice do you roll to hit
hit         <- 3   ### what do you hit on
wound       <- 3   ### what do you wound on  / set to 0 to get back number of hits (and set save to 7)
save        <- 3   ### what is their armour save / set to 7 to get back number of wounds 
ap          <- 0  ### what is the ap
damage      <- 1  ### how many damage per wound / if random use expected damage (e.g., d6 damage is 3.5 d3 damage is 2)
rerollhit1   <- T   ### do you re roll hits of 1 F for no T for yes
rerollwound1 <- F   ### do you re roll wounds of 1 F for no T for yes
rerollhit   <- F   ### do you re roll hits F for no T for yes
rerollwound <- F   ### do you re roll wounds F for no T for yes

rolldice   <- function(hitdice,hit,wound,save,ap,rerollhit,rerollwound){
  savemod    <- save-ap
  
  # Roll for hits
  hitresult  <- sample(c(1:6),hitdice,replace=T)
  hitresult2 <- 0
  if(rerollhit1 == T) {
    rerolls1    <- length(which(hitresult<hit))
    if(rerolls1 > 0) {
      hitresult2 <- sample(c(1:6),rerolls1,replace=T)
    }
  }
  if(rerollhit == T) {
    rerolls    <- length(which(hitresult<hit))
    if(rerolls > 0) {
      hitresult2 <- sample(c(1:6),rerolls,replace=T)
    }
  }
    
  hits       <- sum(c(hitresult >= hit, hitresult2 >= hit))
  
  # Roll for wounds
  woundresult  <- sample(c(1:6),hits,replace=T)
  woundresult2 <- 0
  if(rerollwound1 == T) {
    rerolls1    <- length(which(woundresult<wound))
    if(rerolls1 > 0) {
      woundresult2 <- sample(c(1:6),rerolls1,replace=T)
    }
  }
  if(rerollwound == T) {
    rerolls    <- length(which(woundresult<wound))
    if(rerolls > 0) {
      woundresult2 <- sample(c(1:6),rerolls,replace=T)
    }
  }
    
    wounds     <- sum(c(woundresult >= wound, woundresult2 >= wound))
    saves      <- sum(sample(c(1:6),wounds,replace=T) >= savemod)
    dead       <- wounds - saves
    dead
  }
  
  tmp <- NULL
  for(i in 1:20000) {
    tmp[i] <- rolldice(hitdice,hit,wound,save,ap,rerollhit,rerollwound)
  }
  paste("Number of Wounds Lost:",mean(tmp)*damage, sep = " ")

[ + ] Show input
Absolute running time: 3.04 sec, cpu time: 4.99 sec, memory peak: 37 Mb, absolute service time: 3,05 sec