Commit 9016bb36 authored by HEMERY Mathieu's avatar HEMERY Mathieu
Browse files

Avoid Error and Warning in independant use of numerical_simulation

parent 11dc1360
......@@ -379,12 +379,12 @@ list_hybrid_ode(Stream, Hybrid_original_list) :- % Only MA are converted, need t
/* Can try to get the kinetics function in BIOCHAM 4 directly and translate, but i didn't know how
Also can direclty adapt the kinetics with the hybrid species list, but need to translate the kinetics first if needed */
%! convert_MA()
%! convert_MA(+Reaction, +List_Reactions, +Hybrid_original_list, -NewKinetics)
%
%
convert_MA(Reaction, List_Reactions, Hybrid_original_list, NewKinetics) :-
reaction(Reaction, [kinetics: Kinetics,reactants: Reactants,inhibitors: _, products: _]),
reaction(Reaction, [kinetics: Kinetics, reactants: Reactants, inhibitors: _, products: _]),
Kinetics =.. ['MA',K],
kinetics:eval_kinetics(Reactants, _, product(in('S'*'M', [reactants]), 'M'^'S'), Value),
convert_into_hybrid(Value, Hybrid_original_list, std, NewValue),!,
......@@ -402,9 +402,9 @@ convert_MA(Reaction, List_Reactions, Hybrid_original_list, NewKinetics) :-
).
convert_MA2(Reaction, ODESpecies, NewKinetics, Reactants, Products) :-
reaction(Reaction, [kinetics: Kinetics,reactants: Reactants,inhibitors: _, products: Products]),
reaction(Reaction, [kinetics: Kinetics, reactants: Reactants, inhibitors: _, products: Products]),
Kinetics =.. ['MA',K],
kinetics:eval_kinetics(Reactants, _, product('S'*'M' in [reactants], 'M'^'S'), Value),
kinetics:eval_kinetics(Reactants, _, product(in('S'*'M', [reactants]), 'M'^'S'), Value),
convert_into_hybrid(Value, ODESpecies, stoch, NewValue),!,
NewKinetics = K * NewValue,
(
......@@ -496,6 +496,7 @@ build_stoechiometry(Reactant_list, Product_list, [Sp|TailSp], [Stoech|TailStoech
),
build_stoechiometry(Reactant_list, Product_list, TailSp, TailStoech).
% build_reaction_data(+List_of_reactions, +Species_list, -New_reaction_list)
%
% format the reaction data to the format: [Stoechiometry, Rate]
......
......@@ -18,6 +18,7 @@
:- dynamic(rosenbrock_running2/0).
% Only for separate compilation/linting
:- use_module(biocham).
:- use_module(doc).
:- use_module(objects).
......
Markdown is supported
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