Non-Negative Tensor Factorization Accelerated Using GPGPU

This article presents an optimized algorithm for Non-Negative Tensor Factorization (NTF), implemented in the CUDA (Compute Uniform Device Architecture) framework, that runs on contemporary graphics processors and exploits their massive parallelism. The NTF implementation is primarily targeted for analysis of high-dimensional spectral images, including dimensionality reduction, feature extraction, and other tasks related to spectral imaging; however, the algorithm and its implementation are not limited to spectral imaging. The speed-ups measured on real spectral images are around 60-100x compared to a traditional  C implementation compiled with an optimizing compiler.  Since common problems in the field of spectral imaging may take hours on a state-of-the-art CPU, the speed-up achieved using a graphics card is attractive.  The implementation is publicly available in the form of a dynamically linked library, including an interface to MATLAB, and thus may be of help to researchers and engineers using NTF on large problems.