test_that("AEL outputs right output", {
    # -----------------------------
    # Initialise Variables
    # -----------------------------
    # Generate toy variables
    set.seed(1)
    x    <- runif(30, min = -5, max = 5)
    elip <- rnorm(30, mean = 0, sd = 1)
    y    <- 0.75 - x + elip
    T    <- 10
    
    # Set initial values for AEL computation
    lam0 <- matrix(c(0,0), nrow = 2)
    th   <- matrix(c(0.8277, -1.0050), nrow = 2)
    a    <- 0.00001
    
    # Define Dataset and h-function
    z    <- cbind(x, y)
    h    <- function(z, th) {
        xi <- z[1]
        yi <- z[2]
        h_zith <- c(yi - th[1] - th[2] * xi, xi*(yi - th[1] - th[2] * xi))
        matrix(h_zith, nrow = 2)
    }
    
    # -----------------------------
    # Main
    # -----------------------------
    ansAELRcppprez  <- compute_AEL(th, h, lam0, a, z, useR_forz = TRUE)
    ansAELRcpp      <- compute_AEL(th, h, lam0, a, z, useR_forz = FALSE)
    
    expect_length(ansAELRcppprez, 1)
    expect_length(ansAELRcpp, 1)
    
    # Testing for discrepencies, FALSE if the same (might be floating point errors even with rounding)
    expect_equal(ansAELRcppprez, -106.45360424034257107)
    expect_equal(ansAELRcpp,     -106.45360424034257107)
    # expect_equal(round(ansAELRcpp,5), round(ansAELRcppprez,5))
    
    set.seed(NULL) # Reset seed
})