modulo-vki 2.1.0__py3-none-any.whl → 2.1.1__py3-none-any.whl

{modulo_vki-2.1.0.dist-info → modulo_vki-2.1.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: modulo_vki
-Version: 2.1.0
+Version: 2.1.1
 Summary: MODULO (MODal mULtiscale pOd) is a software developed at the von Karman Institute to perform Multiscale Modal Analysis of numerical and experimental data.
 Home-page: https://github.com/mendezVKI/MODULO/tree/master/modulo_python_package/
 Author: ['R. Poletti', 'L. Schena', 'D. Ninni', 'M. A. Mendez']
@@ -36,7 +36,7 @@ Dynamic: summary
 
 
 MODULO: a python toolbox for data-driven modal decomposition
------------------------------------------------------------
+------------------------------------------------------------
 
 .. image:: https://readthedocs.org/projects/modulo/badge/?version=latest
     :target: https://modulo.readthedocs.io/en/latest/?badge=latest
@@ -49,12 +49,11 @@ MODULO: a python toolbox for data-driven modal decomposition
 
 .. |PyPI| image:: https://img.shields.io/pypi/v/modulo_vki
     :target: https://pypi.org/project/modulo_vki/
-
-.. raw:: html
-
-   <div style="text-align: center;">
-       <img src="https://modulo.readthedocs.io/en/latest/_images/modulo_logo.png" alt="Modulo Logo" width="500"/>
-   </div>
+	
+.. image:: https://modulo.readthedocs.io/en/latest/_images/modulo_logo.png
+   :alt: MODULO logo
+   :width: 500px
+   :align: center
 
 
 **MODULO** is a modal decomposition package developed at the von Karman Institute for Fluid Dynamics (VKI). It offers a wide range of decomposition techniques, allowing users to choose the most appropriate method for their specific problem. MODULO can efficiently handle large datasets natively, thanks to a memory-saving feature that partitions the data and processes the decomposition in chunks (ninni2020modulo). Moreover, it supports non-uniform meshes through a weighted inner product formulation. Currently, MODULO heavily relies on NumPy routines and does not offer additional parallel computing capabilities beyond those naturally provided by NumPy.
@@ -70,65 +69,47 @@ These are available at the `MODULO YouTube channel <https://www.youtube.com/@mod
 Modal decompositions
 --------------------
 Modal decompositions aim to describe the data as a linear combination of *modes*, obtained by projecting the data 
-onto a suitable set of basis. For instance, consider a matrix $D(x, t)$, where $x$ and $t$ are the spatial and temporal
-coordinates, respectively, the modal decomposition can be written as:
-
-$D(x_i, t_k) = \\phi(x_i) \\Sigma \\psi(t_k)^T$
-
-where $\\phi(x_i)$ is the spatial basis, $\\psi(t_k)$ is the temporal basis, and $\\Sigma$ is the modal coefficients. 
-Different decompositions employ different bases, such as prescribed Fourier basis ($\\psi_\\mathcal{F}$) for 
+onto a suitable set of basis. Different decompositions employ different bases, such as prescribed Fourier basis for 
 the Discrete Fourier Transform (DFT), or data-driven basis, i.e. tailored on the dataset at hand, 
-for the Proper Orthogonal Decomposition (POD). 
-
-We refer to (mendez2022statistical, mendez2022generalizedmultiscalemodalanalysis, Mendez_2023) for an introduction to the topic.
-
+for the Proper Orthogonal Decomposition (POD). We refer to `Mendez (2022) <https://doi.org/10.48550/arXiv.2201.03847>`_, `Mendez (2023) <https://doi.org/10.1017/9781108896214.013>`_, and `Mendez et al. (2023) <https://arxiv.org/abs/2208.07746>`_,  for an introduction to the topic. 
 MODULO currently features the following decompositions: 
 
-- Discrete Fourier Transform (DFT) (briggs1995dft)
-- Proper Orthogonal Decomposition (POD) (sirovich1987turbulence, berkooz1993proper)
-- Multi-Scale Proper Orthogonal Decomposition (mPOD) (mendez2019multi)
-- Dynamic Mode Decomposition (DMD) (schmid2010dynamic)
-- Spectral Proper Orthogonal Decomposition (SPOD) (csieber2016spectral, towne2018spectral), 
+- Discrete Fourier Transform (DFT) (`Briggs et al. (1997) <https://epubs.siam.org/doi/book/10.1137/1.9781611971514>`_)
+- Proper Orthogonal Decomposition (POD) (`Sirovich et al. (1987) <https://www.ams.org/journals/qam/1987-45-03/S0033-569X-1987-0910464-1/S0033-569X-1987-0910464-1.pdf>`_ , `Berkooz et al. (1993) <https://doi.org/10.1146/annurev.fl.25.010193.002543>`_)
+- Multi-Scale Proper Orthogonal Decomposition (mPOD) (`Mendez et al. (2019) <https://arxiv.org/abs/1804.09646>`_)
+- Dynamic Mode Decomposition (DMD) (`Schmid (2010) <https://doi.org/10.1017/S0022112010001217>`_)
+- Spectral Proper Orthogonal Decomposition (SPOD) (`Sieber et al. (2016) <https://doi.org/10.48550/arXiv.1508.04642>`_, `Towne et al. (2018) <https://doi.org/10.48550/arXiv.1708.04393>`_), 
   note that the two are different formulations, and both are available in MODULO.
-- Kernel Proper Orthogonal Decomposition (KPOD) (mika1998kernel)
+- Kernel Proper Orthogonal Decomposition (KPOD) (`Mika et al. (1998) <https://proceedings.neurips.cc/paper_files/paper/1998/file/226d1f15ecd35f784d2a20c3ecf56d7f-Paper.pdf>`_)
 
 We remind the curious reader to the respective references for a detailed description of each decomposition, and to the
 documentation for a practical guide on how to use them in MODULO.
 
-
 Release Notes
 -------------
-The latest version of MODULO (v2.0) includes the following updates:
-
-1. **Faster EIG/SVD algorithms**, using powerful randomized svd solvers from scikit_learn 
-    (see `here <https://scikit-learn.org/stable/modules/generated/sklearn.decomposition.TruncatedSVD.html>`_ 
-    and `here <https://scikit-learn.org/stable/modules/generated/sklearn.utils.extmath.randomized_svd.html>`_.) 
-    It is now possible to select various options as "eig_solver" and "svd_solver", 
-    offering different trade-offs in terms of accuracy vs computational time.
-
-2. **Computation the POD directly via SVD**, using any of the four "svd_solver" options.
-This is generally faster but requires more memory.
-
-3. **Faster subscale estimators for the mPOD:** the previous version used the rank of the correlation matrix in each scale to define the number of modes to be computed in each portion of the splitting vector before assembling the full basis. This is computationally very demanding. This estimation has been replaced by a 
-frequency-based threshold (i.e. based on the frequency bins within each portion) since one can show that the 
-frequency-based estimator is always more "conservative" than the rank-based estimator.
+This version of MODULO includes the following updates:
 
-4. **Major improvement on the memory saving option** : the previous version of modulo always required in input the matrix D. 
-Then, if the memory saving option was active, the matrix was partitioned and stored locally to free the RAM before computing the 
-correlation matrix (see `this tutorial by D. Ninni <https://www.youtube.com/watch?v=LclxO1WTuao>`_). 
-In the new version, it is possible to initialize a modulo object *without* the matrix D (see exercise 5 in the examples). 
-Instead, one can create the partitions without loading the matrix D.
+1. **mPOD bug fix:** the previous version of mPOD was skipping the last scale of the frequency splitting vector. Fixed in this version.
 
-5. **Implementation of Dynamic Mode Decomposition (DMD)** from (Schmid, P.J 2010)
+2. **SPOD parallelisation:** CSD - SPOD can now be parallelized, leveraging `joblib`. The user needs just to pass the argument `n_processes` for the computation to be
+   split between different workers.
+   
+3. **Simplified decomposition interface:** the interface of the decomposition methods has been simplified to improve user experience. 
 
-6. **Implementation of the two Spectral POD formulations**, namely the one from (Sieber et al 2016), 
-   and the one from (Towne et al 2018).
+4. **Enhanced POD selection:** the POD function has been redesigned, allowing users to easily choose between different POD methods.  
+   
+5. **Improved computational efficiency:** the code of the decomposition functions has been optimised, resulting in reduced computation time. mPOD now includes two additional optional arguments to enable faster filtering and to avoid recomputing the Sigmas after QR polishing. 
+	
+6. **Extended documentation:** the documentation has been significantly enriched, now including theoretical foundations for all the supported modal decomposition techniques. 
 
-7. **Implementation of a kernel version of the POD**, in which the correlation matrix is replaced by a kernel matrix. This is described in Lecture 15 of the course `Hands on Machine Learning for Fluid dynamics 2023 <https://www.vki.ac.be/index.php/events-ls/events/eventdetail/552/-/online-on-site-hands-on-machine-learning-for-fluid-dynamics-2023>`_. We refer also to: `Mendez, 2022 <https://arxiv.org/abs/2208.07746>`_. 
-
-8. **Implementation of a formulation for non-uniform meshes**, using a weighted matrix for all the relevant inner products. This is currently available only for POD and mPOD but allows for handling data produced from CFD simulation without resampling on a uniform grid (see exercise 4). 
-It can be used both with and without the memory-saving option.
+	
+Documentation
+-------------
 
+The documentation of MODULO is available `here <https://lorenzoschena.github.io/MODULO/intro.html>`_. It 
+contains a comprehensive guide on how to install and use the package, as well as a detailed description of the
+decompositions required inputs and outputs. A `list of YouTube videos <https://www.youtube.com/@modulompod5682>`_ 
+is also available to guide the introduce the user to modal decomposition and MODULO.
 
 Installation
 -------------
@@ -162,14 +143,6 @@ or, if you have pip installed in your environment,
     $ pip install .
 
 
-Documentation
--------------
-
-The documentation of MODULO is available `here <https://modulo.readthedocs.io/en/latest/intro.html>`_. It 
-contains a comprehensive guide on how to install and use the package, as well as a detailed description of the
-decompositions required inputs and outputs. A `list of YouTube videos <https://www.youtube.com/@modulompod5682>`_ 
-is also available to guide the introduce the user to modal decomposition and MODULO.
-
 Example 
 -------------
 
@@ -242,60 +215,7 @@ you can use the weighted inner product formulation (refer to `examples/ex_05_non
 
     # Compute the POD decomposition
     phi_POD, Sigma_POD, psi_POD = m.POD()
-
-Computational Cost Estimates
-----------------------------
-We here provide a rough estimate of the amoung of RAM required to decompose a test case with and without the memory saving option. 
-This option reduces the memory usage at the cost of increasing the computational time (see https://www.youtube.com/watch?v=LclxO1WTuao)
-
-Given a dataset $D \\in \\mathbb{R}^{n_s \\times n_t}$, we consider the computation of $n_r$ modes. When using the memory saving option, we refer to 
-$n_t' = n_t / n_p$ as the number of time steps in each partition, and to $n_s' = n_s / n_p$ as the number of spatial points in each partition.
-
-.. list-table::
-   :header-rows: 1
-
-   * - 
-     - Phase 1: $D$
-     - Phase 2: $K$
-     - Phase 3: $\\Psi$
-     - Phase 4: $\\Phi$
-   * - No Memory Saving
-     - $n_s \\times n_t$
-     - $n_t^2$
-     - $n_t^2 + n_t \\times n_r$
-     - $n_s \\times n_t + n_t \\times n_r + n_s \\times n_r$
-   * - Memory Saving
-     - /
-     - $n_s \\times n_t' + n_t' \\times n_t'$
-     - $n_t^2 + n_t \\times n_r$
-     - $n_s \\times n_t' + n_s' \\times n_t + n_s \\times n_r$
-
-If the memory saving option is active, the memory requirements are mostly linked to the storage of the correlation matrix $K$ in Phase 2. 
-This table can be used to estimate if a dataset is too large for the available RAM, recalling that data in single precision requires 4 bytes (or 32 bits).
-
-For example, for a dataset with n_s=1 000 000 and n_t = 5000  the following table estimates the RAM required in the two cases, considering n_b=10 partitions in the case of memory saving:
-
-.. list-table::
-   :header-rows: 1
-
-   * - 
-     - Phase 1: $D$
-     - Phase 2: $K$
-     - Phase 3: $\\Psi$
-     - Phase 4: $\\Phi$
-   * - No Memory Saving
-     - 18.6 GB
-     - 0.093 GB
-     - ≈0.112 GB
-     - ≈22.39 GB
-   * - Memory Saving
-     - /
-     - ≈1.86 GB
-     - ≈ 0.0026 GB
-     - ≈ 4.20 GB
-
-
-
+	
 Community guidelines
 ---------------------
 
@@ -326,31 +246,36 @@ Ask for help
 If you have troubles using MODULO, or you need help with a specific decomposition, please open an issue on the MODULO GitHub repository.
 
 Citation
----------
-If you use MODULO in your research, please cite it as follows:
-
-``Poletti, R., Schena, L., Ninni, D. Mendez, M. A. (2024). MODULO: A Python toolbox for data-driven modal decomposition. Journal of Open Source Software, 9(102), 6753, https://doi.org/10.21105/joss.06753
-
+--------
 
-.. code-block:: text 
-
-  @article{Poletti2024, 
-    doi = {10.21105/joss.06753}, 
-    url = {https://doi.org/10.21105/joss.06753}, 
-    year = {2024},
-    publisher = {The Open Journal}, 
-    volume = {9}, 
-    number = {102}, 
-    pages = {6753},
-    author = {R. Poletti and L. Schena and D. Ninni and M. A. Mendez}, 
-    title = {MODULO: A Python toolbox for data-driven modal decomposition},
-    journal = {Journal of Open Source Software} }
+If you use MODULO in your research, please cite it as follows:
 
+Poletti, R., Schena, L., Ninni, D., Mendez, M. A. (2024).  
+*MODULO: A Python toolbox for data-driven modal decomposition*.  
+Journal of Open Source Software, 9(102), 6753. https://doi.org/10.21105/joss.06753
+
+.. code-block:: text
+
+    @article{Poletti2024, 
+      doi = {10.21105/joss.06753}, 
+      url = {https://doi.org/10.21105/joss.06753}, 
+      year = {2024},
+      publisher = {The Open Journal}, 
+      volume = {9}, 
+      number = {102}, 
+      pages = {6753},
+      author = {R. Poletti and L. Schena and D. Ninni and M. A. Mendez}, 
+      title = {MODULO: A Python toolbox for data-driven modal decomposition},
+      journal = {Journal of Open Source Software}
+    }
+	
 and 
 
-``Ninni, D., & Mendez, M. A. (2020). MODULO: A software for Multiscale Proper Orthogonal Decomposition of data. SoftwareX, 12, 100622.``
+Ninni, D., & Mendez, M. A. (2020).  
+*MODULO: A software for Multiscale Proper Orthogonal Decomposition of data*.  
+SoftwareX, 12, 100622. https://doi.org/10.1016/j.softx.2020.100622
 
-.. code-block:: text 
+.. code-block:: text
 
     @article{ninni2020modulo,
         title={MODULO: A software for Multiscale Proper Orthogonal Decomposition of data},
@@ -363,7 +288,6 @@ and
     }
 
 
-
 References
 ----------
 
@@ -395,3 +319,8 @@ There are also decomposition-specific projects, some of which are listed below:
 - Ichinaga, Andreuzzi, Demo, Tezzele, Lapo, Rozza, Brunton, Kutz. PyDMD: A Python package for robust dynamic mode decomposition. arXiv preprint, 2024.
 - Rogowski, Marcin, et al. "Unlocking massively parallel spectral proper orthogonal decompositions in the PySPOD package." Computer Physics Communications 302 (2024): 109246.
 
+
+
+
+
+

{modulo_vki-2.1.0.dist-info → modulo_vki-2.1.1.dist-info}/RECORD

@@ -19,8 +19,8 @@ modulo_vki/utils/_plots.py,sha256=m43t08cVq-TY0BW0YPqT71hN-54hBphIYKZEn8Kw16E,14
 modulo_vki/utils/_utils.py,sha256=WFD7nwjSzVHpevVwTEvMdjAmcbeqwoXT9M48tIIniJw,14355
 modulo_vki/utils/others.py,sha256=26ES5EmsLhwkvcXTwNhDMkblGrxoWepX5c9TXeLTRWg,17336
 modulo_vki/utils/read_db.py,sha256=lJFauxJxS0_mYoxrbn-43UqZjOkr-qb9f6RTUq4IxZU,15149
-modulo_vki-2.1.0.dist-info/licenses/LICENSE,sha256=5TivriXFErrYrJgBq3M72kHNHqtSiCft3xESM1zHc0k,1091
-modulo_vki-2.1.0.dist-info/METADATA,sha256=DwFf-uEygrIts5jdauD3tnjELaecgAKLsDr4uPO_haM,18201
-modulo_vki-2.1.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-modulo_vki-2.1.0.dist-info/top_level.txt,sha256=4PA4AmafKU6M7us7gvt_Q976Khx3qjNUEThRRM5zxeA,11
-modulo_vki-2.1.0.dist-info/RECORD,,
+modulo_vki-2.1.1.dist-info/licenses/LICENSE,sha256=5TivriXFErrYrJgBq3M72kHNHqtSiCft3xESM1zHc0k,1091
+modulo_vki-2.1.1.dist-info/METADATA,sha256=MGJRX_sOUGm4_18AnPXzz3lT8PKi7KxPhbHZ9yiQEVQ,14933
+modulo_vki-2.1.1.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+modulo_vki-2.1.1.dist-info/top_level.txt,sha256=4PA4AmafKU6M7us7gvt_Q976Khx3qjNUEThRRM5zxeA,11
+modulo_vki-2.1.1.dist-info/RECORD,,