Now 80% faster! Replaced naive solve with a schur complement one... did some presolving to combine inequality constraints and second order matrix to reduced KKT matrix size.. still no cuda version yet