elif constraint.lower_bound == constraint.upper_bound:
                A_eq = numpy.concatenate((A_eq, A_row))
                b_eq = numpy.concatenate((b_eq, [float(constraint.lower_bound)]))