tests/testthat/test_makeTs.R

In rte-antares-rpackage/antaresFlowbased: Run FlowBased Simulations on Antares

context("createFBTS")


test_that("make ts", {

  op5 <- antaresRead::setSimulationPath(testStudy2)
  
  matProb <- readRDS(system.file("testdata/proba.RDS", package = "antaresFlowbased"))
  
  setnames(matProb[[1]],"FR_load", "fr@load" )
  setnames(matProb[[2]],"FR_load", "fr@load" )
  
  setnames(matProb[[1]],"DE_wind", "de@wind" )
  setnames(matProb[[2]],"DE_wind", "de@wind" )
  
  setnames(matProb[[1]],"DE_solar", "be@wind" )
  setnames(matProb[[2]],"DE_solar", "be@wind" )
  
  multiplier <- data.frame(variable = c("fr@load", "de@wind", "be@wind"),
                           coef = c(1, 352250, 246403))
  firstDay <- suppressWarnings(identifyFirstDay(op5, firstArea = "fr", secondArea = NULL))
  
  
  interSeasonBegin <- as.Date(c("2017-09-03", "2018-02-02"))
  interSeasonEnd <- as.Date(c("2017-10-04", "2018-05-02"))
  
  firstF <- secondF <-  NULL
  for(k in 1:10)
  {
    ts <- suppressWarnings(createFBTS(opts = op5, probabilityMatrix = matProb, multiplier = multiplier,
                     interSeasonBegin = interSeasonBegin, interSeasonEnd = interSeasonEnd,
                     firstDay = firstDay, seed = k, silent = TRUE, outputPath =  tempdir()))
    
    
    frLoad <- suppressWarnings(readInputTS(load = "fr", timeStep = "daily", showProgress = FALSE))
    windbe <- suppressWarnings(readInputTS(wind = c("be"), timeStep = "daily", showProgress = FALSE))
    windde <- suppressWarnings(readInputTS(wind = c("de"), timeStep = "daily", showProgress = FALSE))
    allDta <- data.table(frLoad, be = windbe[["wind"]],de = windde[["wind"]])
    allDta <- allDta[tsId == 1]
    
    dates <- allDta$time
    calendar <- .getVirtualCalendar(dates, interSeasonBegin, interSeasonEnd, firstDay)
    
    data1 <- allDta[180]
   
    firstF <- c(firstF, ts[ts$Date == data1$time]$`1`)
    
    data2 <- allDta[7]
    
    secondF <- c(secondF, ts[ts$Date == data2$time]$`1`)
    
    
  }
  expect_true(2%in%firstF)
  expect_true(1%in%firstF)
  expect_true(all(firstF%in%c(1, 2)))
  expect_true(all(secondF == 6))
  
  
  
  
  
})