CP decomposition of semi-nonnegative semisymmetric tensors based on QR matrix factorization

International audienceThe problem of Canonical Polyadic (CP) decomposition of semi-nonnegative semi-symmetric three-way arrays is often encountered in Independent Component Analysis (ICA), where the cumulant of a nonnegative mixing process is frequently involved, such as the Magnetic Resonance Spectroscopy (MRS). We propose a new method, called JD+QR, to solve such a problem. The nonnegativity constraint is imposed by means of a square change of variable. Then the high-dimensional optimization problem is decomposed into several sequential rational subproblems using QR matrix factorization. A numerical experiment on simulated arrays emphasizes its good performance. A BSS application on MRS data confirms the validity and improvement of the proposed method