To shorten calculation time in finite element simulation by using multithreading computer, a parallel algorithm for stiffness matrix decomposition based on threadpool method is proposed. Firstly, a decomposition method of applicability to parallel computation is put forward by transferring the Cholesky's LLT method. Then, the threadplool is employed to generate multithreading for repeating use and the optimization is conducted considering load-balancing of each thread. Finally, numerical tests by using proposed algorithm in decomposition of one-dimensional array stored stiffness matrix are carried out on different calculation platforms with multi-processors. It is shown that the parallel algorithm can overcome the limitations of OpenMP when being applied in nested loops and is of high efficiency on stiffness matrix decomposition with low platform demands. The algorithm has explicit concept and minor programming difficulty and is applicable to solve problems caused by limitation of OpenMP in particular.
A novel network regularized matrix decomposition method to detect mutated cancer genes in tumour samples with inter-patient heterogeneity