LU factorization indexing needs to be perfomed back-to-front in order to get it in the right order.  It is not intuitive, also the solve_LU_fast is breaking at size>2
added main() source file Test.cpp