capytaine 2.3__cp38-cp38-macosx_13_0_x86_64.whl → 2.3.1__cp38-cp38-macosx_13_0_x86_64.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/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.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,5 +1,5 @@
 capytaine/__init__.py,sha256=ffYPfnzVMU03arvItDaf5VoDa68KEGP4J9aE1MQlmQg,1822
-capytaine/__about__.py,sha256=di-zISnbYXwp9FpVIxlF_mgV6Ind8W8YBuMiy82ypCA,413
+capytaine/__about__.py,sha256=14fwPZt3SRhKmIUSRPB6gpUnkNBLGgZLmq6_R2Xvy_k,415
 capytaine/ui/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 capytaine/ui/rich.py,sha256=ZfMHUpCSVztNseIGpN2I69ts_Aw_HJ5RnZc0f73av5Q,233
 capytaine/ui/cli.py,sha256=o3MH2OCHXFww1kb2zM9R3qs4jjfkCeeZFRAgnFnb2NE,797
@@ -9,7 +9,7 @@ capytaine/ui/vtk/animation.py,sha256=FYH-YAu0LY2zr7fTl36aOSkM8z74_gxQRKOJFhUSluo
 capytaine/ui/vtk/helpers.py,sha256=sOK_a679a7n9ozhcEZY3IhTWmPyctoLkFwKTyzDHrm8,2640
 capytaine/ui/vtk/body_viewer.py,sha256=gGIla-gADWiYuiNdS2io-0whpO6AaH8A8KkEmqX4OSY,1032
 capytaine/green_functions/__init__.py,sha256=527Pd0SIqwR64-d4sT-_jyE89XU-PLgFL0et96ChAIo,107
-capytaine/green_functions/delhommeau.py,sha256=Z3lXv4MP1MPPJKmyISMVDqTQIpEHdroAl-OUjYJVE2w,23357
+capytaine/green_functions/delhommeau.py,sha256=y9nkyVKj9ES0nL2jDLdx_5ejQ5nbMSKV0qYUnPnfTNg,23550
 capytaine/green_functions/hams.py,sha256=yxT0ysNFxsm4VbRtywHnhtWnHrxkU8q6_DL8eEwZPTs,8252
 capytaine/green_functions/abstract_green_function.py,sha256=6_YjjTKbYNgteWw4FWKQQI9nfpPzEETQxD3WOWFqosw,2834
 capytaine/green_functions/FinGreen3D/test_program.f90,sha256=BlQv3O0z0tYXejz261S8l7x5H01wHAdn9-YfsRurW1w,926
@@ -20,15 +20,15 @@ capytaine/green_functions/FinGreen3D/FinGreen3D.f90,sha256=tBsSc6BfQJpaFQ6j6tEYW
 capytaine/green_functions/FinGreen3D/.gitignore,sha256=GBMUBl3y8v2vkgsai1MR2qIWotZImgatpbScxRTYlBc,7
 capytaine/green_functions/LiangWuNoblesse/test_program.f90,sha256=cettkRXfJDmhigo7oUACjU3qCNFLJVDU0vI6Og4UFs8,494
 capytaine/green_functions/LiangWuNoblesse/LICENSE,sha256=IMF9i4xIpgCADf0U-V1cuf9HBmqWQd3qtI3FSuyW4zE,26526
-capytaine/green_functions/LiangWuNoblesse/Makefile,sha256=S0zjKvMQ9z2nWHEKNe-ES32LWyNI1_ogcJ6s9qVFEho,413
+capytaine/green_functions/LiangWuNoblesse/Makefile,sha256=txvf2ppXhnT9OPDcE8zQcEq5jQcQ3B_A0ge1vdR6huA,411
 capytaine/green_functions/LiangWuNoblesse/README.md,sha256=0pQbZxIAWTc-9eQ2Qb62SfM5Eh3dNx0bPSznuOCx8xU,102
 capytaine/green_functions/LiangWuNoblesse/.gitignore,sha256=GBMUBl3y8v2vkgsai1MR2qIWotZImgatpbScxRTYlBc,7
 capytaine/green_functions/LiangWuNoblesse/LiangWuNoblesseWaveTerm.f90,sha256=BYaVrslAbglcI-7-kOz7qH3AQKLY4rvCmXaa7zfDCaM,22552
