radia 1.3.10__tar.gz → 1.3.13__tar.gz

{radia-1.3.10/src/radia.egg-info → radia-1.3.13}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: radia
-Version: 1.3.10
+Version: 1.3.13
 Summary: Radia 3D Magnetostatics with NGSolve Integration and OpenMP Parallelization
 Home-page: https://github.com/ksugahar/Radia_NGSolve
 Author: Pascal Elleaume

{radia-1.3.10 → radia-1.3.13}/docs/MESH_MSC_API_DESIGN.md

@@ -27,16 +27,17 @@ Understanding the differences is critical for choosing the right approach.
 
 ### Method Comparison Table
 
-| Property | ObjRecMag (Analytical) | ObjPolyhdr Hex (MSC) | ObjPolyhdr Tet (MSC) |
-|----------|----------------------|---------------------|---------------------|
-| **Element Type** | Axis-aligned rectangular | General hexahedron | Tetrahedron |
-| **Implementation** | 8-vertex atan formula | 6-face surface integral | 4-face surface integral |
+| Property | ObjRecMag (Analytical) | ObjPolyhdr Hex (6-face MSC) | ObjPolyhdr Tet (4-face MSC) |
+|----------|----------------------|------------------------------|------------------------------|
+| **Element Type** | Axis-aligned rectangular | General hexahedron (6 quad faces) | Tetrahedron (4 triangular faces) |
+| **Implementation** | 8-vertex atan formula | 6-quad -> 12-tri surface integral | 4-face surface integral |
 | **Source File** | rad_rectangular_block.cpp | rad_polyhedron.cpp | rad_polyhedron.cpp |
 | **Geometry Constraint** | Axis-aligned only | Any convex hexahedron | Any tetrahedron |
 | **Speed** | Fastest | Medium | Medium |
-| **Accuracy** | Exact (analytical) | High (numerical) | High (numerical) |
-| **Mesh Import** | ObjDivMag only | External mesh | External mesh |
+| **Accuracy** | Exact (analytical) | High (MSC matching ELF_MAGIC) | High (MSC matching ELF_MAGIC) |
+| **Mesh Import** | ObjDivMag only | External mesh (6-face format) | External mesh (4-face format) |
 | **DOF per Element** | 3 (Mx, My, Mz) | 3 (Mx, My, Mz) | 3 (Mx, My, Mz) |
+| **Radia Face Format** | N/A | 1-indexed quad faces | 1-indexed triangular faces |
 
 ### Detailed Method Descriptions
 
@@ -657,5 +658,6 @@ The NGSolve integration benefits especially from:
 ---
 
 **Created**: 2025-12-05
+**Updated**: 2025-12-11 (6-face hexahedral MSC implementation)
 **Author**: Claude Code
 **Project**: Radia Magnetic Field Computation

{radia-1.3.10 → radia-1.3.13}/docs/MMM_MSC_IMPLEMENTATION.md

@@ -109,7 +109,7 @@ void radTPolygon::B_comp(radTField* FieldPtr)
 | **Method** | MMM with MSC | MMM with MSC |
 | **Rectangular elements** | Analytical formula (8-vertex atan) | MSC (6 quad faces -> 12 triangles) |
 | **Tetrahedral elements** | MSC (4 triangular faces) | MSC (4 triangular faces) |
-| **Hexahedral elements** | Analytical formula (radTRecMag) | MSC (6 quad faces) |
+| **Hexahedral elements** | Analytical (radTRecMag) or MSC (radTPolyhedron) | MSC (6 quad faces) |
 | **General polyhedra** | MSC (N polygon faces) | Not directly supported |
 | **Polygon field formula** | log + atan (same as ELF_MAGIC) | log + atan |
 | **Self-demagnetization** | Correct (N=1/3 for cube) | Correct (N=1/3 for cube) |
@@ -170,6 +170,6 @@ The performance difference comes from:
 
 ---
 
-**Last Updated**: 2025-12-05
+**Last Updated**: 2025-12-11 (6-face hexahedral MSC added)
 **Author**: Claude Code
 **Project**: Radia Magnetic Field Computation

{radia-1.3.10 → radia-1.3.13}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "radia"
-version = "1.3.10"
+version = "1.3.13"
 description = "Radia 3D Magnetostatics with NGSolve Integration and OpenMP Parallelization"
 readme = "README.md"
 requires-python = ">=3.12"

{radia-1.3.10 → radia-1.3.13}/src/radia/__init__.py

@@ -2,7 +2,7 @@
 # This module re-exports all symbols from the C++ extension module (radia.pyd)
 # so that 'import radia' works correctly when installed via pip
 
-__version__ = "1.3.10"
+__version__ = "1.3.13"
 
 # Import all symbols from the C++ extension module
 try:

