capytaine 2.3.1__cp39-cp39-macosx_14_0_arm64.whl → 3.0.0a1__cp39-cp39-macosx_14_0_arm64.whl

capytaine/meshes/symmetric_meshes.py

@@ -0,0 +1,641 @@
+# Copyright 2025 Mews Labs
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+#
+#     http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+from __future__ import annotations
+
+import logging
+from typing import Optional, Union, Dict, Literal
+from functools import cached_property, lru_cache
+from itertools import chain
+
+import numpy as np
+
+from .abstract_meshes import AbstractMesh
+from .meshes import Mesh
+
+LOG = logging.getLogger(__name__)
+
+
+class ReflectionSymmetricMesh(AbstractMesh):
+    """A mesh with reflection symmetry across a plane.
+
+    This class represents a mesh that has reflection symmetry across either
+    the xOz plane (y=0) or yOz plane (x=0). Only half of the mesh is stored,
+    and the full mesh can be reconstructed by reflecting across the symmetry plane.
+
+    Supports nested symmetries: if the half mesh is itself a ReflectionSymmetricMesh,
+    this represents a quarter mesh with symmetries across both planes.
+
+    Attributes
+    ----------
+    half: AbstractMesh
+        The half mesh
+    plane: str
+        The symmetry plane, either "xOz" or "yOz"
+    faces_metadata: Dict[str, np.ndarray], optional
+        Some arrays with the same first dimension (should be the number
+        of faces of the whole mesh) storing some fields defined on all the
+        faces of the mesh.
+    name: str, optional
+        Name for the mesh
+
+    Examples
+    --------
+    >>> # Create a mesh with xOz symmetry (y=0 plane)
+    >>> half_mesh = Mesh(vertices=..., faces=...)
+    >>> symmetric_mesh = ReflectionSymmetricMesh(half=half_mesh, plane="xOz")
+    >>>
+    >>> # Create a mesh with both xOz and yOz symmetries (quarter mesh)
+    >>> quarter_mesh = Mesh(vertices=..., faces=...)
+    >>> sym_xOz = ReflectionSymmetricMesh(half=quarter_mesh, plane="xOz")
+    >>> sym_both = ReflectionSymmetricMesh(half=sym_xOz, plane="yOz")
+    >>>
+    >>> # Get the full merged mesh
+    >>> full_mesh = symmetric_mesh.merged()
+    """
+
+    def __init__(
+            self,
+            half: AbstractMesh, *,
+            plane: str,
+            faces_metadata: Optional[Dict[str, np.ndarray]] = None,
+            name: Optional[str] = None
+            ):
+
+        if plane not in ["xOz", "yOz"]:
+            raise ValueError(f"Plane must be 'xOz' or 'yOz', got '{plane}'")
+
+        self.half = half
+        self.plane = plane
+        if self.half.nb_faces > 0:
+            self.other_half = self.half.mirrored(plane)
+        else:
+            self.other_half = half  # Degenerate case without any face...
+
+        self.faces_metadata = {k: np.concatenate([v, v]) for k, v in half.faces_metadata.items()}
+        if faces_metadata is not None:
+            self.faces_metadata.update(**{k: np.asarray(faces_metadata[k]) for k in faces_metadata})
+
+        for m in self.faces_metadata:
+            assert self.faces_metadata[m].shape[0] == self.nb_faces
+
+        self.name = str(name) if name is not None else None
+
+    @cached_property
+    def nb_vertices(self) -> int:
+        return 2*self.half.nb_vertices
+
+    @cached_property
+    def nb_faces(self) -> int:
+        return 2*self.half.nb_faces
+
+    @cached_property
+    def vertices(self) -> np.ndarray:
+        return np.concatenate([self.half.vertices, self.other_half.vertices])
+
+    @cached_property
+    def faces(self) -> np.ndarray:
+        return np.concatenate([self.half.faces, self.other_half.faces])
+
+    @cached_property
+    def faces_normals(self) -> np.ndarray:
+        return np.concatenate([self.half.faces_normals, self.other_half.faces_normals])
+
+    @cached_property
+    def faces_areas(self) -> np.ndarray:
+        return np.concatenate([self.half.faces_areas, self.other_half.faces_areas])
+
+    @cached_property
+    def faces_centers(self) -> np.ndarray:
+        return np.concatenate([self.half.faces_centers, self.other_half.faces_centers])
+
+    @cached_property
+    def faces_radiuses(self) -> np.ndarray:
+        return np.concatenate([self.half.faces_radiuses, self.other_half.faces_radiuses])
+
+    @cached_property
+    def quadrature_points(self) -> np.ndarray:
+        return (
+                np.concatenate([self.half.quadrature_points[0], self.other_half.quadrature_points[0]]),
+                np.concatenate([self.half.quadrature_points[1], self.other_half.quadrature_points[1]]),
+                )
+
+    def with_quadrature(self, quadrature_method):
+        return ReflectionSymmetricMesh(
+                self.half.with_quadrature(quadrature_method),
+                plane=self.plane,
+                faces_metadata=self.faces_metadata,
+                name=self.name,
+                )
+
+    def __str__(self) -> str:
+        return (f"ReflectionSymmetricMesh(half={str(self.half)}"
+                + f", plane='{self.plane}'"
+                + (f", name=\"{self.name}\")" if self.name is not None else ")"))
+
+    def __short_str__(self) -> str:
+        return self.__str__()
+
+    def __repr__(self) -> str:
+        return self.__str__()
+
+    def _repr_pretty_(self, p, cycle):
+        p.text(self.__str__())
+
+    def extract_faces(self, faces_id, *, name=None) -> Mesh:
+        return self.merged().extract_faces(faces_id, name=name)
+
+    def translated(self, shift, *, name=None) -> Union[ReflectionSymmetricMesh, Mesh]:
+        if ((self.plane == 'xOz' and abs(shift[1]) < 1e-12)
+            or(self.plane == 'yOz' and abs(shift[0]) < 1e-12)):
+            return ReflectionSymmetricMesh(
+                half=self.half.translated(shift),
+                plane=self.plane,
+                faces_metadata=self.faces_metadata,
+                name=name
+            )
+        else:
+            return self.merged().translated(shift, name=name)
+
+    def rotated_with_matrix(self, R, *, name=None) -> Mesh:
+        return self.merged().rotated_with_matrix(R, name=name)
+
+    def mirrored(self, plane: Literal['xOz', 'yOz'], *, name=None) -> ReflectionSymmetricMesh:
+        return ReflectionSymmetricMesh(
+            half=self.half.mirrored(plane),
+            plane=self.plane,
+            faces_metadata=self.faces_metadata,
+            name=name
+        )
+
+    def join_meshes(self, *meshes, return_masks=False, name=None) -> Union[ReflectionSymmetricMesh, Mesh]:
+        if (all(isinstance(m, ReflectionSymmetricMesh) for m in meshes) and
+            all(m.plane == self.plane for m in meshes)):
+                if return_masks:
+                    joined_halves, half_masks = self.half.join_meshes(
+                        *[m.half for m in meshes],
+                        return_masks=True
+                    )
+                    masks = [np.concatenate([half_mask, half_mask]) for half_mask in half_masks]
+                else:
+                    joined_halves = self.half.join_meshes(
+                        *[m.half for m in meshes],
+                        return_masks=False
+                    )
+
+                faces_metadata = {k: np.concatenate(
+                    [m.faces_metadata[k][:m.nb_faces//2] for m in chain([self], meshes)]
+                    + [m.faces_metadata[k][m.nb_faces//2:] for m in chain([self], meshes)],
+                    axis=0
+                    )
+                                  for k in AbstractMesh._common_metadata_keys(*meshes)}
+
+                joined_mesh = ReflectionSymmetricMesh(
+                        half=joined_halves,
+                        plane=self.plane,
+                        faces_metadata=faces_metadata,
+                        name=name,
+                        )
+
+                if return_masks:
+                    return joined_mesh, masks
+                else:
+                    return joined_mesh
+
+        else:
+            return Mesh.join_meshes(
+                self.merged(),
+                *[m.merged() for m in meshes],
+                return_masks=return_masks,
+                name=name
+            )
+
+    def generate_lid(self, z=0.0, faces_max_radius=None, name=None):
+        return ReflectionSymmetricMesh(
+                self.half.generate_lid(z=z, faces_max_radius=faces_max_radius),
+                plane=self.plane,
+                name=name
+                )
+
+    def extract_lid(self, z=0.0):
+        half_hull, half_lid = self.half.extract_lid(z=z)
+        return (
+                ReflectionSymmetricMesh(half_hull, plane=self.plane),
+                ReflectionSymmetricMesh(half_lid, plane=self.plane),
+                )
+
+    def with_normal_vector_going_down(self, **kwargs) -> ReflectionSymmetricMesh:
+        return ReflectionSymmetricMesh(
+                half=self.half.with_normal_vector_going_down(),
+                plane=self.plane,
+                faces_metadata=self.faces_metadata,
+                name=self.name)
+
+    def copy(self, *, faces_metadata=None, name=None) -> ReflectionSymmetricMesh:
+        if faces_metadata is None:
+            faces_metadata = self.faces_metadata.copy()
+        if name is None:
+            name = self.name
+        return ReflectionSymmetricMesh(
+                half=self.half.copy(),
+                plane=self.plane,
+                faces_metadata=faces_metadata,
+                name=self.name)
+
+    @lru_cache
+    def merged(self, name=None) -> Mesh:
+        return Mesh.join_meshes(
+            self.half.merged(),
+            self.other_half.merged()
+        ).with_metadata(
+            **self.faces_metadata
+        )
+
+    def clipped(self, *, origin, normal, name=None) -> Union[ReflectionSymmetricMesh, Mesh]:
+        if ((self.plane == 'xOz' and abs(normal[0]) < 1e-12)
+            or(self.plane == 'yOz' and abs(normal[1]) < 1e-12)):
+            clipped_half, indices = (
+                    self.half
+                    .with_metadata(index=np.arange(self.half.nb_faces))
+                    .clipped(origin=origin, normal=normal)
+                    .pop_metadata("index")
+                    )
+            all_indices = np.concatenate([indices, indices + self.half.nb_faces])
+            metadata = {k: self.faces_metadata[k][all_indices] for k in self.faces_metadata.keys()}
+            return ReflectionSymmetricMesh(
+                    half=clipped_half,
+                    plane=self.plane,
+                    faces_metadata=metadata,
+                    name=name)
+        else:
+            LOG.warning("Dropping mesh reflection symmetry with respect to "
+                        f"{self.plane} when clipping with respect to plane "
+                        f"with origin {origin} and normal {normal}")
+            return self.merged().clipped(origin=origin, normal=normal, name=name)
+
+    def show(self, *, backend=None, ghost_meshes=None, **kwargs):
+        if ghost_meshes is None:
+            ghost_meshes = []
+        ghost_meshes = ghost_meshes + [self.other_half.merged()]
+        return self.half.show(backend=backend, ghost_meshes=ghost_meshes, **kwargs)
+
+    def export(self, format, **kwargs):
+        LOG.warning(f"Losing symmetric structure when exporting {self} to {format}")
+        return self.merged().export(format, **kwargs)
+
+
+# For some backward compatibility:
+yOz_Plane = "yOz"
+xOz_Plane = "xOz"
+
+
+class RotationSymmetricMesh(AbstractMesh):
+    """A mesh with rotation symmetry around the Oz axis.
+
+    This class represents a mesh that has n-fold rotational symmetry about
+    the z-axis. Only a wedge (1/n of the full mesh) is stored, and the full
+    mesh can be reconstructed by rotating the wedge n times.
+
+    Supports nested symmetries: the wedge mesh can be a ReflectionSymmetricMesh
+    for dihedral symmetry.
+
+    Attributes
+    ----------
+    wedge: AbstractMesh
+        The wedge mesh (1/n of the full mesh)
+    n: int
+        The rotation order (number of rotations to complete full circle)
+    axis: either 'z+' or 'z-'
+        Only the z-axis is supported, but two possible orientations can be used.
+        Both are equivalent, except for the ordering of the sub-meshes.
+    faces_metadata: Dict[str, np.ndarray], optional
+        Some arrays with the same first dimension (should be the number
+        of faces of the whole mesh) storing some fields defined on all the
+        faces of the mesh.
+    name: str, optional
+        Name for the mesh
+
+    Examples
+    --------
+    >>> # Create a mesh with 3-fold rotation symmetry about z-axis
+    >>> wedge_mesh = Mesh(vertices=..., faces=...)
+    >>> symmetric_mesh = RotationSymmetricMesh(wedge=wedge_mesh, n=3)
+    >>>
+    >>> # Get the full merged mesh
+    >>> full_mesh = symmetric_mesh.merged()
+    """
+
+    def __init__(
+            self,
+            wedge: AbstractMesh,
+            n: int, *,
+            axis: Literal['z+', 'z-'] = 'z+',
+            faces_metadata: Optional[Dict[str, np.ndarray]] = None,
+            name: Optional[str] = None
+            ):
+
+        if isinstance(wedge, ReflectionSymmetricMesh) and n > 4:
+            LOG.warning("RotationSymmetricMesh with n>4 and containing ReflectionSymmetricMesh are not fully supported at the moment. "
+                        "You might prefer to define instead a ReflectionSymmetricMesh of a RotationSymmetricMesh.")
+
+        if n < 2:
+            raise ValueError(f"Rotation order must be >= 2, got: {n}")
+
+        self.wedge = wedge
+        self.n = n
+        self.axis = axis
+        if self.axis == 'z+':
+            self.all_wedges = [self.wedge] + [self.wedge.rotated_z(2*i*np.pi/n) for i in range(1, n)]
+        elif self.axis == 'z-':
+            self.all_wedges = [self.wedge] + [self.wedge.rotated_z(-2*i*np.pi/n) for i in range(1, n)]
+        else:
+            raise ValueError(f"Unsupported axis for RotationSymmetricMesh: {axis}")
+
+        self.faces_metadata = {k: np.concatenate([v]*n) for k, v in wedge.faces_metadata.items()}
+        if faces_metadata is not None:
+            self.faces_metadata.update(**{k: np.asarray(faces_metadata[k]) for k in faces_metadata})
+
+        for m in self.faces_metadata:
+            assert self.faces_metadata[m].shape[0] == self.nb_faces
+
+        self.name = str(name) if name is not None else None
+
+    @classmethod
+    def from_profile_points(cls, points: np.ndarray, n: int, *, faces_metadata=None, name=None):
+        """Return the mesh defined by the set of `points` repeated `n` times around the z-axis.
+
+        Points will be sorted by increasing z-coordinate before making a mesh,
+        in order to ensure that the normal vector are outwards.
+
+        Parameters
+        ---------
+        points: array of shape (..., 3)
+            A list of points in 3D.
+        n: int
+            The rotation order (number of rotations to complete full circle)
+        faces_metadata: Dict[str, np.ndarray], optional
+            Some arrays with the same first dimension (should be the number
+            of faces of the whole mesh) storing some fields defined on all the
+            faces of the mesh.
+        name: str, optional
+            Name for the mesh
+
+        Example
+        -------
+        >>> meridian_points = np.array([(np.sqrt(1-z**2), 0.0, z) for z in np.linspace(-1.0, 1.0, 10)])
+        >>> sphere = RotationSymmetricMesh.from_profile_points(meridian_points, n=10)
+        """
+        c, s = np.cos(2*np.pi/n), np.sin(2*np.pi/n)
+        rotation_matrix = np.array([[c, -s, 0], [s, c, 0], [0, 0, 1]])
+
+        points = np.asarray(sorted(list(points), key=lambda p: p[2]))  # Sort by increasing z
+        vertices = np.concatenate([points, points @ rotation_matrix.T])
+        faces = np.array([(i, i+len(points), i+len(points)+1, i+1) for i in range(len(points)-1)])
+        wedge = Mesh(vertices=vertices, faces=faces)
+
+        return RotationSymmetricMesh(
+            wedge=wedge,
+            n=n,
+            faces_metadata=faces_metadata,
+            name=name
+        )
+
+    @cached_property
+    def nb_vertices(self) -> int:
+        return self.n * self.wedge.nb_vertices
+
+    @cached_property
+    def nb_faces(self) -> int:
+        return self.n * self.wedge.nb_faces
+
+    @cached_property
+    def vertices(self) -> np.ndarray:
+        return np.concatenate([w.vertices for w in self.all_wedges])
+
+    @cached_property
+    def faces(self) -> np.ndarray:
+        return np.concatenate([w.faces for w in self.all_wedges])
+
+    @cached_property
+    def faces_normals(self) -> np.ndarray:
+        return np.concatenate([w.faces_normals for w in self.all_wedges])
+
+    @cached_property
+    def faces_areas(self) -> np.ndarray:
+        return np.concatenate([w.faces_areas for w in self.all_wedges])
+
+    @cached_property
+    def faces_centers(self) -> np.ndarray:
+        return np.concatenate([w.faces_centers for w in self.all_wedges])
+
+    @cached_property
+    def faces_radiuses(self) -> np.ndarray:
+        return np.concatenate([w.faces_radiuses for w in self.all_wedges])
+
+    @cached_property
+    def quadrature_points(self) -> np.ndarray:
+        return (
+                np.concatenate([w.quadrature_points[0] for w in self.all_wedges]),
+                np.concatenate([w.quadrature_points[1] for w in self.all_wedges]),
+                )
+
+    def with_quadrature(self, quadrature_method):
+        return RotationSymmetricMesh(
+                self.wedge.with_quadrature(quadrature_method),
+                n=self.n,
+                axis=self.axis,
+                faces_metadata=self.faces_metadata,
+                name=self.name,
+                )
+
+    def __str__(self) -> str:
+        return (f"RotationSymmetricMesh(wedge={str(self.wedge)}"
+                + f", n={self.n}"
+                + (f", name=\"{self.name}\")" if self.name is not None else ")"))
+
+    def __short_str__(self) -> str:
+        return self.__str__()
+
+    def __repr__(self) -> str:
+        return self.__str__()
+
+    def _repr_pretty_(self, p, cycle):
+        p.text(self.__str__())
+
+    def extract_faces(self, faces_id, *, name=None) -> Mesh:
+        return self.merged().extract_faces(faces_id, name=name)
+
+    def translated(self, shift, *, name=None) -> Union[RotationSymmetricMesh, Mesh]:
+        if (abs(shift[0]) < 1e-12 and abs(shift[1] < 1e-12)):
+            # Vertical translation
+            return RotationSymmetricMesh(
+                    self.wedge.translated_z(shift[2]),
+                    n=self.n,
+                    axis=self.axis,
+                    faces_metadata=self.faces_metadata,
+                    name=name)
+        else:
+            return self.merged().translated(shift, name=name)
+
+    def rotated_with_matrix(self, R, *, name=None) -> Union[RotationSymmetricMesh, Mesh]:
+        if (np.allclose(R[:, 2], [0.0, 0.0, 1.0])
+            and np.allclose(R[2, :], [0.0, 0.0, 1.0])):
+            # Rotation around the z-axis: we keep the symmetry
+            return RotationSymmetricMesh(
+                    self.wedge.rotated_with_matrix(R),
+                    n=self.n,
+                    axis=self.axis,
+                    faces_metadata=self.faces_metadata,
+                    name=name,
+                    )
+        else:
+            return self.merged().rotated_with_matrix(R, name=name)
+
+    @property
+    def _opposite_axis(self):
+        if self.axis == 'z+':
+            return 'z-'
+        else:
+            return 'z+'
+
+
+    def mirrored(self, plane: Literal['xOz', 'yOz'], *, name=None) -> RotationSymmetricMesh:
+        return RotationSymmetricMesh(
+            wedge=self.wedge.mirrored(plane),
+            n=self.n,
+            axis=self._opposite_axis,
+            faces_metadata=self.faces_metadata,
+            name=name
+        )
+
+    def _metadata_of_wedge(self, k, i):
+        return self.faces_metadata[k][i*self.wedge.nb_faces:(i+1)*self.wedge.nb_faces]
+
+    def join_meshes(self, *meshes, return_masks=False, name=None) -> Union[RotationSymmetricMesh, Mesh]:
+        if (all(isinstance(m, RotationSymmetricMesh) for m in meshes) and
+            all(m.n == self.n for m in meshes)):
+                if return_masks:
+                    joined_wegdes, wedges_masks = self.wedge.join_meshes(
+                        *[m.wedge for m in meshes],
+                        return_masks=True
+                    )
+                    masks = [np.concatenate([w_mesh]*self.n) for w_mesh in wedges_masks]
+                else:
+                    joined_wegdes = self.wedge.join_meshes(
+                        *[m.wedge for m in meshes],
+                        return_masks=False
+                    )
+
+                faces_metadata = {k: np.concatenate(
+                    [m._metadata_of_wedge(k, i) for i in range(self.n) for m in chain([self], meshes)] ,
+                    axis=0
+                    )
+                                  for k in AbstractMesh._common_metadata_keys(*meshes)}
+
+                joined_mesh = RotationSymmetricMesh(
+                        wedge=joined_wegdes,
+                        n=self.n,
+                        axis=self.axis,
+                        faces_metadata=faces_metadata,
+                        name=name,
+                        )
+
+                if return_masks:
+                    return joined_mesh, masks
+                else:
+                    return joined_mesh
+
+        else:
+            return Mesh.join_meshes(
+                self.merged(),
+                *[m.merged() for m in meshes],
+                return_masks=return_masks,
+                name=name
+            )
+
+    def generate_lid(self, z=0.0, faces_max_radius=None, name=None):
+        return RotationSymmetricMesh(
+                self.wedge.generate_lid(z=z, faces_max_radius=faces_max_radius),
+                axis=self.axis,
+                n=self.n,
+                name=name
+                )
+
+    def extract_lid(self, z=0.0):
+        wedge_hull, wedge_lid = self.wedge.extract_lid(z=z)
+        return (
+                RotationSymmetricMesh(wedge_hull, axis=self.axis, n=self.n),
+                RotationSymmetricMesh(wedge_lid, axis=self.axis, n=self.n),
+                )
+
+    def with_normal_vector_going_down(self, **kwargs) -> RotationSymmetricMesh:
+        return RotationSymmetricMesh(
+                wedge=self.wedge.with_normal_vector_going_down(),
+                n=self.n,
+                axis=self.axis,
+                faces_metadata=self.faces_metadata,
+                name=self.name)
+
+    def copy(self, *, faces_metadata=None, name=None) -> RotationSymmetricMesh:
+        if faces_metadata is None:
+            faces_metadata = self.faces_metadata.copy()
+        if name is None:
+            name = self.name
+        return RotationSymmetricMesh(
+                wedge=self.wedge.copy(),
+                n=self.n,
+                axis=self.axis,
+                faces_metadata=faces_metadata,
+                name=self.name)
+
+    @lru_cache
+    def merged(self, name=None) -> Mesh:
+        return Mesh.join_meshes(
+                *[w.merged() for w in self.all_wedges]
+        ).with_metadata(
+            **self.faces_metadata
+        )
+
+    def clipped(self, *, origin, normal, name=None) -> Union[RotationSymmetricMesh, Mesh]:
+        if (abs(normal[0]) < 1e-12 and abs(normal[1]) < 1e-12):
+            # Horizontal plane
+            clipped_wedge, indices = (
+                    self.wedge
+                    .with_metadata(index=np.arange(self.wedge.nb_faces))
+                    .clipped(origin=origin, normal=normal)
+                    .pop_metadata("index")
+                    )
+            all_indices = np.concatenate([indices + i*self.wedge.nb_faces for i in range(self.n)])
+            metadata = {k: self.faces_metadata[k][all_indices] for k in self.faces_metadata.keys()}
+            return RotationSymmetricMesh(
+                    wedge=clipped_wedge,
+                    n=self.n,
+                    axis=self.axis,
+                    faces_metadata=metadata,
+                    name=name)
+        else:
+            LOG.warning("Dropping mesh rotation symmetry with respect to "
+                        f"z-axis when clipping with respect to plane "
+                        f"with origin {origin} and normal {normal}")
+            return self.merged().clipped(origin=origin, normal=normal, name=name)
+
+    def show(self, *, backend=None, ghost_meshes=None, **kwargs):
+        if ghost_meshes is None:
+            ghost_meshes = []
+        ghost_meshes = ghost_meshes + [w.merged() for w in self.all_wedges[1:]]
+        return self.wedge.show(backend=backend, ghost_meshes=ghost_meshes, **kwargs)
+
+    def export(self, format, **kwargs):
+        LOG.warning(f"Losing symmetric structure when exporting {self} to {format}")
+        return self.merged().export(format, **kwargs)