function [c, ceq] = columnConstraintsWithDisc(x)

x = columnMapVariables(x);


[c, ceq] = columnConstraints(x);


