cad-to-dagmc 0.7.7__tar.gz → 0.8.1__tar.gz

{cad_to_dagmc-0.7.7 → cad_to_dagmc-0.8.1}/.github/workflows/ci_with_conda_install.yml

@@ -72,3 +72,6 @@ jobs:
           python examples/surface_mesh/single_stp_file.py
           python examples/unstrucutred_volume_mesh/curved_cadquery_object_to_dagmc_volume_mesh.py
           python examples/unstrucutred_volume_mesh/simulate_unstrucutred_volume_mesh_with_openmc.py
+          python examples/unstrucutred_volume_mesh/different_resolution_meshes.py
+          python examples/surface_mesh/different_resolution_meshes.py
+          python examples/surface_mesh/cadquery_assembly_with_scaled_geometry.py

{cad_to_dagmc-0.7.7 → cad_to_dagmc-0.8.1}/.github/workflows/ci_with_pip_install.yml

@@ -58,4 +58,6 @@ jobs:
           python examples/surface_mesh/from_gmsh_mesh_file.py
           python examples/unstrucutred_volume_mesh/curved_cadquery_object_to_dagmc_volume_mesh.py
           python examples/unstrucutred_volume_mesh/simulate_unstrucutred_volume_mesh_with_openmc.py
+          python examples/unstrucutred_volume_mesh/different_resolution_meshes.py
+          python examples/surface_mesh/different_resolution_meshes.py
           python examples/surface_mesh/cadquery_assembly_with_scaled_geometry.py

