make crout work in place, also speed up U a little bit
move denominator computation out of main loop
1.0.1