bsplyne 1.0.0__py3-none-any.whl

bsplyne/save_utils.py

@@ -0,0 +1,141 @@
+import numpy as np
+import meshio as io
+from typing import Iterable
+from functools import reduce
+import os
+import xml
+
+
+def writePVD(fileName: str, groups: dict[str, dict]):
+    """
+    Writes a Paraview Data (PVD) file that references multiple VTK files.
+
+    Creates an XML-based PVD file that collects and organizes multiple VTK files into groups,
+    allowing visualization of multi-part, time-series data in Paraview.
+
+    Parameters
+    ----------
+    fileName : str
+        Base name for the mesh files (without numbers and extension)
+    groups : dict[str, dict]
+        Nested dictionary specifying file groups with format:
+        {"group_name": {"ext": "file_extension", "npart": num_parts, "nstep": num_timesteps}}
+
+    Notes
+    -----
+    VTK files must follow naming pattern: {fileName}_{group}_{part}_{timestep}.{ext}
+    Example: for fileName="mesh", group="fluid", part=1, timestep=5, ext="vtu":
+            mesh_fluid_1_5.vtu
+
+    Returns
+    -------
+    None
+    """
+    rep, fname = os.path.split(fileName)
+    pvd = xml.dom.minidom.Document()
+    pvd_root = pvd.createElementNS("VTK", "VTKFile")
+    pvd_root.setAttribute("type", "Collection")
+    pvd_root.setAttribute("version", "0.1")
+    pvd_root.setAttribute("byte_order", "LittleEndian")
+    pvd.appendChild(pvd_root)
+    collection = pvd.createElementNS("VTK", "Collection")
+    pvd_root.appendChild(collection)
+    for name, grp in groups.items():
+        for jp in range(grp["npart"]):
+            for js in range(grp["nstep"]):
+                dataSet = pvd.createElementNS("VTK", "DataSet")
+                dataSet.setAttribute("timestep", str(js))
+                dataSet.setAttribute("group", name)
+                dataSet.setAttribute("part", str(jp))
+                dataSet.setAttribute("file", f"{fname}_{name}_{jp}_{js}.{grp['ext']}")
+                dataSet.setAttribute("name", f"{name}_{jp}")
+                collection.appendChild(dataSet)
+    outFile = open(fileName + ".pvd", "w")
+    pvd.writexml(outFile, newl="\n")
+    print("VTK: " + fileName + ".pvd written")
+    outFile.close()
+
+
+def merge_meshes(meshes: Iterable[io.Mesh]) -> io.Mesh:
+    """
+    Merges multiple meshio.Mesh objects into a single mesh.
+
+    Parameters
+    ----------
+    meshes : Iterable[io.Mesh]
+        An iterable of meshio.Mesh objects to merge.
+
+    Returns
+    -------
+    io.Mesh
+        A single meshio.Mesh object containing the merged meshes with combined
+        vertices, cells and point data.
+    """
+    vertices = np.vstack([m.points for m in meshes])
+    all_cell_types = reduce(lambda a, b: a & b, [m.cells_dict.keys() for m in meshes])
+    cells = {}
+    for cell_type in all_cell_types:
+        counter = 0
+        cell = []
+        for m in meshes:
+            if cell_type in m.cells_dict:
+                cell.append(m.cells_dict[cell_type] + counter)
+            counter += m.points.shape[0]
+        cells[cell_type] = np.vstack(cell)
+    point_data_names = reduce(lambda a, b: a & b, [m.point_data.keys() for m in meshes])
+    point_datas = {}
+    for point_data_name in point_data_names:
+        for m in meshes:
+            if point_data_name in m.point_data:
+                shape = m.point_data[point_data_name].shape[1:]  # type: ignore
+                break
+        point_data = []
+        for m in meshes:
+            if point_data_name in m.point_data:
+                point_data.append(m.point_data[point_data_name])
+            else:
+                point_data.append(np.zeros((m.points.shape[0], *shape)))
+        point_datas[point_data_name] = np.concatenate(point_data, axis=0)
+    return io.Mesh(vertices, cells, point_data=point_datas)
+
+
+def merge_saves(
+    path: str, name: str, nb_patchs: int, nb_steps: int, group_names: list[str]
+) -> None:
+    """
+    Merge multiple mesh files and save the merged results.
+
+    This function reads multiple mesh files for each group and time step,
+    merges them into a single mesh, and writes the merged mesh to a new file.
+    It also generates a PVD file to describe the collection of merged meshes.
+
+    Parameters
+    ----------
+    path : str
+        The directory path where the mesh files are located.
+    name : str
+        The base name of the mesh files.
+    nb_patchs : int
+        The number of patches to merge for each group and time step.
+    nb_steps : int
+        The number of time steps for which meshes are available.
+    group_names : list[str]
+        A list of group names to process.
+
+    Returns
+    -------
+    None
+    """
+    filename = os.path.join(path, name)
+    filename_out = filename + "_merged"
+    groups = {}
+    for group_name in group_names:
+        groups[group_name] = {"ext": "vtu", "npart": 1, "nstep": nb_steps}
+        for step in range(nb_steps):
+            merge_meshes(
+                [
+                    io.read(f"{filename}_{group_name}_{p}_{step}.vtu")
+                    for p in range(nb_patchs)
+                ]
+            ).write(f"{filename_out}_{group_name}_0_{step}.vtu")
+    writePVD(filename_out, groups)

