LU solve routines changed to use permutation vector and solve using L,U overwritten in A (no longer form L,U)