It uses U0 but never explains what it is. It also calls a the "shift" but
never explains what it is. Perhaps some notes explaining this stuff would make
the solver more transparent. And a reference to further material that uses this
approach would be great if it exists.Barry
