Remove error bounds calculation from triangular solver

Thanks to @ivanslapnicar for pointing out that out triangular solver was very slow. This pull request makes a 1000x1000 problem with 1000 right hand sides about 70 times faster. The reason was that an error bounds calculation was done inside the solver. It could right as well be done outside the solver.