Skip to content

Commit

Permalink
remove unsupported/unused code
Browse files Browse the repository at this point in the history
  • Loading branch information
exaexa committed Jan 8, 2025
1 parent f00e48c commit b3acc60
Show file tree
Hide file tree
Showing 2 changed files with 18 additions and 5 deletions.
14 changes: 12 additions & 2 deletions docs/src/examples/07a-srba.jl
Original file line number Diff line number Diff line change
Expand Up @@ -42,7 +42,7 @@ download_model(
import AbstractFBCModels as A
import JSONFBCModels
import ConstraintTrees as C
import SCIP
import HiGHS

# ## Collect data for RBA model
# RBA models require a lot of data. Below we have processed the data into a
Expand Down Expand Up @@ -493,6 +493,11 @@ end

mus = range(0.1, 0.95, 10) # simulate at these growth rates

settings = [
set_optimizer_attribute("solver", "simplex"),
set_optimizer_attribute("simplex_strategy", 4),
]

res = screen(mus) do mu
@info "sRBA step" mu
rbat = with_srba_constraints(ct, mu)
Expand All @@ -512,8 +517,9 @@ res = screen(mus) do mu
rbat;
settings = [silence],
objective = rbat.lp_objective.value,
optimizer = SCIP.Optimizer,
optimizer = HiGHS.Optimizer,
sense = Minimal,
settings
)
isnothing(sol) && return nothing

Expand All @@ -529,6 +535,8 @@ end

# finally, we can plot the data, to see if we can recapitulate known phenomena

#= TODO
using CairoMakie
# load measured ribosome protein mass fractions
Expand All @@ -554,3 +562,5 @@ lines!(ax2, mus, [abs(r.glc_flux) for r in res], label = "Glucose")
lines!(ax2, mus, [abs(r.o2_flux) for r in res], label = "Oxygen")
axislegend(ax2, position = :lt)
fig
=#
9 changes: 6 additions & 3 deletions docs/src/examples/07b-community-ecfba.jl
Original file line number Diff line number Diff line change
Expand Up @@ -45,7 +45,7 @@ download_model(
import AbstractFBCModels as A
import JSONFBCModels
import ConstraintTrees as C
import SCIP
import HiGHS

# ## Collect data for enzyme constrained models
# Like in the previous example, we will spend some time gathering the necessary
Expand Down Expand Up @@ -356,10 +356,9 @@ function auxotrophe_fba(wt, aa_ko; fbc_only = false)

sol = optimized_values(
x;
optimizer = SCIP.Optimizer,
optimizer = HiGHS.Optimizer,
objective = x.objective.value,
sense = Maximal,
settings = [silence],
)
isnothing(sol) && return nothing

Expand Down Expand Up @@ -444,6 +443,8 @@ eccfba_res = screen(specs) do spec
auxotrophe_fba(wt, spec)
end

#=
# ## Plot cFBA vs. ec-cFBA vs. experimental data
# Now that we have simulated the models, we need to compare them to data. The
# data from the paper is processed below, and then plotted.
Expand Down Expand Up @@ -513,3 +514,5 @@ stem!(
xlims!(ax2, 0, 1)
ylims!(ax2, 0, 1)
fig
=#

0 comments on commit b3acc60

Please sign in to comment.