Remove error bounds calculation from triangular solver and make errorbounds a separate function. Add some tests, adjust the test code and add some spaces.