% We don't support initials here yet
% x0 = []; % ZL
% To simplify code, we currently normalize everything to z'*z<z0^2
% This done by writting [c^Tx+d;Ax+b] in cone as
% [c^Tx+d;Ax+b]==[z0;z], [z0;z] in cone
nNew = sum(K.q);
x0 = [x0; NaN(nNew,1)]; % ZL