verified type check changes

peterdsharpe · peterdsharpe · commit 4e78784ae422 · 2025-12-24T16:53:41.000+01:00
diff --git a/aerosandbox/aerodynamics/aero_3D/aero_buildup.py b/aerosandbox/aerodynamics/aero_3D/aero_buildup.py
@@ -3,6 +3,7 @@
 from aerosandbox.performance import OperatingPoint
 import aerosandbox.library.aerodynamics as aero
 import aerosandbox.numpy as np
+from aerosandbox.numpy.typing import Vectorizable
 from aerosandbox.aerodynamics.aero_3D.aero_buildup_submodels.fuselage_aerodynamics_utilities import *
 from aerosandbox.library.aerodynamics import transonic
 import aerosandbox.library.aerodynamics as aerolib
@@ -137,78 +138,78 @@ def __repr__(self):
             )
 
         @property
-        def F_b(self) -> list[float | np.ndarray]:
+        def F_b(self) -> tuple[Vectorizable, Vectorizable, Vectorizable]:
             """
-            An [x, y, z] list of forces in body axes [N]
+            An (x, y, z) tuple of forces in body axes [N]
             """
             return self.op_point.convert_axes(
                 *self.F_g, from_axes="geometry", to_axes="body"
             )
 
         @property
-        def F_w(self) -> list[float | np.ndarray]:
+        def F_w(self) -> tuple[Vectorizable, Vectorizable, Vectorizable]:
             """
-            An [x, y, z] list of forces in wind axes [N]
+            An (x, y, z) tuple of forces in wind axes [N]
             """
             return self.op_point.convert_axes(
                 *self.F_g, from_axes="geometry", to_axes="wind"
             )
 
         @property
-        def M_b(self) -> list[float | np.ndarray]:
+        def M_b(self) -> tuple[Vectorizable, Vectorizable, Vectorizable]:
             """
-            An [x, y, z] list of moments about body axes [Nm]
+            An (x, y, z) tuple of moments about body axes [Nm]
             """
             return self.op_point.convert_axes(
                 *self.M_g, from_axes="geometry", to_axes="body"
             )
 
         @property
-        def M_w(self) -> list[float | np.ndarray]:
+        def M_w(self) -> tuple[Vectorizable, Vectorizable, Vectorizable]:
             """
-            An [x, y, z] list of moments about wind axes [Nm]
+            An (x, y, z) tuple of moments about wind axes [Nm]
             """
             return self.op_point.convert_axes(
                 *self.M_g, from_axes="geometry", to_axes="wind"
             )
 
         @property
-        def L(self) -> float | np.ndarray:
+        def L(self) -> Vectorizable:
             """
             The lift force [N]. Definitionally, this is in wind axes.
             """
             return -self.F_w[2]
 
         @property
-        def Y(self) -> float | np.ndarray:
+        def Y(self) -> Vectorizable:
             """
             The side force [N]. Definitionally, this is in wind axes.
             """
             return self.F_w[1]
 
         @property
-        def D(self) -> float | np.ndarray:
+        def D(self) -> Vectorizable:
             """
             The drag force [N]. Definitionally, this is in wind axes.
             """
             return -self.F_w[0]
 
         @property
-        def l_b(self) -> float | np.ndarray:
+        def l_b(self) -> Vectorizable:
             """
             The rolling moment [Nm] in body axes. Positive is roll-right.
             """
             return self.M_b[0]
 
         @property
-        def m_b(self) -> float | np.ndarray:
+        def m_b(self) -> Vectorizable:
             """
             The pitching moment [Nm] in body axes. Positive is nose-up.
             """
             return self.M_b[1]
 
         @property
-        def n_b(self) -> float | np.ndarray:
+        def n_b(self) -> Vectorizable:
             """
             The yawing moment [Nm] in body axes. Positive is nose-right.
             """
diff --git a/aerosandbox/aerodynamics/aero_3D/lifting_line.py b/aerosandbox/aerodynamics/aero_3D/lifting_line.py
@@ -9,6 +9,7 @@
 )
 from typing import Callable
 from aerosandbox.aerodynamics.aero_3D.aero_buildup import AeroBuildup
+from aerosandbox.numpy.typing import Vectorizable
 from dataclasses import dataclass
 
 
@@ -157,78 +158,78 @@ def __repr__(self):
             )
 
         @property
-        def F_b(self) -> list[float | np.ndarray]:
+        def F_b(self) -> tuple[Vectorizable, Vectorizable, Vectorizable]:
             """
-            An [x, y, z] list of forces in body axes [N]
+            An (x, y, z) tuple of forces in body axes [N]
             """
             return self.op_point.convert_axes(
                 *self.F_g, from_axes="geometry", to_axes="body"
             )
 
         @property
