[Help Wanted] Combining stage 2 and 3, trouble with optimization problem

[chusloj]

In my branch stage_2+3, I've been trying to combine stages 2 and 3 of the optimization process for PUF weights. I've expanded the variable for interest income (INTS) to 12 new variables, each of them pertaining to the 12 income bins specified in the SOI update instructions. . Running stage1_just_INT.py from my branch to produce Stage 2 targets based on these new variables and then running stage2.py, I get the following error from Julia's JuMP solver:

A similar error was thrown by the CVXOPT optimization package used previously for Stage 2. CVXOPT's error said the coefficient matrix produced in dataprep.py (I use dataprep_INT.py in my branch) does not have full rank.

I'm trying to identify why the optimization problem is infeasible, and I've had no luck so far.