-capytaine/green_functions/libs/Delhommeau_float64.cpython-38-darwin.so,sha256=20MVKC6HB72LdbrEUE0ok0SGQGUBsTCIcTJLgMswOsc,465904
-capytaine/green_functions/libs/Delhommeau_float32.cpython-38-darwin.so,sha256=1XYJ6lO6Kr5TzB5huYLPYWZXYwdEESju1gFzurSFRNU,412800
+capytaine/green_functions/libs/Delhommeau_float64.cpython-38-darwin.so,sha256=ig5Ftsg-9km_wABbmDelvFjrZ4sBq73LA5f5-pce4NI,465904
+capytaine/green_functions/libs/Delhommeau_float32.cpython-38-darwin.so,sha256=UXpvw50jmadyl7QIQW4t5-xgqxsQFvQab2RISqUMuZg,412800
 capytaine/green_functions/libs/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 capytaine/bodies/dofs.py,sha256=M5nWTeVJxmxKl2BowOnrE3plv1rrqDHNb76D_qHroS8,598
-capytaine/bodies/bodies.py,sha256=w_2kx-8Hm9qDWI9Zl-RsQQQ9RHAqAZ0lldxSseL4OFw,52243
+capytaine/bodies/bodies.py,sha256=KOBeLwpXDBQhZAjLePqXMQf8LcAl2TdOtYS_4rzIjMs,52708
 capytaine/bodies/__init__.py,sha256=wIYFs14NAWqaHLQWvDGUEZTBWp_kdthKkbOlXHup6LA,157
 capytaine/bodies/predefined/rectangles.py,sha256=PtNK7_RUjRIi1oJNlsBhOosC3Ib2cXX1SnW9llTmD44,4500
 capytaine/bodies/predefined/spheres.py,sha256=8lHMqCQBQneN8EqTSXQanumt6j5TeR70CR4ckrnDfo4,2598
@@ -36,17 +36,17 @@ capytaine/bodies/predefined/__init__.py,sha256=IA1y5XDGmyzNDpgQvz8kjyGajBCYoClAn
 capytaine/bodies/predefined/cylinders.py,sha256=Zfp3ecpKD72dfzv5Dd2y87EvDBBRsOJPq6GR6B76r5Y,5596
 capytaine/tools/lru_cache.py,sha256=ssPo7sZMMe15GIxPiQ7rA0GrjwYWsE66x3A-d-4FQTI,1709
 capytaine/tools/cache_on_disk.py,sha256=Cj2J5pNKw9PbIvxUn53swSmuGJoMibZzIvB5Skt52ig,912
-capytaine/tools/timer.py,sha256=3-PXNn6lMfbPP6z0HNCE9YcnIdrepzMHFd_CZ1F4aiM,1464
+capytaine/tools/timer.py,sha256=0kEIvYFnwa0F4ayecxLl_VtiCrVL-8i9_NrAKrlSVmo,1520
 capytaine/tools/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 capytaine/tools/deprecation_handling.py,sha256=_QFNrkreLerm3-LzRZ3l6B_O249PM-fCCAXMsHA7kIU,809
 capytaine/tools/optional_imports.py,sha256=HNjtvNZGxsw-eYmr4ejq001KhLevrvqq7ASAqTOw8W8,860
-capytaine/tools/symbolic_multiplication.py,sha256=Ef7oJftyAyw7kbPl5F4xD3YGDTWQFVNrdYCBii5e9po,4425
+capytaine/tools/symbolic_multiplication.py,sha256=DCmVlUwGtyUBU7VZZ0xxNnIjNSC0rQbgSuHRnlFQ8VU,4842
 capytaine/tools/lists_of_points.py,sha256=6DyL9kr8jbkvB2IuCq288eSd4Xq6pwa8wiQCYZm9M1Q,1828
 capytaine/tools/prony_decomposition.py,sha256=Eg8CcedPuLKVUA-wpjwXBJsnXjBGZruHVy4jMU2cCSM,5147