bsplyne-1.0.0.dist-info/METADATA

@@ -0,0 +1,91 @@
+Metadata-Version: 2.4
+Name: bsplyne
+Version: 1.0.0
+Summary: A package for N dimensions B-splines.
+Home-page: https://github.com/Dorian210/bsplyne
+Author: Dorian Bichet
+Author-email: dbichet@insa-toulouse.fr
+Classifier: Programming Language :: Python :: 3
+Classifier: Operating System :: OS Independent
+Classifier: License :: OSI Approved :: CEA CNRS Inria Logiciel Libre License, version 2.1 (CeCILL-2.1)
+Description-Content-Type: text/markdown
+License-File: LICENSE.txt
+Requires-Dist: numpy
+Requires-Dist: numba
+Requires-Dist: scipy
+Requires-Dist: matplotlib
+Requires-Dist: meshio
+Requires-Dist: tqdm
+Dynamic: author
+Dynamic: author-email
+Dynamic: classifier
+Dynamic: description
+Dynamic: description-content-type
+Dynamic: home-page
+Dynamic: license-file
+Dynamic: requires-dist
+Dynamic: summary
+
+# bsplyne
+
+<p align="center">
+  <img src=docs/logo.png width="500" />
+</p>
+
+**bsplyne** is a Python library for working with N-dimensional B-splines. It implements the Cox-de Boor algorithm for basis evaluation, order elevation, knot insertion, and provides a connectivity class for multi-patch structures. Additionally, it includes visualization tools with export capabilities to Paraview.
+
+## Installation
+
+Since **bsplyne** is not yet on PyPI, you can install it locally as follows:
+
+```bash
+git clone https://github.com/Dorian210/bsplyne
+cd bsplyne
+pip install -e .
+```
+
+### Dependencies
+Make sure you have the following dependencies installed:
+- `numpy`
+- `numba`
+- `scipy`
+- `matplotlib`
+- `meshio`
+- `tqdm`
+
+## Main Modules
+
+- **BSplineBasis**  
+  Implements B-spline basis function evaluation using the Cox-de Boor recursion formula.
+  
+- **BSpline**  
+  Provides methods for creating and manipulating N-dimensional B-splines, including order elevation and knot insertion.
+  
+- **MultiPatchBSplineConnectivity**  
+  Manages the connectivity between multiple N-dimensional B-spline patches.
+  
+- **CouplesBSplineBorder** (less documented)  
+  Handles coupling between B-spline borders.
+
+## Examples
+
+Several example scripts demonstrating the usage of **bsplyne** can be found in the `examples/` directory. These scripts cover:
+- Basis evaluation on a curved line
+- Plotting with Matplotlib
+- Order elevation
+- Knot insertion
+- Surface examples
+- Exporting to Paraview
+
+## Documentation
+
+The full API documentation is available in the `docs/` directory of the project or via the [online documentation portal](https://dorian210.github.io/bsplyne/).
+
+## Contributions
+
+At the moment, I am not actively reviewing contributions. However, if you encounter issues or have suggestions, feel free to open an issue.
+
+## License
+
+This project is licensed under the [CeCILL License](LICENSE.txt).
+

bsplyne-1.0.0.dist-info/RECORD

@@ -0,0 +1,13 @@
+bsplyne/__init__.py,sha256=sajGBdm1-pu6iMgD-G_nimXrbXZsXwrJT6bjA_JyLO8,1899
+bsplyne/b_spline.py,sha256=CJYIS_kAJFOy1ay7kIOACyX1j20BkT9TRn1t1j8WkGU,110105
+bsplyne/b_spline_basis.py,sha256=J4vrGx-WcW7suZDJzvLr3zDQyGCCLAZlMbzuGeQ1fZ4,35315
+bsplyne/geometries_in_3D.py,sha256=krQPNZEg-PZCTnlhd1YrJjwCwv7R9T7oAxxlZJms_Zs,36519
+bsplyne/multi_patch_b_spline.py,sha256=4lsxS3m7EV_igfzOBjCKaB7Pd8pxI-ZHyLmnQRP4NhQ,75925
+bsplyne/my_wide_product.py,sha256=JqFy2s3ks0wKf4iq22NyzrZo-3Uo7ita-r6fGcGEfEE,5678
+bsplyne/parallel_utils.py,sha256=M0Tuy8wA3AYltsnr_iYOnREWRGZnPxavq9cpWaEC_lU,15326
+bsplyne/save_utils.py,sha256=p-gDHTAU4Fr8iI_yPfdrMHtEF6bOlgkm8b4yHxA5qtw,5023
+bsplyne-1.0.0.dist-info/licenses/LICENSE.txt,sha256=NuJeFieCkXCWPQlUUkzn5TrngzZmXp09-Q8fk15_r6A,3196
+bsplyne-1.0.0.dist-info/METADATA,sha256=5Gb889Dpk_bkCi3Pgds_5EintBhxWTNCREE7C4bNjkY,2690
+bsplyne-1.0.0.dist-info/WHEEL,sha256=wUyA8OaulRlbfwMtmQsvNngGrxQHAvkKcvRmdizlJi0,92
+bsplyne-1.0.0.dist-info/top_level.txt,sha256=Mwj9pTQb41GC_guOzmpvMcdkXeeAIWF6VtZcXghgx0c,8
+bsplyne-1.0.0.dist-info/RECORD,,

bsplyne-1.0.0.dist-info/WHEEL

@@ -0,0 +1,5 @@
+Wheel-Version: 1.0
+Generator: setuptools (80.10.2)
+Root-Is-Purelib: true
+Tag: py3-none-any
+

bsplyne-1.0.0.dist-info/licenses/LICENSE.txt

@@ -0,0 +1,70 @@
+
+  CeCILL FREE SOFTWARE LICENSE AGREEMENT
+
+Version 2.1 dated 2013-06-21
+
+
+    Notice
+
+This Agreement is a Free Software license agreement that is the result
+of discussions between its authors in order to ensure compliance with
+the two main principles guiding its drafting:
+
+  * firstly, compliance with the principles governing the distribution
+    of Free Software: access to source code, broad rights granted to users,
+  * secondly, the election of a governing law, French law, with which it
+    is conformant, both as regards the law of torts and intellectual
+    property law, and the protection that it offers to both authors and
+    holders of the economic rights over software.
+
+The authors of the CeCILL (for Ce[a] C[nrs] I[nria] L[ogiciel] L[ibre]) 
+license are: 
+
+Commissariat à l'énergie atomique et aux énergies alternatives - CEA, a
+public scientific, technical and industrial research establishment,
+having its principal place of business at 25 rue Leblanc, immeuble Le
+Ponant D, 75015 Paris, France.
+
+Centre National de la Recherche Scientifique - CNRS, a public scientific
+and technological establishment, having its principal place of business
+at 3 rue Michel-Ange, 75794 Paris cedex 16, France.
+
+Institut National de Recherche en Informatique et en Automatique -
+Inria, a public scientific and technological establishment, having its
+principal place of business at Domaine de Voluceau, Rocquencourt, BP
+105, 78153 Le Chesnay cedex, France.
+
+
+    Preamble
+
+The purpose of this Free Software license agreement is to grant users
+the right to modify and redistribute the software governed by this
+license within the framework of an open source distribution model.
+
+The exercising of this right is conditional upon certain obligations for
+users so as to preserve this status for all subsequent redistributions.
+
+In consideration of access to the source code and the rights to copy,
+modify and redistribute granted by the license, users are provided only
+with a limited warranty and the software's author, the holder of the
+economic rights, and the successive licensors only have limited liability.
+
+In this respect, the risks associated with loading, using, modifying
+and/or developing or reproducing the software by the user are brought to
+the user's attention, given its Free Software status, which may make it
+complicated to use, with the result that its use is reserved for
+developers and experienced professionals having in-depth computer
+knowledge. Users are therefore encouraged to load and test the
+suitability of the software as regards their requirements in conditions
+enabling the security of their systems and/or data to be ensured and,
+more generally, to use and operate it in the same conditions of
+security. This Agreement may be freely reproduced and published,
+provided it is not altered, and that no provisions are either added or
+removed herefrom.
+
+This Agreement may apply to any or all software for which the holder of
+the economic rights decides to submit the use thereof to its provisions.
+
+Frequently asked questions can be found on the official website of the
+CeCILL licenses family (http://www.cecill.info/index.en.html) for any 
+necessary clarification.

bsplyne-1.0.0.dist-info/top_level.txt

@@ -0,0 +1 @@
+bsplyne