def test_point_calculation(self):
expected = [point(i, 0, 0) for i in [
Dl("0"), Dl("29.41800004"), Dl("56.38450008"),
Dl("80.89950012"), Dl("102.96300016"), Dl("122.5750002"),
Dl("139.73550024"), Dl("154.44450028"), Dl("166.70200032"),
Dl("176.50800036"), Dl("183.8625004"), Dl("188.76550044"),
Dl("191.21700048"), Dl("191.21700052"), Dl("188.76550056"),
Dl("183.8625006"), Dl("176.50800064"), Dl("166.70200068"),
Dl("154.44450072"), Dl("139.73550076"), Dl("122.5750008"),
Dl("102.96300084"), Dl("80.89950088"), Dl("56.38450092"),
Dl("29.41800096"), Dl("0.000001")]
]
dt = Dl("0.5")
vel = vector(Dl("61.28750008"), 0, 0)
accn = vector(Dl("-9.806"), 0, 0)
proc = trajectory()
proc.send(None)
for n, x in enumerate(expected):
if n == 0:
imp = proc.send(Impulse(
Dl(0), Dl("0.5"), accn, point(0, 0, 0)))
elif n == 1:
p = (
point(0, 0, 0) + vel * dt +
Dl("0.5") * accn * dt * dt
)
imp = proc.send(Impulse(
imp.tEnd, imp.tEnd + dt, accn, p))
else:
imp = proc.send(Impulse(
imp.tEnd, imp.tEnd + dt, accn, imp.pos))
with self.subTest(n=n):
self.assertEqual(x, imp.pos)
def test_point_calculation(self):
expected = [point(i, 0, 0) for i in [
Dl("0"), Dl("0.502656"), Dl("0.824448"), Dl("0.986112"),
Dl("1.008384"), Dl("0.912"), Dl("0.717696"),
Dl("0.446208"), Dl("0.118272"), Dl("-0.245376"),
Dl("-0.624"), Dl("-0.996864"), Dl("-1.343232"),
Dl("-1.642368"), Dl("-1.873536"), Dl("-2.016"),
Dl("-2.049024"), Dl("-1.951872"), Dl("-1.703808"),
Dl("-1.284096"), Dl("-0.672"), Dl("0.153216"),
Dl("1.212288"), Dl("2.525952"), Dl("4.114944"), Dl("6")]
]
accns = [vector(i, 0, 0) for i in [
Dl("-14"), Dl("-12.56"), Dl("-11.12"), Dl("-9.68"),
Dl("-8.24"), Dl("-6.8"), Dl("-5.36"), Dl("-3.92"),
Dl("-2.48"), Dl("-1.04"), Dl("0.4"), Dl("1.84"),
Dl("3.28"), Dl("4.72"), Dl("6.16"), Dl("7.6"), Dl("9.04"),
Dl("10.48"), Dl("11.92"), Dl("13.36"), Dl("14.8"),
Dl("16.24"), Dl("17.68"), Dl("19.12"), Dl("20.56"),
Dl("22")]
]
dt = Dl("0.12")
proc = trajectory()
proc.send(None)
for n, x in enumerate(expected):
accn = accns[n]
if n == 0:
p = point(Dl(0), Dl(0), Dl(0))
imp = proc.send(Impulse(Dl(0), dt, accn, p))
elif n == 1:
p = point(Dl("0.502656"), Dl(0), Dl(0))
imp = proc.send(Impulse(
imp.tEnd, imp.tEnd + dt, accn, p))
else:
imp = proc.send(Impulse(
imp.tEnd, imp.tEnd + dt, accn, imp.pos))
with self.subTest(n=n):
self.assertEqual(x, imp.pos)