I'm keeping this one as a nonincremental solve and throwing it away; I need to get the `boundaryConstraints` into the solver after this solve. I could make this solve incremental, but then I would have to complicate the logic in the double-while loop to add the `boundaryConstraints` only for the first time.