LU Decomposition Part 2

This method uses a single array to represent an nxn matrix instead of a matrix. This should make the code more robust since C naturally reads arrays across columns. For example:

Code 1: 

for j 1 to 10

   for i 1 to 10

      a[i][j] += 2;

   end

end

Code 2: 

for j 1 to 10

   for i 1 to 10

      a[j][i] += 2;

   end

end

Code 1 wont be as fast as code 2 but in fortran the opposite is true. 

Here's a link to my code.