context("check if the jacobian matrix contains NA") load("groups.level.RData") test_that('all values defined',{ mat.fluxes = fluxing(groups.level$mat, groups.level$biomasses, 0.71*groups.level$bodymasses, groups.level$efficiencies) met.types = rep("animal", nrow(mat.fluxes)) met.types[groups.level$efficiencies == 0.545] = "plant" met.types[1] = "detritus" # not true, only for testing purpose jacob = create.jacob(mat.fluxes, groups.level$biomasses, groups.level$efficiencies, met.types) expect_equal(0, sum(is.na(jacob))) })