radia 1.3.3__py3-none-any.whl → 1.3.5__py3-none-any.whl

radia/radia_vtk_export.py

@@ -0,0 +1,221 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+"""
+Radia VTK Export Utilities
+
+Functions for exporting Radia geometry to VTK format for visualization
+in ParaView and other VTK-compatible tools.
+"""
+
+import csv
+from itertools import accumulate
+
+def _get_radia_length_unit():
+	"""
+	Get current Radia length unit by querying rad.FldUnits().
+
+	Returns:
+		tuple: (unit_name, scale_to_meters)
+		       ('mm', 0.001) if Radia is using millimeters
+		       ('m', 1.0) if Radia is using meters
+
+	Raises:
+		ValueError: If length unit cannot be determined
+	"""
+	import radia as rad
+
+	units_str = rad.FldUnits()
+
+	if 'Length:  mm' in units_str:
+		return ('mm', 0.001)
+	elif 'Length:  m' in units_str:
+		return ('m', 1.0)
+	else:
+		raise ValueError(f"Cannot determine Radia length unit from: {units_str}")
+
+def chunks(lst, n):
+	"""
+	Yield successive n-sized chunks from a list.
+
+	Args:
+		lst: List to be chunked
+		n: Chunk size
+
+	Yields:
+		Chunks of size n from the input list
+	"""
+	for i in range(0, len(lst), n):
+		yield lst[i:i + n]
+
+def exportGeometryToVTK(obj, fileName='radia_Geometry'):
+	"""
+	Export Radia object geometry to VTK Legacy format file.
+
+	Writes the geometry of a Radia object to a .vtk file for visualization
+	in ParaView. The format is VTK Legacy (ASCII), consisting of polygons only.
+
+	Args:
+		obj: Radia object ID (integer)
+		fileName: Output filename without extension (default: 'radia_Geometry')
+
+	Output:
+		Creates fileName.vtk in the current directory
+
+	Example:
+		>>> import radia as rad
+		>>> from radia_vtk_export import exportGeometryToVTK
+		>>> mag = rad.ObjRecMag([0,0,0], [10,10,10], [0,0,1])
+		>>> exportGeometryToVTK(mag, 'my_magnet')
+	"""
+	import radia as rad
+
+	vtkData = rad.ObjDrwVTK(obj, 'Axes->False')
+
+	lengths = vtkData['polygons']['lengths']
+	nPoly = len(lengths)
+	offsets = list(accumulate(lengths))
+	offsets.insert(0, 0) # prepend list with a zero
+	points = vtkData['polygons']['vertices']
+	nPnts = int(len(points)/3)
+
+	# Get Radia's current length unit and convert to meters for VTK standard
+	unit_name, scale_to_meters = _get_radia_length_unit()
+	points = [round(num * scale_to_meters, 8) for num in points]
+
+	# define colours array
+	colors = vtkData['polygons']['colors']
+
+	# pre-process the output lists to have chunkLength items per line
+	chunkLength = 9 # this writes 9 numbers per line
+	points = list(chunks(points, chunkLength))
+	colors = list(chunks(colors, chunkLength))
+
+	# write the data to file (VTK Legacy format - most compatible)
+	with open(fileName + ".vtk", "w", newline="") as f:
+		f.write('# vtk DataFile Version 3.0\n')
+		f.write('vtk output\nASCII\nDATASET POLYDATA\n')
+		f.write('POINTS ' + str(nPnts) + ' float\n')
+		writer = csv.writer(f, delimiter=" ")
+		writer.writerows(points)
+		f.write('\n')
+
+		# POLYGONS in classic format (most compatible with all ParaView versions)
+		# Format: nPoly totalSize
+		# Each line: nVertices v1 v2 v3 ...
+		total_size = sum(lengths) + nPoly  # sum of (nVertices + nVertices) for each polygon
+		f.write('POLYGONS ' + str(nPoly) + ' ' + str(total_size) + '\n')
+		for i in range(nPoly):
+			n_vertices = lengths[i]
+			start = offsets[i]
+			end = offsets[i+1]
+			f.write(str(n_vertices))
+			for j in range(start, end):
+				f.write(' ' + str(j))
+			f.write('\n')
+
+		f.write('\n')
+		f.write('CELL_DATA ' + str(nPoly) + '\n')
+		f.write('COLOR_SCALARS Radia_colours 3\n')
+		writer.writerows(colors)
+
+	print(f"VTK file exported: {fileName}.vtk")
+	print(f"  Polygons: {nPoly}")
+	print(f"  Points: {nPnts}")
+
+
+def exportFieldToVTK(points, field_data, fileName='field_distribution', field_name='B_field'):
+	"""
+	Export magnetic field distribution to VTK Legacy format file.
+
+	Args:
+		points: List of observation points [[x, y, z], ...] in Radia's current units
+		field_data: List of field vectors [[Bx, By, Bz], ...] in Tesla
+		fileName: Output file name (without .vtk extension)
+		field_name: Name for the vector field in VTK file
+
+	Note:
+		Points are automatically converted to meters (VTK standard) based on
+		Radia's current unit setting (queried via rad.FldUnits()).
+	"""
+	import numpy as np
+
+	points = np.array(points)
+	field_data = np.array(field_data)
+
+	if points.shape[0] != field_data.shape[0]:
+		raise ValueError(f"Points ({points.shape[0]}) and field data ({field_data.shape[0]}) must have same length")
+
+	if points.shape[1] != 3 or field_data.shape[1] != 3:
+		raise ValueError("Points and field data must be Nx3 arrays")
+
+	# Get Radia's current length unit
+	unit_name, scale_to_meters = _get_radia_length_unit()
+
+	# Ensure .vtk extension
+	if not fileName.endswith('.vtk'):
+		fileName += '.vtk'
+
+	with open(fileName, 'w') as f:
+		# Header
+		f.write('# vtk DataFile Version 3.0\n')
+		f.write('Magnetic Field Distribution\n')
+		f.write('ASCII\n')
+		f.write('DATASET POLYDATA\n')
+
+		# Points (convert to meters for VTK standard)
+		f.write(f'POINTS {len(points)} float\n')
+		for pt in points:
+			# Convert to meters using detected unit scale
+			f.write(f'{pt[0]*scale_to_meters} {pt[1]*scale_to_meters} {pt[2]*scale_to_meters}\n')
+
+		# Point data (vector field)
+		f.write(f'\nPOINT_DATA {len(points)}\n')
+		f.write(f'VECTORS {field_name} float\n')
+		for B in field_data:
+			f.write(f'{B[0]} {B[1]} {B[2]}\n')
+
+	print(f"VTK field file exported: {fileName}")
+	print(f"  Points: {len(points)}")
+	print(f"  Field: {field_name}")
+
+
+if __name__ == '__main__':
+	"""
+	Demo: Export a simple Radia geometry to VTK format
+	"""
+	import sys
+	import os
+
+	# Add build directory to path
+	sys.path.insert(0, os.path.join(os.path.dirname(__file__), '..', '..', 'build', 'Release'))
+
+	import radia as rad
+
+	print("=" * 60)
+	print("Radia VTK Export Demo")
+	print("=" * 60)
+
+	# Create a simple test geometry
+	print("\nCreating test geometry...")
+
+	# Rectangular magnet
+	mag = rad.ObjRecMag([0, 0, 0], [30, 30, 10], [0, 0, 1])
+
+	# Cylindrical magnet
+	cyl = rad.ObjCylMag([50, 0, 0], 15, 20, 16, 'z', [0, 0, 1])
+
+	# Container
+	container = rad.ObjCnt([mag, cyl])
+
+	# Export to VTK
+	output_file = 'radia_demo_geometry'
+	print(f"\nExporting geometry to {output_file}.vtk...")
+	exportGeometryToVTK(container, output_file)
+
+	print("\n" + "=" * 60)
+	print("Export complete!")
+	print("\nTo view in ParaView:")
+	print(f"  1. Open ParaView")
+	print(f"  2. File → Open → {output_file}.vtk")
+	print(f"  3. Click 'Apply' in the Properties panel")
+	print("=" * 60)

