test_that("Test getRmw function", {
  expect_equal(getRmw(seq(0, 80, 5), -15), c(
    60, 55, 51, 47, 44, 41, 38, 35, 32,
    30, 28, 25, 24, 22, 20, 19, 17
  ))
})





test_that("Test Willougbhy / Holland / Boose functions", {
  expect_equal(
    willoughby(seq(1, 100), 20, 50, -15),
    c(
      1.813, 3.905, 6.118, 8.414, 10.772, 13.182, 15.635, 18.127, 20.652, 23.208,
      25.791, 28.399, 31.032, 33.685, 36.360, 39.053, 41.765, 44.494, 47.239, 50.000,
      49.327, 48.676, 48.044, 47.433, 46.839, 46.264, 45.706, 45.165, 44.640, 44.131,
      43.636, 43.155, 42.688, 42.235, 41.794, 41.366, 40.949, 40.543, 40.149, 39.765,
      39.391, 39.027, 38.672, 38.327, 37.990, 37.661, 37.341, 37.028, 36.723, 36.425,
      36.134, 35.850, 35.572, 35.300, 35.035, 34.775, 34.521, 34.272, 34.028, 33.789,
      33.555, 33.326, 33.102, 32.881, 32.665, 32.453, 32.245, 32.040, 31.839, 31.642,
      31.448, 31.257, 31.070, 30.885, 30.704, 30.525, 30.349, 30.176, 30.005, 29.837,
      29.671, 29.508, 29.347, 29.188, 29.031, 28.876, 28.724, 28.573, 28.424, 28.277,
      28.131, 27.988, 27.846, 27.706, 27.567, 27.430, 27.294, 27.159, 27.026, 26.895
    )
  )

  expect_equal(
    holland(seq(1, 100), 20, 50, 1020 * 100, 915 * 100, -15),
    c(
      25.621, 31.978, 36.024, 38.946, 41.180, 42.944, 44.365, 45.524, 46.474, 47.257,
      47.901, 48.428, 48.856, 49.199, 49.469, 49.675, 49.825, 49.926, 49.983, 50.001,
      49.984, 49.937, 49.861, 49.761, 49.638, 49.495, 49.333, 49.155, 48.962, 48.755,
      48.536, 48.306, 48.066, 47.817, 47.559, 47.295, 47.023, 46.746, 46.463, 46.175,
      45.883, 45.587, 45.288, 44.985, 44.681, 44.373, 44.064, 43.753, 43.441, 43.128,
      42.814, 42.499, 42.183, 41.868, 41.552, 41.236, 40.921, 40.606, 40.291, 39.977,
      39.664, 39.352, 39.040, 38.729, 38.420, 38.112, 37.805, 37.499, 37.195, 36.892,
      36.591, 36.291, 35.993, 35.696, 35.401, 35.108, 34.817, 34.527, 34.239, 33.953,
      33.669, 33.386, 33.106, 32.827, 32.550, 32.275, 32.003, 31.732, 31.463, 31.196,
      30.931, 30.667, 30.406, 30.147, 29.890, 29.635, 29.381, 29.130, 28.881, 28.633
    )
  )

  expect_equal(
    boose(seq(1, 100), 20, 50, 1020 * 100, 915 * 100, 0.1, 0.1, 0.25, 0.25, 15, 1, -15),
    c(
      19.786, 23.774, 26.206, 27.913, 29.192, 30.186, 30.975, 31.611, 32.129, 32.551,
      32.895, 33.176, 33.402, 33.583, 33.724, 33.831, 33.909, 33.961, 33.991, 34.000,
      33.991, 33.967, 33.928, 33.877, 33.814, 33.741, 33.659, 33.569, 33.471, 33.366,
      33.255, 33.138, 33.017, 32.891, 32.760, 32.626, 32.489, 32.349, 32.205, 32.060,
      31.912, 31.762, 31.610, 31.457, 31.302, 31.146, 30.989, 30.831, 30.672, 30.513,
      30.353, 30.192, 30.031, 29.870, 29.708, 29.546, 29.385, 29.223, 29.061, 28.900,
      28.738, 28.577, 28.416, 28.256, 28.095, 27.935, 27.776, 27.617, 27.459, 27.301,
      27.143, 26.986, 26.830, 26.674, 26.519, 26.365, 26.211, 26.058, 25.905, 25.754,
      25.603, 25.453, 25.303, 25.154, 25.006, 24.859, 24.712, 24.567, 24.422, 24.278,
      24.134, 23.992, 23.850, 23.709, 23.569, 23.429, 23.291, 23.153, 23.016, 22.880
    )
  )
})

test_that("Test spatialBehaviour function", {
  suppressWarnings(sds <- defStormsDataset(verbose = 0))
  pam <- defStormsList(sds = sds, loi = "Vanuatu", names = "PAM", verbose = 0)
  sb <- spatialBehaviour(pam, method = "Holland", verbose = 0)

  expect_equal(
    terra::minmax(sb[["PAM_MSW"]]),
    matrix(
      c(6.846, 76.733942), 
      dimnames =list(c('min', 'max'), c('PAM_MSW'))
    )
  )

  sb <- spatialBehaviour(pam, method = "Boose", verbose = 0)
  expect_equal(
    terra::minmax(sb[["PAM_MSW"]]),
    matrix(
      c(8.086, 73.804678), 
      dimnames =list(c('min', 'max'), c('PAM_MSW'))
    )
  )
})



