radia 1.3.2__tar.gz → 1.3.4__tar.gz

{radia-1.3.2 → radia-1.3.4}/LICENSE

@@ -1,93 +1,93 @@
-RADIA SOFTWARE LICENSE
-======================
-
-This software contains multiple components with different licenses:
-
-1. RADIA Core (BSD-style License)
-2. HACApK_LH-Cimplm Library (MIT License)
-
-All licenses are permissive open-source licenses allowing redistribution
-and modification.
-
-================================================================================
-PART 1: RADIA CORE
-================================================================================
-
-Copyright © 1997 - 2018, European Synchrotron Radiation Facility
-
-The RADIA software is distributed subject to the following license conditions:
-
-Redistribution and use in source and binary forms, with or without
-modification, are permitted provided that the following conditions are met:
-
-1. Redistributions of source code must retain the above copyright notice,
-   this list of conditions and the following disclaimer.
-
-2. Redistributions in binary form must reproduce the above copyright notice,
-   this list of conditions and the following disclaimer in the documentation
-   and/or other materials provided with the distribution.
-
-DISCLAIMER: THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS "AS IS"
-AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
-IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
-ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE FOR
-ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
-DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
-SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
-CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
-OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
-OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
-
-The views and conclusions contained in the software and documentation are
-those of the authors and should not be interpreted as representing official
-policies, either expressed or implied, of the RADIA project.
-
-================================================================================
-PART 2: HACApK_LH-Cimplm (H-matrix Library)
-================================================================================
-
-Software Name : HACApK
-Version : 1.3.0
-Location: src/ext/HACApK_LH-Cimplm/
-
-Copyright (c) 2015 Akihiro Ida and Takeshi Iwashita
-
-MIT License
-
-Permission is hereby granted, free of charge, to any person obtaining a copy
-of this software and associated documentation files (the "Software"), to deal
-in the Software without restriction, including without limitation the rights
-to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
-copies of the Software, and to permit persons to whom the Software is
-furnished to do so, subject to the following conditions:
-
-The above copyright notice and this permission notice shall be included in all
-copies or substantial portions of the Software.
-
-THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
-IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
-FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
-AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
-LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
-OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
-SOFTWARE.
-
-ppOpen-HPC project:
-  Open Source Infrastructure for Development and Execution of Large-Scale
-  Scientific Applications on Post-Peta-Scale Supercomputers with Automatic
-  Tuning (AT).
-
-Sponsorship:
-  Japan Science and Technology Agency (JST), Basic Research Programs: CREST,
-  Development of System Software Technologies for post-Peta Scale High
-  Performance Computing.
-
-================================================================================
-ADDITIONAL COMPONENTS
-================================================================================
-
-This software may include additional third-party libraries and components,
-each governed by their respective licenses. Please refer to the individual
-component directories for specific license information.
-
-================================================================================
+RADIA SOFTWARE LICENSE
+======================
+
+This software contains multiple components with different licenses:
+
+1. RADIA Core (BSD-style License)
+2. HACApK_LH-Cimplm Library (MIT License)
+
+All licenses are permissive open-source licenses allowing redistribution
+and modification.
+
+================================================================================
+PART 1: RADIA CORE
+================================================================================
+
+Copyright © 1997 - 2018, European Synchrotron Radiation Facility
+
+The RADIA software is distributed subject to the following license conditions:
+
+Redistribution and use in source and binary forms, with or without
+modification, are permitted provided that the following conditions are met:
+
+1. Redistributions of source code must retain the above copyright notice,
+   this list of conditions and the following disclaimer.
+
+2. Redistributions in binary form must reproduce the above copyright notice,
+   this list of conditions and the following disclaimer in the documentation
+   and/or other materials provided with the distribution.
+
+DISCLAIMER: THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS "AS IS"
+AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE FOR
+ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
+SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
+CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
+OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+
+The views and conclusions contained in the software and documentation are
+those of the authors and should not be interpreted as representing official
+policies, either expressed or implied, of the RADIA project.
+
+================================================================================
+PART 2: HACApK_LH-Cimplm (H-matrix Library)
+================================================================================
+
+Software Name : HACApK
+Version : 1.3.0
+Location: src/ext/HACApK_LH-Cimplm/
+
+Copyright (c) 2015 Akihiro Ida and Takeshi Iwashita
+
+MIT License
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.
+
+ppOpen-HPC project:
+  Open Source Infrastructure for Development and Execution of Large-Scale
+  Scientific Applications on Post-Peta-Scale Supercomputers with Automatic
+  Tuning (AT).
+
+Sponsorship:
+  Japan Science and Technology Agency (JST), Basic Research Programs: CREST,
+  Development of System Software Technologies for post-Peta Scale High
+  Performance Computing.
+
+================================================================================
+ADDITIONAL COMPONENTS
+================================================================================
+
+This software may include additional third-party libraries and components,
+each governed by their respective licenses. Please refer to the individual
+component directories for specific license information.
+
+================================================================================

