def test_reaction_rate_str():
try:
vp1 = np.array([[0.5,-2],['j3',0],[0,1]])
vpp1 = np.array([[0,0],[1,1],[2,1]])
x1 = np.array([[4],[2],[1]])
reaction_rate(vp1,vpp1,x1,2)
except TypeError as err:
assert(type(err)==TypeError)
def test_reaction_rate_neg():
vp1 = np.array([[0.5,-2],[3,0],[0,1]])
vpp1 = np.array([[0,0],[1,1],[2,1]])
x1 = np.array([[4],[2],[1]])
assert(reaction_rate(vp1,vpp1,x1,2)==[-15.75,-63.875,64.0])
vp2 = np.array([[1,2],[2,1],[0,2]])
vpp2 = np.array([[0,0],[0,1],[2,1]])
x2 = np.array([[4],[-2],[1]])
assert(reaction_rate(vp2,vpp2,x2,2.5)==[120.0,-80.0,160.0])
def test_reaction_rate_len():
try:
vp1 = np.array([[0.5,-2,7],[3,0,1],[0,1,2]])
vpp1 = np.array([[0,0],[1,1],[2,1]])
x1 = np.array([[4],[2],[1]])
reaction_rate(vp1,vpp1,x1,2)
except ValueError as err:
assert(type(err)==ValueError)
try:
vp2 = np.array([[0.5,-2],[3,0],[0,1]])
vpp2 = np.array([[0,0],[1,1],[2,1]])
x2 = np.array([[4,1],[2,2],[1,3]])
reaction_rate(vp2,vpp2,x2,2)
except ValueError as err:
assert(type(err)==ValueError)
try:
vp3 = np.array([[0.5,-2],[3,0],[0,1]])
vpp3 = np.array([[0,0],[1,1],[2,1]])
x3 = np.array([[4],[2],[1],[2]])
reaction_rate(vp3,vpp3,x3,2)
except ValueError as err:
assert(type(err)==ValueError)
try:
vp4 = np.array([[0.5,-2],[3,0],[0,1],[1,2]])
vpp4 = np.array([[0,0],[1,1],[2,1],[0,3]])
x4 = np.array([[4],[2],[1]])
reaction_rate(vp4,vpp4,x4,2)
except ValueError as err:
assert(type(err)==ValueError)
def test_reaction_rate_int():
vp1 = np.array([[0.5,2],[3,0],[0,1]])
vpp1 = np.array([[0,0],[1,1],[2,1]])
x1 = np.array([[1],[2],[1]])
assert(reaction_rate(vp1,vpp1,x1,5)==[-30.0,-75.0,80.0])
vp2 = np.array([[0.5,0.5],[2,1],[0,2]])
vpp2 = np.array([[0,0],[0,1],[2,1]])
x2 = np.array([[4],[2],[1]])
assert(reaction_rate(vp2,vpp2,x2,2.5)==[-15.0,-40.0,30.0])
vp3 = np.array([[1,2],[1,2],[0,1]])
vpp3 = np.array([[0,0],[0,1],[2,1]])
x3 = np.array([[0.5],[4],[2]])
assert(reaction_rate(vp3,vpp3,x3,2)==[-36.0,-20.0,8.0])
def test_reaction_rate():
assert chemkin.reaction_rate([[3.0, 1.0, 1.0]], [[1.0, 2.0, 1.0]],
[1.0, 2.0, 3.0], 10) == [-10.0]
def reaction_rate():
try:
chemkin.reaction_rate(np.array([[1.0, 2.0, 0.0],[2.0, 0.0, 2.0]]),np.array([[0.0, 0.0, 2.0], [0.0, 1.0, 1.0]]),\
["", "green", "hi"], [10,10])
except TypeError as err:
assert(type(err) == TypeError)