-capytaine/bem/solver.py,sha256=KvDzuf4fOnUiklnBMutSWS-IoEd5ST-9BEXkVDdGiD0,28795
+capytaine/bem/solver.py,sha256=j2ZoVGA8x46xQe61s1Dh0zCFVXnTDuehAKa3f3-bZ0o,28790
 capytaine/bem/engines.py,sha256=YauXkOHKLbrDbzM1qisKiJ1D3LZ0-ALOWIIRNzNchUU,18070
 capytaine/bem/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-capytaine/bem/problems_and_results.py,sha256=X5va5A8dMCpXeuhAdDr_OWLYkI_KcerhLq99z-qkS20,27040
+capytaine/bem/problems_and_results.py,sha256=4pdwswnrbp7D1v1xOnBA3aYSCxm5qiaFZDJt_u5BxaA,27430
 capytaine/bem/airy_waves.py,sha256=GAhX5xNEgjoOo2uEHXfKCxb2TnmbPBCqW6VWfvzJTM8,3498
 capytaine/io/wamit.py,sha256=dvmthwIZnzOM8L1kq814pilCpQPLeTHvjjgg21hV2po,16028
 capytaine/io/legacy.py,sha256=lARoaPGjTsUQ4YiLmD8ZYMSCU_UDmnOQx4_JppaCcRc,14730
@@ -55,18 +55,18 @@ capytaine/io/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 capytaine/io/mesh_loaders.py,sha256=NjSWrCR20ozIFu94B7ZveiSS1kzt-IWKoFnkvLHcNXU,31771
 capytaine/io/meshio.py,sha256=QmuHGJt_a1VDrCpA_OKVWftgHPdZyjNkao38QpGuPT8,1469
 capytaine/io/mesh_writers.py,sha256=zLHRdzTpbpXyPsf7c09s6Nk2ZXIAQCCfHYIF3vKKDBo,21084
-capytaine/io/xarray.py,sha256=LdbWiLAYqN9uz-TaYSqhVfSAKUnFnt2NmbJA4XUoO3Y,26909
-capytaine/.dylibs/libgfortran.5.dylib,sha256=vrn9LxAxYhy74BVfINa1ANs4zkvZq3LJR583KM3eBV4,3497360
-capytaine/.dylibs/libquadmath.0.dylib,sha256=yKg9ePtGRSTCvOoq-MKeXyL-KwJ5fFc7ber_FvOWU9I,381088
-capytaine/.dylibs/libgcc_s.1.1.dylib,sha256=TkxUgIZRQtgSBW5wG3sHiwurvUOrRHsQElu6_LTXlDI,244384
+capytaine/io/xarray.py,sha256=IF-l5O9gGeBTRpcJontvvi11sYfgqwfK7wp1_MAtFHY,27057
+capytaine/.dylibs/libgfortran.5.dylib,sha256=58lGzOQAJ8maiYrLh0oaYocIfJMrLQhHzo0HaPuJp54,3497360
+capytaine/.dylibs/libquadmath.0.dylib,sha256=aU-tHcHpGTPL1MkPenAw_94eJx2POBThXvwMvrIaubE,381088
+capytaine/.dylibs/libgcc_s.1.1.dylib,sha256=8NDI5rBKk_mhvebcgJkLy3WpXbzgBdVPMiqlgN4RrYo,244384
 capytaine/meshes/clipper.py,sha256=CB38_cmha8q6hnaxKukKQPtB7sXx34fU2HBujF8hD24,18677
 capytaine/meshes/properties.py,sha256=eJHhTtXzW-aiokuvluqLSfumAGMFZGmI3ZLP-pfVh2I,10790