{radia-1.3.2/src/radia.egg-info → radia-1.3.4}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: radia
-Version: 1.3.2
+Version: 1.3.4
 Summary: Radia 3D Magnetostatics with NGSolve Integration and OpenMP Parallelization
 Home-page: https://github.com/ksugahar/Radia_NGSolve
 Author: Pascal Elleaume
@@ -88,12 +88,12 @@ pip install radia-ngsolve
 # 1. Build radia.pyd (core module)
 .\Build.ps1
 
-# 2. Build rad_ngsolve.pyd (optional, for NGSolve integration)
+# 2. Build radia_ngsolve.pyd (optional, for NGSolve integration)
 .\Build_NGSolve.ps1
 
 # Outputs:
 # - dist/radia.pyd
-# - build/Release/rad_ngsolve.pyd
+# - build/Release/radia_ngsolve.pyd
 ```
 
 See [README_BUILD.md](README_BUILD.md) for detailed build instructions.
@@ -117,12 +117,12 @@ print(f"Field: {field} T")
 
 ### NGSolve Integration
 
-The `rad_ngsolve` module provides a C++ CoefficientFunction interface for using Radia magnetic fields in NGSolve FEM analysis.
+The `radia_ngsolve` module provides a C++ CoefficientFunction interface for using Radia magnetic fields in NGSolve FEM analysis.
 
 **Function Specification:**
 
 ```python
-rad_ngsolve.RadiaField(radia_obj, field_type='b')
+radia_ngsolve.RadiaField(radia_obj, field_type='b')
 ```
 
 **Parameters:**
@@ -149,7 +149,7 @@ import ngsolve
 from ngsolve import Mesh, H1, GridFunction
 
 import radia as rad
-import rad_ngsolve
+import radia_ngsolve
 
 # Create Radia magnet
 magnet = rad.ObjRecMag([0,0,0], [20,20,20], [0,0,1.2])  # mm units
@@ -157,10 +157,10 @@ rad.MatApl(magnet, rad.MatPM(1.2, 900000, [0,0,1]))  # NdFeB permanent magnet
 rad.Solve(magnet, 0.0001, 10000)
 
 # Create NGSolve CoefficientFunction for different field types
-B_field = rad_ngsolve.RadiaField(magnet, 'b')  # Flux density [T]
-H_field = rad_ngsolve.RadiaField(magnet, 'h')  # Magnetic field [A/m]
-A_field = rad_ngsolve.RadiaField(magnet, 'a')  # Vector potential [T·m]
-M_field = rad_ngsolve.RadiaField(magnet, 'm')  # Magnetization [A/m]
+B_field = radia_ngsolve.RadiaField(magnet, 'b')  # Flux density [T]
+H_field = radia_ngsolve.RadiaField(magnet, 'h')  # Magnetic field [A/m]
+A_field = radia_ngsolve.RadiaField(magnet, 'a')  # Vector potential [T·m]
+M_field = radia_ngsolve.RadiaField(magnet, 'm')  # Magnetization [A/m]
 
 # Use in FEM analysis (NGSolve mesh in meters)
 gf = GridFunction(fes)
@@ -177,13 +177,10 @@ See [examples/Radia_to_NGSolve_CoefficientFunction/](examples/Radia_to_NGSolve_C
 - [docs/DIRECTORY_STRUCTURE.md](docs/DIRECTORY_STRUCTURE.md) - Project structure
 
 ### NGSolve Integration
-- [RAD_NGSOLVE_BUILD_SUCCESS.md](RAD_NGSOLVE_BUILD_SUCCESS.md) - Complete rad_ngsolve documentation
+- [RAD_NGSOLVE_BUILD_SUCCESS.md](RAD_NGSOLVE_BUILD_SUCCESS.md) - Complete radia_ngsolve documentation
 - [examples/Radia_to_NGSolve_CoefficientFunction/README.md](examples/Radia_to_NGSolve_CoefficientFunction/README.md) - NGSolve examples overview
 - [examples/Radia_to_NGSolve_CoefficientFunction/EXAMPLES_GUIDE.md](examples/Radia_to_NGSolve_CoefficientFunction/EXAMPLES_GUIDE.md) - Detailed usage guide
-- [tests/test_rad_ngsolve.py](tests/test_rad_ngsolve.py) - Integration tests
-
-### Visualization
-- [docs/PYVISTA_VIEWER.md](examples/2024_02_03_振分電磁石/PYVISTA_VIEWER.md) - PyVista viewer guide
+- [tests/test_radia_ngsolve.py](tests/test_radia_ngsolve.py) - Integration tests
 
 ### Development
 - [docs/TAB_CONVERSION_REPORT.md](docs/TAB_CONVERSION_REPORT.md) - Code style conversion
@@ -239,7 +236,7 @@ python tests/test_radia.py
 python tests/test_advanced.py
 
 # NGSolve integration test
-python tests/test_rad_ngsolve.py
+python tests/test_radia_ngsolve.py
 
 # OpenMP performance test
 python tests/test_parallel_performance.py
@@ -248,7 +245,7 @@ python tests/test_parallel_performance.py
 pytest tests/
 
 # Run specific test suite
-pytest tests/test_rad_ngsolve.py -v
+pytest tests/test_radia_ngsolve.py -v
 ```
 
 See [tests/README.md](tests/README.md) for detailed testing documentation.
