interval_union    812 src/Solve.cpp              result = interval_union(ia, ib);
interval_union    830 src/Solve.cpp              result = interval_union(ia, ib);