test_that("Geometric mean calculation", {

    x <- c(1,2,3,4,5,6,7,8,9)
    expect_equal(QFASA:::mean_geometric(x), exp(mean(log(x))))
})