We study the low rank phase retrieval (LRPR) problem defined as follows:
recover an $n \times q$ matrix $X^*$ of rank $r$ from a different and
independent set of $m$ phaseless (magnitude-only) linear projections of each of
its columns. To be precise, we need to recover $X^*$ from $y_k