Multidimensional empirical mode decomposition: Difference between revisions

Revision as of 14:14, 12 February 2025

Latest revision as of 14:16, 12 February 2025

(3 intermediate revisions by the same user not shown)

Line 181:

The EEMDs comprising MEEMD are assigned to independent threads for parallel execution, relying on the OpenMP runtime to resolve any load imbalance issues. Stride memory accesses of high-dimensional data are eliminated by transposing these data to lower dimensions, resulting in better utilization of cache lines. The partial results of each EEMD are made thread-private for correct functionality. Memory requirements depend on the number of OpenMP threads and are managed by OpenMP runtime.<ref name=":8" />

=== CUDA implementation<ref name=":8" /> ===

In the GPU CUDA implementation, each EMD, is mapped to a thread. The memory layout, especially of high-dimensional data, is rearranged to meet memory coalescing requirements and fit into the 128-byte cache lines. The data is first loaded along the lowest dimension and then consumed along a higher dimension. This step is performed when the Gaussian noise is added to form the ensemble data. In the new memory layout, the ensemble dimension is added to the lowest dimension to reduce possible branch divergence. The impact of the unavoidable branch divergence from data irregularity, caused by the noise, is minimized via a regularization technique using the on-chip memory. Moreover, the cache memory is utilized to amortize unavoidable uncoalesced memory accesses.<ref name=":8" />

Line 191:

Based on the algorithm of BEMD, the implementation method of FABEMD is really similar to BEMD, but the FABEMD approach just changes the interpolation step into a direct envelope estimation method and restricts the number of iterations for every BIMF to one. As a result, two order statistics, including MAX and MIN, will be used for approximating the upper and lower envelope. The size of the filter will depend on the maxima and minima maps obtained from the input. The steps of the FABEMD algorithm are listed below.

===FABEMD algorithm===
Source:<ref name=":0" />===

;Step 1 – Determine and detect local maximum and minimum

Line 259 ⟶ 260:

Here, we take 2-D PDE-based EMD as an example.

===PDE-based BEMD algorithm===
Source:<ref name=":2" />===

;Step 1 – Extend super diffusion model from 1-D to 2-D

Line 318 ⟶ 320:

There are some problems in BEMD and boundary extending implementation in the iterative sifting process, including time-consuming, shape and continuity of the edges, decomposition results comparison and so on. In order to fix these problems, the '''Boundary Processing in Bidimensional Empirical Decomposition (BPBEMD)''' method was created. The main points of the new method algorithm will be described next.

===BPBEMD algorithm<ref name=":3" />===

The few core steps for BPBEMD algorithm are:<ref name=":3" />

;Step 1