qnty 0.0.7__py3-none-any.whl → 0.0.9__py3-none-any.whl

qnty-0.0.9.dist-info/METADATA

@@ -0,0 +1,199 @@
+Metadata-Version: 2.3
+Name: qnty
+Version: 0.0.9
+Summary: High-performance unit system library for Python with dimensional safety and fast unit conversions
+License: Apache-2.0
+Keywords: units,dimensional analysis,engineering,physics,quantities,measurements
+Author: tn3wman
+Requires-Python: >=3.11, <3.14
+Classifier: Development Status :: 4 - Beta
+Classifier: Intended Audience :: Developers
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: Apache Software License
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Programming Language :: Python :: 3.13
+Classifier: Topic :: Scientific/Engineering
+Classifier: Topic :: Scientific/Engineering :: Physics
+Classifier: Topic :: Software Development :: Libraries :: Python Modules
+Provides-Extra: benchmark
+Provides-Extra: dev
+Requires-Dist: Pint (>=0.24.4) ; extra == "benchmark"
+Requires-Dist: numpy (>=2.3.2)
+Requires-Dist: pytest (>=8.4.1) ; extra == "dev"
+Requires-Dist: ruff (>=0.1.0) ; extra == "dev"
+Project-URL: Bug Tracker, https://github.com/tn3wman/qnty/issues
+Project-URL: Documentation, https://github.com/tn3wman/qnty#readme
+Project-URL: Homepage, https://github.com/tn3wman/qnty
+Project-URL: Repository, https://github.com/tn3wman/qnty
+Description-Content-Type: text/markdown
+
+# Qnty
+
+**High-performance unit system library for Python with dimensional safety and fast unit conversions for engineering calculations.**
+
+[![Python 3.11+](https://img.shields.io/badge/python-3.11+-blue.svg)](https://www.python.org/downloads/)
+[![License: Apache-2.0](https://img.shields.io/badge/License-Apache%202.0-blue.svg)](https://opensource.org/licenses/Apache-2.0)
+[![Development Status](https://img.shields.io/badge/status-beta-orange.svg)](https://pypi.org/project/qnty/)
+
+## ⚠️ Important Disclaimer
+
+**🚧 Work in Progress**: Qnty is currently in active development and has not been thoroughly vetted for production engineering calculations. While we strive for accuracy, this library should not be used for critical engineering applications without independent verification.
+
+**📐 Accuracy Notice**: The authors are not responsible or liable for incorrect results, calculation errors, or any consequences arising from the use of this library. Always validate calculations independently using established engineering tools and practices.
+
+*Use Qnty to help prevent unit errors, but always verify critical calculations through multiple methods.*
+
+---
+
+## ✨ Key Features
+
+- **🚀 Ultra-Fast Performance**: Prime number encoding and pre-computed conversion tables
+- **🛡️ Type Safety**: Compile-time dimensional analysis prevents unit errors
+- **⚡ Zero-Cost Abstractions**: Optimized operations with `__slots__` and caching
+- **🔗 Fluent API**: Intuitive method chaining for readable code
+- **🧮 Engineering-Focused**: Built for real-world engineering calculations
+- **🧬 Mathematical System**: Built-in equation solving and expression trees
+- **📊 Comprehensive Testing**: 187 tests with performance benchmarks
+
+## 🚀 Quick Start
+
+### Installation
+
+```bash
+pip install qnty
+```
+
+### Basic Usage
+
+```python
+from qnty import Length, Pressure, Area
+
+# Create quantities with dimensional safety
+width = Length(3, "meter", "Width")
+height = Length(2, "meter", "Height")
+
+# Solve mathematical expressions
+area = Area("area", is_known=False)
+area.solve_from(width * height)
+print(f"Area: {area}")  # Area: 6.0 m²
+```
+
+### Engineering Example
+
+```python
+from qnty import Problem, Length, Pressure
+
+class PipeThickness(Problem):
+    """Calculate pipe wall thickness"""
+    
+    # Known parameters
+    pressure = Pressure(150, "pound_force_per_square_inch", "Internal Pressure")
+    diameter = Length(6, "inch", "Pipe Diameter")
+    allowable_stress = Pressure(20000, "pound_force_per_square_inch", "Allowable Stress")
+    
+    # Unknown to solve for
+    thickness = Length("thickness", is_known=False)
+    
+    # Engineering equation: t = (P × D) / (2 × S)
+    equation = thickness.equals((pressure * diameter) / (2 * allowable_stress))
+
+# Solve the problem
+problem = PipeThickness()
+problem.solve()
+print(f"Required thickness: {problem.thickness}")
+```
+
+### Mathematical Operations
+
+```python
+from qnty import Length, sqrt, Area
+
+# Dimensional analysis with mathematical functions
+area = Area(25, "square_meter", "Square Area")
+side = Length("side", is_known=False)
+side.solve_from(sqrt(area))  # Returns Length, not Area!
+print(f"Side length: {side}")  # Side length: 5.0 m
+```
+
+## 📚 Documentation
+
+- **[📖 Tutorial](docs/TUTORIAL.md)** - Step-by-step learning guide
+- **[📋 API Reference](docs/API_REFERENCE.md)** - Complete API documentation
+- **[🏗️ Examples](examples/)** - Real-world engineering examples
+- **[📁 Full Documentation](docs/)** - Complete documentation index
+
+## 🚀 Performance
+
+Qnty significantly outperforms other unit libraries with **18.9x average speedup** over Pint:
+
+| Operation | Qnty | Pint | **Speedup** |
+|-----------|------|------|-------------|
+| Mixed Unit Addition | 0.76 μs | 17.52 μs | **23.1x** |
+| Complex ASME Equation | 4.07 μs | 106.17 μs | **26.1x** 🚀 |
+| Type-Safe Variables | 0.98 μs | 9.65 μs | **9.8x** |
+| **AVERAGE** | **1.89 μs** | **35.83 μs** | **18.9x** 🏆 |
+
+*Run `pytest tests/test_benchmark.py -v -s` to verify on your system.*
+
+## 🧮 100+ Engineering Quantities
+
+Qnty provides comprehensive coverage of engineering domains:
+
+```python
+from qnty import (
+    # Mechanical
+    Length, Area, Volume, Mass, Force, Pressure, Temperature,
+    # Electrical  
+    ElectricPotential, ElectricCurrentIntensity, ElectricResistance,
+    # Thermal
+    ThermalConductivity, HeatTransferCoefficient,
+    # Fluid Dynamics
+    ViscosityDynamic, MassFlowRate, VolumetricFlowRate,
+    # And 80+ more...
+)
+```
+
+## 🔧 Development
+
+```bash
+# Install dependencies
+pip install -r requirements.txt
+
+# Run tests
+pytest
+
+# Run specific test
+pytest tests/test_dimension.py -v
+
+# Run benchmarks
+python tests/test_benchmark.py
+
+# Lint code
+ruff check src/ tests/
+ruff format src/ tests/
+```
+
+## 📄 License
+
+This project is licensed under the Apache License 2.0 - see the [LICENSE](LICENSE) file for details.
+
+## 🤝 Contributing
+
+We welcome contributions! Please see [CLAUDE.md](CLAUDE.md) for development guidelines and:
+
+1. Fork the repository
+2. Create a feature branch
+3. Add tests for new functionality
+4. Ensure all tests pass: `pytest`
+5. Submit a pull request
+
+---
+
+**Ready to supercharge your engineering calculations?** 🚀
+
+- Start with the **[Tutorial](docs/TUTORIAL.md)**
+- Browse the **[API Reference](docs/API_REFERENCE.md)**
+- Try the **[Examples](examples/)**
+

qnty-0.0.9.dist-info/RECORD

@@ -0,0 +1,63 @@
+qnty/__init__.py,sha256=nixvrkyF4sIwtLCPVIfADSEKkQoKlFZs_7TbpuBFWXE,5563
+qnty/_backup/problem_original.py,sha256=M58iUTiWBO1sS3rM3aQSlpIwxbSRiW2LICYK-Lnoz90,56698
+qnty/_backup/quantity.py,sha256=oD-TNDkxK1Tfa6S_sxNnDWyv1_SuW0TrNKz01e45fxg,1509
+qnty/codegen/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+qnty/codegen/cli.py,sha256=C66IQ-je69UP0llVy7kgSrfUYM3PHuantcefYvFBDHw,4092
+qnty/codegen/generators/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+qnty/codegen/generators/data/unit_data.json,sha256=KbJHRv5U2bF_uS_Az9G-uFSjwvxCcvmn1ZCxR5H0uGI,253446
+qnty/codegen/generators/data_processor.py,sha256=3Cibqkegi48KuaD9I9i1QinP2h7Fi76EOJ2hwmZq6ug,13741
+qnty/codegen/generators/dimensions_gen.py,sha256=03Ryq8CPIU3gMSOZNnOjc2JMR3Vi-lgVmBYFnOliL8w,16379
+qnty/codegen/generators/doc_generator.py,sha256=zmhDWwxODII5d-_FlKskhMTHHBmS8NRZOzOyPzPMyTs,6184
+qnty/codegen/generators/out/dimension_mapping.json,sha256=_YRPCzMvT9RW_YhURbuHH8gKhYtdu8gaYSqJ-bZxSAE,23864
+qnty/codegen/generators/out/dimension_metadata.json,sha256=-2uXcB8grYMZew9A6mDSajrJ9yCoegrTx23g70QTK_Y,2929
+qnty/codegen/generators/out/units_metadata.json,sha256=5MMDsT6Y99tthus8YT5Xu6iaG4cn99sP-JomKYZ084k,5890
+qnty/codegen/generators/quantities_gen.py,sha256=9Nf_3PUpTf0A1yLe8L4CodrpsD4Oz7DTXX3E-sssj3I,6271
+qnty/codegen/generators/setters_gen.py,sha256=u76OwdH7Sl6XsIDGZn6ID8Xp0VEptBXO1QJ4h6sne1c,7145
+qnty/codegen/generators/stubs_gen.py,sha256=F6ZEk48W9C8nRKP6_pKoGPO3gnDtGfvn-42XgrhA4uk,6659
+qnty/codegen/generators/units_gen.py,sha256=iESIFLCVNWi_4Ftss2_s2QTGTOTutInAeltLQdMjrI8,10592
+qnty/codegen/generators/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+qnty/equations/__init__.py,sha256=ngficq_uwosz7O2BaqDH-G6dqp1burjfkcwE4MnGyIo,108
+qnty/equations/equation.py,sha256=dv527D-fE_0jLmK6BJI5SNCjIgzql1JpBcxZkej06-A,10848
+qnty/equations/system.py,sha256=mgslAh_YgEU1V_y48Yy1_eHIFO37vDsSG4Gb3zeuPxc,5115
+qnty/expressions/__init__.py,sha256=qo5I-Z3O1hP9Kk1Jgf5r7x4tAzMLS-a96rY5cqRNVKc,917
+qnty/expressions/cache.py,sha256=ahJ-AyxJJQyGAVQWSEcQAng1TdQ2qysBPiujk9hhink,3413
+qnty/expressions/functions.py,sha256=X7Cfchmc4qY3JE3NnacI5DomYE0jPbAAQHOsrxG3vqk,3618
+qnty/expressions/nodes.py,sha256=bTmpJ36sjnlxgKb2npJfh56K1arN-FOxoiI30Gx31-g,23394
+qnty/generated/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+qnty/generated/dimensions.py,sha256=WHVr28Ou2DB2QuZCDcxLY2AcauqdE4QGrNFNUJ2EDXY,20661
+qnty/generated/quantities.py,sha256=1FC1P_hIcXotUgDiJXy28P9OCq2jdWEmO5aVZJepNpY,245345
+qnty/generated/quantities.pyi,sha256=URYr8K2JVeXlwbSnUvOLoLYS9beVhGrL33A7Y1hs1Hs,151386
+qnty/generated/setters.py,sha256=Go62i_TlfcBztI8rI7h2FMWGmtoVPGc3JTx6ZmJMrWI,446272
+qnty/generated/units.py,sha256=WXkSL_G4hKzbTiGmN3jreWCuUyutyGb50tRwfT0Xje4,285245
+qnty/problem/__init__.py,sha256=9bML5bjM-uPAHet8n98_OwQLivEs2mzaHjZhENbA7sE,3140
+qnty/problem/base.py,sha256=jepN3vmMyh9eXptJrAQfPICHbDPLQsD2g9k4EY7IqrA,5680
+qnty/problem/composition.py,sha256=DUdJ7oMl1d4ETqgR4rDflqfBOdE2ykr-f6qONOErfZQ,13894
+qnty/problem/composition_mixin.py,sha256=hH9bEZyWc5vz5LfLb7GM2K6WDP7SUgHWUu_Vx4BtEM0,17603
+qnty/problem/equations.py,sha256=QefhpWBQt9OD0vUeYRZs54iv7kQ9xkO0XBtYBtOjfPE,19334
+qnty/problem/metaclass.py,sha256=00T4fg3ACeTo7w1njRgZyOGP2jZm2P2_yxO0rZlnwhE,11944
+qnty/problem/reconstruction.py,sha256=hE7L5Eh2OH3nR_3PFKsSiNRHDSbgSumrhXxPRESmJXI,40042
+qnty/problem/solving.py,sha256=DJ_Dj63vn2ykdFUp77QNxEGOWmPrAzpLe7A_mnCkMzs,6657
+qnty/problem/validation.py,sha256=1Kbebn4JvYO0qw2VOhfeyzzkuSmnTpOANpvVmsXcNuU,2157
+qnty/problem/variables.py,sha256=Fmhc-GmW-tHMQk1YVtHgOQlfvRJxH1le94Nc0HmeyBY,9911
+qnty/quantities/__init__.py,sha256=C8Z032FGpMVt9-PvtIQLNZkxmXjbHUm2OHoytebUL1g,107
+qnty/quantities/expression_quantity.py,sha256=khJ5izSaj497-rd6Ci1p8XyBidwAoXFlqpBiOqVbi3A,13278
+qnty/quantities/quantity.py,sha256=PK1vGMZv55fy87ysvakVipO0lWK_XySKZCobJfKMjTQ,18855
+qnty/quantities/typed_quantity.py,sha256=-PL_6qhJvOLdGXoxJjCHiCx3ZB5Oao9l-e2zrON3cug,9681
+qnty/solving/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+qnty/solving/manager.py,sha256=FfG3DOrwAKeR0vHll46RtkLB9ZuWN43mthYHnr5AW98,3594
+qnty/solving/order.py,sha256=Z7I_7QonCaGC47DthC2O0m77On34GcaD7nn2oLIFmqs,14411
+qnty/solving/solvers/__init__.py,sha256=CBI2PGLhnQ3AVsfJPfBB7b0ff_BlPCRamGUIHjEMStM,501
+qnty/solving/solvers/base.py,sha256=KYACgawfo8GoRaXEduZRQjJeo948QNFU76N2VrCwcA4,3041
+qnty/solving/solvers/iterative.py,sha256=5SxIGMdP1HkAKsT6cXZR7Lxv3LIgM7dXY3PsTn0w_d4,8139
+qnty/solving/solvers/simultaneous.py,sha256=pqYUwh067z26kP49Nxq8HcCxsxB0dirPaVHMWG4QMTE,23225
+qnty/units/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+qnty/units/prefixes.py,sha256=U-dGj5ZdseErwFIzvz_RC3fID0D-DazE6ozD6vn4DSs,7442
+qnty/units/registry.py,sha256=rnKnEaG8rZg56ejK1nGl9xpe9hm1WmWAtD8tJ0A0MvA,7524
+qnty/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+qnty/utils/logging.py,sha256=QSYgZg6w1n9L-lIj5c3Ufee84E8VvRZOUgoW84vTUbA,1159
+qnty/validation/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+qnty/validation/registry.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+qnty/validation/rules.py,sha256=4X3pR6bl9YA5skT2vRmJFB78FR588Rz_cSTTJXD0BDI,6285
+qnty-0.0.9.dist-info/METADATA,sha256=AKPCK06zvVod4ss5nmYHmCgActFnehcUGRF8b7cIfLo,6761
+qnty-0.0.9.dist-info/WHEEL,sha256=b4K_helf-jlQoXBBETfwnf4B04YC67LOev0jo4fX5m8,88
+qnty-0.0.9.dist-info/RECORD,,

qnty/dimension.py

@@ -1,186 +0,0 @@
-"""
-Dimension System
-================
-
-Compile-time dimensional analysis using type system for ultra-fast operations.
-"""
-
-from dataclasses import dataclass
-from enum import IntEnum
-from typing import final
-
-
-class BaseDimension(IntEnum):
-    """Base dimensions as prime numbers for efficient bit operations."""
-    LENGTH = 2
-    MASS = 3
-    TIME = 5
-    CURRENT = 7
-    TEMPERATURE = 11
-    AMOUNT = 13
-    LUMINOSITY = 17
-    DIMENSIONLESS = 1  # Must be 1 to act as multiplicative identity
-
-
-@final
-@dataclass(frozen=True)
-class DimensionSignature:
-    """Immutable dimension signature for zero-cost dimensional analysis."""
-    
-    # Store as bit pattern for ultra-fast comparison
-    _signature: int | float = 1
-    
-    @classmethod
-    def create(cls, length=0, mass=0, time=0, current=0, temp=0, amount=0, luminosity=0):
-        """Create dimension from exponents."""
-        signature = 1
-        if length != 0:
-            signature *= BaseDimension.LENGTH ** length
-        if mass != 0:
-            signature *= BaseDimension.MASS ** mass
-        if time != 0:
-            signature *= BaseDimension.TIME ** time
-        if current != 0:
-            signature *= BaseDimension.CURRENT ** current
-        if temp != 0:
-            signature *= BaseDimension.TEMPERATURE ** temp
-        if amount != 0:
-            signature *= BaseDimension.AMOUNT ** amount
-        if luminosity != 0:
-            signature *= BaseDimension.LUMINOSITY ** luminosity
-        
-        return cls(signature)
-    
-    def __mul__(self, other):
-        return DimensionSignature(self._signature * other._signature)
-    
-    def __truediv__(self, other):
-        return DimensionSignature(self._signature / other._signature)
-    
-    def __pow__(self, power):
-        return DimensionSignature(self._signature ** power)
-    
-    def is_compatible(self, other):
-        """Ultra-fast dimensional compatibility check."""
-        return self._signature == other._signature
-    
-    def __eq__(self, other):
-        """Fast equality check for dimensions."""
-        return isinstance(other, DimensionSignature) and self._signature == other._signature
-    
-    def __hash__(self):
-        """Enable dimensions as dictionary keys."""
-        return hash(self._signature)
-
-
-# Pre-defined dimension constants (alphabetically ordered)
-ABSORBED_DOSE = DimensionSignature.create(length=2, time=-2)  # L^2 T^-2
-ACCELERATION = DimensionSignature.create(length=1, time=-2)  # L T^-2
-ACTIVATION_ENERGY = DimensionSignature.create(amount=-1, length=2, time=-2)  # N^-1 L^2 T^-2
-AMOUNT = DimensionSignature.create(amount=1)  # N
-AMOUNT_OF_SUBSTANCE = DimensionSignature.create(amount=1)  # N
-ANGLE_PLANE = DimensionSignature(BaseDimension.DIMENSIONLESS)  # Dimensionless
-ANGLE_SOLID = DimensionSignature(BaseDimension.DIMENSIONLESS)  # Dimensionless
-ANGULAR_ACCELERATION = DimensionSignature.create(time=-2)  # T^-2
-ANGULAR_MOMENTUM = DimensionSignature.create(length=2, mass=1, time=-1)  # L^2 M T^-1
-AREA = DimensionSignature.create(length=2)  # L^2
-AREA_PER_UNIT_VOLUME = DimensionSignature.create(length=-1)  # L^-1
-ATOMIC_WEIGHT = DimensionSignature.create(amount=-1, mass=1)  # N^-1 M
-CONCENTRATION = DimensionSignature.create(length=-3, mass=1)  # L^-3 M
-CURRENT = DimensionSignature.create(current=1)  # A
-DIMENSIONLESS = DimensionSignature(BaseDimension.DIMENSIONLESS)  # Dimensionless
-DYNAMIC_FLUIDITY = DimensionSignature.create(length=1, mass=-1, time=1)  # L M^-1 T
-ELECTRICAL_CONDUCTANCE = DimensionSignature.create(current=2, length=-2, mass=-1, time=3)  # A^2 L^-2 M^-1 T^3
-ELECTRICAL_PERMITTIVITY = DimensionSignature.create(current=2, length=-3, mass=-1, time=4)  # A^2 L^-3 M^-1 T^4
-ELECTRICAL_RESISTIVITY = DimensionSignature.create(current=-2, length=3, mass=1, time=-3)  # A^-2 L^3 M T^-3
-ELECTRIC_CAPACITANCE = DimensionSignature.create(current=2, length=-2, mass=-1, time=4)  # A^2 L^-2 M^-1 T^4
-ELECTRIC_CHARGE = DimensionSignature.create(amount=-1, current=1, time=1)  # N^-1 A T
-ELECTRIC_CURRENT_INTENSITY = DimensionSignature.create(current=1)  # A
-ELECTRIC_DIPOLE_MOMENT = DimensionSignature.create(current=1, length=1, time=1)  # A L T
-ELECTRIC_FIELD_STRENGTH = DimensionSignature.create(current=-1, length=1, mass=1, time=-3)  # A^-1 L M T^-3
-ELECTRIC_INDUCTANCE = DimensionSignature.create(current=-2, length=2, mass=1, time=-2)  # A^-2 L^2 M T^-2
-ELECTRIC_POTENTIAL = DimensionSignature.create(current=-1, length=2, mass=1, time=-3)  # A^-1 L^2 M T^-3
-ELECTRIC_RESISTANCE = DimensionSignature.create(current=-2, length=2, mass=1, time=-3)  # A^-2 L^2 M T^-3
-ENERGY_FLUX = DimensionSignature.create(mass=1, time=-3)  # M T^-3
-ENERGY_HEAT_WORK = DimensionSignature.create(length=2, mass=1, time=-2)  # L^2 M T^-2
-ENERGY_PER_UNIT_AREA = DimensionSignature.create(mass=1, time=-2)  # M T^-2
-FORCE = DimensionSignature.create(length=1, mass=1, time=-2)  # L M T^-2
-FORCE_BODY = DimensionSignature.create(length=-2, mass=1, time=-2)  # L^-2 M T^-2
-FORCE_PER_UNIT_MASS = DimensionSignature.create(length=1, time=-2)  # L T^-2
-FREQUENCY_VOLTAGE_RATIO = DimensionSignature.create(current=1, length=-2, mass=-1, time=3)  # A L^-2 M^-1 T^3
-FUEL_CONSUMPTION = DimensionSignature.create(length=-2)  # L^-2
-HEAT_OF_COMBUSTION = DimensionSignature.create(length=2, time=-2)  # L^2 T^-2
-HEAT_OF_FUSION = DimensionSignature.create(length=2, time=-2)  # L^2 T^-2
-HEAT_OF_VAPORIZATION = DimensionSignature.create(length=2, time=-2)  # L^2 T^-2
-HEAT_TRANSFER_COEFFICIENT = DimensionSignature.create(mass=1, temp=-1, time=-3)  # M Θ^-1 T^-3
-ILLUMINANCE = DimensionSignature.create(length=-2, luminosity=1)  # L^-2 J
-KINETIC_ENERGY_OF_TURBULENCE = DimensionSignature.create(length=2, time=-2)  # L^2 T^-2
-LENGTH = DimensionSignature.create(length=1)  # L
-LINEAR_MASS_DENSITY = DimensionSignature.create(length=-1, mass=1)  # L^-1 M
-LINEAR_MOMENTUM = DimensionSignature.create(length=1, mass=1, time=-1)  # L M T^-1
-LUMINANCE_SELF = DimensionSignature.create(length=-2, luminosity=1)  # L^-2 J
-LUMINOSITY = DimensionSignature.create(luminosity=1)  # J
-LUMINOUS_FLUX = DimensionSignature.create(luminosity=1)  # J
-LUMINOUS_INTENSITY = DimensionSignature.create(luminosity=1)  # J
-MAGNETIC_FIELD = DimensionSignature.create(current=1, length=-1)  # A L^-1
-MAGNETIC_FLUX = DimensionSignature.create(current=-1, length=2, mass=1, time=-2)  # A^-1 L^2 M T^-2
-MAGNETIC_INDUCTION_FIELD_STRENGTH = DimensionSignature.create(current=-1, mass=1, time=-2)  # A^-1 M T^-2
-MAGNETIC_MOMENT = DimensionSignature.create(current=1, length=2)  # A L^2
-MAGNETIC_PERMEABILITY = DimensionSignature.create(current=-2, length=2, mass=1, time=-2)  # A^-2 L^2 M T^-2
-MAGNETOMOTIVE_FORCE = DimensionSignature.create(current=1)  # A
-MASS = DimensionSignature.create(mass=1)  # M
-MASS_DENSITY = DimensionSignature.create(length=-3, mass=1)  # L^-3 M
-MASS_FLOW_RATE = DimensionSignature.create(mass=1, time=-1)  # M T^-1
-MASS_FLUX = DimensionSignature.create(length=-2, mass=1, time=-1)  # L^-2 M T^-1
-MASS_FRACTION_OF_I = DimensionSignature(BaseDimension.DIMENSIONLESS)  # Dimensionless
-MASS_TRANSFER_COEFFICIENT = DimensionSignature.create(length=-2, mass=1, time=-1)  # L^-2 M T^-1
-MOLALITY_OF_SOLUTE_I = DimensionSignature.create(amount=1, mass=-1)  # N M^-1
-MOLARITY_OF_I = DimensionSignature.create(amount=1, length=-3)  # N L^-3
-MOLAR_CONCENTRATION_BY_MASS = DimensionSignature.create(amount=1)  # N
-MOLAR_FLOW_RATE = DimensionSignature.create(amount=1, time=-1)  # N T^-1
-MOLAR_FLUX = DimensionSignature.create(amount=1, length=-2, time=-1)  # N L^-2 T^-1
-MOLAR_HEAT_CAPACITY = DimensionSignature.create(amount=-1, length=2, temp=-1, time=-2)  # N^-1 L^2 Θ^-1 T^-2
-MOLE_FRACTION_OF_I = DimensionSignature(BaseDimension.DIMENSIONLESS)  # Dimensionless
-MOMENTUM_FLOW_RATE = DimensionSignature.create(length=1, mass=1, time=-2)  # L M T^-2
-MOMENTUM_FLUX = DimensionSignature.create(length=-1, mass=1, time=-2)  # L^-1 M T^-2
-MOMENT_OF_INERTIA = DimensionSignature.create(length=2, mass=1)  # L^2 M
-NORMALITY_OF_SOLUTION = DimensionSignature.create(amount=1, length=-3)  # N L^-3
-PARTICLE_DENSITY = DimensionSignature.create(length=-3)  # L^-3
-PERCENT = DimensionSignature(BaseDimension.DIMENSIONLESS)  # Dimensionless
-PERMEABILITY = DimensionSignature.create(length=2)  # L^2
-PHOTON_EMISSION_RATE = DimensionSignature.create(length=-2, time=-1)  # L^-2 T^-1
-POWER_PER_UNIT_MASS = DimensionSignature.create(length=2, time=-3)  # L^2 T^-3
-POWER_PER_UNIT_VOLUME = DimensionSignature.create(length=-1, mass=1, time=-3)  # L^-1 M T^-3
-POWER_THERMAL_DUTY = DimensionSignature.create(length=2, mass=1, time=-3)  # L^2 M T^-3
-PRESSURE = DimensionSignature.create(length=-1, mass=1, time=-2)  # L^-1 M T^-2
-RADIATION_DOSE_EQUIVALENT = DimensionSignature.create(length=2, time=-2)  # L^2 T^-2
-RADIATION_EXPOSURE = DimensionSignature.create(current=1, mass=-1, time=1)  # A M^-1 T
-RADIOACTIVITY = DimensionSignature.create(time=-1)  # T^-1
-SECOND_MOMENT_OF_AREA = DimensionSignature.create(length=4)  # L^4
-SECOND_RADIATION_CONSTANT_PLANCK = DimensionSignature.create(length=1, temp=1)  # L Θ
-SPECIFIC_ENTHALPY = DimensionSignature.create(length=2, time=-2)  # L^2 T^-2
-SPECIFIC_GRAVITY = DimensionSignature(BaseDimension.DIMENSIONLESS)  # Dimensionless
-SPECIFIC_HEAT_CAPACITY_CONSTANT_PRESSURE = DimensionSignature.create(length=2, mass=1, temp=-1, time=-2)  # L^2 M Θ^-1 T^-2
-SPECIFIC_LENGTH = DimensionSignature.create(length=1, mass=-1)  # L M^-1
-SPECIFIC_SURFACE = DimensionSignature.create(length=2, mass=-1)  # L^2 M^-1
-SPECIFIC_VOLUME = DimensionSignature.create(length=3, mass=-1)  # L^3 M^-1
-STRESS = DimensionSignature.create(length=-1, mass=1, time=-2)  # L^-1 M T^-2
-SURFACE_MASS_DENSITY = DimensionSignature.create(length=-2, mass=1)  # L^-2 M
-SURFACE_TENSION = DimensionSignature.create(mass=1, time=-2)  # M T^-2
-TEMPERATURE = DimensionSignature.create(temp=1)  # Θ
-THERMAL_CONDUCTIVITY = DimensionSignature.create(length=1, mass=1, temp=1, time=-3)  # L M Θ T^-3
-TIME = DimensionSignature.create(time=1)  # T
-TORQUE = DimensionSignature.create(length=2, mass=1, time=-2)  # L^2 M T^-2
-TURBULENCE_ENERGY_DISSIPATION_RATE = DimensionSignature.create(length=2, time=-3)  # L^2 T^-3
-VELOCITY_ANGULAR = DimensionSignature.create(time=-1)  # T^-1
-VELOCITY_LINEAR = DimensionSignature.create(length=1, time=-1)  # L T^-1
-VISCOSITY_DYNAMIC = DimensionSignature.create(length=-1, mass=1, time=-1)  # L^-1 M T^-1
-VISCOSITY_KINEMATIC = DimensionSignature.create(length=2, time=-1)  # L^2 T^-1
-VOLUME = DimensionSignature.create(length=3)  # L^3
-VOLUMETRIC_CALORIFIC_HEATING_VALUE = DimensionSignature.create(length=-1, mass=1, time=-2)  # L^-1 M T^-2
-VOLUMETRIC_COEFFICIENT_OF_EXPANSION = DimensionSignature.create(length=-3, mass=1, temp=-1)  # L^-3 M Θ^-1
-VOLUMETRIC_FLOW_RATE = DimensionSignature.create(length=3, time=-1)  # L^3 T^-1
-VOLUMETRIC_FLUX = DimensionSignature.create(length=1, time=-1)  # L T^-1
-VOLUMETRIC_MASS_FLOW_RATE = DimensionSignature.create(length=-3, mass=1, time=-1)  # L^-3 M T^-1
-VOLUME_FRACTION_OF_I = DimensionSignature(BaseDimension.DIMENSIONLESS)  # Dimensionless
-WAVENUMBER = DimensionSignature.create(length=-1)  # L^-1

qnty/equation.py

@@ -1,216 +0,0 @@
-"""
-Equation System
-===============
-
-Mathematical equations for qnty variables with solving capabilities.
-"""
-
-from __future__ import annotations
-
-from typing import cast
-
-from .expression import Expression, VariableReference
-from .variable import TypeSafeVariable
-
-
-class Equation:
-    """Represents a mathematical equation with left-hand side equal to right-hand side."""
-    
-    def __init__(self, name: str, lhs: TypeSafeVariable | Expression, rhs: Expression):
-        self.name = name
-        
-        # Convert Variable to VariableReference if needed
-        # Use duck typing to avoid circular import
-        if hasattr(lhs, 'name') and hasattr(lhs, 'quantity') and hasattr(lhs, 'is_known'):
-            # It's a TypeSafeVariable-like object
-            self.lhs = VariableReference(cast('TypeSafeVariable', lhs))
-        else:
-            # It's already an Expression
-            self.lhs = cast(Expression, lhs)
-            
-        self.rhs = rhs
-        self.variables = self.get_all_variables()
-    
-    def get_all_variables(self) -> set[str]:
-        """Get all variable names used in this equation."""
-        # Both lhs and rhs should be Expressions after __init__ conversion
-        lhs_vars = self.lhs.get_variables()
-        rhs_vars = self.rhs.get_variables()
-        return lhs_vars | rhs_vars
-    
-    def get_unknown_variables(self, known_vars: set[str]) -> set[str]:
-        """Get variables that are unknown (not in known_vars set)."""
-        return self.variables - known_vars
-    
-    def get_known_variables(self, known_vars: set[str]) -> set[str]:
-        """Get variables that are known (in known_vars set)."""
-        return self.variables & known_vars
-    
-    def can_solve_for(self, target_var: str, known_vars: set[str]) -> bool:
-        """Check if this equation can solve for target_var given known_vars."""
-        if target_var not in self.variables:
-            return False
-        # Direct assignment case: lhs is the variable
-        if isinstance(self.lhs, VariableReference) and self.lhs.name == target_var:
-            rhs_vars = self.rhs.get_variables()
-            return rhs_vars.issubset(known_vars)
-        unknown_vars = self.get_unknown_variables(known_vars)
-        # Can solve if target_var is the only unknown
-        return unknown_vars == {target_var}
-    
-    def solve_for(self, target_var: str, variable_values: dict[str, TypeSafeVariable]) -> TypeSafeVariable:
-        """
-        Solve the equation for target_var.
-        Returns the target variable with updated quantity.
-        """
-        if target_var not in self.variables:
-            raise ValueError(f"Variable '{target_var}' not found in equation")
-        
-        # Handle direct assignment: target = expression
-        if isinstance(self.lhs, VariableReference) and self.lhs.name == target_var:
-            # Direct assignment: target_var = rhs
-            result_qty = self.rhs.evaluate(variable_values)
-            
-            # Update existing variable object to preserve references
-            var_obj = variable_values.get(target_var)
-            if var_obj is not None:
-                # Convert result to the target variable's original unit if it had one
-                if var_obj.quantity is not None and var_obj.quantity.unit is not None:
-                    # Convert to the target variable's defined unit
-                    try:
-                        result_qty = result_qty.to(var_obj.quantity.unit)
-                    except Exception:
-                        # If conversion fails, keep the calculated unit
-                        pass
-                
-                var_obj.quantity = result_qty
-                var_obj.is_known = True
-                return var_obj
-            
-            # Create new variable if not found - this shouldn't happen in normal usage
-            raise ValueError(f"Variable '{target_var}' not found in variable_values")
-        
-        # For more complex equations, we would need algebraic manipulation
-        # Currently focusing on direct assignment which covers most engineering cases
-        raise NotImplementedError(f"Cannot solve for {target_var} in equation {self}. "
-                                f"Only direct assignment equations (var = expression) are supported.")
-    
-    def check_residual(self, variable_values: dict[str, TypeSafeVariable], tolerance: float = 1e-10) -> bool:
-        """
-        Check if equation is satisfied by evaluating residual (LHS - RHS).
-        Returns True if |residual| < tolerance, accounting for units.
-        """
-        try:
-            # Both lhs and rhs should be Expressions after __init__ conversion
-            lhs_value = self.lhs.evaluate(variable_values)
-            rhs_value = self.rhs.evaluate(variable_values)
-            
-            # Check dimensional compatibility
-            if lhs_value._dimension_sig != rhs_value._dimension_sig:
-                return False
-            
-            # Convert to same units for comparison
-            rhs_converted = rhs_value.to(lhs_value.unit)
-            residual = abs(lhs_value.value - rhs_converted.value)
-            
-            return residual < tolerance
-        except Exception:
-            return False
-    
-    def __str__(self) -> str:
-        return f"{self.lhs} = {self.rhs}"
-    
-    def __repr__(self) -> str:
-        return f"Equation(name='{self.name}', lhs={self.lhs!r}, rhs={self.rhs!r})"
-
-
-class EquationSystem:
-    """System of equations that can be solved together."""
-    
-    def __init__(self, equations: list[Equation] | None = None):
-        self.equations = equations or []
-        self.variables = {}  # Dict[str, TypeSafeVariable]
-    
-    def add_equation(self, equation: Equation):
-        """Add an equation to the system."""
-        self.equations.append(equation)
-    
-    def add_variable(self, variable: TypeSafeVariable):
-        """Add a variable to the system."""
-        self.variables[variable.name] = variable
-    
-    def get_known_variables(self) -> set[str]:
-        """Get names of all known variables."""
-        return {name for name, var in self.variables.items() if var.is_known and var.quantity is not None}
-    
-    def get_unknown_variables(self) -> set[str]:
-        """Get names of all unknown variables."""
-        return {name for name, var in self.variables.items() if not var.is_known or var.quantity is None}
-    
-    def can_solve_any(self) -> bool:
-        """Check if any equation can be solved with current known variables."""
-        known_vars = self.get_known_variables()
-        unknown_vars = self.get_unknown_variables()
-        
-        for equation in self.equations:
-            for unknown_var in unknown_vars:
-                if equation.can_solve_for(unknown_var, known_vars):
-                    return True
-        return False
-    
-    def solve_step(self) -> bool:
-        """Solve one step - find and solve one equation. Returns True if progress made."""
-        known_vars = self.get_known_variables()
-        unknown_vars = self.get_unknown_variables()
-        
-        # Find an equation that can be solved
-        for equation in self.equations:
-            for unknown_var in unknown_vars:
-                if equation.can_solve_for(unknown_var, known_vars):
-                    # Solve for this variable
-                    equation.solve_for(unknown_var, self.variables)
-                    return True  # Progress made
-        
-        return False  # No progress possible
-    
-    def solve(self, max_iterations: int = 100) -> bool:
-        """Solve the system iteratively. Returns True if fully solved."""
-        for _ in range(max_iterations):
-            if not self.can_solve_any():
-                break
-            if not self.solve_step():
-                break
-        
-        # Check if all variables are known
-        unknown_vars = self.get_unknown_variables()
-        return len(unknown_vars) == 0
-    
-    def get_solving_order(self) -> list[str]:
-        """Get the order in which variables can be solved."""
-        order = []
-        temp_system = EquationSystem(self.equations.copy())
-        temp_system.variables = self.variables.copy()
-        
-        while temp_system.can_solve_any():
-            known_vars = temp_system.get_known_variables()
-            unknown_vars = temp_system.get_unknown_variables()
-            
-            # Find next solvable variable
-            for equation in temp_system.equations:
-                for unknown_var in unknown_vars:
-                    if equation.can_solve_for(unknown_var, known_vars):
-                        order.append(unknown_var)
-                        # Mark as known for next iteration
-                        temp_system.variables[unknown_var].is_known = True
-                        break
-                else:
-                    continue
-                break
-        
-        return order
-    
-    def __str__(self) -> str:
-        return f"EquationSystem({len(self.equations)} equations, {len(self.variables)} variables)"
-    
-    def __repr__(self) -> str:
-        return f"EquationSystem(equations={self.equations!r})"