From 11a6549896d977a2af06fe763d75042a865c1ba8 Mon Sep 17 00:00:00 2001 From: tfaour Date: Mon, 2 Jun 2025 12:28:39 -0400 Subject: [PATCH] Delete orbiter/units.py --- orbiter/units.py | 74 ------------------------------------------------ 1 file changed, 74 deletions(-) delete mode 100644 orbiter/units.py diff --git a/orbiter/units.py b/orbiter/units.py deleted file mode 100644 index be121ae..0000000 --- a/orbiter/units.py +++ /dev/null @@ -1,74 +0,0 @@ -import numpy as np -from decimal import Decimal - -#Declare these as units to make code clearer -class HighPrecisionVector(np.ndarray): - def __new__(cls, input_array, *args, **kwargs): - decimal_array = [Decimal(i) for i in input_array] - obj = np.asarray(decimal_array).view(cls) - return obj - -class HighPrecisionMatrix(np.ndarray): - def __new__(cls, dim1, dim2, *args, **kwargs): - decimal_array = [Decimal(0) for _ in range(dim1)] - decimal_matrix = [decimal_array for _ in range(dim2)] - obj = np.asarray(decimal_matrix).view(cls) - return obj - -class Position(HighPrecisionVector): - pass -class Velocity(HighPrecisionVector): - pass -class Acceleration(HighPrecisionVector): - pass - -Mass = Decimal - -#################### -# USEFUL CONSTANTS -#################### -EARTH_MASS = Decimal(5972 * 10**21) #kg -EARTH_RADIUS = Decimal(6378 * 10**3) #meters -EARTH_ORBITAL_VELOCITY = Decimal(29780) # m/s -AU = Decimal(149_597_870_700) #meters - -MOON_MASS = Decimal(734767309 * 10**14) -MOON_ORBITAL_VELOCITY = Decimal(1022) #m/s relative to earth - -SUN_MASS = Decimal(1989 * 10**27) #kg -SUN_RADIUS = Decimal(6957 * 10**5) #meters - -pi_approx = Decimal("3.14159265358979323846264338327950288419716939937510") - - -#NORMALIZING CONSTANTS -G = Decimal(6.67430e-11) -r_0 = Decimal(EARTH_RADIUS) #1.496e11 -m_0 = Decimal(5.972e24) -t_0 = np.sqrt((r_0**3) / (G*m_0)) - -def norm_pos(pos): - return Decimal(pos) / Decimal(r_0) - -def real_pos(pos): - return pos * Decimal(r_0) - -def norm_mass(mass): - return Decimal(mass) / Decimal(m_0) - -def real_mass(mass): - return Decimal(mass) * Decimal(m_0) - -def norm_time(time): - return Decimal(time) / Decimal(t_0) - -def real_time(time): - return Decimal(time) * Decimal(t_0) - -def norm_vel(vel): - return vel / Decimal(r_0/t_0) - -def real_vel(vel): - return vel * Decimal(r_0/t_0) - -