capytaine 2.3__cp312-cp312-win_amd64.whl → 2.3.1__cp312-cp312-win_amd64.whl

capytaine/__about__.py

@@ -5,7 +5,7 @@ __all__ = ["__title__", "__description__", "__version__", "__author__", "__uri__
 __title__ = "capytaine"
 __description__ = """Python BEM solver for linear potential flow, based on Nemoh"""
 
-__version__ = "2.3"
+__version__ = "2.3.1"
 
 __author__ = "Matthieu Ancellin"
 __uri__ = "https://github.com/capytaine/capytaine"

capytaine/__init__.py

@@ -3,14 +3,14 @@
 
 
 # start delvewheel patch
-def _delvewheel_patch_1_10_1():
+def _delvewheel_patch_1_11_2():
     import os
     if os.path.isdir(libs_dir := os.path.abspath(os.path.join(os.path.dirname(__file__), os.pardir, 'capytaine.libs'))):
         os.add_dll_directory(libs_dir)
 
 
-_delvewheel_patch_1_10_1()
-del _delvewheel_patch_1_10_1
+_delvewheel_patch_1_11_2()
+del _delvewheel_patch_1_11_2
 # end delvewheel patch
 
 from .__about__ import (

capytaine/bem/problems_and_results.py

@@ -360,16 +360,18 @@ class DiffractionProblem(LinearPotentialFlowProblem):
 
         if self.body is not None:
 
-            self.boundary_condition = -(
+            self.boundary_condition = np.zeros(
+                shape=(self.body.mesh_including_lid.nb_faces,),
+                dtype=np.complex128
+            )
+
+            self.boundary_condition[self.body.hull_mask] = -(
                     airy_waves_velocity(self.body.mesh.faces_centers, self)
                     * self.body.mesh.faces_normals
             ).sum(axis=1)
             # Note that even with forward speed, this is computed based on the
             # frequency and not the encounter frequency.
 
-            if self.body.lid_mesh is not None:
-                self.boundary_condition = np.concatenate([self.boundary_condition, np.zeros(self.body.lid_mesh.nb_faces)])
-
             if len(self.body.dofs) == 0:
                 LOG.warning(f"The body {self.body.name} used in diffraction problem has no dofs!")
 
@@ -421,11 +423,16 @@ class RadiationProblem(LinearPotentialFlowProblem):
 
             dof = self.body.dofs[self.radiating_dof]
 
-            displacement_on_face = np.sum(dof * self.body.mesh.faces_normals, axis=1)  # This is a dot product on each face
-            if self.body.lid_mesh is not None:
-                displacement_on_face = np.concatenate([displacement_on_face, np.zeros(self.body.lid_mesh.nb_faces)])
+            self.boundary_condition = self.encounter_omega * np.zeros(
+                shape=(self.body.mesh_including_lid.nb_faces,),
+                dtype=np.complex128
+            )
+            # The multiplication by encounter_omega is just a programming trick to ensure that boundary_condition
+            # is implemented with the correct type (for zero and infinite frequencies), it does not affect the value.
+            # Below the value is update on the hull. It remains zero on the lid.
 
-            self.boundary_condition = -1j * self.encounter_omega * displacement_on_face
+            displacement_on_face = np.sum(dof * self.body.mesh.faces_normals, axis=1)  # This is a dot product on each face
+            self.boundary_condition[self.body.hull_mask] = -1j * self.encounter_omega * displacement_on_face
 
             if self.forward_speed != 0.0:
                 if self.radiating_dof.lower() == "pitch":
@@ -440,7 +447,8 @@ class RadiationProblem(LinearPotentialFlowProblem):
                             "Only radiating dofs with name in {'Surge', 'Sway', 'Heave', 'Roll', 'Pitch', 'Yaw'} are supported.\n"
                             f"Got instead `radiating_dof={self.radiating_dof}`"
                             )
-                self.boundary_condition += self.forward_speed * ddofdx_dot_n
+
+                self.boundary_condition[self.body.hull_mask] += self.forward_speed * ddofdx_dot_n
 
 
 

capytaine/bem/solver.py

@@ -183,7 +183,7 @@ class BEMSolver:
                 nabla_phi = self._compute_potential_gradient(problem.body.mesh_including_lid, result)
                 pressure += problem.rho * problem.forward_speed * nabla_phi[:, 0]
 
-        pressure_on_hull = pressure[:problem.body.mesh.nb_faces]  # Discards pressure on lid if any
+        pressure_on_hull = pressure[problem.body.hull_mask]  # Discards pressure on lid if any
         forces = problem.body.integrate_pressure(pressure_on_hull)
 
         if not keep_details:

capytaine/bodies/bodies.py

@@ -114,10 +114,21 @@ class FloatingBody(ClippableMixin, Abstract3DObject):
         else:
             self.dofs = dofs
 
+        self._evaluate_full_mesh()
+
         LOG.info(f"New floating body: {self.__str__()}.")
 
         self._check_dofs_shape_consistency()
 
+    def _evaluate_full_mesh(self):
+        """Merge the mesh and lid_mesh, while keeping track of where each panel came from."""
+        if self.lid_mesh is None:
+            self.mesh_including_lid = self.mesh
+            self.hull_mask = np.full((self.mesh.nb_faces,), True)
+        else:
+            self.mesh_including_lid, masks = self.mesh.join_meshes(self.lid_mesh, return_masks=True)
+            self.hull_mask = masks[0]
+
     @staticmethod
     def from_meshio(mesh, name=None) -> 'FloatingBody':
         """Create a FloatingBody from a meshio mesh object.
@@ -138,13 +149,6 @@ class FloatingBody(ClippableMixin, Abstract3DObject):
         """Arbitrary order. The point is to sort together the problems involving the same body."""
         return self.name < other.name
 
-    @cached_property
-    def mesh_including_lid(self):
-        if self.lid_mesh is not None:
-            return self.mesh.join_meshes(self.lid_mesh)
-        else:
-            return self.mesh
-
     ##########
     #  Dofs  #
     ##########
@@ -570,7 +574,7 @@ class FloatingBody(ClippableMixin, Abstract3DObject):
                 )
             for radiating_dof_name in self.dofs
             for influenced_dof_name in self.dofs
-            ])
+            ], compat='no_conflicts', join="outer")
 
         # Reorder dofs
         K = hs_set.hydrostatic_stiffness.sel(influenced_dof=list(self.dofs.keys()), radiating_dof=list(self.dofs.keys()))
@@ -660,7 +664,7 @@ class FloatingBody(ClippableMixin, Abstract3DObject):
                                             'radiating_dof': body_dof_names},
                                     name="inertia_matrix")
 
-        total_mass_xr = xr.merge([rigid_inertia_matrix_xr, other_dofs_inertia_matrix_xr], compat="override").inertia_matrix
+        total_mass_xr = xr.merge([rigid_inertia_matrix_xr, other_dofs_inertia_matrix_xr], compat="override", join="outer").inertia_matrix
 
         non_rigid_dofs = set(body_dof_names) - set(rigid_dof_names)
 
@@ -980,6 +984,7 @@ respective inertia coefficients are assigned as NaN.")
         self.mesh.mirror(plane)
         if self.lid_mesh is not None:
             self.lid_mesh.mirror(plane)
+        self._evaluate_full_mesh()
         for dof in self.dofs:
             self.dofs[dof] -= 2 * np.outer(np.dot(self.dofs[dof], plane.normal), plane.normal)
         for point_attr in ('geometric_center', 'rotation_center', 'center_of_mass'):
@@ -994,6 +999,7 @@ respective inertia coefficients are assigned as NaN.")
         self.mesh.translate(vector, *args, **kwargs)
         if self.lid_mesh is not None:
             self.lid_mesh.translate(vector, *args, **kwargs)
+        self._evaluate_full_mesh()
         for point_attr in ('geometric_center', 'rotation_center', 'center_of_mass'):
             if point_attr in self.__dict__ and self.__dict__[point_attr] is not None:
                 self.__dict__[point_attr] = np.array(self.__dict__[point_attr]) + vector
@@ -1004,6 +1010,7 @@ respective inertia coefficients are assigned as NaN.")
         self.mesh.rotate(axis, angle)
         if self.lid_mesh is not None:
             self.lid_mesh.rotate(axis, angle)
+        self._evaluate_full_mesh()
         for point_attr in ('geometric_center', 'rotation_center', 'center_of_mass'):
             if point_attr in self.__dict__ and self.__dict__[point_attr] is not None:
                 self.__dict__[point_attr] = axis.rotate_points([self.__dict__[point_attr]], angle)[0, :]
@@ -1023,6 +1030,7 @@ respective inertia coefficients are assigned as NaN.")
             if self.lid_mesh.nb_faces == 0:
                 LOG.warning("Lid mesh %s is empty after clipping. The lid mesh is removed.", self.lid_mesh)
                 self.lid_mesh = None
+        self._evaluate_full_mesh()
 
         # Clip dofs
         ids = self.mesh._clipping_data['faces_ids']

capytaine/green_functions/LiangWuNoblesse/Makefile

@@ -14,5 +14,3 @@ $(BUILD_DIR)/LiangWuNoblesseWaveTerm.o: LiangWuNoblesseWaveTerm.f90
 clean:
 	rm -rf $(BUILD_DIR)
 .PHONY: run_test clean
-
-

capytaine/green_functions/delhommeau.py

@@ -215,21 +215,25 @@ class Delhommeau(AbstractGreenFunction):
         filepath = os.path.join(tabulation_cache_dir, filename)
 
         if os.path.exists(filepath):
-            LOG.info("Loading tabulation from %s", filepath)
-            loaded_arrays = np.load(filepath)
-            self.tabulated_r_range = loaded_arrays["r_range"]
-            self.tabulated_z_range = loaded_arrays["z_range"]
-            self.tabulated_integrals = loaded_arrays["values"]
-
-        else:
-            self._create_tabulation(tabulation_nr, tabulation_rmax,
-                                    tabulation_nz, tabulation_zmin,
-                                    tabulation_nb_integration_points)
-            LOG.debug("Saving tabulation in %s", filepath)
-            np.savez_compressed(
-                filepath, r_range=self.tabulated_r_range, z_range=self.tabulated_z_range,
-                values=self.tabulated_integrals
-            )
+            try:
+                LOG.info("Loading tabulation from %s", filepath)
+                loaded_arrays = np.load(filepath)
+                self.tabulated_r_range = loaded_arrays["r_range"]
+                self.tabulated_z_range = loaded_arrays["z_range"]
+                self.tabulated_integrals = loaded_arrays["values"]
+                return filename
+            except (EOFError, FileNotFoundError, KeyError, ValueError):
+                LOG.warning("Error loading tabulation from %s", filepath)
+
+        self._create_tabulation(tabulation_nr, tabulation_rmax,
+                                tabulation_nz, tabulation_zmin,
+                                tabulation_nb_integration_points)
+        LOG.debug("Saving tabulation in %s", filepath)
+        np.savez_compressed(
+            filepath, r_range=self.tabulated_r_range, z_range=self.tabulated_z_range,
+            values=self.tabulated_integrals
+        )
+        return filename
 
     def _create_tabulation(self, tabulation_nr, tabulation_rmax,
                                    tabulation_nz, tabulation_zmin,

capytaine/io/xarray.py

@@ -257,7 +257,7 @@ def hydrostatics_dataset(bodies: Sequence[FloatingBody]) -> xr.Dataset:
         if len(bodies_properties) > 0:
             bodies_properties = xr.concat(bodies_properties.values(), pd.Index(bodies_properties.keys(), name='body_name'))
             bodies_properties = _squeeze_dimensions(bodies_properties, dimensions=['body_name'])
-            dataset = xr.merge([dataset, {body_property: bodies_properties}])
+            dataset = xr.merge([dataset, {body_property: bodies_properties}], compat="no_conflicts", join="outer")
     return dataset
 
 
@@ -416,7 +416,7 @@ def assemble_dataset(results,
             variables=['added_mass', 'radiation_damping'],
             dimensions=[main_freq_type, 'radiating_dof', 'influenced_dof'],
             optional_dims=optional_dims + ['wave_direction'])
-        dataset = xr.merge([dataset, radiation_cases])
+        dataset = xr.merge([dataset, radiation_cases], compat="no_conflicts", join="outer")
 
     # DIFFRACTION RESULTS
     if "DiffractionResult" in kinds_of_results:
@@ -425,7 +425,7 @@ def assemble_dataset(results,
             variables=['diffraction_force', 'Froude_Krylov_force'],
             dimensions=[main_freq_type, 'wave_direction', 'influenced_dof'],
             optional_dims=optional_dims)
-        dataset = xr.merge([dataset, diffraction_cases])
+        dataset = xr.merge([dataset, diffraction_cases], compat="no_conflicts", join="outer")
         dataset['excitation_force'] = dataset['Froude_Krylov_force'] + dataset['diffraction_force']
 
     # OTHER FREQUENCIES TYPES
@@ -519,7 +519,7 @@ def assemble_dataset(results,
             LOG.warning('Bemio data import being used, hydrostatics=True is ignored.')
         else:
             bodies = list({result.body for result in results})
-            dataset = xr.merge([dataset, hydrostatics_dataset(bodies)])
+            dataset = xr.merge([dataset, hydrostatics_dataset(bodies)], compat="no_conflicts", join="outer")
 
     for var in set(dataset) | set(dataset.coords):
         if var in VARIABLES_ATTRIBUTES:

capytaine/meshes/collections.py

@@ -223,8 +223,16 @@ class CollectionOfMeshes(ClippableMixin, SurfaceIntegralsMixin, Abstract3DObject
     # Transformation #
     ##################
 
-    def join_meshes(*meshes, name=None):
-        return CollectionOfMeshes(meshes, name=name)
+    def join_meshes(*meshes, name=None, return_masks=False):
+        coll = CollectionOfMeshes(meshes, name=name)
+        if return_masks:
+            masks = []
+            for i_mesh in range(len(meshes)):
+                mask = np.full((coll.nb_faces,), False)
+                mask[coll.indices_of_mesh(i_mesh)] = True
+                masks.append(mask)
+            return coll, masks
+        return coll
 
     def __add__(self, mesh_to_add):
         return self.join_meshes(mesh_to_add)

capytaine/meshes/mesh_like_protocol.py

@@ -30,7 +30,7 @@ class MeshLike(Protocol):
     def extract_faces(self, faces_id):
         ...
 
-    def join_meshes(*meshes):
+    def join_meshes(*meshes, return_mask):
         ...
 
     def with_normal_vector_going_down(self, **kwargs):

capytaine/meshes/meshes.py

@@ -644,9 +644,18 @@ class Mesh(ClippableMixin, SurfaceIntegralsMixin, Abstract3DObject):
     #  Combine meshes  #
     ####################
 
-    def join_meshes(*meshes, name=None):
+    def join_meshes(*meshes, name=None, return_masks=False):
         from capytaine.meshes.collections import CollectionOfMeshes
-        return CollectionOfMeshes(meshes, name=name).merged()
+        coll = CollectionOfMeshes(meshes, name=name)
+        if return_masks:
+            masks = []
+            for i_mesh in range(len(meshes)):
+                mask = np.full((coll.nb_faces,), False)
+                mask[coll.indices_of_mesh(i_mesh)] = True
+                masks.append(mask)
+            return coll.merged(), masks
+        else:
+            return coll.merged()
 
     def __add__(self, mesh_to_add) -> 'Mesh':
         return self.join_meshes(mesh_to_add)

capytaine/meshes/symmetric.py

@@ -84,13 +84,27 @@ class ReflectionSymmetricMesh(SymmetricMesh):
     def __deepcopy__(self, *args):
         return ReflectionSymmetricMesh(self.half.copy(), self.plane, name=self.name)
 
-    def join_meshes(*meshes, name=None):
+    def join_meshes(*meshes, name=None, return_masks=False):
         assert all(isinstance(mesh, ReflectionSymmetricMesh) for mesh in meshes), \
             "Only meshes with the same symmetry can be joined together."
         assert all(meshes[0].plane == mesh.plane for mesh in meshes), \
             "Only reflection symmetric meshes with the same reflection plane can be joined together."
-        half_mesh = meshes[0].half.join_meshes(*(mesh.half for mesh in meshes[1:]), name=f"half_of_{name}" if name is not None else None)
-        return ReflectionSymmetricMesh(half_mesh, plane=meshes[0].plane, name=name)
+        if not return_masks:
+            name = name=f"half_of_{name}" if name is not None else None
+            half_mesh = meshes[0].half.join_meshes(
+                *(mesh.half for mesh in meshes[1:]),
+                name=name, return_masks=False
+            )
+            return ReflectionSymmetricMesh(half_mesh, plane=meshes[0].plane, name=name)
+        else:
+            name = name=f"half_of_{name}" if name is not None else None
+            half_mesh, half_masks = meshes[0].half.join_meshes(
+                *(mesh.half for mesh in meshes[1:]),
+                name=name, return_masks=True
+            )
+            masks = [np.concatenate([half_mask, half_mask]) for half_mask in half_masks]
+            joined = ReflectionSymmetricMesh(half_mesh, plane=meshes[0].plane, name=name)
+            return joined, masks
 
     @inplace_transformation
     def translate(self, vector):
@@ -203,15 +217,43 @@ class TranslationalSymmetricMesh(SymmetricMesh):
         CollectionOfMeshes.mirror(self, plane)
         return self
 
-    def join_meshes(*meshes, name=None):
+    def join_meshes(*meshes, name=None, return_masks=False):
         assert all(isinstance(mesh, TranslationalSymmetricMesh) for mesh in meshes), \
             "Only meshes with the same symmetry can be joined together."
         assert all(np.allclose(meshes[0].translation, mesh.translation) for mesh in meshes), \
             "Only translation symmetric meshes with the same translation vector can be joined together."
         assert all(len(meshes[0]) == len(mesh) for mesh in meshes), \
             "Only symmetric meshes with the same number of elements can be joined together."
-        mesh_strip = CollectionOfMeshes([mesh.first_slice for mesh in meshes], name=f"strip_of_{name}" if name is not None else None)
-        return TranslationalSymmetricMesh(mesh_strip, translation=meshes[0].translation, nb_repetitions=len(meshes[0]) - 1, name=name)
+        if not return_masks:
+            strip_name = f"strip_of_{name}" if name is not None else None
+            mesh_strip = meshes[0].first_slice.join_meshes(
+                *(mesh.first_slice for mesh in meshes[1:]),
+                name=strip_name,
+                return_masks=False
+            )
+            return TranslationalSymmetricMesh(
+                mesh_strip,
+                translation=meshes[0].translation,
+                nb_repetitions=len(meshes[0]) - 1,
+                name=name
+            )
+        else:
+            strip_name = f"strip_of_{name}" if name is not None else None
+            mesh_strip, strip_masks = meshes[0].first_slice.join_meshes(
+                *(mesh.first_slice for mesh in meshes[1:]),
+                name=strip_name,
+                return_masks=True
+            )
+            joined = TranslationalSymmetricMesh(
+                mesh_strip,
+                translation=meshes[0].translation,
+                nb_repetitions=len(meshes[0]) - 1,
+                name=name
+            )
+            masks = [np.concatenate([
+                strip_mask for _ in range(len(meshes[0]))
+            ]) for strip_mask in strip_masks]
+            return joined, masks
 
 
 def build_regular_array_of_meshes(base_mesh, distance, nb_bodies):
@@ -363,15 +405,43 @@ class AxialSymmetricMesh(SymmetricMesh):
     def __deepcopy__(self, *args):
         return AxialSymmetricMesh(self.first_slice.copy(), axis=self.axis.copy(), nb_repetitions=len(self) - 1, name=self.name)
 
-    def join_meshes(*meshes, name=None):
+    def join_meshes(*meshes, name=None, return_masks=False):
         assert all(isinstance(mesh, AxialSymmetricMesh) for mesh in meshes), \
             "Only meshes with the same symmetry can be joined together."
         assert all(meshes[0].axis == mesh.axis for mesh in meshes), \
             "Only axisymmetric meshes with the same symmetry axis can be joined together."
         assert all(len(meshes[0]) == len(mesh) for mesh in meshes), \
             "Only axisymmetric meshes with the same number of elements can be joined together."
-        mesh_slice = CollectionOfMeshes([mesh.first_slice for mesh in meshes], name=f"slice_of_{name}" if name is not None else None)
-        return AxialSymmetricMesh(mesh_slice, axis=meshes[0].axis, nb_repetitions=len(meshes[0]) - 1, name=name)
+        if not return_masks:
+            slice_name = f"slice_of_{name}" if name is not None else None
+            mesh_slice = meshes[0].first_slice.join_meshes(
+                *(mesh.first_slice for mesh in meshes[1:]),
+                name=slice_name,
+                return_masks=False
+            )
+            return AxialSymmetricMesh(
+                mesh_slice,
+                axis=meshes[0].axis,
+                nb_repetitions=len(meshes[0]) - 1,
+                name=name
+            )
+        else:
+            slice_name = f"slice_of_{name}" if name is not None else None
+            mesh_slice, slice_masks = meshes[0].first_slice.join_meshes(
+                *(mesh.first_slice for mesh in meshes[1:]),
+                name=slice_name,
+                return_masks=True
+            )
+            joined = AxialSymmetricMesh(
+                mesh_slice,
+                axis=meshes[0].axis,
+                nb_repetitions=len(meshes[0]) - 1,
+                name=name
+            )
+            masks = [np.concatenate([
+                slice_mask for _ in range(len(meshes[0]))
+            ]) for slice_mask in slice_masks]
+            return joined, masks
 
     @inplace_transformation
     def translate(self, vector):

capytaine/tools/symbolic_multiplication.py

@@ -11,7 +11,6 @@ output of the form `SymbolicMultiplication("0", np.array(...))`
 import numpy as np
 from functools import wraps, total_ordering
 
-@total_ordering
 class SymbolicMultiplication:
     def __init__(self, symbol, value=1.0):
         self.symbol = symbol
@@ -82,11 +81,26 @@ class SymbolicMultiplication:
     def __getitem__(self, item):
         return SymbolicMultiplication(self.symbol, self.value[item])
 
-    def __eq__(self, x):
-        return float(self) == x
+    def __setitem__(self, item, val):
+        if isinstance(val, SymbolicMultiplication) and self.symbol == val.symbol:
+            self.value.__setitem__(item, val.value)
+        else:
+            raise NotImplementedError
 
     def __lt__(self, x):
-        return float(self) < x
+        return self._concretize() < x
+
+    def __le__(self, x):
+        return self._concretize() <= x
+
+    def __eq__(self, x):
+        return self._concretize() == x
+
+    def __ge__(self, x):
+        return self._concretize() >= x
+
+    def __gt__(self, x):
+        return self._concretize() > x
 
     def __hash__(self):
         return hash((self.symbol, self.value))

capytaine/tools/timer.py

@@ -8,22 +8,24 @@ class Timer:
 
     Example
     -------
-    timer = Timer()
-    with timer:
-        sleep(1.0)
+    ::
 
-    print(timer.total)  # 1.0...
+        timer = Timer()
+        with timer:
+            sleep(1.0)
 
-    @timer.wraps_function
-    def my_function():
-        sleep(0.5)
+        print(timer.total)  # 1.0...
 
-    my_function()
-    print(timer.total)  # 1.5...
-    my_function()
-    print(timer.total)  # 2.0...
+        @timer.wraps_function
+        def my_function():
+            sleep(0.5)
 
-    print(timer.timings)  # [1.0, 0.5, 0.5]
+        my_function()
+        print(timer.total)  # 1.5...
+        my_function()
+        print(timer.total)  # 2.0...
+
+        print(timer.timings)  # [1.0, 0.5, 0.5]
     """
 
     def __init__(self, timings=None):

capytaine-2.3.1.dist-info/DELVEWHEEL

@@ -0,0 +1,2 @@
+Version: 1.11.2
+Arguments: ['C:\\Users\\runneradmin\\AppData\\Local\\Temp\\cibw-run-f4_ko6cf\\cp312-win_amd64\\build\\venv\\Scripts\\delvewheel', 'repair', '-w', 'C:\\Users\\runneradmin\\AppData\\Local\\Temp\\cibw-run-f4_ko6cf\\cp312-win_amd64\\repaired_wheel', 'C:\\Users\\runneradmin\\AppData\\Local\\Temp\\cibw-run-f4_ko6cf\\cp312-win_amd64\\built_wheel\\capytaine-2.3.1-cp312-cp312-win_amd64.whl', '--no-mangle-all']

{capytaine-2.3.dist-info → capytaine-2.3.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: capytaine
-Version: 2.3
+Version: 2.3.1
 Summary: Python BEM solver for linear potential flow, based on Nemoh
 Author-Email: Matthieu Ancellin <matthieu.ancellin@mews-labs.com>
 License:                     GNU GENERAL PUBLIC LICENSE
@@ -689,7 +689,7 @@ Requires-Python: >=3.8
 Requires-Dist: numpy>=1.20; python_version >= "3.9"
 Requires-Dist: numpy>=1.24; python_version == "3.8"
 Requires-Dist: scipy
-Requires-Dist: pandas>=1.3
+Requires-Dist: pandas<3,>=1.3
 Requires-Dist: xarray
 Requires-Dist: rich
 Provides-Extra: optional
@@ -698,17 +698,6 @@ Requires-Dist: joblib>=1.3; extra == "optional"
 Requires-Dist: meshio; extra == "optional"
 Requires-Dist: netcdf4; extra == "optional"
 Requires-Dist: vtk; extra == "optional"
-Provides-Extra: more-optional
-Requires-Dist: pygmsh; extra == "more-optional"
-Requires-Dist: gmsh; extra == "more-optional"
-Provides-Extra: test
-Requires-Dist: pytest; extra == "test"
-Requires-Dist: capytaine[optional]; extra == "test"
-Provides-Extra: docs
-Requires-Dist: sphinx; extra == "docs"
-Requires-Dist: sphinx-toolbox; extra == "docs"
-Requires-Dist: sphinxcontrib-proof; extra == "docs"
-Requires-Dist: sphinxcontrib-mermaid; extra == "docs"
 Description-Content-Type: text/markdown
 
 # Capytaine: a linear potential flow BEM solver with Python.

{capytaine-2.3.dist-info → capytaine-2.3.1.dist-info}/RECORD

@@ -1,11 +1,11 @@
-capytaine/__about__.py,sha256=LRZ6r4hFLAiFQOnekKww3qXw0AwGCtRHyfduv0eZKcc,429
-capytaine/__init__.py,sha256=BHnatA7zRAuTQpRuNgHp9cx0WwXmHV1qkT9ev98eUzk,2184
+capytaine/__about__.py,sha256=OCkVmS8l1azFK-SbX6-WSZiZEozXQVZyU5MkvgrZiq0,431
+capytaine/__init__.py,sha256=zj_J7weaxC7WkSUov6RcQftY-JLZUu9iaZ5cdWBibGk,2184
 capytaine/bem/airy_waves.py,sha256=vap1Qmdo-m1lZ-WLj_5TJY2v3_bNJNtJ_04OkJX83fA,3609
 capytaine/bem/engines.py,sha256=8ShZTaYTYUKSQwt_sHwWH8BfBlt1G0RkDM3mSfuyeE4,18511
-capytaine/bem/problems_and_results.py,sha256=j3xdiKZ77fvTMQ8l_JyArTPtcZ284IX1dxgciVbuDZg,27632
-capytaine/bem/solver.py,sha256=e1oHAC8c6paOMpeClg34zA3II90Ao_Rzdx-uHIozr1E,29389
+capytaine/bem/problems_and_results.py,sha256=_DirKZ-Wo6-Bs99Q9YuOB6W77eH1CaEzu7HiRli3-Ok,28030
+capytaine/bem/solver.py,sha256=ha4psz3_B2s2RPNG3F0LQ4nEhdtaY__7Pf3SCxLgFak,29384
 capytaine/bem/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-capytaine/bodies/bodies.py,sha256=kTPL3cyk47erZXaVAi9Z1Oqfz0Cfv_GjS9vZO0KqplI,53456
+capytaine/bodies/bodies.py,sha256=Nm8pJ7qC6EWJWr0MgCCZ10N2Jox_u7P_WM8n4uy3igM,53929
 capytaine/bodies/dofs.py,sha256=CCGNYqilJE6eppyKxiF_712Vc5iCkbbbrRFCZwvKLdY,617
 capytaine/bodies/__init__.py,sha256=1157LRyNMo7TFZOMvlsC7NPPSMflsZVPfMxeLyNUQkU,161
 capytaine/bodies/predefined/cylinders.py,sha256=1F9aLXp5ZzUcdo0_xncIliNIK3KAIFMc2XspjfWPh-k,5747
@@ -13,7 +13,7 @@ capytaine/bodies/predefined/rectangles.py,sha256=4iHUNJO_5oX8Xfhj2zTg4mIbOu7Myab
 capytaine/bodies/predefined/spheres.py,sha256=vHCuB9wq7Fkgx62tYSSpxqo8mgNna49nAQ_bS4GcfAI,2668
 capytaine/bodies/predefined/__init__.py,sha256=oXMF4lq-wIWO5IWtUt4ZywVAGIi926zYTUz2j_iogT4,344
 capytaine/green_functions/abstract_green_function.py,sha256=PjV2qANb71JEK1YI1GPCF4_2UmBZLN7YCKllvMP6ELY,2898
-capytaine/green_functions/delhommeau.py,sha256=CwLaP2r-qaEZeAr0LaLJL7Hdt4cCGLF_QITsQ27jjcA,23860
+capytaine/green_functions/delhommeau.py,sha256=OhVs6LujvB3ohvhxXCPdr91byDP7mU75XYGoVpTjdwE,24057
 capytaine/green_functions/hams.py,sha256=MS04A6ucbxk7tqGcgrihthbz1Uj4jpL1B4LvdN0C-jc,8456
 capytaine/green_functions/__init__.py,sha256=Rvl4Xv0MbD9NOpcnYkiSpeD097l_SQX8C5r-9CvdIZs,109
 capytaine/green_functions/FinGreen3D/.gitignore,sha256=IxYDvaMqbeEuzSKLq9nRgJIH2AtjPp-ugUuPMhlNirE,8
@@ -25,13 +25,13 @@ capytaine/green_functions/FinGreen3D/test_program.f90,sha256=JYDnew3326C_hCZISpW
 capytaine/green_functions/LiangWuNoblesse/.gitignore,sha256=IxYDvaMqbeEuzSKLq9nRgJIH2AtjPp-ugUuPMhlNirE,8
 capytaine/green_functions/LiangWuNoblesse/LiangWuNoblesseWaveTerm.f90,sha256=BYaVrslAbglcI-7-kOz7qH3AQKLY4rvCmXaa7zfDCaM,22552
 capytaine/green_functions/LiangWuNoblesse/LICENSE,sha256=f_4ZVFh8d9-6HPjrmy6nQ2cfpuY_nnovJYEZ1C4U7v4,27030
-capytaine/green_functions/LiangWuNoblesse/Makefile,sha256=25uL3q1VNE9ladcJf9CWMHJ1j2nm8klpKQ-uFEUWctQ,431
+capytaine/green_functions/LiangWuNoblesse/Makefile,sha256=tnPiOa2ZWD5db2cb9fJA8sqB1tcAtLN4UiV7ElPFpOQ,427
 capytaine/green_functions/LiangWuNoblesse/README.md,sha256=mrO9uR0UcNbDOXRjNwahU7AIto-DXHXmNiA5gJ1jKpk,104
 capytaine/green_functions/LiangWuNoblesse/test_program.f90,sha256=xBkCHsS17Mik7IKgTKRTZNz4EwQcPgnJM1WSmopEWqQ,522
-capytaine/green_functions/libs/Delhommeau_float32.cp312-win_amd64.dll.a,sha256=r4wFk6iDdip4yxlHygfmJqe6p0guEyOaQlgwjTOI0eo,1886
-capytaine/green_functions/libs/Delhommeau_float32.cp312-win_amd64.pyd,sha256=zIjB4QCK0bYLHphVVU8SPimtJ2qSowmtj2abcG-DYHc,424875
-capytaine/green_functions/libs/Delhommeau_float64.cp312-win_amd64.dll.a,sha256=42MxooO3YY0nvXEcAegQec1jbTU71E5tnqWgOcTXg1g,1886
-capytaine/green_functions/libs/Delhommeau_float64.cp312-win_amd64.pyd,sha256=Wxprjr_sm1_CQIhJiS-17wYlEsDVgiaQgq5W6I1KDEU,483290
+capytaine/green_functions/libs/Delhommeau_float32.cp312-win_amd64.dll.a,sha256=V9sEqHh7XTVEOlflKiG_0n8ehCcUyox9bCUXxWrU89w,1886
+capytaine/green_functions/libs/Delhommeau_float32.cp312-win_amd64.pyd,sha256=2LgPWI_EDaRGUETyJDZiN8yWmJdRaZKL-2rMP2jT2F8,424875
+capytaine/green_functions/libs/Delhommeau_float64.cp312-win_amd64.dll.a,sha256=kFZZDMyRRLLoB3Vl7TvYpzrJkDXrUdUilcLox2etZ0o,1886
+capytaine/green_functions/libs/Delhommeau_float64.cp312-win_amd64.pyd,sha256=G2CdpA96yg4S-4TFEevcpFTMFxxUSceqzXBes8CmS4w,483290
 capytaine/green_functions/libs/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 capytaine/io/bemio.py,sha256=_UKSNaI-uuT3ICSEpBRSHhWqN3cFjuBab2t-EJLRtxU,7487
 capytaine/io/legacy.py,sha256=YwVxJKzofVyTshOcb84wr1l0Q2W4USeJNKau5PSBx5o,15058
@@ -39,7 +39,7 @@ capytaine/io/meshio.py,sha256=bu575xBYbl_7VOUwJz7cvtEPe1twYud2hZDWlbX3jjg,1507
 capytaine/io/mesh_loaders.py,sha256=5Zb33dgfEZ2B-3SzbMML__N8yGfpbXDDgwuh6Fd561M,32857
 capytaine/io/mesh_writers.py,sha256=pk88khE-8uXFJj1jmy80C6-saVn8-B7GohSmCWtco9k,21776
 capytaine/io/wamit.py,sha256=ZCYQBPa2vBWLadtZP9Gkvrq_GpYdhxs4ekPLFDHTHRQ,16507
-capytaine/io/xarray.py,sha256=g0wNvPqvfP6ZqPscKmqDCmyafBdZg5HapXsfEMa67MA,27577
+capytaine/io/xarray.py,sha256=9vXThds5aWpfdwBsbXJyXe7-pC-ghmNUYkHlhJGKq_U,27725
 capytaine/io/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 capytaine/matrices/block.py,sha256=Kqsyx1NkFyvwGuNB5Ws-UW5tcMa6pAYDfs6BfiVO0Tw,26078
 capytaine/matrices/block_toeplitz.py,sha256=2G6Hg3cY6nhAsGPpz84utaXKPd0rALdRtGH9MQyd7LQ,13255
@@ -48,15 +48,15 @@ capytaine/matrices/linear_solvers.py,sha256=yobnf7afu7a5Sgby9z8ZcHAjSRlQxVQj7a4Q
 capytaine/matrices/low_rank.py,sha256=ftWVq35dP-9fNZcqMU6FOXyQcrg_3-0u3ba7g1oW0Oo,16391
 capytaine/matrices/__init__.py,sha256=BYCiuVbu4sp3Z6mdbK5uTphUTalvaLhHX1aPLljlsRA,720
 capytaine/meshes/clipper.py,sha256=3C1vxe6mlAYdv0uvtvc3ZA_Lb43ks7UZmDZ1mE9ImRc,19142
-capytaine/meshes/collections.py,sha256=8El0zPaUIWbGu_hoqo_CNIv_rbuLfO77AaxNmTr3R7c,11423
+capytaine/meshes/collections.py,sha256=qD_guE5cTAQKvF6wBbtVyyatVa3wwWvf8mJI2qVFThA,11745
 capytaine/meshes/geometry.py,sha256=QSowVTd4dva0pe4osfFZj0giaL-7Dj-T3wdVZ4k-Clg,14516
-capytaine/meshes/meshes.py,sha256=-avX1RpGhCp2rLyNtuC7kSBSivjpaguV_K_9x4LfUN8,33165
-capytaine/meshes/mesh_like_protocol.py,sha256=ACIwdvJO_qK6wwvbXUT8s5EoHAIUWWJNoBpDrR_3WCY,1066
+capytaine/meshes/meshes.py,sha256=irSdVC2cG68d-DcM56cFN3BmX4-k-4hcpCqTj7kO9aw,33515
+capytaine/meshes/mesh_like_protocol.py,sha256=QdL4uYjfq4-DckrJfE_8PiEhTUAjilJwtW1IeaDKBGg,1079
 capytaine/meshes/properties.py,sha256=zNnQCTAMReiJDMguU9pPwca_Zj_NDTTRqT-uVgABI0E,11066
 capytaine/meshes/quadratures.py,sha256=j6waiQlEoRfzUsnTL2dQg9Ynq5WYkwZT7vV8KOCwN4E,3307
 capytaine/meshes/quality.py,sha256=YVS294DJCTEyB9qGXEP56h9yqxjWIl6hg8P0QI2msIo,14831
 capytaine/meshes/surface_integrals.py,sha256=7s_LOtyz2I1SFFsUE1YMoeSR-02MsocRb8AJS6hTEjc,2326
-capytaine/meshes/symmetric.py,sha256=_yQ8iRiCooFtvH6ryoOUYv0nIbutYiz5atU2rD8USIM,16089
+capytaine/meshes/symmetric.py,sha256=fFN0Pv8BxixSJRatl7KJBSE2RN7i8GAFugvRyQ79vUM,18701
 capytaine/meshes/__init__.py,sha256=FKeQT_I3TAb9iHNqJi0qf9HPo_O001P4fXhZQZ8Oj50,367
 capytaine/meshes/predefined/cylinders.py,sha256=r2UC8uMb60QyyUBaMEvCEkJ9EmazZvDnO-p0GsS2Us8,15160
 capytaine/meshes/predefined/rectangles.py,sha256=wX0vtY8YLMmJ3UZk2oi7q5eKakDnULCXxMhO1KC8eC8,12555
@@ -73,8 +73,8 @@ capytaine/tools/lists_of_points.py,sha256=OOCZkFBsPZj25KGYAvX2BhVLNnxbR-1lT3BSoE
 capytaine/tools/lru_cache.py,sha256=OH1J6XhXyKQGKfDhngBne_D_C2Jub3g92Y5JWWcN4ik,1758
 capytaine/tools/optional_imports.py,sha256=L2_6h1My2jjl7GHL0_hmHsjOE8hRWFrs0eiq_OIDTH8,887
 capytaine/tools/prony_decomposition.py,sha256=_iqHEhop8oBREz0zY35WXV-UkPw1QHbQhP6L_BrtCQM,5297
-capytaine/tools/symbolic_multiplication.py,sha256=MQLsbAfQcWEn7YgcSNRhBmHMRRXY_8wA9Okt4IGaDmI,4560
-capytaine/tools/timer.py,sha256=Ch47sRiDcOF-EcVrdM7Eyv4Lch2wqf2eQpIXNgHENMY,1528
+capytaine/tools/symbolic_multiplication.py,sha256=Aloloo9KEaPaQgmqUEFcYcTHQqy_gQSZHQsMJAMR5LA,4991
+capytaine/tools/timer.py,sha256=VOte4SjqgvjUGiwhi3YalifFclZe2X_6VLvuecKofD4,1586
 capytaine/tools/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 capytaine/ui/cli.py,sha256=DLgsJdmMASM70BqN7bbNQxBkiHKReJDJqgf4E9mr_yk,825
 capytaine/ui/rich.py,sha256=kxI4CVa4ycJBloEHvLpLqrEArkLT4wAdcPKhbEw2hcY,238
@@ -84,12 +84,12 @@ capytaine/ui/vtk/body_viewer.py,sha256=1WGMjHBRhyGkCFxgvJKs_dY-eL2YCCKhqqehKgIww
 capytaine/ui/vtk/helpers.py,sha256=V1_DeB6svtTlzdNzdcL7-GGj5QvCW6wpVzhXVmpVlwI,2722
 capytaine/ui/vtk/mesh_viewer.py,sha256=V6ReLbsArkVAZRX6QZdqHhyo1dStT2ql8T2AzpP8HbA,15211
 capytaine/ui/vtk/__init__.py,sha256=mtXE0si8Fc8QGneQIQNdGPPzhRLNwi4bPMUmQ9Vs1io,164
-capytaine-2.3.dist-info/DELVEWHEEL,sha256=wLqEmWBECLItOZ_kx9BJQv6mBGXyIiOFTVa8M6AbvU0,419
-capytaine-2.3.dist-info/entry_points.txt,sha256=R72-je8lc6ELm8ftt7lJ7f1aalnQs5BWYrGDBMexHII,53
-capytaine-2.3.dist-info/LICENSE,sha256=Czg9WmPaZE9ijZnDOXbqZIftiaqlnwsyV5kt6sEXHms,35821
-capytaine-2.3.dist-info/METADATA,sha256=NSFNwYHZoOfMqn7if-xXY-r9qOLAGrlbUSG8qtlahy8,45516
-capytaine-2.3.dist-info/RECORD,,
-capytaine-2.3.dist-info/WHEEL,sha256=50PeAbplA6PkI0hYOYoeacB9US1R6EguyfOnsccH0WU,85
+capytaine-2.3.1.dist-info/DELVEWHEEL,sha256=BIx_Bwuob1wkJAAid4pQ5MgSYbgcHYpURgKHT2fXRCE,421
+capytaine-2.3.1.dist-info/entry_points.txt,sha256=R72-je8lc6ELm8ftt7lJ7f1aalnQs5BWYrGDBMexHII,53
+capytaine-2.3.1.dist-info/LICENSE,sha256=Czg9WmPaZE9ijZnDOXbqZIftiaqlnwsyV5kt6sEXHms,35821
+capytaine-2.3.1.dist-info/METADATA,sha256=ftEZnZzzOAlQ85oSpm66yuOls3Sh6ypK-bEYP3EcKME,45072
+capytaine-2.3.1.dist-info/RECORD,,
+capytaine-2.3.1.dist-info/WHEEL,sha256=50PeAbplA6PkI0hYOYoeacB9US1R6EguyfOnsccH0WU,85
 capytaine.libs/libgcc_s_seh-1.dll,sha256=q2m46g2pAq7jvn1d0qbVW6WqyOpxQ6KtnAdh5r3gCxY,111616
 capytaine.libs/libgfortran-5.dll,sha256=wLLASxjhj_hbogJghmWeKSyi6w4L9h8xR04-iS1Ejq4,3118080
 capytaine.libs/libgomp-1.dll,sha256=yCJD_z8GQ7Hi7iXgwG2hkjFAOqf-KeE45YnavaKCpf8,280576

capytaine-2.3.dist-info/DELVEWHEEL

@@ -1,2 +0,0 @@
-Version: 1.10.1
-Arguments: ['C:\\Users\\runneradmin\\AppData\\Local\\Temp\\cibw-run-9rerwt60\\cp312-win_amd64\\build\\venv\\Scripts\\delvewheel', 'repair', '-w', 'C:\\Users\\runneradmin\\AppData\\Local\\Temp\\cibw-run-9rerwt60\\cp312-win_amd64\\repaired_wheel', 'C:\\Users\\runneradmin\\AppData\\Local\\Temp\\cibw-run-9rerwt60\\cp312-win_amd64\\built_wheel\\capytaine-2.3-cp312-cp312-win_amd64.whl', '--no-mangle-all']