Commit 1c7529d2 authored by Miroslav Kratochvil's avatar Miroslav Kratochvil
Browse files

the usual set of fixes

parent 35a73fab
......@@ -16,7 +16,7 @@ function atom_exchange(model::MetabolicModel, flux_dict::Dict{String,Float64})
atoms = metabolite_formula(model, mids[midx])
isnothing(atoms) && continue
for (atom, abundance) in atoms
atom_flux[atom] = get(atom_flux, atom, 0.0) + flux * mstoi * abundance
atom_flux[atom] = get(atom_flux, atom, 0.0) + rflux * mstoi * abundance
end
end
end
......
......@@ -17,7 +17,7 @@
)
# bounds setting
cbm = make_optimization_model(model, optimizer)
cbm = make_optimization_model(model, Tulip.Optimizer)
ubs = cbm[:ubs]
lbs = cbm[:lbs]
glucose_index = first(indexin(["EX_glc__D_e"], reactions(model)))
......
......@@ -8,7 +8,7 @@
)
# atom tracker
atom_fluxes = atom_exchange(fluxes, model)
atom_fluxes = atom_exchange(model, fluxes)
@test isapprox(atom_fluxes["C"], 37.19016648975907; atol = TEST_TOLERANCE)
@test isapprox(atom_fluxes["N"], 37.19016648975907; atol = TEST_TOLERANCE)
end
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment