We implemented all methods in Python with SciPy [58] and JAX [9] and executed them on a computer with Apple M1 Chip (8 cores, 3.2 GHz) and 16 GB RAM. 7.1 Problem setting We consider three types of instances: (i) compressed sensing with quadratic measurement, (ii) nonnegative matrix ...