Store CVXPY variables in lists/numpy arrays

Question

As a part of my optimization formulation, I need to implement the multinomial logit. A very simplified version of my problem is like this:

I want to store all