@@ -296,7 +293,7 @@ export_geometry_to_vtk(mag, 'geometry.vtk')
 ### Runtime Requirements
 - Python 3.12
 - NumPy
-- NGSolve (optional, for FEM coupling via rad_ngsolve)
+- NGSolve (optional, for FEM coupling via radia_ngsolve)
 - PyVista (optional, for 3D visualization)
 
 ## Changes from Original Radia
@@ -348,7 +345,7 @@ See:
 
 **This Fork**:
 - OpenMP parallelization
-- NGSolve C++ integration (rad_ngsolve)
+- NGSolve C++ integration (radia_ngsolve)
 - Python 3.12 optimization
 - Build system modernization
 - PyVista integration

{radia-1.3.2 → radia-1.3.4}/README.md

@@ -46,12 +46,12 @@ pip install radia-ngsolve
 # 1. Build radia.pyd (core module)
 .\Build.ps1
 
-# 2. Build rad_ngsolve.pyd (optional, for NGSolve integration)
+# 2. Build radia_ngsolve.pyd (optional, for NGSolve integration)
 .\Build_NGSolve.ps1
 
 # Outputs:
 # - dist/radia.pyd
-# - build/Release/rad_ngsolve.pyd
+# - build/Release/radia_ngsolve.pyd
 ```
 
 See [README_BUILD.md](README_BUILD.md) for detailed build instructions.
@@ -75,12 +75,12 @@ print(f"Field: {field} T")
 
 ### NGSolve Integration
 
-The `rad_ngsolve` module provides a C++ CoefficientFunction interface for using Radia magnetic fields in NGSolve FEM analysis.
+The `radia_ngsolve` module provides a C++ CoefficientFunction interface for using Radia magnetic fields in NGSolve FEM analysis.
 
 **Function Specification:**
 
 ```python
