Multivariate Empirical Mode Decomposition (MEMD) is a powerful tool to analyze nonlinear properties, however, it is very computationally time-consuming. The purpose of this study is to parallelize and accelerate MEMD, and to apply the accelerated MEMD to the analysis of electromagnetic wave propagation in photonic crystals, to improve the efficiency of generation of a demultiplexer. In the proposed method, part of the procedures that do not depend on channels of input signal were parallelized. Speedup which is based on 24 CPU threads of 5.12 times is achieved when 192 threads are adopted.