diff --git a/src/base/types/FluxSummary.jl b/src/base/types/FluxSummary.jl
index 470fe75b928ba0c5b18f8731dbe409d68b1019fd..3f99390003b37c9a97d7ac48b6da24dfd636c0a5 100644
--- a/src/base/types/FluxSummary.jl
+++ b/src/base/types/FluxSummary.jl
@@ -19,7 +19,6 @@ end
                 exclude_biomass = false,
                 small_flux_bound = 1.0/_constants.default_reaction_bound^2,
                 large_flux_bound = _constants.default_reaction_bound,
-                round_digits = 3,
                 keep_unbounded = false,
                 )::FluxSummary
 
diff --git a/test/base/types/FluxSummary.jl b/test/base/types/FluxSummary.jl
new file mode 100644
index 0000000000000000000000000000000000000000..7b4e40ed8f71a0600af439f84e0d3313780addbf
--- /dev/null
+++ b/test/base/types/FluxSummary.jl
@@ -0,0 +1,18 @@
+@testset "Flux summary" begin
+    model = load_model(model_paths["e_coli_core.json"])
+
+    sol = flux_balance_analysis_dict(
+        model,
+        Tulip.Optimizer;
+        modifications = [
+            change_optimizer_attribute("IPM_IterationsLimit", 200),
+        ],
+    )
+
+    fr = flux_summary(sol; keep_unbounded=true, large_flux_bound=25)
+
+    @test isapprox(fr.biomass_fluxes["BIOMASS_Ecoli_core_w_GAM"], 0.8739215022690006; atol=TEST_TOLERANCE)
+    @test isapprox(fr.export_fluxes["EX_co2_e"], 22.80983339307183; atol=TEST_TOLERANCE)
+    @test isapprox(fr.import_fluxes["EX_o2_e"], -21.799492758430517; atol=TEST_TOLERANCE)
+    @test isapprox(fr.unbounded_fluxes["EX_h2o_e"], 29.175827202663395; atol=TEST_TOLERANCE) 
+end
diff --git a/test/base/types/FluxVariabilitySummary.jl b/test/base/types/FluxVariabilitySummary.jl
new file mode 100644
index 0000000000000000000000000000000000000000..24409e3bca50cbf813831298a098aaafecd7b574
--- /dev/null
+++ b/test/base/types/FluxVariabilitySummary.jl
@@ -0,0 +1,13 @@
+@testset "Flux variability summary" begin
+    model = load_model(model_paths["e_coli_core.json"])
+
+    sol = flux_balance_analysis_dict(
+        model,
+        Tulip.Optimizer;
+        modifications = [
+            change_optimizer_attribute("IPM_IterationsLimit", 200),
+        ],
+    )
+
+    fr = flux_variability_summary(sol)
+end