{radia-1.3.3.dist-info → radia-1.3.5.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: radia
-Version: 1.3.3
+Version: 1.3.5
 Summary: Radia 3D Magnetostatics with NGSolve Integration and OpenMP Parallelization
 Home-page: https://github.com/ksugahar/Radia_NGSolve
 Author: Pascal Elleaume
@@ -182,9 +182,6 @@ See [examples/Radia_to_NGSolve_CoefficientFunction/](examples/Radia_to_NGSolve_C
 - [examples/Radia_to_NGSolve_CoefficientFunction/EXAMPLES_GUIDE.md](examples/Radia_to_NGSolve_CoefficientFunction/EXAMPLES_GUIDE.md) - Detailed usage guide
 - [tests/test_radia_ngsolve.py](tests/test_radia_ngsolve.py) - Integration tests
 
-### Visualization
-- [docs/PYVISTA_VIEWER.md](examples/2024_02_03_振分電磁石/PYVISTA_VIEWER.md) - PyVista viewer guide
-
 ### Development
 - [docs/TAB_CONVERSION_REPORT.md](docs/TAB_CONVERSION_REPORT.md) - Code style conversion
 - [docs/CLAUDE.md](docs/CLAUDE.md) - Development notes

radia-1.3.5.dist-info/RECORD

@@ -0,0 +1,18 @@
+radia/__init__.py,sha256=edxD1jM_-2LddwoxqtgATa_vTaruXk5WZsntpFJ9qtM,47
+radia/nastran_mesh_import.py,sha256=CSoVhZCXa85lPiTF2hlspE2clBKOD7-_sCp1bxu_IK0,18147
+radia/netgen_mesh_import.py,sha256=UopXk-5bbfj1j9_hyiq8jbjb4SQXnWaeuaC7TDf17wA,19872
+radia/rad_ngsolve.pyd,sha256=s5ZvpTIrjYeV_B3vQv5HrtP4u_RaPhvxuZf5aKv3oC0,566272
+radia/rad_ngsolve_fast.py,sha256=GkC7ruKy3MESHsO0iRSWsrgLU4-DPPgctOi6pPpsWg4,5675
+radia/radia_coil_builder.py,sha256=nQkiAbfhueNvvxUARHdPD0C68ImidHmUQv_q4RsImeY,11253
+radia/radia_field_cached.py,sha256=Bjw3ecNe3u7AAXnLob5m_tjYIY7HwB9DpgFg9a-aP_8,9509
+radia/radia_ngsolve.pyd,sha256=s5ZvpTIrjYeV_B3vQv5HrtP4u_RaPhvxuZf5aKv3oC0,566272
+radia/radia_ngsolve_field.py,sha256=suJr4wacfYFKOkyV-5AQuHWnW5rtUMb0gSSjq8VRSXc,10166
+radia/radia_ngsolve_utils.py,sha256=xZCR9DOIKMwdEjmC28rOXVZiWFY5BQYH2VfopfuVBps,8406
+radia/radia_pyvista_viewer.py,sha256=JS33Mx4azGI7hUX0bzefc6zJfhv6qfRjM3Kl1bE9Mjs,4275
+radia/radia_tetra_mesh.py,sha256=nEpNpOBdJ-0DjS7WUx_bRHskrDKSKH3rSEWr-CwyQRM,21147
+radia/radia_vtk_export.py,sha256=I8Vyyt9tky78Qw1xPru9f0Rii6QEmdEgTFjQtamyooc,6540
+radia-1.3.5.dist-info/licenses/LICENSE,sha256=yaWxyzG9DpJ44dDNdGni4nukwiZ8pU-r_aW-1tYNAjk,4374
+radia-1.3.5.dist-info/METADATA,sha256=LuqDACqJI2J9ucBkZe1UQFfh77ZuYnZvZIwI9JwRpX0,13202
+radia-1.3.5.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+radia-1.3.5.dist-info/top_level.txt,sha256=kNDonE5X3Q2xnmOsWleQnGKQobuFsKM6Px5_Ta1I6x8,6
+radia-1.3.5.dist-info/RECORD,,

{radia-1.3.3.dist-info → radia-1.3.5.dist-info}/licenses/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.5.dist-info/top_level.txt

@@ -0,0 +1 @@
+radia

python/__init__.py

@@ -1,2 +0,0 @@
-# Radia Python package
-__version__ = "1.1.0"