from parameters import *
from subsystem import *
from geometry import *
# ---- instantiation ----
# ---- release sites ----
rel_A = m.ReleaseSite(name='rel_A',
complex=m.Complex('A', orientation=m.Orientation.UP),
region=up,
number_to_release=100)
rel_B = m.ReleaseSite(name='rel_B',
complex=m.Complex('B', orientation=m.Orientation.UP),
region=bottom,
number_to_release=100)
# ---- surface classes assignment ----
# ---- compartments assignment ----
# ---- create instantiation object and add components ----
instantiation = m.Instantiation()
instantiation.add_geometry_object(up)
instantiation.add_geometry_object(bottom)
instantiation.add_release_site(rel_A)
instantiation.add_release_site(rel_B)
import os
MCELL_PATH = os.environ.get('MCELL_PATH', '')
if MCELL_PATH:
sys.path.append(os.path.join(MCELL_PATH, 'lib'))
else:
print("Error: variable MCELL_PATH that is used to find the mcell library was not set.")
sys.exit(1)
import mcell as m
# single property
box1 = m.geometry_utils.create_box('b1', 1)
box2 = m.geometry_utils.create_box('b1', 1.1)
i1 = m.Instantiation()
i1.add_geometry_object(box1)
# error
try:
i1.add_geometry_object(box2)
assert False
except ValueError as err:
print(err)
box3 = m.geometry_utils.create_box('b2', 1)
box4 = m.geometry_utils.create_box('b2', 1)
i2 = m.Instantiation()
i2.add_geometry_object(box3)
# warning
