/** rPri := Ax - b */ @Override protected DoubleMatrix1D rPri(DoubleMatrix1D X) { if (getMeq() == 0) { return F1.make(0); } return ColtUtils.zMult(getA(), X, getB(), -1); }