-rad_ngsolve.RadiaField(radia_obj, field_type='b')
+radia_ngsolve.RadiaField(radia_obj, field_type='b')
 ```
 
 **Parameters:**
@@ -107,7 +107,7 @@ import ngsolve
 from ngsolve import Mesh, H1, GridFunction
 
 import radia as rad
-import rad_ngsolve
+import radia_ngsolve
 
 # Create Radia magnet
 magnet = rad.ObjRecMag([0,0,0], [20,20,20], [0,0,1.2])  # mm units
@@ -115,10 +115,10 @@ rad.MatApl(magnet, rad.MatPM(1.2, 900000, [0,0,1]))  # NdFeB permanent magnet
 rad.Solve(magnet, 0.0001, 10000)
 
 # Create NGSolve CoefficientFunction for different field types
-B_field = rad_ngsolve.RadiaField(magnet, 'b')  # Flux density [T]
-H_field = rad_ngsolve.RadiaField(magnet, 'h')  # Magnetic field [A/m]
-A_field = rad_ngsolve.RadiaField(magnet, 'a')  # Vector potential [T·m]
-M_field = rad_ngsolve.RadiaField(magnet, 'm')  # Magnetization [A/m]
+B_field = radia_ngsolve.RadiaField(magnet, 'b')  # Flux density [T]
+H_field = radia_ngsolve.RadiaField(magnet, 'h')  # Magnetic field [A/m]
+A_field = radia_ngsolve.RadiaField(magnet, 'a')  # Vector potential [T·m]
+M_field = radia_ngsolve.RadiaField(magnet, 'm')  # Magnetization [A/m]
 
 # Use in FEM analysis (NGSolve mesh in meters)
 gf = GridFunction(fes)
@@ -135,13 +135,10 @@ See [examples/Radia_to_NGSolve_CoefficientFunction/](examples/Radia_to_NGSolve_C
 - [docs/DIRECTORY_STRUCTURE.md](docs/DIRECTORY_STRUCTURE.md) - Project structure
 
 ### NGSolve Integration
-- [RAD_NGSOLVE_BUILD_SUCCESS.md](RAD_NGSOLVE_BUILD_SUCCESS.md) - Complete rad_ngsolve documentation
+- [RAD_NGSOLVE_BUILD_SUCCESS.md](RAD_NGSOLVE_BUILD_SUCCESS.md) - Complete radia_ngsolve documentation
 - [examples/Radia_to_NGSolve_CoefficientFunction/README.md](examples/Radia_to_NGSolve_CoefficientFunction/README.md) - NGSolve examples overview
 - [examples/Radia_to_NGSolve_CoefficientFunction/EXAMPLES_GUIDE.md](examples/Radia_to_NGSolve_CoefficientFunction/EXAMPLES_GUIDE.md) - Detailed usage guide
-- [tests/test_rad_ngsolve.py](tests/test_rad_ngsolve.py) - Integration tests
-
-### Visualization
-- [docs/PYVISTA_VIEWER.md](examples/2024_02_03_振分電磁石/PYVISTA_VIEWER.md) - PyVista viewer guide
+- [tests/test_radia_ngsolve.py](tests/test_radia_ngsolve.py) - Integration tests
 
 ### Development
 - [docs/TAB_CONVERSION_REPORT.md](docs/TAB_CONVERSION_REPORT.md) - Code style conversion
@@ -197,7 +194,7 @@ python tests/test_radia.py
 python tests/test_advanced.py
 
 # NGSolve integration test
-python tests/test_rad_ngsolve.py
+python tests/test_radia_ngsolve.py
 
 # OpenMP performance test
 python tests/test_parallel_performance.py
@@ -206,7 +203,7 @@ python tests/test_parallel_performance.py
 pytest tests/
 
 # Run specific test suite
-pytest tests/test_rad_ngsolve.py -v
+pytest tests/test_radia_ngsolve.py -v
 ```
 
 See [tests/README.md](tests/README.md) for detailed testing documentation.
@@ -254,7 +251,7 @@ export_geometry_to_vtk(mag, 'geometry.vtk')
 ### Runtime Requirements
 - Python 3.12
 - NumPy
-- NGSolve (optional, for FEM coupling via rad_ngsolve)
+- NGSolve (optional, for FEM coupling via radia_ngsolve)
 - PyVista (optional, for 3D visualization)
 
 ## Changes from Original Radia
@@ -306,7 +303,7 @@ See:
 
 **This Fork**:
 - OpenMP parallelization
-- NGSolve C++ integration (rad_ngsolve)
+- NGSolve C++ integration (radia_ngsolve)
 - Python 3.12 optimization
 - Build system modernization
 - PyVista integration

{radia-1.3.2 → radia-1.3.4}/docs/API_EXTENSIONS.md

@@ -4,7 +4,7 @@ This document describes custom extensions to the original ESRF Radia API.
 
 **Original Documentation**: https://www.esrf.fr/home/Accelerators/instrumentation--equipment/Software/Radia/Documentation/ReferenceGuide.html
 
-**Date**: 2025-11-08
+**Date**: 2025-11-27
 **Status**: Active development
 
 ---