-        def F_w(self) -> list[float | np.ndarray]:
+        def F_w(self) -> tuple[Vectorizable, Vectorizable, Vectorizable]:
             """
-            An [x, y, z] list of forces in wind axes [N]
+            An (x, y, z) tuple of forces in wind axes [N]
             """
             return self.op_point.convert_axes(
                 *self.F_g, from_axes="geometry", to_axes="wind"
             )
 
         @property
-        def M_b(self) -> list[float | np.ndarray]:
+        def M_b(self) -> tuple[Vectorizable, Vectorizable, Vectorizable]:
             """
-            An [x, y, z] list of moments about body axes [Nm]
+            An (x, y, z) tuple of moments about body axes [Nm]
             """
             return self.op_point.convert_axes(
                 *self.M_g, from_axes="geometry", to_axes="body"
             )
 
         @property
-        def M_w(self) -> list[float | np.ndarray]:
+        def M_w(self) -> tuple[Vectorizable, Vectorizable, Vectorizable]:
             """
-            An [x, y, z] list of moments about wind axes [Nm]
+            An (x, y, z) tuple of moments about wind axes [Nm]
             """
             return self.op_point.convert_axes(
                 *self.M_g, from_axes="geometry", to_axes="wind"
             )
 
         @property
-        def L(self) -> float | np.ndarray:
+        def L(self) -> Vectorizable:
             """
             The lift force [N]. Definitionally, this is in wind axes.
             """
             return -self.F_w[2]
 
         @property
-        def Y(self) -> float | np.ndarray:
+        def Y(self) -> Vectorizable:
             """
             The side force [N]. Definitionally, this is in wind axes.
             """
             return self.F_w[1]
 
         @property
-        def D(self) -> float | np.ndarray:
+        def D(self) -> Vectorizable:
             """
             The drag force [N]. Definitionally, this is in wind axes.
             """
             return -self.F_w[0]
 
         @property
-        def l_b(self) -> float | np.ndarray:
+        def l_b(self) -> Vectorizable:
             """
             The rolling moment [Nm] in body axes. Positive is roll-right.
             """
             return self.M_b[0]
 
         @property
-        def m_b(self) -> float | np.ndarray:
+        def m_b(self) -> Vectorizable:
             """
             The pitching moment [Nm] in body axes. Positive is nose-up.
             """
             return self.M_b[1]
 
         @property
-        def n_b(self) -> float | np.ndarray:
+        def n_b(self) -> Vectorizable:
             """
             The yawing moment [Nm] in body axes. Positive is nose-right.
             """
