# from roadrunner import RoadRunner
rr_to_1 = RoadRunnerPy(sbml=rr_from.getCurrentSBML())
rr_to_1.stepSize = step_size
print
for i in range(20):
rr_from.timestep()
rr_to_1.timestep()
print 'STEP=',i
print 'orig S1=',rr_from.model['S1'],' new S1=',rr_to_1.model['S1']
print 'orig S2=', rr_from.model['S2'], ' new S2=', rr_to_1.model['S2']
rr_to = RoadRunnerPy(_path=path_from)
rr_to.stepSize = step_size
rr_to.loadSBML(_externalPath=path_from)
# # setting initial conditions - this has to be done after loadingSBML
# for name, value in initial_conditions_from.iteritems():
# if name.startswith('init('): continue
# try: # have to catch exceptions in case initial conditions contain "unsettable" entries such as reaction rate etc...
# rr_to.model[name] = value
# except:
# pass
#
independent_floating_species = rr.getIndependentFloatingSpeciesIds()
print independent_floating_species
boundary_species = rr.getIndependentFloatingSpeciesIds()
