How can I solve this kind of large scale problem?

Try using CVXQUAD per my April post in CVX reports unbounded status but with viable solution . You will have to rewrite log in terms of explicit exponential cone cone construct, as mentioned there (also see mvg’s answer at Solve optimization problems of exp function for reformulation).