petit/demo/lu.t

   1 real a(100,100),b(100)
   2 integer n
   3 !
   4 ! For LU decomposition, we can also
   5 ! obtain all 6 permutations. Normal
   6 ! loop interchange techniques can't
   7 ! produce the ijk ordering
   8 !
   9 !
  10 for k = 1 to n do
  11  for i = k+1 to n do
  12   a(i,k) = a(i,k) / a(k,k)
  13   for j = k+1 to n do
  14    a(i,j) = a(i,j) - a(i,k)*a(k,j)
  15   endfor
  16  endfor
  17 endfor