{cad_to_dagmc-0.7.7 → cad_to_dagmc-0.8.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: cad_to_dagmc
-Version: 0.7.7
+Version: 0.8.1
 Summary: Converts CAD files to a DAGMC h5m file
 Author-email: Jonathan Shimwell <mail@jshimwell.com>
 Project-URL: Homepage, https://github.com/fusion-energy/cad_to_dagmc
@@ -35,7 +35,7 @@ Requires-Dist: vtk; extra == "tests"
 [![PyPI](https://img.shields.io/pypi/v/cad_to_dagmc?color=brightgreen&label=pypi&logo=grebrightgreenen&logoColor=green)](https://pypi.org/project/cad_to_dagmc/)
 
 
-A minimal package that converts CAD geometry to [DAGMC](https://github.com/svalinn/DAGMC/) compatible meshes
+A minimal package that converts CAD geometry to [DAGMC](https://github.com/svalinn/DAGMC/) (h5m) files, [unstructured mesh](https://docs.openmc.org/en/latest/pythonapi/generated/openmc.UnstructuredMesh.html) files (vtk) and Gmsh (msh) files ready for use in neutronics simulations.
 
 cad-to-dagmc can create DAGMC compatible:
 - surface meshes / faceted geometry / triangular meshes
@@ -43,23 +43,23 @@ cad-to-dagmc can create DAGMC compatible:
 
 cad-to-dagmc can convert the following in to DAGMC compatible meshes:
 - STEP files
-- CadQuery objects
+- CadQuery objects (in memory)
 - Gmsh meshes
 
-Options include
-- Ability to scale the output mesh (e.g from meters to cm)
-- Ability to skip imprinting stage
-- Specify meshing parameters for finer or coarser meshes or change the mesh algorithm
-- Ability to specify a material tag for the DAGMC implicit complement
-- Pass objects from CadQuery in memory (without writting to disk)
+Cad-to-dagmc offers a wide range of features including.
+- Geometry scaling with ```scale_factor``` argument
+- Model wide mesh size parameters with ```min_mesh_size``` and ```max_mesh_size``` arguments
+- Volume specific mesh sizing parameters with the ```set_size``` argument
+- Parallel meshing to quickly mesh the geometry using multiple CPU cores
+- Imprint and merging of CAD geometry, or disable with the ```imprint``` argument
+- Add geometry from multiple sources ([STEP](http://www.steptools.com/stds/step/) files, [CadQuery](https://cadquery.readthedocs.io) objects and [Gmsh](https://gmsh.info/) meshes)
+- Ability to tag the DAGMC implicit complement material using the ```implicit_complement_material_tag``` argument
+- Selected different Gmesh mesh algorithms (defaults to 1) using the ```mesh_algorithm``` argument
+- Pass CadQuery objects in memory for fast transfer of geometry using the ```method``` argument
+- Easy to install with [pip](https://pypi.org/project/cad-to-dagmc/) and [Conda/Mamba](https://anaconda.org/conda-forge/cad_to_dagmc)
+- Well tested both with [CI unit tests](https://github.com/fusion-energy/cad_to_dagmc/tree/main/tests), integration tests and the CSG [Model Benchmark Zoo](https://github.com/fusion-energy/model_benchmark_zoo).
+- Compatible with [Paramak](https://github.com/fusion-energy/paramak) geometry for fusion simulations.
 
-cad-to-dagmc aims to produce DAGMC compatible h5m and vtk files from CAD geometry is intended to convert [STEP](http://www.steptools.com/stds/step/) files, [CadQuery](https://cadquery.readthedocs.io) or [Gmsh]([https://cadquery.readthedocs.io](https://gmsh.info/)) meshes objects to a [DAGMC](https://github.com/svalinn/DAGMC/) compatible h5m file or a DAGMC compatible Unstruuctre vtk file.
-
-The resulting DAGMC geometry can then be used for simulations in [OpenMC](https://github.com/openmc-dev/openmc/) or [other supported codes](https://svalinn.github.io/DAGMC/).
-
-This package is tested with [pytest tests](https://github.com/fusion-energy/cad_to_dagmc/tree/main/tests) and also the DAGMC geometry made with this package is compared to simulation carried out with native constructive solid geometry, see [Model Benchmark Zoo](https://github.com/fusion-energy/model_benchmark_zoo) for more details.
-
-Also checkout these other software projects that also create DAGMC geometry [CAD-to-OpenMC](https://github.com/openmsr/CAD_to_OpenMC), [Stellarmesh](https://github.com/Thea-Energy/stellarmesh) and [Coreform Cubit](https://coreform.com/products/coreform-cubit/)
 
 # Installation options
 
@@ -176,3 +176,7 @@ Installing one of the package dependancies (Gmsh) with pip appears to result in
 For examples showing creation of DAGMC h5m files, vtk files and usage within OpenMC transport code see the [examples folder](https://github.com/fusion-energy/cad_to_dagmc/tree/main/examples)
 
 For more examples see the CAD tasks in the [neutronics-workshop](https://github.com/fusion-energy/neutronics-workshop) and [model benchmark zoo](https://github.com/fusion-energy/model_benchmark_zoo)
+
+# Related software
+
+Also checkout these other software projects that also create DAGMC geometry [CAD-to-OpenMC](https://github.com/openmsr/CAD_to_OpenMC), [Stellarmesh](https://github.com/Thea-Energy/stellarmesh) and [Coreform Cubit](https://coreform.com/products/coreform-cubit/).

{cad_to_dagmc-0.7.7 → cad_to_dagmc-0.8.1}/README.md

@@ -12,7 +12,7 @@
 [![PyPI](https://img.shields.io/pypi/v/cad_to_dagmc?color=brightgreen&label=pypi&logo=grebrightgreenen&logoColor=green)](https://pypi.org/project/cad_to_dagmc/)
 
 
-A minimal package that converts CAD geometry to [DAGMC](https://github.com/svalinn/DAGMC/) compatible meshes
+A minimal package that converts CAD geometry to [DAGMC](https://github.com/svalinn/DAGMC/) (h5m) files, [unstructured mesh](https://docs.openmc.org/en/latest/pythonapi/generated/openmc.UnstructuredMesh.html) files (vtk) and Gmsh (msh) files ready for use in neutronics simulations.
 
 cad-to-dagmc can create DAGMC compatible:
 - surface meshes / faceted geometry / triangular meshes
@@ -20,23 +20,23 @@ cad-to-dagmc can create DAGMC compatible:
 
 cad-to-dagmc can convert the following in to DAGMC compatible meshes:
 - STEP files
-- CadQuery objects
+- CadQuery objects (in memory)
 - Gmsh meshes
 
-Options include
-- Ability to scale the output mesh (e.g from meters to cm)
-- Ability to skip imprinting stage
-- Specify meshing parameters for finer or coarser meshes or change the mesh algorithm
-- Ability to specify a material tag for the DAGMC implicit complement
-- Pass objects from CadQuery in memory (without writting to disk)
+Cad-to-dagmc offers a wide range of features including.
+- Geometry scaling with ```scale_factor``` argument
+- Model wide mesh size parameters with ```min_mesh_size``` and ```max_mesh_size``` arguments
+- Volume specific mesh sizing parameters with the ```set_size``` argument
+- Parallel meshing to quickly mesh the geometry using multiple CPU cores
+- Imprint and merging of CAD geometry, or disable with the ```imprint``` argument
+- Add geometry from multiple sources ([STEP](http://www.steptools.com/stds/step/) files, [CadQuery](https://cadquery.readthedocs.io) objects and [Gmsh](https://gmsh.info/) meshes)
+- Ability to tag the DAGMC implicit complement material using the ```implicit_complement_material_tag``` argument
+- Selected different Gmesh mesh algorithms (defaults to 1) using the ```mesh_algorithm``` argument
+- Pass CadQuery objects in memory for fast transfer of geometry using the ```method``` argument
+- Easy to install with [pip](https://pypi.org/project/cad-to-dagmc/) and [Conda/Mamba](https://anaconda.org/conda-forge/cad_to_dagmc)
+- Well tested both with [CI unit tests](https://github.com/fusion-energy/cad_to_dagmc/tree/main/tests), integration tests and the CSG [Model Benchmark Zoo](https://github.com/fusion-energy/model_benchmark_zoo).
+- Compatible with [Paramak](https://github.com/fusion-energy/paramak) geometry for fusion simulations.
 
-cad-to-dagmc aims to produce DAGMC compatible h5m and vtk files from CAD geometry is intended to convert [STEP](http://www.steptools.com/stds/step/) files, [CadQuery](https://cadquery.readthedocs.io) or [Gmsh]([https://cadquery.readthedocs.io](https://gmsh.info/)) meshes objects to a [DAGMC](https://github.com/svalinn/DAGMC/) compatible h5m file or a DAGMC compatible Unstruuctre vtk file.
-
-The resulting DAGMC geometry can then be used for simulations in [OpenMC](https://github.com/openmc-dev/openmc/) or [other supported codes](https://svalinn.github.io/DAGMC/).
-
-This package is tested with [pytest tests](https://github.com/fusion-energy/cad_to_dagmc/tree/main/tests) and also the DAGMC geometry made with this package is compared to simulation carried out with native constructive solid geometry, see [Model Benchmark Zoo](https://github.com/fusion-energy/model_benchmark_zoo) for more details.
-
-Also checkout these other software projects that also create DAGMC geometry [CAD-to-OpenMC](https://github.com/openmsr/CAD_to_OpenMC), [Stellarmesh](https://github.com/Thea-Energy/stellarmesh) and [Coreform Cubit](https://coreform.com/products/coreform-cubit/)
 
 # Installation options
 
@@ -153,3 +153,7 @@ Installing one of the package dependancies (Gmsh) with pip appears to result in
 For examples showing creation of DAGMC h5m files, vtk files and usage within OpenMC transport code see the [examples folder](https://github.com/fusion-energy/cad_to_dagmc/tree/main/examples)
 
 For more examples see the CAD tasks in the [neutronics-workshop](https://github.com/fusion-energy/neutronics-workshop) and [model benchmark zoo](https://github.com/fusion-energy/model_benchmark_zoo)
+
+# Related software
+
+Also checkout these other software projects that also create DAGMC geometry [CAD-to-OpenMC](https://github.com/openmsr/CAD_to_OpenMC), [Stellarmesh](https://github.com/Thea-Energy/stellarmesh) and [Coreform Cubit](https://coreform.com/products/coreform-cubit/).

cad_to_dagmc-0.8.1/examples/surface_mesh/different_resolution_meshes.py

@@ -0,0 +1,44 @@
+# This example makes 3 CAD boxes
+# Meshes the 3 volumes with different resolutions
+# exports the mesh to a DAGMC h5m file and GMsh msh file
+import cadquery as cq
+from cad_to_dagmc import CadToDagmc
+
+box_set_size_course_mesh = cq.Workplane().box(1, 1, 2)
+box_set_size_fine_mesh = cq.Workplane().moveTo(1, 0.5).box(1, 1, 1.5)
+box_set_global_mesh = cq.Workplane().moveTo(2, 1).box(1, 1, 1)
+
+assembly = cq.Assembly()
+assembly.add(box_set_size_course_mesh, color=cq.Color(0, 0, 1))
+assembly.add(box_set_size_fine_mesh, color=cq.Color(0, 1, 0))
+assembly.add(box_set_global_mesh, color=cq.Color(1, 0, 0))
+
+assembly.export("different_resolution_meshes.step")
+
+# uncomment to see the assembly in a pop up vtk viewer
+# from cadquery import vis
+# vis.show(assembly)
+
+model = CadToDagmc()
+model.add_cadquery_object(assembly, material_tags=["mat1", "mat2", "mat3"])
+
+model.export_dagmc_h5m_file(
+    filename="different_resolution_meshes.h5m",
+    min_mesh_size=0.01,
+    max_mesh_size=10,
+    set_size={
+        1: 0.9,
+        2: 0.1,
+    },  # not volume 3 is not specified in the set_size so it uses only the min max mesh sizes
+)
+
+model.export_gmsh_mesh_file(
+    filename="different_resolution_meshes.msh",
+    dimensions=2,
+    min_mesh_size=0.01,
+    max_mesh_size=10,
+    set_size={
+        1: 0.9,
+        2: 0.1,
+    },  # not volume 3 is not specified in the set_size so it uses only the min max mesh sizes
+)

cad_to_dagmc-0.8.1/examples/unstrucutred_volume_mesh/different_resolution_meshes.py

@@ -0,0 +1,44 @@
+# This example makes 3 CAD boxes
+# Meshes the 3 volumes with different resolutions
+# exports the mesh to a DAGMC unstructured VTK file and Gmsh msh file
+import cadquery as cq
+from cad_to_dagmc import CadToDagmc
+
+box_set_size_course_mesh = cq.Workplane().box(1, 1, 2)
+box_set_size_fine_mesh = cq.Workplane().moveTo(1, 0.5).box(1, 1, 1.5)
+box_set_global_mesh = cq.Workplane().moveTo(2, 1).box(1, 1, 1)
+
+assembly = cq.Assembly()
+assembly.add(box_set_size_course_mesh, color=cq.Color(0, 0, 1))
+assembly.add(box_set_size_fine_mesh, color=cq.Color(0, 1, 0))
+assembly.add(box_set_global_mesh, color=cq.Color(1, 0, 0))
+
+assembly.export("different_resolution_meshes.step")
+
+# uncomment to see the assembly in a pop up vtk viewer
+# from cadquery import vis
+# vis.show(assembly)
+
+model = CadToDagmc()
+model.add_cadquery_object(assembly, material_tags=["mat1", "mat2", "mat3"])
+
+model.export_gmsh_mesh_file(
+    filename="different_resolution_meshes.msh",
+    dimensions=3,
+    min_mesh_size=0.01,
+    max_mesh_size=10,
+    set_size={
+        1: 0.9,
+        2: 0.1,
+    },  # not volume 3 is not specified in the set_size so it uses only the min max mesh sizes
+)
+
+model.export_unstructured_mesh_file(
+    filename="different_resolution_meshes.vtk",
+    min_mesh_size=0.01,
+    max_mesh_size=10,
+    set_size={
+        1: 0.9,
+        2: 0.1,
+    },  # not volume 3 is not specified in the set_size so it uses only the min max mesh sizes
+)

{cad_to_dagmc-0.7.7 → cad_to_dagmc-0.8.1}/src/_version.py

@@ -12,5 +12,5 @@ __version__: str
 __version_tuple__: VERSION_TUPLE
 version_tuple: VERSION_TUPLE
 
-__version__ = version = '0.7.7'
-__version_tuple__ = version_tuple = (0, 7, 7)
+__version__ = version = '0.8.1'
+__version_tuple__ = version_tuple = (0, 8, 1)

{cad_to_dagmc-0.7.7 → cad_to_dagmc-0.8.1}/src/cad_to_dagmc/core.py

@@ -235,10 +235,11 @@ def init_gmsh():
 
 def _mesh_brep(
     gmsh,
-    min_mesh_size: float = 1,
-    max_mesh_size: float = 10,
+    min_mesh_size: float | None = None,
+    max_mesh_size: float | None = None,
     mesh_algorithm: int = 1,
     dimensions: int = 2,
+    set_size: dict[int, float] | None = None,
 ):
     """Creates a conformal surface meshes of the volumes in a Brep file using Gmsh.
 
@@ -252,15 +253,49 @@ def _mesh_brep(
             gmsh.option.setNumber("Mesh.Algorithm", mesh_algorithm)
         dimensions: The number of dimensions, 2 for a surface mesh 3 for a
             volume mesh. Passed to gmsh.model.mesh.generate()
+        set_size: a dictionary of volume ids (int) and target mesh sizes
+            (floats) to set for each volume, passed to gmsh.model.mesh.setSize.
 
     Returns:
         The resulting gmsh object and volumes
     """
+    if min_mesh_size and max_mesh_size:
+        if min_mesh_size > max_mesh_size:
+            raise ValueError(
+                f"min_mesh_size must be less than or equal to max_mesh_size. Currently min_mesh_size is set to {min_mesh_size} and max_mesh_size is set to {max_mesh_size}"
+            )
+
+    if min_mesh_size:
+        gmsh.option.setNumber("Mesh.MeshSizeMin", min_mesh_size)
+
+    if max_mesh_size:
+        gmsh.option.setNumber("Mesh.MeshSizeMax", max_mesh_size)
 
     gmsh.option.setNumber("Mesh.Algorithm", mesh_algorithm)
-    gmsh.option.setNumber("Mesh.MeshSizeMin", min_mesh_size)
-    gmsh.option.setNumber("Mesh.MeshSizeMax", max_mesh_size)
     gmsh.option.setNumber("General.NumThreads", 0)  # Use all available cores
+
+    if set_size:
+        volumes = gmsh.model.getEntities(3)
+        available_volumes = [volume[1] for volume in volumes]
+        print("volumes", volumes)
+        for volume_id, size in set_size.items():
+            if volume_id in available_volumes:
+                size = set_size[volume_id]
+                if isinstance(size, dict):
+                    # TODO face specific mesh sizes
+                    pass
+                else:
+                    boundaries = gmsh.model.getBoundary(
+                        [(3, volume_id)], recursive=True
+                    )  # dim must be set to 3
+                    print("boundaries", boundaries)
+                    gmsh.model.mesh.setSize(boundaries, size)
+                    print(f"Set size of {size} for volume {volume_id}")
+            else:
+                raise ValueError(
+                    f"volume ID of {volume_id} set in set_sizes but not found in available volumes {volumes}"
+                )
+
     gmsh.model.mesh.generate(dimensions)
 
     return gmsh
@@ -468,6 +503,7 @@ class CadToDagmc:
             cq.assembly.Assembly | cq.occ_impl.shapes.Compound | cq.occ_impl.shapes.Solid
         ),
         material_tags: list[str] | None,
+        scale_factor: float = 1.0,
     ) -> int:
         """Loads the parts from CadQuery object into the model.
 
@@ -479,6 +515,10 @@ class CadToDagmc:
                 same order as the volumes in the geometry added (STP file and
                 CadQuery objects) and match the material tags used in the
                 neutronics code (e.g. OpenMC).
+            scale_factor: a scaling factor to apply to the geometry that can be
+                used to increase the size or decrease the size of the geometry.
+                Useful when converting the geometry to cm for use in neutronics
+                simulations.
 
         Returns:
             int: number of volumes in the stp file.
@@ -492,12 +532,17 @@ class CadToDagmc:
         else:
             iterable_solids = cadquery_object.val().Solids()
 
-        _check_material_tags(material_tags, iterable_solids)
+        if scale_factor == 1.0:
+            scaled_iterable_solids = iterable_solids
+        else:
+            scaled_iterable_solids = [part.scale(scale_factor) for part in iterable_solids]
+
+        _check_material_tags(material_tags, scaled_iterable_solids)
         if material_tags:
             self.material_tags = self.material_tags + material_tags
-        self.parts = self.parts + iterable_solids
+        self.parts = self.parts + scaled_iterable_solids
 
-        return len(iterable_solids)
+        return len(scaled_iterable_solids)
 
     def export_unstructured_mesh_file(
         self,
@@ -508,6 +553,7 @@ class CadToDagmc:
         method: str = "file",
         scale_factor: float = 1.0,
         imprint: bool = True,
+        set_size: dict[int, float] | None = None,
     ):
         """
         Exports an unstructured mesh file in VTK format for use with openmc.UnstructuredMesh.
@@ -539,6 +585,9 @@ class CadToDagmc:
                 normally needed to ensure the geometry is meshed correctly. However if
                 you know your geometry does not need imprinting you can set this to False
                 and this can save time.
+            set_size: a dictionary of volume ids (int) and target mesh sizes
+                (floats) to set for each volume, passed to gmsh.model.mesh.setSize.
+
 
         Returns:
         --------
@@ -570,6 +619,7 @@ class CadToDagmc:
             max_mesh_size=max_mesh_size,
             mesh_algorithm=mesh_algorithm,
             dimensions=3,
+            set_size=set_size,
         )
 
         # makes the folder if it does not exist
@@ -589,13 +639,14 @@ class CadToDagmc:
     def export_gmsh_mesh_file(
         self,
         filename: str = "mesh.msh",
-        min_mesh_size: float = 1,
-        max_mesh_size: float = 5,
+        min_mesh_size: float | None = None,
+        max_mesh_size: float | None = None,
         mesh_algorithm: int = 1,
         dimensions: int = 2,
         method: str = "file",
         scale_factor: float = 1.0,
         imprint: bool = True,
+        set_size: dict[int, float] | None = None,
     ):
         """Saves a GMesh msh file of the geometry in either 2D surface mesh or
         3D volume mesh.
@@ -622,6 +673,8 @@ class CadToDagmc:
                 normally needed to ensure the geometry is meshed correctly. However if
                 you know your geometry does not need imprinting you can set this to False
                 and this can save time.
+            set_size: a dictionary of volume ids (int) and target mesh sizes
+                (floats) to set for each volume, passed to gmsh.model.mesh.setSize.
         """
 
         assembly = cq.Assembly()
@@ -643,6 +696,7 @@ class CadToDagmc:
             max_mesh_size=max_mesh_size,
             mesh_algorithm=mesh_algorithm,
             dimensions=dimensions,
+            set_size=set_size,
         )
 
         # makes the folder if it does not exist
@@ -662,13 +716,14 @@ class CadToDagmc:
     def export_dagmc_h5m_file(
         self,
         filename: str = "dagmc.h5m",
-        min_mesh_size: float = 1,
-        max_mesh_size: float = 5,
+        min_mesh_size: float | None = None,
+        max_mesh_size: float | None = None,
         mesh_algorithm: int = 1,
         implicit_complement_material_tag: str | None = None,
         method: str = "file",
         scale_factor: float = 1.0,
         imprint: bool = True,
+        set_size: dict[int, float] | None = None,
     ) -> str:
         """Saves a DAGMC h5m file of the geometry
 
@@ -695,6 +750,8 @@ class CadToDagmc:
                 normally needed to ensure the geometry is meshed correctly. However if
                 you know your geometry does not need imprinting you can set this to False
                 and this can save time.
+            set_size: a dictionary of volume ids (int) and target mesh sizes
+                (floats) to set for each volume, passed to gmsh.model.mesh.setSize.
 
         Returns:
             str: the DAGMC filename saved
@@ -739,6 +796,7 @@ class CadToDagmc:
             min_mesh_size=min_mesh_size,
             max_mesh_size=max_mesh_size,
             mesh_algorithm=mesh_algorithm,
+            set_size=set_size,
         )
 
         dims_and_vol_ids = volumes

{cad_to_dagmc-0.7.7 → cad_to_dagmc-0.8.1}/src/cad_to_dagmc.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: cad_to_dagmc
-Version: 0.7.7
+Version: 0.8.1
 Summary: Converts CAD files to a DAGMC h5m file
 Author-email: Jonathan Shimwell <mail@jshimwell.com>
 Project-URL: Homepage, https://github.com/fusion-energy/cad_to_dagmc
@@ -35,7 +35,7 @@ Requires-Dist: vtk; extra == "tests"
 [![PyPI](https://img.shields.io/pypi/v/cad_to_dagmc?color=brightgreen&label=pypi&logo=grebrightgreenen&logoColor=green)](https://pypi.org/project/cad_to_dagmc/)
 
 
-A minimal package that converts CAD geometry to [DAGMC](https://github.com/svalinn/DAGMC/) compatible meshes
+A minimal package that converts CAD geometry to [DAGMC](https://github.com/svalinn/DAGMC/) (h5m) files, [unstructured mesh](https://docs.openmc.org/en/latest/pythonapi/generated/openmc.UnstructuredMesh.html) files (vtk) and Gmsh (msh) files ready for use in neutronics simulations.
 
 cad-to-dagmc can create DAGMC compatible:
 - surface meshes / faceted geometry / triangular meshes
@@ -43,23 +43,23 @@ cad-to-dagmc can create DAGMC compatible:
 
 cad-to-dagmc can convert the following in to DAGMC compatible meshes:
 - STEP files
-- CadQuery objects
+- CadQuery objects (in memory)
 - Gmsh meshes
 
-Options include
-- Ability to scale the output mesh (e.g from meters to cm)
-- Ability to skip imprinting stage
-- Specify meshing parameters for finer or coarser meshes or change the mesh algorithm
-- Ability to specify a material tag for the DAGMC implicit complement
-- Pass objects from CadQuery in memory (without writting to disk)
+Cad-to-dagmc offers a wide range of features including.
+- Geometry scaling with ```scale_factor``` argument
+- Model wide mesh size parameters with ```min_mesh_size``` and ```max_mesh_size``` arguments
+- Volume specific mesh sizing parameters with the ```set_size``` argument
+- Parallel meshing to quickly mesh the geometry using multiple CPU cores
+- Imprint and merging of CAD geometry, or disable with the ```imprint``` argument
+- Add geometry from multiple sources ([STEP](http://www.steptools.com/stds/step/) files, [CadQuery](https://cadquery.readthedocs.io) objects and [Gmsh](https://gmsh.info/) meshes)
+- Ability to tag the DAGMC implicit complement material using the ```implicit_complement_material_tag``` argument
+- Selected different Gmesh mesh algorithms (defaults to 1) using the ```mesh_algorithm``` argument
+- Pass CadQuery objects in memory for fast transfer of geometry using the ```method``` argument
+- Easy to install with [pip](https://pypi.org/project/cad-to-dagmc/) and [Conda/Mamba](https://anaconda.org/conda-forge/cad_to_dagmc)
+- Well tested both with [CI unit tests](https://github.com/fusion-energy/cad_to_dagmc/tree/main/tests), integration tests and the CSG [Model Benchmark Zoo](https://github.com/fusion-energy/model_benchmark_zoo).
+- Compatible with [Paramak](https://github.com/fusion-energy/paramak) geometry for fusion simulations.
 
-cad-to-dagmc aims to produce DAGMC compatible h5m and vtk files from CAD geometry is intended to convert [STEP](http://www.steptools.com/stds/step/) files, [CadQuery](https://cadquery.readthedocs.io) or [Gmsh]([https://cadquery.readthedocs.io](https://gmsh.info/)) meshes objects to a [DAGMC](https://github.com/svalinn/DAGMC/) compatible h5m file or a DAGMC compatible Unstruuctre vtk file.
-
-The resulting DAGMC geometry can then be used for simulations in [OpenMC](https://github.com/openmc-dev/openmc/) or [other supported codes](https://svalinn.github.io/DAGMC/).
-
-This package is tested with [pytest tests](https://github.com/fusion-energy/cad_to_dagmc/tree/main/tests) and also the DAGMC geometry made with this package is compared to simulation carried out with native constructive solid geometry, see [Model Benchmark Zoo](https://github.com/fusion-energy/model_benchmark_zoo) for more details.
-
-Also checkout these other software projects that also create DAGMC geometry [CAD-to-OpenMC](https://github.com/openmsr/CAD_to_OpenMC), [Stellarmesh](https://github.com/Thea-Energy/stellarmesh) and [Coreform Cubit](https://coreform.com/products/coreform-cubit/)
 
 # Installation options
 
@@ -176,3 +176,7 @@ Installing one of the package dependancies (Gmsh) with pip appears to result in
 For examples showing creation of DAGMC h5m files, vtk files and usage within OpenMC transport code see the [examples folder](https://github.com/fusion-energy/cad_to_dagmc/tree/main/examples)
 
 For more examples see the CAD tasks in the [neutronics-workshop](https://github.com/fusion-energy/neutronics-workshop) and [model benchmark zoo](https://github.com/fusion-energy/model_benchmark_zoo)
+
+# Related software
+
+Also checkout these other software projects that also create DAGMC geometry [CAD-to-OpenMC](https://github.com/openmsr/CAD_to_OpenMC), [Stellarmesh](https://github.com/Thea-Energy/stellarmesh) and [Coreform Cubit](https://coreform.com/products/coreform-cubit/).

{cad_to_dagmc-0.7.7 → cad_to_dagmc-0.8.1}/src/cad_to_dagmc.egg-info/SOURCES.txt

@@ -14,6 +14,7 @@ examples/surface_mesh/cadquery_compound.py
 examples/surface_mesh/cadquery_object_and_stp_file.py
 examples/surface_mesh/cadquery_text.py
 examples/surface_mesh/curved_cadquery_object_to_dagmc_surface_mesh.py
+examples/surface_mesh/different_resolution_meshes.py
 examples/surface_mesh/from_gmsh_mesh_file.py
 examples/surface_mesh/multiple_cadquery_objects.py
 examples/surface_mesh/multiple_stp_files.py
@@ -21,6 +22,7 @@ examples/surface_mesh/single_cadquery_object.py
 examples/surface_mesh/single_stp_file.py
 examples/surface_mesh/single_stp_file_multiple_volumes.py
 examples/unstrucutred_volume_mesh/curved_cadquery_object_to_dagmc_volume_mesh.py
+examples/unstrucutred_volume_mesh/different_resolution_meshes.py
 examples/unstrucutred_volume_mesh/simulate_unstrucutred_volume_mesh_with_openmc.py
 src/_version.py
 src/cad_to_dagmc/__init__.py

{cad_to_dagmc-0.7.7 → cad_to_dagmc-0.8.1}/tests/test_python_api.py

@@ -1,16 +1,16 @@
 import os
-from pathlib import Path
-import pytest
-import cadquery as cq
-from cad_to_dagmc import CadToDagmc
 import warnings
-from cad_to_dagmc.core import _check_material_tags
-
 from pathlib import Path
 
+import cadquery as cq
+import gmsh
 import pymoab as mb
+import pytest
 from pymoab import core, types
 
+from cad_to_dagmc import CadToDagmc
+from cad_to_dagmc.core import _check_material_tags
+
 
 def get_volumes_and_materials_from_h5m(filename: str) -> dict:
     """Reads in a DAGMC h5m file and uses PyMoab to find the volume ids with
@@ -219,3 +219,117 @@ def test_check_material_tags_too_long():
         assert issubclass(w[-1].category, UserWarning)
         assert "Material tag" in str(w[-1].message)
         assert "a" * 29 in str(w[-1].message)
+
+
+@pytest.mark.parametrize(
+    "scale_factor, expected_width",
+    [
+        (1, 10.0),
+        (2, 20.0),
+        (10, 100.0),
+    ],
+)
+def test_scaling_factor_when_adding_stp(scale_factor, expected_width):
+
+    c2d = CadToDagmc()
+    c2d.add_stp_file("tests/single_cube.stp", scale_factor=scale_factor)
+    c2d.export_gmsh_mesh_file(f"st_test_scaling_factor_{scale_factor}.msh")
+
+    gmsh.initialize()
+    gmsh.open(f"st_test_scaling_factor_{scale_factor}.msh")
+    _, node_coords, _ = gmsh.model.mesh.getNodes()
+
+    # Reshape the node coordinates into a 2D array
+    node_coords = node_coords.reshape(-1, 3)
+
+    # Calculate the bounding box
+    min_coords = node_coords.min(axis=0)
+    max_coords = node_coords.max(axis=0)
+
+    width_x = max_coords[0] - min_coords[0]
+    width_y = max_coords[1] - min_coords[1]
+    width_z = max_coords[2] - min_coords[2]
+
+    gmsh.finalize()
+
+    assert width_x == expected_width
+    assert width_y == expected_width
+    assert width_z == expected_width
+
+
+@pytest.mark.parametrize(
+    "scale_factor, expected_width",
+    [
+        (1, 10.0),
+        (2, 20.0),
+        (10, 100.0),
+    ],
+)
+def test_scaling_factor_when_adding_cq_object(scale_factor, expected_width):
+
+    box = cq.Workplane("XY").box(10, 10, 10)
+    c2d = CadToDagmc()
+    c2d.add_cadquery_object(box, scale_factor=scale_factor, material_tags=["mat1"])
+    c2d.export_gmsh_mesh_file(f"cq_test_scaling_factor_{scale_factor}.msh")
+
+    gmsh.initialize()
+    gmsh.open(f"cq_test_scaling_factor_{scale_factor}.msh")
+    _, node_coords, _ = gmsh.model.mesh.getNodes()
+
+    # Reshape the node coordinates into a 2D array
+    node_coords = node_coords.reshape(-1, 3)
+
+    # Calculate the bounding box
+    min_coords = node_coords.min(axis=0)
+    max_coords = node_coords.max(axis=0)
+
+    width_x = max_coords[0] - min_coords[0]
+    width_y = max_coords[1] - min_coords[1]
+    width_z = max_coords[2] - min_coords[2]
+
+    gmsh.finalize()
+
+    assert width_x == expected_width
+    assert width_y == expected_width
+    assert width_z == expected_width
+
+
+@pytest.mark.parametrize(
+    "scale_factor, expected_x_width, expected_y_width, expected_z_width",
+    [
+        (1, 20.0, 10.0, 10.0),
+        (2, 40.0, 20.0, 20.0),
+        (10, 200.0, 100.0, 100.0),
+    ],
+)
+def test_two_box_scaling_factor_when_adding_cq_object(
+    scale_factor, expected_x_width, expected_y_width, expected_z_width
+):
+
+    box = cq.Workplane("XY").box(10, 10, 10)
+    box2 = cq.Workplane("XY").moveTo(10, 0).box(10, 10, 10)
+    c2d = CadToDagmc()
+    c2d.add_cadquery_object(box, scale_factor=scale_factor, material_tags=["mat1"])
+    c2d.add_cadquery_object(box2, scale_factor=scale_factor, material_tags=["mat1"])
+    c2d.export_gmsh_mesh_file(f"cq_test_2_box_scaling_factor_{scale_factor}.msh")
+
+    gmsh.initialize()
+    gmsh.open(f"cq_test_2_box_scaling_factor_{scale_factor}.msh")
+    _, node_coords, _ = gmsh.model.mesh.getNodes()
+
+    # Reshape the node coordinates into a numpy 2D array
+    node_coords = node_coords.reshape(-1, 3)
+
+    # Calculate the bounding box
+    min_coords = node_coords.min(axis=0)
+    max_coords = node_coords.max(axis=0)
+
+    width_x = max_coords[0] - min_coords[0]
+    width_y = max_coords[1] - min_coords[1]
+    width_z = max_coords[2] - min_coords[2]
+
+    gmsh.finalize()
+
+    assert width_x == expected_x_width
+    assert width_y == expected_y_width
+    assert width_z == expected_z_width