simulateNetworkBasic <- function(
  initDifficulty,
  initReward,
  numCompetitions,
  idealBlockTime,
  idealNetworkPower,
  priceList = NULL,
  nodePowers = NULL, 
  nodeCFs = NULL, 
  difficultyAdaptationRate = 0.1,
  rewardAdaptationRate = 0.1,
  adaptationStrategy = "adaptive"){

  numNodes = length(nodePowers)
  
  #
  # Results List
  #
  difficultyList = NULL
  rewardList = NULL
  blockTimeList = NULL
  powerList = NULL
  participationList = NULL
  
  # play is as long as the number of nodes and flags which node will play the next round
  play = numeric()
  
  reward = initReward
  difficulty = initDifficulty
  networkPower = sum(nodePowers)
  
  print(paste0("Running ", numCompetitions, " competitions for ", numNodes))
  pb <- txtProgressBar(1, numCompetitions, style=3)
  
  for (i in 1:numCompetitions){
    
    #
    # Each node decides whether to participate.
    #
    for (j in 1:numNodes){
      bkEvenR = getBreakEvenReward (d = difficulty,
                                    x = priceList[i],
                                    c = nodeCFs[j],
                                    N = (networkPower-nodePowers[j])/nodePowers[j]
      )

    if (reward > bkEvenR){
        play[j] = 1
      } else {
        play[j] = -1
      }
    }

        
    numParticipants = sum(play[play==1])
    networkPower = sum(nodePowers[play==1]) # Full observability
    
    #
    # If the network power is zero then the network has simply collapsed.
    # [we choose to end simulation for that]
    #
    if (numParticipants == 0) {
      print("Network Collapsed.")
      return(data.frame(Difficulty = difficultyList,
                        BlockTime = blockTimeList,
                        Reward = rewardList,
                        Price = priceList[1:length(difficultyList)],
                        NetPower = powerList,
                        Participation = participationList
      ))
    }
    
    
    blockTime = difficulty/networkPower
    
    # Capture Statistics
    difficultyList = c(difficultyList,difficulty)
    rewardList = c(rewardList,reward)
    powerList = c(powerList,networkPower)
    blockTimeList = c(blockTimeList,blockTime)
    participationList = c(participationList,numParticipants)
    
    
    #
    # Netowrk Adapts Reward and Difficulty.
    #
    
    # Adjust d
    if (adaptationStrategy == "fixed"){
      # Do nothing
    } else {
      newDiff = idealBlockTime * networkPower
      difficulty = difficulty + (newDiff - difficulty)*difficultyAdaptationRate
    }
    
    # Adjust r
    if (adaptationStrategy == "fixed"){
      # Do nothing
    } else {
      newR = getOptimalReward(difficulty = difficulty,
                              price = priceList[i],
                              nodePowers = nodePowers,
                              nodeCFs = nodeCFs,
                              networkPower = networkPower,
                              idealPower = idealNetworkPower)
      reward = reward + (newR - reward)*rewardAdaptationRate          
    }
    setTxtProgressBar(pb,i)
  }
  
  return(data.frame(Difficulty = difficultyList,
                    BlockTime = blockTimeList,
                    Reward = rewardList,
                    Price = priceList,
                    NetPower = powerList,
                    Participation = participationList
  ))
  
}