Merge branch 'blitz' into feat/download_stride

Author: b-matteo

doc/changelog.d/1675.added.md

@@ -0,0 +1 @@
+implementation of NURBS curves

doc/changelog.d/1684.added.md

@@ -0,0 +1 @@
+get assigned material

doc/changelog.d/1685.fixed.md

@@ -0,0 +1 @@
+filter set export id to only CoreService based backends

doc/changelog.d/1689.added.md

@@ -0,0 +1 @@
+matrix rotation and translation

doc/changelog.d/1690.fixed.md

@@ -0,0 +1 @@
+cleanup unsupported module

pyproject.toml

@@ -29,6 +29,7 @@ dependencies = [
     "ansys-tools-visualization-interface>=0.2.6,<1",
     "attrs!=24.3.0",
     "beartype>=0.11.0,<0.20",
+    "geomdl>=5,<6",
     "grpcio>=1.35.0,<1.68",
     "grpcio-health-checking>=1.45.0,<1.68",
     "numpy>=1.20.3,<3",
@@ -54,6 +55,7 @@ tests = [
     "ansys-tools-visualization-interface==0.8.1",
     "beartype==0.19.0",
     "docker==7.1.0",
+    "geomdl==5.3.1",
     "grpcio==1.67.1",
     "grpcio-health-checking==1.67.1",
     "numpy==2.2.1",
@@ -82,6 +84,7 @@ doc = [
     "ansys-tools-visualization-interface==0.8.1",
     "beartype==0.19.0",
     "docker==7.1.0",
+    "geomdl==5.3.1",
     "grpcio==1.67.1",
     "grpcio-health-checking==1.67.1",
     "ipyvtklink==0.2.3",

src/ansys/geometry/core/connection/conversions.py

@@ -23,6 +23,8 @@
 
 from typing import TYPE_CHECKING
 
+from pint import Quantity, UndefinedUnitError
+
 from ansys.api.geometry.v0.models_pb2 import (
     Arc as GRPCArc,
     Circle as GRPCCircle,
@@ -32,6 +34,8 @@
     Frame as GRPCFrame,
     Geometries as GRPCGeometries,
     Line as GRPCLine,
+    Material as GRPCMaterial,
+    MaterialProperty as GRPCMaterialProperty,
     Matrix as GRPCMatrix,
     Plane as GRPCPlane,
     Point as GRPCPoint,
@@ -42,12 +46,17 @@
     TrimmedCurve as GRPCTrimmedCurve,
     TrimmedSurface as GRPCTrimmedSurface,
 )
+from ansys.geometry.core.materials.material import (
+    Material,
+    MaterialProperty,
+    MaterialPropertyType,
+)
 from ansys.geometry.core.math.frame import Frame
 from ansys.geometry.core.math.matrix import Matrix44
 from ansys.geometry.core.math.plane import Plane
 from ansys.geometry.core.math.point import Point2D, Point3D
 from ansys.geometry.core.math.vector import UnitVector3D
-from ansys.geometry.core.misc.measurements import DEFAULT_UNITS
+from ansys.geometry.core.misc.measurements import DEFAULT_UNITS, UNITS
 from ansys.geometry.core.shapes.curves.circle import Circle
 from ansys.geometry.core.shapes.curves.curve import Curve
 from ansys.geometry.core.shapes.curves.ellipse import Ellipse
@@ -680,3 +689,37 @@ def line_to_grpc_line(line: Line) -> GRPCLine:
     start = line.origin
     end = line.origin + line.direction
     return GRPCLine(start=point3d_to_grpc_point(start), end=point3d_to_grpc_point(end))
+
+
+def grpc_material_to_material(material: GRPCMaterial) -> Material:
+    """Convert a material gRPC message to a ``Material`` class."""
+    properties = []
+    density = Quantity(0, UNITS.kg / UNITS.m**3)
+    for property in material.material_properties:
+        mp = grpc_material_property_to_material_property(property)
+        properties.append(mp)
+        if mp.type == MaterialPropertyType.DENSITY:
+            density = mp.quantity
+
+    return Material(material.name, density, properties)
+
+
+def grpc_material_property_to_material_property(
+    material_property: GRPCMaterialProperty,
+) -> MaterialProperty:
+    """Convert a material property gRPC message to a ``MaterialProperty`` class."""
+    try:
+        mp_type = MaterialPropertyType.from_id(material_property.id)
+    except ValueError:
+        mp_type = material_property.id
+
+    try:
+        mp_quantity = Quantity(material_property.value, material_property.units)
+    except (
+        UndefinedUnitError,
+        TypeError,
+    ):
+        mp_quantity = material_property.value
+
+    mp = MaterialProperty(mp_type, material_property.display_name, mp_quantity)
+    return mp

src/ansys/geometry/core/designer/body.py

@@ -60,6 +60,7 @@
 from ansys.geometry.core.connection.client import GrpcClient
 from ansys.geometry.core.connection.conversions import (
     frame_to_grpc_frame,
+    grpc_material_to_material,
     plane_to_grpc_plane,
     point3d_to_grpc_point,
     sketch_shapes_to_grpc_geometries,
@@ -254,6 +255,17 @@ def assign_material(self, material: Material) -> None:
         """
         return
 
+    @abstractmethod
+    def get_assigned_material(self) -> Material:
+        """Get the assigned material of the body.
+
+        Returns
+        -------
+        Material
+            Material assigned to the body.
+        """
+        return
+
     @abstractmethod
     def add_midsurface_thickness(self, thickness: Quantity) -> None:
         """Add a mid-surface thickness to a surface body.
@@ -908,6 +920,13 @@ def assign_material(self, material: Material) -> None:  # noqa: D102
             SetAssignedMaterialRequest(id=self._id, material=material.name)
         )
 
+    @property
+    @protect_grpc
+    def get_assigned_material(self) -> Material:  # noqa: D102
+        self._grpc_client.log.debug(f"Retrieving assigned material for body {self.id}.")
+        material_response = self._bodies_stub.GetAssignedMaterial(self._grpc_id)
+        return grpc_material_to_material(material_response)
+
     @protect_grpc
     @check_input_types
     def add_midsurface_thickness(self, thickness: Quantity) -> None:  # noqa: D102
@@ -1463,6 +1482,11 @@ def volume(self) -> Quantity:  # noqa: D102
     def assign_material(self, material: Material) -> None:  # noqa: D102
         self._template.assign_material(material)
 
+    @property
+    @ensure_design_is_active
+    def get_assigned_material(self) -> Material:  # noqa: D102
+        return self._template.get_assigned_material
+
     @ensure_design_is_active
     def add_midsurface_thickness(self, thickness: Quantity) -> None:  # noqa: D102
         self._template.add_midsurface_thickness(thickness)

src/ansys/geometry/core/math/matrix.py

@@ -21,7 +21,10 @@
 # SOFTWARE.
 """Provides matrix primitive representations."""
 
-from typing import Union
+from typing import TYPE_CHECKING, Union
+
+if TYPE_CHECKING:
+    from ansys.geometry.core.math.vector import Vector3D  # For type hints
 
 from beartype import beartype as check_input_types
 import numpy as np
@@ -129,3 +132,165 @@ def __new__(cls, input: np.ndarray | RealSequence | Matrix = DEFAULT_MATRIX44):
             raise ValueError("Matrix44 should only be a 2D array of shape (4,4).")
 
         return obj
+
+    @classmethod
+    def create_translation(cls, translation: "Vector3D") -> "Matrix44":
+        """Create a matrix representing the specified translation.
+
+        Parameters
+        ----------
+        translation : Vector3D
+            The translation vector representing the translation. The components of the vector
+            should be in meters.
+
+        Returns
+        -------
+        Matrix44
+            A 4x4 matrix representing the translation.
+
+        Examples
+        --------
+        >>> translation_vector = Vector3D(1.0, 2.0, 3.0)
+        >>> translation_matrix = Matrix44.create_translation(translation_vector)
+        >>> print(translation_matrix)
+        [[1. 0. 0. 1.]
+         [0. 1. 0. 2.]
+         [0. 0. 1. 3.]
+         [0. 0. 0. 1.]]
+        """
+        matrix = cls(
+            [
+                [1, 0, 0, translation.x],
+                [0, 1, 0, translation.y],
+                [0, 0, 1, translation.z],
+                [0, 0, 0, 1],
+            ]
+        )
+        return matrix
+
+    def is_translation(self, including_identity=False):
+        """Check if the matrix represents a translation.
+
+        This method checks if the matrix represents a translation transformation.
+        A translation matrix has the following form:
+            [1 0 0 tx]
+            [0 1 0 ty]
+            [0 0 1 tz]
+            [0 0 0  1]
+
+        Parameters
+        ----------
+        including_identity : bool, optional
+            If True, the method will return True for the identity matrix as well.
+            If False, the method will return False for the identity matrix.
+
+        Returns
+        -------
+        bool
+            ``True`` if the matrix represents a translation, ``False`` otherwise.
+
+        Examples
+        --------
+        >>> matrix = Matrix44([[1, 0, 0, 5], [0, 1, 0, 3], [0, 0, 1, 2], [0, 0, 0, 1]])
+        >>> matrix.is_translation()
+        True
+        >>> identity_matrix = Matrix44([[1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 1, 0], [0, 0, 0, 1]])
+        >>> identity_matrix.is_translation()
+        False
+        >>> identity_matrix.is_translation(including_identity=True)
+        True
+        """
+        if not (
+            self.__is_close(self[0][0], 1)
+            and self.__is_close(self[0][1], 0)
+            and self.__is_close(self[0][2], 0)
+        ):
+            return False
+        if not (
+            self.__is_close(self[1][0], 0)
+            and self.__is_close(self[1][1], 1)
+            and self.__is_close(self[1][2], 0)
+        ):
+            return False
+        if not (
+            self.__is_close(self[2][0], 0)
+            and self.__is_close(self[2][1], 0)
+            and self.__is_close(self[2][2], 1)
+        ):
+            return False
+        if not self.__is_close(self[2][2], 1):
+            return False
+
+        if (
+            not including_identity
+            and self.__is_close(self[0][3], 0)
+            and self.__is_close(self[1][3], 0)
+            and self.__is_close(self[2][3], 0)
+        ):
+            return False
+
+        return True
+
+    def __is_close(self, a, b, tol=1e-9):
+        """Check if two values are close to each other within a tolerance."""
+        return np.isclose(a, b, atol=tol)
+
+    @classmethod
+    def create_rotation(
+        cls, direction_x: "Vector3D", direction_y: "Vector3D", direction_z: "Vector3D" = None
+    ) -> "Matrix44":
+        """Create a matrix representing the specified rotation.
+
+        Parameters
+        ----------
+        direction_x : Vector3D
+            The X direction vector.
+        direction_y : Vector3D
+            The Y direction vector.
+        direction_z : Vector3D, optional
+            The Z direction vector. If not provided, it will be calculated
+            as the cross product of direction_x and direction_y.
+
+        Returns
+        -------
+        Matrix44
+            A 4x4 matrix representing the rotation.
+
+        Examples
+        --------
+        >>> direction_x = Vector3D(1.0, 0.0, 0.0)
+        >>> direction_y = Vector3D(0.0, 1.0, 0.0)
+        >>> rotation_matrix = Matrix44.create_rotation(direction_x, direction_y)
+        >>> print(rotation_matrix)
+        [[1. 0. 0. 0.]
+        [0. 1. 0. 0.]
+        [0. 0. 1. 0.]
+        [0. 0. 0. 1.]]
+        """
+        if not direction_x.is_perpendicular_to(direction_y):
+            raise ValueError("The provided direction vectors are not orthogonal.")
+
+        # Normalize the vectors
+        direction_x = direction_x.normalize()
+        direction_y = direction_y.normalize()
+
+        # Calculate the third direction vector if not provided
+        if direction_z is None:
+            direction_z = direction_x.cross(direction_y)
+        else:
+            if not (
+                direction_x.is_perpendicular_to(direction_z)
+                and direction_y.is_perpendicular_to(direction_z)
+            ):
+                raise ValueError("The provided direction vectors are not orthogonal.")
+            direction_z = direction_z.normalize()
+
+        matrix = cls(
+            [
+                [direction_x.x, direction_y.x, direction_z.x, 0],
+                [direction_x.y, direction_y.y, direction_z.y, 0],
+                [direction_x.z, direction_y.z, direction_z.z, 0],
+                [0, 0, 0, 1],
+            ]
+        )
+        return matrix

src/ansys/geometry/core/modeler.py

@@ -40,10 +40,10 @@
 from ansys.geometry.core.logger import LOG
 from ansys.geometry.core.misc.checks import check_type, min_backend_version
 from ansys.geometry.core.misc.options import ImportOptions
-from ansys.geometry.core.misc.unsupported import UnsupportedCommands
 from ansys.geometry.core.tools.measurement_tools import MeasurementTools
 from ansys.geometry.core.tools.prepare_tools import PrepareTools
 from ansys.geometry.core.tools.repair_tools import RepairTools
+from ansys.geometry.core.tools.unsupported import UnsupportedCommands
 from ansys.geometry.core.typing import Real
 
 if TYPE_CHECKING:  # pragma: no cover
@@ -120,6 +120,9 @@ def __init__(
             backend_type=backend_type,
         )
 
+        # Maintaining references to all designs within the modeler workspace
+        self._designs: dict[str, "Design"] = {}
+
         # Initialize the RepairTools - Not available on Linux
         # TODO: delete "if" when Linux service is able to use repair tools
         # https://github.com/ansys/pyansys-geometry/issues/1319
@@ -135,9 +138,6 @@ def __init__(
         self._geometry_commands = GeometryCommands(self._grpc_client)
         self._unsupported = UnsupportedCommands(self._grpc_client, self)
 
-        # Maintaining references to all designs within the modeler workspace
-        self._designs: dict[str, "Design"] = {}
-
         # Check if the backend allows for multiple designs and throw warning if needed
         if not self.client.multiple_designs_allowed:
             LOG.warning(