Alternating maximization: unifying framework for 8 sparse PCA formulations and efficient parallel codes

Peter Richtarik, Majid Jahani, Selin Damla Ahipaşaoğlu, Martin Takáč

Research output: Contribution to journalArticlepeer-review

3 Scopus citations


Given a multivariate data set, sparse principal component analysis (SPCA) aims to extract several linear combinations of the variables that together explain the variance in the data as much as possible, while controlling the number of nonzero loadings in these combinations. In this paper we consider 8 different optimization formulations for computing a single sparse loading vector: we employ two norms for measuring variance (L2, L1) and two sparsity-inducing norms (L0, L1), which are used in two ways (constraint, penalty). Three of our formulations, notably the one with L0 constraint and L1 variance, have not been considered in the literature. We give a unifying reformulation which we propose to solve via the alternating maximization (AM) method. We show that AM is equivalent to GPower for all formulations. Besides this, we provide 24 efficient parallel SPCA implementations: 3 codes (multi-core, GPU and cluster) for each of the 8 problems. Parallelism in the methods is aimed at (1) speeding up computations (our GPU code can be 100 times faster than an efficient serial code written in C++), (2) obtaining solutions explaining more variance and (3) dealing with big data problems (our cluster code can solve a 357 GB problem in a minute).
Original languageEnglish (US)
JournalOptimization and Engineering
StatePublished - Sep 22 2020

Bibliographical note

KAUST Repository Item: Exported on 2020-10-01


Dive into the research topics of 'Alternating maximization: unifying framework for 8 sparse PCA formulations and efficient parallel codes'. Together they form a unique fingerprint.

Cite this