-capytaine/meshes/symmetric.py,sha256=kK6eELnByoZZmkpziiJ_EJsOFc89MMQqAIG810V5VoI,15697
-capytaine/meshes/mesh_like_protocol.py,sha256=9D0PiJHPThc-eezgtJ2FBhxkw50d2sMIHLowA7LVUew,1029
+capytaine/meshes/symmetric.py,sha256=CHBusSwTufKiGYOqVrUurQojYqlLsSWqxB-Z2kjVEYE,18239
+capytaine/meshes/mesh_like_protocol.py,sha256=BxFLfJsHZFobyN5L0siuPRPq0DFpJbGsQw1yuFN06GQ,1042
 capytaine/meshes/surface_integrals.py,sha256=ihxbuiSsgep4rbrBQlq2q_2p_uDXDkRPyz0S6OVofc8,2263
 capytaine/meshes/__init__.py,sha256=5-AE32KOmOcq_F09uZanLafy7SzOe0-Cl6PaDbbOLMQ,361
-capytaine/meshes/meshes.py,sha256=GgxXIeVDxrZJ1G9CQqPgGpNToXWs4MJ9z_A_wSZGHrk,32284
-capytaine/meshes/collections.py,sha256=KAnGlw5VQZtAlKj15O_8jw0WUlrYmKiFgUkFVN8HYn8,11089
+capytaine/meshes/meshes.py,sha256=_b0DtxR_UDmV0UOuB3bZ47ClGkIKJynTWESJxiIsbio,32625
+capytaine/meshes/collections.py,sha256=uTXQovnoPsKZfl3HqUveHc94SyscyfELILikRuAcpG0,11403
 capytaine/meshes/geometry.py,sha256=ffeHvFzke3sXHzGexivRmSBQIOCk-pCMoveQF4eYyGM,14107
 capytaine/meshes/quality.py,sha256=OTZvqyDSJGKIHfFWwn-QROH1Je4PjgFxbwp3Ce-nzzI,14383
 capytaine/meshes/quadratures.py,sha256=eT5tsfonBJEDkrCnoUJbEPLX34I1rCkgzD_E6XhDu-E,3227
@@ -85,8 +85,8 @@ capytaine/matrices/linear_solvers.py,sha256=_5R-CR7eJ07Z3GcOd1YiX2bRyukuXQnHE0Dg
 capytaine/matrices/block_toeplitz.py,sha256=yL4GFPghfdVZTKtden0RKVDhpiBZz3D1oZ-pAVjDBjg,12930
 capytaine/matrices/__init__.py,sha256=_b868RYRI8ANHHf2PUVI-2_c0RdBfRDHLvJW5Uq1qys,704
 capytaine/matrices/block.py,sha256=9Gv800-regvNmD62-h9MvbW7vtPb4N-_3o8Y8HPnxgc,25486
-capytaine-2.3.dist-info/RECORD,,
-capytaine-2.3.dist-info/LICENSE,sha256=jOtLnuWt7d5Hsx6XXB2QxzrSe2sWWh3NgMfFRetluQM,35147
-capytaine-2.3.dist-info/WHEEL,sha256=GvSWRiUfFrukiVvOyV_xTozA93yOU6whvfgQHhrp_b0,92
-capytaine-2.3.dist-info/entry_points.txt,sha256=R72-je8lc6ELm8ftt7lJ7f1aalnQs5BWYrGDBMexHII,53
-capytaine-2.3.dist-info/METADATA,sha256=NSFNwYHZoOfMqn7if-xXY-r9qOLAGrlbUSG8qtlahy8,45516
+capytaine-2.3.1.dist-info/RECORD,,
+capytaine-2.3.1.dist-info/LICENSE,sha256=jOtLnuWt7d5Hsx6XXB2QxzrSe2sWWh3NgMfFRetluQM,35147
+capytaine-2.3.1.dist-info/WHEEL,sha256=GvSWRiUfFrukiVvOyV_xTozA93yOU6whvfgQHhrp_b0,92
+capytaine-2.3.1.dist-info/entry_points.txt,sha256=R72-je8lc6ELm8ftt7lJ7f1aalnQs5BWYrGDBMexHII,53
+capytaine-2.3.1.dist-info/METADATA,sha256=ftEZnZzzOAlQ85oSpm66yuOls3Sh6ypK-bEYP3EcKME,45072