def steady_state_function(self, function, parameter_values): if isinstance(function, Expression): expr = function else: expr = Expression(function) p_vals = parameter_values.copy() p_vals.update(self.steady_state(parameter_values, log_out=False)) p_vals.update(self.steady_state_flux(parameter_values, log_out=False)) for i in self.dependent_variables: for j in self.independent_variables: p_vals["$L_" + i + "_" + j] = self.log_gain(i, j) value = expr.eval_with_values(p_vals=p_vals) return value
def steady_state_function(self, function, parameter_values): if isinstance(function, Expression): expr = function else: expr = Expression(function) p_values = VariablePool(names=self.independent_variables) for key in p_values: p_values[key] = parameter_values[key] ss = self.steady_state(p_values, log_out=False) flux = self.steady_state_flux(p_values, log_out=False) values = dict() for subcase in ss: p_values.update(ss[subcase]) p_values.update(flux[subcase]) value = expr.eval_with_values(p_vals=p_values) values[subcase] = value return values