bool proj_sub(const IntervalMatrix& y, IntervalMatrix& x1, IntervalMatrix& x2) { x1 &= y+x2; x2 &= x1-y; return !x1.is_empty() && !x2.is_empty(); }
bool proj_add(const IntervalMatrix& y, IntervalMatrix& x1, IntervalMatrix& x2) { x1 &= y-x2; x2 &= y-x1; return !x1.is_empty() && !x2.is_empty(); }