@@ -21,7 +21,7 @@ This document describes custom extensions to the original ESRF Radia API.
 - [Performance Features](#performance-features)
   - [SolverHMatrixDisable/Enable](#solverhmatrixdisableenable)
 - [NGSolve Integration](#ngsolve-integration)
-  - [rad_ngsolve.RadiaField](#rad_ngsolveradiafield)
+  - [radia_ngsolve.RadiaField](#radia_ngsolveradiafield)
 
 ---
 
@@ -106,7 +106,7 @@ For use with NGSolve's CoefficientFunction:
 
 ```python
 import ngsolve as ngs
-import rad_ngsolve
+import radia_ngsolve
 
 # Create NGSolve coefficient function (in meters)
 def ngsolve_field_wrapper(x, y, z):
@@ -116,7 +116,7 @@ def ngsolve_field_wrapper(x, y, z):
 	return B
 
 # Create background field
-bg_cf = rad_ngsolve.BackgroundFieldCF(ngsolve_field_wrapper)
+bg_cf = radia_ngsolve.BackgroundFieldCF(ngsolve_field_wrapper)
 ```
 
 **Notes**:
@@ -416,8 +416,97 @@ rad.SolverHMatrixEnable()
 
 ---
 
+### SolverTetraMethod
+
+**Purpose**: Set the field computation method for tetrahedral mesh elements.
+
+**Syntax**:
+```python
+rad.SolverTetraMethod(method)
+```
+
+**Parameters**:
+- `method`: Integer selecting the computation method
+  - `0`: Original Radia polygon-based method (default)
+  - `1`: Analytical method using magnetic charge formulation
+
+**Returns**: None
+
+**Background**:
+
+Tetrahedral elements in Radia can compute their magnetic field contribution using two different algorithms:
+
+| Method | Description | Best For |
+|--------|-------------|----------|
+| 0 | Original Radia polygon method | General purpose, standard accuracy |
+| 1 | Analytical method | High-permeability materials (mu_r > 100) |
+
+**Method 0 (Original)**:
+- Uses Radia's standard polygon-based field calculation
+- Treats each tetrahedral face as a magnetic charge sheet
+- Well-tested and reliable for most applications
+
+**Method 1 (Analytical)**:
+- Uses analytical formulas for magnetic field from surface charges
+- May provide better accuracy for high-permeability soft magnetic materials
+- Mathematically equivalent to Method 0 but with different numerical implementation
+
+**Usage Example**:
+
+```python
+import radia as rad
+from netgen_mesh_import import netgen_mesh_to_radia
+
+# Use analytical method for high-permeability material
+rad.SolverTetraMethod(1)
+
+# Import tetrahedral mesh with soft iron (mu_r = 4000)
+mesh_data = netgen_mesh_to_radia(
+    'soft_iron.vol',
+    material={'magnetization': [0, 0, 0]},  # Zero initial M
+    units='m'
+)
+
+# Apply linear material
+mat = rad.MatLin(3999)  # ksi = mu_r - 1
+rad.MatApl(mesh_data, mat)
+
+# Solve
+rad.Solve(mesh_data, 0.0001, 10000)
+
+# Compute field
+B = rad.Fld(mesh_data, 'b', [0.1, 0.1, 0.1])
+```
+
+**Migration from Environment Variable**:
+
+Previously, the tetrahedral method was controlled via the `RADIA_TETRA_METHOD` environment variable. This has been replaced by the `SolverTetraMethod()` API for better discoverability and control.
+
+```python
+# Old way (deprecated):
+import os
+os.environ['RADIA_TETRA_METHOD'] = '1'  # No longer works
+
+# New way (recommended):
+import radia as rad
+rad.SolverTetraMethod(1)
+```
+
+**Important Notes**:
+- Setting is global and persists for the session
+- Must be called **before** `rad.Solve()` or `rad.RlxPre()`
+- Both methods should produce identical results (within numerical precision)
+- Method 0 and Method 1 have been verified to produce 0.00% difference
+
+---
+
 ## Version History
 
+### v1.0.8 (2025-11-27)
+- Added `SolverTetraMethod()` API for tetrahedral mesh computation method control
+- Fixed ANALYTICAL method (Method 1) to match standard polygon method
+- Deprecated `RADIA_TETRA_METHOD` environment variable
+
 ### v1.0.7 (2025-11-08)
 - Added `SetRelaxSubInterval()` for LU decomposition control
 - Extended `RlxMan()` to support Method 5 (LU decomposition)
@@ -436,18 +525,18 @@ rad.SolverHMatrixEnable()
 
 ## NGSolve Integration
 
-### rad_ngsolve.RadiaField
+### radia_ngsolve.RadiaField
 
 **Purpose**: Create NGSolve CoefficientFunction for Radia magnetic field with full control over computation accuracy and performance.
 
-**Module**: `rad_ngsolve` (C++ extension)
+**Module**: `radia_ngsolve` (C++ extension)
 
 **Syntax**:
 ```python
 from ngsolve import *
-import rad_ngsolve
+import radia_ngsolve
 
-cf = rad_ngsolve.RadiaField(
+cf = radia_ngsolve.RadiaField(
     radia_obj,
     field_type='b',
     origin=None,
@@ -512,13 +601,13 @@ Controls accuracy vs speed trade-off. Smaller values = more accurate but slower.
 
 ```python
 # High accuracy (slow)
-B_cf = rad_ngsolve.RadiaField(magnet, 'b', precision=1e-8)
+B_cf = radia_ngsolve.RadiaField(magnet, 'b', precision=1e-8)
 
 # Normal accuracy (default)
-B_cf = rad_ngsolve.RadiaField(magnet, 'b')  # Uses Radia default
+B_cf = radia_ngsolve.RadiaField(magnet, 'b')  # Uses Radia default
 
 # Fast evaluation (less accurate)
-B_cf = rad_ngsolve.RadiaField(magnet, 'b', precision=1e-4)
+B_cf = radia_ngsolve.RadiaField(magnet, 'b', precision=1e-4)
 ```
 
 **Internally calls**: `rad.FldCmpPrc()` with specified precision for B, H, A, M fields.
@@ -531,13 +620,13 @@ Controls H-matrix (hierarchical matrix) acceleration.
 
 ```python
 # Disable H-matrix (accurate for small N)
-B_cf = rad_ngsolve.RadiaField(magnet, 'b', use_hmatrix=False)
+B_cf = radia_ngsolve.RadiaField(magnet, 'b', use_hmatrix=False)
 
 # Enable H-matrix (fast for large N > 1000)
-B_cf = rad_ngsolve.RadiaField(magnet, 'b', use_hmatrix=True)
+B_cf = radia_ngsolve.RadiaField(magnet, 'b', use_hmatrix=True)
 
 # Keep current setting (default)
-B_cf = rad_ngsolve.RadiaField(magnet, 'b')  # use_hmatrix=None
+B_cf = radia_ngsolve.RadiaField(magnet, 'b')  # use_hmatrix=None
 ```
 
 **When to enable H-matrix**:
@@ -558,13 +647,13 @@ B_cf = rad_ngsolve.RadiaField(magnet, 'b')  # use_hmatrix=None
 import radia as rad
 from ngsolve import *
 from netgen.occ import *
-import rad_ngsolve
+import radia_ngsolve
 
 # Create Radia magnet
 magnet = rad.ObjRecMag([0,0,0], [10,10,10], [0,0,1.2])
 
 # Create NGSolve CoefficientFunction
-B_cf = rad_ngsolve.RadiaField(magnet, 'b')
+B_cf = radia_ngsolve.RadiaField(magnet, 'b')
 
 # Use in NGSolve mesh
 mesh = Mesh(...)
@@ -583,7 +672,7 @@ u_axis = [np.cos(angle), np.sin(angle), 0]
 v_axis = [-np.sin(angle), np.cos(angle), 0]
 w_axis = [0, 0, 1]
 
-B_cf = rad_ngsolve.RadiaField(
+B_cf = radia_ngsolve.RadiaField(
     magnet, 'b',
     origin=origin,
     u_axis=u_axis,
@@ -607,7 +696,7 @@ for i in range(20):
 magnet_array = rad.ObjCnt(elements)  # 8000 elements
 
 # Use H-matrix for fast evaluation
-B_cf = rad_ngsolve.RadiaField(magnet_array, 'b', use_hmatrix=True)
+B_cf = radia_ngsolve.RadiaField(magnet_array, 'b', use_hmatrix=True)
 
 # Evaluate on large mesh
 mesh = Mesh(...)  # Large mesh
@@ -618,13 +707,13 @@ gf.Set(B_cf)  # Fast thanks to H-matrix
 **Vector potential with curl**:
 ```python
 # Get vector potential A
-A_cf = rad_ngsolve.RadiaField(magnet, 'a')
+A_cf = radia_ngsolve.RadiaField(magnet, 'a')
 
 # Compute B = curl(A)
 B_from_curl = Curl(A_cf)
 
 # Compare with direct B
-B_direct = rad_ngsolve.RadiaField(magnet, 'b')
+B_direct = radia_ngsolve.RadiaField(magnet, 'b')
 
 # Should match (within numerical precision)
 ```
@@ -684,6 +773,6 @@ See [`examples/Radia_to_NGSolve_CoefficientFunction_A/README.md`](../examples/Ra
 
 ---
 
-**Last Updated**: 2025-11-08
+**Last Updated**: 2025-11-27
 **Maintained By**: Radia Development Team
 **License**: LGPL-2.1 (modifications), BSD-style (original RADIA from ESRF)