test_that("Test checkInputsSpatialBehaviour function", {
  # Checking sts input
  expect_error(checkInputsSpatialBehaviour(
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "Willoughby",
    asymmetry = "Chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  # Checking product input
  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = 1,
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "Willoughby",
    asymmetry = "Chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = "TS",
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "Willoughby",
    asymmetry = "Chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  # Check windThreshold input
  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c("18", "33", "42", "49", "58", "70"),
    method = "Willoughby",
    asymmetry = "Chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(-3, -45),
    method = "Willoughby",
    asymmetry = "Chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  # Checking method input
  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "willoughby",
    asymmetry = "Chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = 1,
    asymmetry = "Chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = TRUE,
    asymmetry = "Chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = c("Willoughby", "Holland"),
    asymmetry = "Chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  # Checking asymmetry input
  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "Willoughby",
    asymmetry = "chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "Willoughby",
    asymmetry = 1,
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "Willoughby",
    asymmetry = TRUE,
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "Willoughby",
    asymmetry = c("Chen", "None"),
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  # Checking empiricalRMW inputs
  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "Willoughby",
    asymmetry = "Chen",
    empiricalRMW = "TRUE",
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "Willoughby",
    asymmetry = "Chen",
    empiricalRMW = 1,
    spaceRes = "2.5min",
    tempRes = 1,
    verbose = 2
  ))

  # Checking tempRes input
  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "Willoughby",
    asymmetry = "Chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = TRUE,
    verbose = 2
  ))

  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "Willoughby",
    asymmetry = "Chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = 2.3,
    verbose = 2
  ))

  expect_error(checkInputsSpatialBehaviour(
    sts = pam,
    product = c("MSW", "PDI", "Exposure"),
    windThreshold = c(18, 33, 42, 49, 58, 70),
    method = "Willoughby",
    asymmetry = "Chen",
    empiricalRMW = FALSE,
    spaceRes = "2.5min",
    tempRes = c(0.25, 0.5),
    verbose = 2
  ))
})





test_that("Test getIndices function", {
  suppressWarnings(sds <- defStormsDataset(verbose = 0))
  pam <- defStormsList(sds = sds, loi = "Vanuatu", names = "PAM", verbose = 0)

  expect_equal(getIndices(pam@data[["PAM"]], 2, "MSW"), seq(26, 49))
  expect_equal(getIndices(pam@data[["PAM"]], 20, "MSW"), seq(8, 57))
  expect_equal(getIndices(pam@data[["PAM"]], 30, "MSW"), seq(1, 57))
  expect_equal(getIndices(pam@data[["PAM"]], 30, "Profiles"), seq(28, 47))
})





test_that("Test getDataInterpolate function", {
  suppressWarnings(sds <- defStormsDataset(verbose = 0))
  pam <- defStormsList(sds = sds, loi = "Vanuatu", names = "PAM", verbose = 0)

  expect_equal(getDataInterpolate(pam@data[["PAM"]], seq(26, 49), 4, 3, FALSE, "Willoughby"), dfGetDataInterpolate)
})





test_that("Test computeDirection function", {
  expect_equal(computeDirection(1, 1, 30), 315)
  expect_equal(computeDirection(1, -1, 30), 45)
  expect_equal(computeDirection(-1, -1, 30), 135)
  expect_equal(computeDirection(-1, 1, 30), 225)

  expect_equal(computeDirection(0, 1, 30), 270)
  expect_equal(computeDirection(1, 0, 30), 0)
  expect_equal(computeDirection(0, -1, 30), 90)
  expect_equal(computeDirection(-1, 0, 30), 180)

  expect_equal(computeDirection(1, 1, -30), 135)
  expect_equal(computeDirection(1, -1, -30), 225)
  expect_equal(computeDirection(-1, -1, -30), 315)
  expect_equal(computeDirection(-1, 1, -30), 45)

  expect_equal(computeDirection(0, 1, -30), 90)
  expect_equal(computeDirection(1, 0, -30), 180)
  expect_equal(computeDirection(0, -1, -30), 270)
  expect_equal(computeDirection(-1, 0, -30), 360)
})


test_that("Test computeDirectionBoose function", {
  expect_equal(computeDirectionBoose(1, 1, 30, 1), 275)
  expect_equal(computeDirectionBoose(1, 1, 30, 0), 295)
  expect_equal(computeDirectionBoose(1, -1, 30, 1), 5)
  expect_equal(computeDirectionBoose(1, -1, 30, 0), 25)
  expect_equal(computeDirectionBoose(-1, 1, 30, 1), 185)
  expect_equal(computeDirectionBoose(-1, 1, 30, 0), 205)

  expect_equal(computeDirectionBoose(1, 1, -30, 1), 175)
  expect_equal(computeDirectionBoose(1, 1, -30, 0), 155)
  expect_equal(computeDirectionBoose(-1, -1, -30, 1), 355)
  expect_equal(computeDirectionBoose(-1, -1, -30, 0), 335)
})