radia 1.3.4__py3-none-any.whl → 1.3.6__py3-none-any.whl

{radia-1.3.4.dist-info → radia-1.3.6.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: radia
-Version: 1.3.4
+Version: 1.3.6
 Summary: Radia 3D Magnetostatics with NGSolve Integration and OpenMP Parallelization
 Home-page: https://github.com/ksugahar/Radia_NGSolve
 Author: Pascal Elleaume
@@ -296,6 +296,82 @@ export_geometry_to_vtk(mag, 'geometry.vtk')
 - NGSolve (optional, for FEM coupling via radia_ngsolve)
 - PyVista (optional, for 3D visualization)
 
+## Linear Material with External Field (ObjBckgCF)
+
+Radia supports soft magnetic materials (linear materials) in external fields using `ObjBckgCF`.
+
+### Working Example
+
+```python
+import radia as rad
+import numpy as np
+
+rad.UtiDelAll()
+rad.FldUnits('m')  # Use SI units
+
+# Create soft iron cube (10cm side)
+cube = rad.ObjRecMag([0, 0, 0], [0.1, 0.1, 0.1], [0, 0, 0])
+
+# Apply isotropic linear material (mu_r = 100)
+chi = 99.0  # chi = mu_r - 1
+mat = rad.MatLin(chi)  # IMPORTANT: Use single argument for isotropic
+rad.MatApl(cube, mat)
+
+# Create external field (1 Tesla in z-direction)
+def uniform_B(pos):
+    return [0, 0, 1.0]  # Returns B in Tesla
+bg = rad.ObjBckgCF(uniform_B)
+
+# Create system and solve
+system = rad.ObjCnt([cube, bg])
+result = rad.Solve(system, 0.0001, 1000)
+
+# Result: M ~ 2.3 MA/m (correct for cube with mu_r=100 in 1T field)
+M = rad.ObjM(cube)
+print(f"Magnetization: {M[1]} A/m")
+```
+
+### Important Notes
+
+1. **Isotropic materials**: Always use `MatLin(chi)` (single argument)
+   - Do NOT use `MatLin(chi, [0, 0, 1e-10])` - this creates a poorly defined anisotropic tensor
+
+2. **Callback returns B in Tesla**: The callback function for `ObjBckgCF` should return magnetic flux density B in Tesla
+
+3. **Units**: When using `rad.FldUnits('m')`, all coordinates are in meters
+
+4. **Tetrahedral elements**: Also works with `ObjPolyhdr` tetrahedral meshes
+
+### Supported Element Types
+
+| Element | API | External Field Support |
+|---------|-----|----------------------|
+| Hexahedron | `ObjRecMag` | Yes |
+| Tetrahedron | `ObjPolyhdr` | Yes |
+| Extruded Polygon | `ObjThckPgn` | Yes |
+
+## Tetrahedral Mesh Import from NGSolve/Netgen
+
+Radia can import tetrahedral meshes from NGSolve/Netgen for complex geometries.
+
+**Features:**
+- Direct Netgen/OCC mesh import to Radia MMM
+- Independent mesh optimization for FEM vs MMM
+- Nastran CTETRA face connectivity standard
+
+**Usage:** See `examples/cube_uniform_field/linear/` for examples.
+
+**Benchmark** (mu_r=100, B0=1T, 100mm cube):
+
+| Mesh | Elements | M_z (MA/m) | Error |
+|------|----------|------------|-------|
+| 1x1x1 hex | 1 | 2.3171 | 0.00% |
+| 2x2x2 hex | 8 | 2.4180 | 4.36% |
+| 5x5x5 hex | 125 | 2.6652 | 15.02% |
+
+Analytical: M = 2.3171 MA/m (N = 1/3 for cube)
+
+
 ## Changes from Original Radia
 
 ### Removed Components
@@ -379,5 +455,5 @@ See [examples/NGSolve_CoefficientFunction_to_Radia_BackgroundField/](examples/NG
 
 ---
 
-**Version**: 1.0 (OpenMP + NGSolve Edition)
-**Last Updated**: 2025-11-02
+**Version**: 1.3.4 (OpenMP + NGSolve Edition)
+**Last Updated**: 2025-11-30

{radia-1.3.4.dist-info → radia-1.3.6.dist-info}/RECORD

@@ -1,17 +1,17 @@
 python/__init__.py,sha256=oUOAjf_vY8DNy5HRU7oArAMic8urvHCR9yHSi4HFkkQ,47
 python/nastran_mesh_import.py,sha256=CSoVhZCXa85lPiTF2hlspE2clBKOD7-_sCp1bxu_IK0,18147
 python/netgen_mesh_import.py,sha256=UopXk-5bbfj1j9_hyiq8jbjb4SQXnWaeuaC7TDf17wA,19872
-python/rad_ngsolve.pyd,sha256=s5ZvpTIrjYeV_B3vQv5HrtP4u_RaPhvxuZf5aKv3oC0,566272
 python/rad_ngsolve_fast.py,sha256=GkC7ruKy3MESHsO0iRSWsrgLU4-DPPgctOi6pPpsWg4,5675
+python/radia.pyd,sha256=K6I8UDH2NI-U857_jTG_jWGmRW8zOvSFgdLRpuIK_Xo,1738752
 python/radia_coil_builder.py,sha256=nQkiAbfhueNvvxUARHdPD0C68ImidHmUQv_q4RsImeY,11253
 python/radia_field_cached.py,sha256=Bjw3ecNe3u7AAXnLob5m_tjYIY7HwB9DpgFg9a-aP_8,9509
-python/radia_ngsolve.pyd,sha256=s5ZvpTIrjYeV_B3vQv5HrtP4u_RaPhvxuZf5aKv3oC0,566272
+python/radia_ngsolve.pyd,sha256=hmPQal8FEqt9Ia5NQ05Vw35RUW5jbWFGQ5PzwXJiLww,566272
 python/radia_ngsolve_field.py,sha256=suJr4wacfYFKOkyV-5AQuHWnW5rtUMb0gSSjq8VRSXc,10166
 python/radia_ngsolve_utils.py,sha256=xZCR9DOIKMwdEjmC28rOXVZiWFY5BQYH2VfopfuVBps,8406
 python/radia_pyvista_viewer.py,sha256=JS33Mx4azGI7hUX0bzefc6zJfhv6qfRjM3Kl1bE9Mjs,4275
 python/radia_vtk_export.py,sha256=I8Vyyt9tky78Qw1xPru9f0Rii6QEmdEgTFjQtamyooc,6540
-radia-1.3.4.dist-info/licenses/LICENSE,sha256=yaWxyzG9DpJ44dDNdGni4nukwiZ8pU-r_aW-1tYNAjk,4374
-radia-1.3.4.dist-info/METADATA,sha256=iLmP3uISANgF67d-fTvmh17jl1YNvuHjiEcBvBRkQm0,13202
-radia-1.3.4.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-radia-1.3.4.dist-info/top_level.txt,sha256=J-z0poNcsv31IHB413--iOY8LoHBKiTHeybHX3abokI,7
-radia-1.3.4.dist-info/RECORD,,
+radia-1.3.6.dist-info/licenses/LICENSE,sha256=yaWxyzG9DpJ44dDNdGni4nukwiZ8pU-r_aW-1tYNAjk,4374
+radia-1.3.6.dist-info/METADATA,sha256=sLwoSNnU0g8jAa_efG0c3xDb2Z_YLIqgGXo55OIiblo,15518
+radia-1.3.6.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+radia-1.3.6.dist-info/top_level.txt,sha256=J-z0poNcsv31IHB413--iOY8LoHBKiTHeybHX3abokI,7
+radia-1.3.6.dist-info/RECORD,,