{radia-1.3.10 → radia-1.3.13/src/radia.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: radia
-Version: 1.3.10
+Version: 1.3.13
 Summary: Radia 3D Magnetostatics with NGSolve Integration and OpenMP Parallelization
 Home-page: https://github.com/ksugahar/Radia_NGSolve
 Author: Pascal Elleaume

{radia-1.3.10 → radia-1.3.13}/src/radia.egg-info/SOURCES.txt

@@ -67,7 +67,6 @@ examples/cube_uniform_field/nonlinear/benchmark_tetra_netgen.py
 examples/cube_uniform_field/nonlinear/benchmark_tetra_unified.py
 examples/cube_uniform_field/nonlinear/compare_radia_elfmagic_field.py
 examples/cube_uniform_field/nonlinear/compare_tetra_methods.py
-examples/cube_uniform_field/nonlinear/test_lu_vs_bicgstab.py
 examples/electromagnet/README.md
 examples/electromagnet/Radia_model.vtk
 examples/electromagnet/York.vtk

radia-1.3.10/examples/cube_uniform_field/nonlinear/test_lu_vs_bicgstab.py

@@ -1,82 +0,0 @@
-"""
-Compare LU vs BiCGSTAB solvers on same tetrahedral mesh.
-The results should match if both converged correctly.
-"""
-import sys
-import os
-sys.path.insert(0, os.path.join(os.path.dirname(__file__), '../../../build/Release'))
-sys.path.insert(0, os.path.join(os.path.dirname(__file__), '../../../src/radia'))
-import radia as rad
-import numpy as np
-
-from ngsolve import Mesh
-from netgen.csg import CSGeometry, OrthoBrick, Pnt
-from netgen_mesh_import import netgen_mesh_to_radia
-
-MU_0 = 4 * np.pi * 1e-7
-H_EXT = 50000.0
-
-# B-H curve
-BH_DATA = [
-    [0.0, 0.0], [100.0, 0.1], [200.0, 0.3], [500.0, 0.8], [1000.0, 1.2],
-    [2000.0, 1.5], [5000.0, 1.7], [10000.0, 1.8], [50000.0, 2.0], [100000.0, 2.1],
-]
-HM_DATA = [[h, b/MU_0 - h] for h, b in BH_DATA]
-
-def run_comparison(maxh):
-    print(f'\n{"="*60}')
-    print(f'maxh = {maxh}m')
-    print("="*60)
-
-    # Create mesh ONCE
-    geo = CSGeometry()
-    box = OrthoBrick(Pnt(-0.5, -0.5, -0.5), Pnt(0.5, 0.5, 0.5))
-    geo.Add(box)
-    ngmesh = geo.GenerateMesh(maxh=maxh)
-    mesh = Mesh(ngmesh)
-
-    n_elem = mesh.ne
-    print(f'Elements: {n_elem}')
-
-    results = {}
-    for method, name in [(0, 'LU'), (1, 'BiCGSTAB')]:
-        rad.FldUnits('m')
-        rad.UtiDelAll()
-
-        # Import SAME mesh
-        mag_obj = netgen_mesh_to_radia(mesh, material={'magnetization': [0, 0, 0]}, units='m')
-
-        # Apply material
-        mat = rad.MatSatIsoTab(HM_DATA)
-        rad.MatApl(mag_obj, mat)
-
-        # External field
-        B_ext = MU_0 * H_EXT
-        ext = rad.ObjBckg([0, 0, B_ext])
-        grp = rad.ObjCnt([mag_obj, ext])
-
-        # Solve
-        result = rad.Solve(grp, 0.001, 1000, method)
-
-        # Get magnetization
-        all_M = rad.ObjM(mag_obj)
-        M_avg_z = np.mean([m[1][2] for m in all_M])
-        n_iter = int(result[3]) if result[3] else 0
-        residual = result[0] if result[0] else 0.0
-
-        print(f'{name:10s}: iter={n_iter:4d}, residual={residual:10.2f}, M_avg_z={M_avg_z:12.0f} A/m')
-        results[name] = {'M_avg_z': M_avg_z, 'iter': n_iter, 'residual': residual}
-
-    diff = abs(results['LU']['M_avg_z'] - results['BiCGSTAB']['M_avg_z'])
-    pct = 100 * diff / max(results['LU']['M_avg_z'], results['BiCGSTAB']['M_avg_z'])
-    print(f'Difference: {diff:.0f} A/m ({pct:.1f}%)')
-
-    return results
-
-# Run for multiple mesh sizes
-if __name__ == '__main__':
-    print('\nComparison of LU vs BiCGSTAB solvers on tetrahedral meshes')
-    print('(Same mesh, same problem - results should be identical)')
-
-    for maxh in [0.4, 0.35, 0.3]:
-        run_comparison(maxh)