diff --git a/aerosandbox/aerodynamics/aero_3D/linear_potential_flow.py b/aerosandbox/aerodynamics/aero_3D/linear_potential_flow.py
@@ -283,7 +283,7 @@ def get_induced_velocity_at_points(
             points: np.ndarray,
             vortex_strengths: np.ndarray,
             sum_across_elements: bool = True,
-        ) -> tuple[np.ndarray]:
+        ) -> tuple[float | np.ndarray, float | np.ndarray, float | np.ndarray]:
             u_induced, v_induced, w_induced = calculate_induced_velocity_horseshoe(
                 x_field=tall(points[:, 0]),
                 y_field=tall(points[:, 1]),
diff --git a/aerosandbox/atmosphere/thermodynamics/gas.py b/aerosandbox/atmosphere/thermodynamics/gas.py
@@ -326,6 +326,7 @@ def process(
         if inplace:
             self.pressure = new_pressure
             self.temperature = new_temperature
+            return self
 
         else:
             return PerfectGas(
diff --git a/aerosandbox/dynamics/point_mass/common_point_mass.py b/aerosandbox/dynamics/point_mass/common_point_mass.py
@@ -111,8 +111,8 @@ def _set_state(
             )  # Pack the iterable into a dict-like, then do the same thing as above.
 
     def unpack_state(
-        self, dict_like_state: dict[str, float | np.ndarray] | None = None
-    ) -> tuple[float | np.ndarray]:
+        self, dict_like_state: dict[str, Vectorizable] | None = None
+    ) -> tuple[Vectorizable, ...]:
         """
         'Unpacks' a Dict-like state into an array-like that represents the state of the dynamical system.
 
diff --git a/aerosandbox/geometry/airfoil/airfoil_families.py b/aerosandbox/geometry/airfoil/airfoil_families.py
@@ -565,24 +565,17 @@ def get_coordinates_from_raw_dat(raw_text: list[str]) -> np.ndarray:
     Returns: A Nx2 ndarray of airfoil coordinates [x, y].
 
     """
-    raw_coordinates = []
-
-    def is_number(s: str) -> bool:
-        # Determines whether a string is representable as a float
-        try:
-            float(s)
-        except ValueError:
-            return False
-        return True
+    raw_coordinates: list[list[float]] = []
 
     def parse_line(line: str) -> list[float] | None:
-        # Given a single line of a `*.dat` file, tries to parse it into a list of two floats [x, y].
-        # If not possible, returns None.
+        """Parse a single line of a *.dat file into [x, y] floats, or None if not parseable."""
         line_split = re.split(r"[;|,|\s|\t]", line)
         line_items = [s for s in line_split if s != ""]
-        if len(line_items) == 2 and all([is_number(item) for item in line_items]):
-            return line_items
-        else:
+        if len(line_items) != 2:
+            return None
+        try:
+            return [float(line_items[0]), float(line_items[1])]
+        except ValueError:
             return None
 
     for line in raw_text:
@@ -593,12 +586,10 @@ def parse_line(line: str) -> list[float] | None:
     if len(raw_coordinates) == 0:
         raise ValueError("Could not read any coordinates from the `raw_text` input!")
 
-    coordinates = np.array(raw_coordinates, dtype=float)
-
-    return coordinates
+    return np.array(raw_coordinates)
 
 
-def get_file_coordinates(filepath: str | os.PathLike):
+def get_file_coordinates(filepath: str | os.PathLike) -> np.ndarray:
     possible_errors = (FileNotFoundError, UnicodeDecodeError)
 
     if isinstance(filepath, np.ndarray):
diff --git a/aerosandbox/library/regulations/far_part_23.py b/aerosandbox/library/regulations/far_part_23.py
@@ -1,4 +1,5 @@
 import aerosandbox.numpy as np
+from aerosandbox.numpy.typing import Vectorizable
 from aerosandbox.tools import units as u
 from typing import Literal
 
@@ -8,9 +9,9 @@
 
 
 def limit_load_factors(
-    design_mass_TOGW: float,
+    design_mass_TOGW: Vectorizable,
     category: Literal["normal", "utility", "acrobatic", "commuter"] = "normal",
-) -> tuple[float, float]:
+) -> tuple[Vectorizable, Vectorizable]:
     """
     Computes the required limit load factors for FAR Part 23 certification.
 
diff --git a/aerosandbox/numpy/calculus.py b/aerosandbox/numpy/calculus.py
@@ -150,7 +150,7 @@ def gradient(
             dx_shape[axis] = shape[axis] - 1
             dx = reshape(dx, dx_shape)
 
-            def get_slice(slice_obj: slice) -> tuple[slice]:
+            def get_slice(slice_obj: slice) -> tuple[slice, ...]:
                 slices = [slice(None)] * len(shape)
                 slices[axis] = slice_obj
                 return tuple(slices)
diff --git a/aerosandbox/numpy/finite_difference_operators.py b/aerosandbox/numpy/finite_difference_operators.py
@@ -53,10 +53,10 @@ def finite_difference_coefficients(
     """
     ### Check inputs
     if derivative_degree < 1:
-        return ValueError("The parameter derivative_degree must be an integer >= 1.")
+        raise ValueError("The parameter derivative_degree must be an integer >= 1.")
     expected_order_of_accuracy = length(x) - derivative_degree
     if expected_order_of_accuracy < 1:
-        return ValueError(
+        raise ValueError(
             "You need to provide at least (derivative_degree+1) grid points in the x vector."
         )
 
diff --git a/aerosandbox/structures/ignore/wing_structure_generator.py b/aerosandbox/structures/ignore/wing_structure_generator.py
@@ -40,7 +40,7 @@ def add_ribs_from_section_span_fractions(
         self,
         section_index: int,
         section_span_fractions: float | int | list[float] | np.ndarray,
-        rib_thickness: float = None,
+        rib_thickness: float | None = None,
     ):
         if rib_thickness is None:
             rib_thickness = self.default_rib_thickness
@@ -81,7 +81,7 @@ def add_ribs_from_section_span_fractions(
             )
 
     def add_ribs_from_xsecs(
-        self, indexes: list[int] = None, rib_thickness: float = None
+        self, indexes: list[int] | None = None, rib_thickness: float | None = None
     ):
         if rib_thickness is None:
             rib_thickness = self.default_rib_thickness
@@ -117,7 +117,7 @@ def add_ribs_from_xsecs(
     def add_ribs_from_span_fractions(
         self,
         span_fractions: float | list[float] | np.ndarray = np.linspace(0, 1, 10),
-        rib_thickness: float = None,
+        rib_thickness: float | None = None,
     ):
         ### Handle span_fractions if it's not an iterable
         try:
@@ -164,7 +164,7 @@ def add_tube_spar(
         y_over_c_location_root=None,
         x_over_c_location_tip=None,
         y_over_c_location_tip=None,
-        diameter_tip: float = None,
+        diameter_tip: float | None = None,
         cut_ribs: bool = True,
     ):
         if diameter_tip is None:
