@@ -65,12 +65,14 @@ Numerical Methods
6565+-------------------+------------------------------------------------------+-------------------------+---------------------------------------------------------------------------------------------------------+
6666| ``interp_method `` | Interpolation method for velocity fields | ``trilinear `` | Simple and fast triliear interpolation method. Very fast on GPUs. |
6767+-------------------+------------------------------------------------------+-------------------------+---------------------------------------------------------------------------------------------------------+
68- | | | (``tricubic ``) | |
68+ | | | (``tricubic ``) | Catmull–Rom 3D convolution kernel. |
6969+-------------------+------------------------------------------------------+-------------------------+---------------------------------------------------------------------------------------------------------+
7070| | | ``hermite `` | |
7171+-------------------+------------------------------------------------------+-------------------------+---------------------------------------------------------------------------------------------------------+
7272| | | ``WENO `` | Weighted Essentially Non-Oscillatory (WENO) scheme with shockwave capture capability. |
7373+-------------------+------------------------------------------------------+-------------------------+---------------------------------------------------------------------------------------------------------+
74+ | | | (``tricubic-FL ``) | ``I.P. `` The high-performance 3D tricubic interpolation method by [Lekien2005 ]_. |
75+ +-------------------+------------------------------------------------------+-------------------------+---------------------------------------------------------------------------------------------------------+
7476| ``grad_order `` | The method (order) of gradient discretization | (``2 ``) | 2nd-order central difference. |
7577+-------------------+------------------------------------------------------+-------------------------+---------------------------------------------------------------------------------------------------------+
7678| | | ``4 `` | ``I.P. `` 4th-order central difference. |
@@ -79,12 +81,14 @@ Numerical Methods
7981+-------------------+------------------------------------------------------+-------------------------+---------------------------------------------------------------------------------------------------------+
8082| | | ``100 `` | ``I.P. `` Global Fast Fourier Transform (FFT) with finite order. Compute extremely intensive! |
8183+-------------------+------------------------------------------------------+-------------------------+---------------------------------------------------------------------------------------------------------+
82- | ``Eigen_method `` | Solver for max eigenvalue of Cauchy–Green tensor | ``eigmax_sym3 `` | Closed-form cubic eigenvalue solver for symmetric 3×3 matrices (noniterative solver). |
84+ | ``eigen_method `` | Solver for max eigenvalue of Cauchy–Green tensor | ``eigmax_sym3 `` | Closed-form cubic eigenvalue solver for symmetric 3×3 matrices (noniterative solver). |
8385+-------------------+------------------------------------------------------+-------------------------+---------------------------------------------------------------------------------------------------------+
8486| | | ``jacobi `` | ``I.P. `` Jacobi eigenvalue algorithm (iterative). |
8587+-------------------+------------------------------------------------------+-------------------------+---------------------------------------------------------------------------------------------------------+
8688
87- As for your reference, and configured as defaults, the *Berkeley LCS Tutorials * used ``RK4 `` for advection, interpolated by ``tricubic ``. Although non detailed, ``grad_order=2 `` was employed by them from the equation, supposing the mesh is sufficiently refined.
89+ As for your reference, and configured as defaults, the *Berkeley LCS Tutorials * used ``RK4 `` for advection.
90+ The velocity fields were interpolated by ``tricubic-FL ``, originating from [Lekien2005 ]_, which has higher performance from
91+ Although non detailed, ``grad_order=2 `` was employed by them from the equation, supposing the mesh is sufficiently refined.
8892Please always notice that, although providing much better numerical precision and looks cool in papers, high-order methods could be resource-consuming, even hundreds of times.
8993
9094Dynamic LCS Window
0 commit comments