executable-engineering 0.1.29__tar.gz → 0.1.32__tar.gz

{executable_engineering-0.1.29 → executable_engineering-0.1.32}/PKG-INFO

@@ -1,7 +1,7 @@
 Metadata-Version: 2.4
 Name: executable_engineering
-Version: 0.1.29
-Summary: Tools for ENGPHYS 3NM4
+Version: 0.1.32
+Summary: Demonstrative tools for the Executable Engineering textbook
 Home-page: https://github.com/themintlab/ExecutableEngineering/tree/main/executable_engineering
 Author: Michael Welland
 License: MIT
@@ -14,13 +14,15 @@ License-File: LICENSE
 Requires-Dist: numpy
 Requires-Dist: matplotlib
 Requires-Dist: networkx
-Requires-Dist: ipywidgets
-Requires-Dist: IPython
-Requires-Dist: pyppeteer
-Requires-Dist: nbconvert
 Requires-Dist: scipy
 Requires-Dist: plotly
 Requires-Dist: pandas
+Provides-Extra: docs
+Requires-Dist: pyppeteer; extra == "docs"
+Requires-Dist: nbconvert; extra == "docs"
+Provides-Extra: jupyter
+Requires-Dist: ipywidgets; extra == "jupyter"
+Requires-Dist: IPython; extra == "jupyter"
 Dynamic: author
 Dynamic: classifier
 Dynamic: description
@@ -28,6 +30,7 @@ Dynamic: description-content-type
 Dynamic: home-page
 Dynamic: license
 Dynamic: license-file
+Dynamic: provides-extra
 Dynamic: requires-dist
 Dynamic: requires-python
 Dynamic: summary

{executable_engineering-0.1.29 → executable_engineering-0.1.32}/executable_engineering/__init__.py

@@ -4,9 +4,11 @@ from .closedrootfinder_module import RootFinderClosed
 from .closedoptimization_module import OptimizerClosed
 from .openoptimization_module import OptimizerOpen
 from .gradientoptimization_module import OptimizerGrad
