def test_merge_incompartible_trajs(self):
""" Test the incompatible trajectories cause an exception to be raised """
net1 = algebraic_net.copy()
net2 = algebraic_net_assignment.copy()
times1 = np.linspace(0, 15, 100)
times2 = np.linspace(15, 25, 100)
traj1 = Dynamics.integrate(net1, times1, fill_traj=False)
traj2 = Dynamics.integrate(net2, times2, fill_traj=False)
self.assertRaises(ValueError, traj1.merge, traj2)
from __future__ import absolute_import
from builtins import range
import unittest
import os
import numpy as np
from SloppyCell.ReactionNetworks import *
from AlgTestNets import algebraic_net_assignment, algebraic_net
base_net = algebraic_net_assignment.copy()
base_net.remove_component('event0')
base_net.set_var_constant('k1', False)
class test_Trajectory(unittest.TestCase):
def test_interpolate_basic(self):
net = base_net.copy('test')
traj = Dynamics.integrate(net, [0, 5.5, 15])
traj.build_interpolated_traj()
tlist_algebraic_net = np.array([0] + [0.8 * x for x in range(1, 15)])
times, interp_traj = traj.evaluate_interpolated_traj(
'X0', tlist_algebraic_net)
self.assertAlmostEqual(interp_traj[1], 0.923116346390738, 5)
self.assertAlmostEqual(interp_traj[4], .726149037073939, 5)
self.assertAlmostEqual(interp_traj[6], 0.618783391805766, 5)
self.assertAlmostEqual(interp_traj[8], 0.527292424043448, 5)
self.assertAlmostEqual(interp_traj[10], 0.449328964119144, 5)
self.assertAlmostEqual(interp_traj[14], 0.326279794615302, 5)