projectile-sim 1.1.0__tar.gz

projectile_sim-1.1.0/PKG-INFO

@@ -0,0 +1,4 @@
+Metadata-Version: 2.4
+Name: projectile-sim
+Version: 1.1.0
+Requires-Dist: numpy

projectile_sim-1.1.0/README.md

@@ -0,0 +1,2 @@
+# projectile-sim
+by inputting various parameters, find the maximum height reached, total horizontal distance traveled, and time it takes to land for an ideal kinematic projectile problem.

projectile_sim-1.1.0/projectile/__init__.py

@@ -0,0 +1 @@
+from .projectile_sim import Projectile

projectile_sim-1.1.0/projectile/projectile_sim.py

@@ -0,0 +1,51 @@
+# OOP projectile motion
+import math
+import numpy as np
+import matplotlib.pyplot as plt
+
+class Particle(object):
+    def __init__(self):
+        self.position = [0.0, 0.0]
+
+class Projectile(Particle):
+    def __init__(self, velocity, angle, initial_height, dt):
+        super().__init__()
+        self.velocity = velocity
+        self.angle = angle
+        self.initial_height = initial_height
+        self.position = [0.0, initial_height]
+        self.time = 0.0
+        self.dt = dt
+        self.g = -9.81
+
+    def update_position(self, dt):
+        self.time += dt
+        self.position[0] = self.velocity * math.cos(math.radians(self.angle)) * self.time
+        self.position[1] = self.initial_height + (self.velocity * math.sin(math.radians(self.angle)) * self.time) + (0.5 * self.g * self.time ** 2)
+
+    def get_position(self):
+        return tuple(self.position)
+    
+    def get_time(self):
+        return self.time
+    
+    def get_max_height(self):
+        return self.initial_height + (self.velocity ** 2 * (math.sin(math.radians(self.angle))) ** 2) / (2 * -self.g)
+
+if __name__ == "__main__":
+    ball = Projectile(velocity=10.0, angle=45.0, initial_height=10.0, dt=0.01)
+    print(f"Initial x, y: {ball.get_position()} (m)")
+
+    x_positions = [ball.get_position()[0]]
+    y_positions = [ball.get_position()[1]]
+
+    while ball.position[1] >= 0:
+        ball.update_position(ball.dt)
+        x_positions.append(ball.get_position()[0])
+        y_positions.append(ball.get_position()[1])
+
+    print(f"Final x, y: ...")
+
+    fig, ax = plt.subplots(figsize=(8, 5))  # <-- moves in here
+    ax.plot(x_positions, y_positions)
+    plt.show()

projectile_sim-1.1.0/projectile_sim.egg-info/PKG-INFO

@@ -0,0 +1,4 @@
+Metadata-Version: 2.4
+Name: projectile-sim
+Version: 1.1.0
+Requires-Dist: numpy

projectile_sim-1.1.0/projectile_sim.egg-info/SOURCES.txt

@@ -0,0 +1,13 @@
+README.md
+pyproject.toml
+projectile/__init__.py
+projectile/projectile_sim.py
+projectile_sim.egg-info/PKG-INFO
+projectile_sim.egg-info/SOURCES.txt
+projectile_sim.egg-info/dependency_links.txt
+projectile_sim.egg-info/requires.txt
+projectile_sim.egg-info/top_level.txt
+tests/test_edgecases.py
+tests/test_outputs.py
+tests/test_params.py
+tests/test_postime.py

projectile_sim-1.1.0/projectile_sim.egg-info/dependency_links.txt

@@ -0,0 +1 @@
+

projectile_sim-1.1.0/projectile_sim.egg-info/requires.txt

@@ -0,0 +1 @@
+numpy

projectile_sim-1.1.0/projectile_sim.egg-info/top_level.txt

@@ -0,0 +1 @@
+projectile

projectile_sim-1.1.0/pyproject.toml

@@ -0,0 +1,12 @@
+[build-system]
+requires = ["setuptools>=61.0"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "projectile-sim"
+version = "1.1.0"
+
+dependencies = [
+    "numpy",
+
+]

projectile_sim-1.1.0/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

projectile_sim-1.1.0/tests/test_edgecases.py

@@ -0,0 +1,15 @@
+from projectile import Projectile
+import math
+def test_90_degree_angle():
+    # straight up — x should stay ~0
+    ball = Projectile(velocity=10.0, angle=90.0, initial_height=0.0, dt=0.01)
+    while ball.position[1] >= 0:
+        ball.update_position(ball.dt)
+    assert abs(ball.get_position()[0]) < 0.01
+
+def test_zero_velocity():
+    # with no velocity, ball should just fall from initial height
+    ball = Projectile(velocity=0.0, angle=45.0, initial_height=10.0, dt=0.01)
+    while ball.position[1] >= 0:
+        ball.update_position(ball.dt)
+    assert abs(ball.get_position()[0]) < 0.01  # no horizontal motion

projectile_sim-1.1.0/tests/test_outputs.py

@@ -0,0 +1,24 @@
+from projectile import Projectile
+import math
+def test_max_height():
+    # analytical formula: h_max = h0 + (v * sin(angle))^2 / (2g)
+    ball = Projectile(velocity=10.0, angle=90.0, initial_height=0.0, dt=0.01)
+    expected = (10.0 ** 2) / (2 * 9.81)  # straight up shot, simplifies nicely
+    assert abs(ball.get_max_height() - expected) < 0.01
+
+def test_horizontal_distance_flat_ground():
+    # at angle=45, v=10, h0=0: range = v^2 * sin(2*angle) / g
+    ball = Projectile(velocity=10.0, angle=45.0, initial_height=0.0, dt=0.01)
+    expected_range = (10.0 ** 2 * math.sin(math.radians(90))) / 9.81
+    while ball.position[1] >= 0:
+        ball.update_position(ball.dt)
+    assert abs(ball.get_position()[0] - expected_range) < 0.1  # tolerance for dt discretization
+
+def test_y_is_always_decreasing_after_peak():
+    ball = Projectile(velocity=10.0, angle=45.0, initial_height=0.0, dt=0.01)
+    positions = []
+    while ball.position[1] >= 0:
+        ball.update_position(ball.dt)
+        positions.append(ball.get_position()[1])
+    peak_idx = positions.index(max(positions))
+    assert all(positions[i] >= positions[i+1] for i in range(peak_idx, len(positions)-1))

projectile_sim-1.1.0/tests/test_params.py

@@ -0,0 +1,9 @@
+from projectile import Projectile
+import math
+def test_initial_position():
+    ball = Projectile(velocity=10.0, angle=45.0, initial_height=5.0, dt=0.01)
+    assert ball.get_position() == (0.0, 5.0)
+
+def test_zero_initial_height():
+    ball = Projectile(velocity=10.0, angle=45.0, initial_height=0.0, dt=0.01)
+    assert ball.get_position()[1] == 0.0

projectile_sim-1.1.0/tests/test_postime.py

@@ -0,0 +1,13 @@
+from projectile import Projectile
+import math
+import pytest
+def test_time_increases():
+    ball = Projectile(velocity=10.0, angle=45.0, initial_height=0.0, dt=0.01)
+    ball.update_position(ball.dt)
+    assert ball.get_time() == pytest.approx(0.01)
+
+def test_time_accumulates():
+    ball = Projectile(velocity=10.0, angle=45.0, initial_height=0.0, dt=0.01)
+    for _ in range(5):
+        ball.update_position(ball.dt)
+    assert ball.get_time() == pytest.approx(0.05)