-
+from .numerical_error_module import python_internal_binary, decimal_to_binary
+from .linearsystems_module import visual_solve_2d
 
 __all__ = [
+  "visual_solve_2d",
   "init_weights",
   "tanh",
   "tanh_derivative",
@@ -29,4 +31,6 @@ __all__ = [
   "OptimizerClosed",
   "OptimizerOpen",
   "OptimizerGrad",
+  "python_internal_binary",
+  "decimal_to_binary",
 ]

executable_engineering-0.1.32/executable_engineering/linearsystems_module.py

@@ -0,0 +1,95 @@
+import numpy as np
+import plotly.graph_objects as go
+import plotly.colors as pcolors
+
+def visual_solve_2d(A, b, x_range=None, title="Linear System", labels=None):
+    """
+    Solves and visualizes a 2D linear system Ax = b.
+    Plots the lines represented by the rows of A and b, and the solution point.
+    Falls back to the pseudoinverse for singular or non-square systems.
+    
+    Args:
+        A: (N, 2) array-like of coefficients.
+        b: (N,) array-like of constants.
+        x_range: Tuple of (xmin, xmax) for plotting. If None, auto-calculated.
+        title: Plot title.
+        labels: Tuple of (x_label, y_label) for axes.
+    """
+    if labels is None:
+        labels = ("x", "y")
+        
+    A = np.array(A, dtype=float)
+    b = np.array(b, dtype=float).flatten()
+    
+    if A.shape[1] != 2:
+        raise ValueError("visual_solve_2d only supports 2D systems (A must have 2 columns).")
+    
+    # 1. Solve the system using the pseudoinverse
+    # For a square, full-rank matrix, pinv is mathematically identical to the standard inverse.
+    # For singular or non-square matrices, it finds the least-squares/minimum-norm solution.
+    sol = np.linalg.pinv(A) @ b
+    
+    # Check if the solution is exact AND unique.
+    # It is exact if A*x = b. It is unique if the rank of A equals the number of variables (columns).
+    is_exact = np.allclose(A @ sol, b)
+    is_unique = np.linalg.matrix_rank(A) == A.shape[1]
+    is_exact_unique = is_exact and is_unique
+        
+    # 2. Determine plotting range
+    if x_range is None:
+        center_x = sol[0]
+        x_range = (center_x - 10, center_x + 10)
+        
+    x_vals = np.linspace(x_range[0], x_range[1], 100)
+    
+    fig = go.Figure()
+    
+    colors = pcolors.qualitative.Plotly
+    
+    # 3. Plot each line
+    for i in range(A.shape[0]):
+        a1, a2 = A[i]
+        bi = b[i]
+        
+        if np.abs(a2) > 1e-10:
+            # y = (bi - a1*x) / a2
+            y_vals = (bi - a1 * x_vals) / a2
+            line_name = f"{a1:g}{labels[0]} + {a2:g}{labels[1]} = {bi:g}"
+            fig.add_trace(go.Scatter(x=x_vals, y=y_vals, mode='lines', name=line_name, line=dict(color=colors[i % len(colors)])))
+        else:
+            # Vertical line: x = bi / a1
+            if np.abs(a1) > 1e-10:
+                x_val = bi / a1
+                line_name = f"{a1:g}{labels[0]} = {bi:g}"
+                fig.add_trace(go.Scatter(x=[x_val, x_val], y=[-100, 100], mode='lines', name=line_name, line=dict(color=colors[i % len(colors)])))
+                
+    # 4. Plot solution point
+    if is_unique:
+        if is_exact:
+            sol_label = "Exact Solution"
+        else:
+            sol_label = "Best Fit Solution"
+            
+        hover_text = f"{labels[0]}={sol[0]:.4g}, {labels[1]}={sol[1]:.4g}"
+        fig.add_trace(go.Scatter(
+            x=[sol[0]], 
+            y=[sol[1]], 
+            mode='markers', 
+            marker=dict(color='black', size=12, symbol='star'),
+            name=sol_label,
+            hoverinfo='text',
+            hovertext=hover_text
+        ))
+    
+    # Auto-adjust y-range nicely around the solution if not a vertical-only mess
+    y_margin = (x_range[1] - x_range[0]) / 2.0
+    
+    fig.update_layout(
+        title=title,
+        xaxis_title=labels[0],
+        yaxis_title=labels[1],
+        yaxis=dict(range=[sol[1] - y_margin, sol[1] + y_margin]),
+        template="plotly_white"
+    )
+    
+    return fig

executable_engineering-0.1.32/executable_engineering/numerical_error_module.py

@@ -0,0 +1,40 @@
+import struct
+
+def python_internal_binary(num):
+    # Pack the float into 8 bytes using IEEE 754 double-precision ('d')
+    # The '!' ensures network byte order (big-endian) so bits read left-to-right
+    packed_bytes = struct.pack('!d', float(num))
+    
+    # Unpack those bytes as a single 64-bit unsigned integer ('Q')
+    packed_int, = struct.unpack('!Q', packed_bytes)
+    
+    # Format the integer as a 64-bit binary string, padded with leading zeros
+    binary_str = f"{packed_int:064b}"
+    
+    # Split into the IEEE 754 components for readability
+    sign = binary_str[0]         # 1 bit
+    exponent = binary_str[1:12]  # 11 bits
+    fraction = binary_str[12:]   # 52 bits
+    
+    return f"{sign} {exponent} {fraction}"
+
+def decimal_to_binary(num, precision=20):
+    sign = "-" if num < 0 else ""
+    num = abs(num)
+    int_part = int(num)
+    frac = num - int_part
+    int_str = bin(int_part)[2:]
+    
+    frac_part = ""
+    for _ in range(precision):
+        if frac == 0: 
+            break
+        frac *= 2
+        bit = int(frac)
+        frac_part += str(bit)
+        frac -= bit
+        
+    if not frac_part:
+        frac_part = "0"
+        
+    return f"{sign}{int_str}.{frac_part}"

{executable_engineering-0.1.29 → executable_engineering-0.1.32}/executable_engineering.egg-info/PKG-INFO

@@ -1,7 +1,7 @@
 Metadata-Version: 2.4
 Name: executable_engineering
-Version: 0.1.29
-Summary: Tools for ENGPHYS 3NM4
+Version: 0.1.32
+Summary: Demonstrative tools for the Executable Engineering textbook
 Home-page: https://github.com/themintlab/ExecutableEngineering/tree/main/executable_engineering
 Author: Michael Welland
 License: MIT
@@ -14,13 +14,15 @@ License-File: LICENSE
 Requires-Dist: numpy
 Requires-Dist: matplotlib
 Requires-Dist: networkx
-Requires-Dist: ipywidgets
-Requires-Dist: IPython
-Requires-Dist: pyppeteer
-Requires-Dist: nbconvert
 Requires-Dist: scipy
 Requires-Dist: plotly
 Requires-Dist: pandas
+Provides-Extra: docs
+Requires-Dist: pyppeteer; extra == "docs"
+Requires-Dist: nbconvert; extra == "docs"
+Provides-Extra: jupyter
+Requires-Dist: ipywidgets; extra == "jupyter"
+Requires-Dist: IPython; extra == "jupyter"
 Dynamic: author
 Dynamic: classifier
 Dynamic: description
@@ -28,6 +30,7 @@ Dynamic: description-content-type
 Dynamic: home-page
 Dynamic: license
 Dynamic: license-file
+Dynamic: provides-extra
 Dynamic: requires-dist
 Dynamic: requires-python
 Dynamic: summary

{executable_engineering-0.1.29 → executable_engineering-0.1.32}/executable_engineering.egg-info/SOURCES.txt

@@ -6,7 +6,9 @@ executable_engineering/__init__.py
 executable_engineering/closedoptimization_module.py
 executable_engineering/closedrootfinder_module.py
 executable_engineering/gradientoptimization_module.py
+executable_engineering/linearsystems_module.py
 executable_engineering/neuralnet_module.py
+executable_engineering/numerical_error_module.py
 executable_engineering/openoptimization_module.py
 executable_engineering/openrootfinder_module.py
 executable_engineering.egg-info/PKG-INFO

{executable_engineering-0.1.29 → executable_engineering-0.1.32}/executable_engineering.egg-info/requires.txt

@@ -1,10 +1,14 @@
 numpy
 matplotlib
 networkx
-ipywidgets
-IPython
-pyppeteer
-nbconvert
 scipy
 plotly
 pandas
+
+[docs]
+pyppeteer
+nbconvert
+
+[jupyter]
+ipywidgets
+IPython

{executable_engineering-0.1.29 → executable_engineering-0.1.32}/setup.py

@@ -6,8 +6,8 @@ long_description = (this_directory / "README.md").read_text(encoding="utf-8")
 
 setup(
     name="executable_engineering",
-    version="0.1.29",
-    description="Tools for ENGPHYS 3NM4",
+    version="0.1.32",
+    description="Demonstrative tools for the Executable Engineering textbook",
     long_description=long_description,
     long_description_content_type="text/markdown",
     url="https://github.com/themintlab/ExecutableEngineering/tree/main/executable_engineering",
@@ -18,14 +18,14 @@ setup(
         "numpy",
         "matplotlib",
         "networkx",
-        "ipywidgets",
-        "IPython",
-        "pyppeteer",
-        "nbconvert",
         "scipy",
         "plotly",
         "pandas",
     ],
+    extras_require={
+        "docs": ["pyppeteer", "nbconvert"],
+        "jupyter": ["ipywidgets", "IPython"],
+    },
     python_requires=">=3.9",
     license="MIT",
     classifiers=[