lup_decomp.f90  n= 4
 A=
 2.00000000000000 3.00000000000000 5.00000000000000 6.00000000000000
 3.00000000000000 5.00000000000000 6.00000000000000 8.00000000000000
 1.00000000000000 4.00000000000000 5.00000000000000 7.00000000000000
 3.00000000000000 4.00000000000000 7.00000000000000 8.00000000000000
 initial U=
 1.00000000000000 0.000000000000000 0.000000000000000 0.000000000000000
 0.000000000000000 1.00000000000000 0.000000000000000 0.000000000000000
 0.000000000000000 0.000000000000000 1.00000000000000 0.000000000000000
 0.000000000000000 0.000000000000000 0.000000000000000 1.00000000000000
 initial L=
 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000
 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000
 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000
 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000
 interchange m= 1 , im= 2
 interchange m= 2 , im= 3
 interchange m= 3 , im= 4
 L( 1 , 1 (= 3.00000000000000
 L( 2 , 1 (= 1.00000000000000
 L( 3 , 1 (= 3.00000000000000
 L( 4 , 1 (= 2.00000000000000
 U( 1 , 2 (= 1.66666666666667
 U( 1 , 3 (= 2.00000000000000
 U( 1 , 4 (= 2.66666666666667
 m= 2 , i= 2 , k= 1 , L(i,m)= 2.33333333333333
 L( 2 , 2 (= 2.33333333333333
 m= 2 , i= 3 , k= 1 , L(i,m)= -1.00000000000000
 L( 3 , 2 (= -1.00000000000000
 m= 2 , i= 4 , k= 1 , L(i,m)= -0.333333333333333
 L( 4 , 2 (= -0.333333333333333
 m= 2 , j= 3 , k= 1 , U(m,j)= 3.00000000000000
 U( 2 , 3 (= 1.28571428571429
 m= 2 , j= 4 , k= 1 , U(m,j)= 4.33333333333333
 U( 2 , 4 (= 1.85714285714286
 m= 3 , i= 3 , k= 1 , L(i,m)= 1.00000000000000
 m= 3 , i= 3 , k= 2 , L(i,m)= 2.28571428571429
 L( 3 , 3 (= 2.28571428571429
 m= 3 , i= 4 , k= 1 , L(i,m)= 1.00000000000000
 m= 3 , i= 4 , k= 2 , L(i,m)= 1.42857142857143
 L( 4 , 3 (= 1.42857142857143
 m= 3 , j= 4 , k= 1 , U(m,j)= 4.44089209850063E-16
 m= 3 , j= 4 , k= 2 , U(m,j)= 1.85714285714286
 U( 3 , 4 (= 0.812500000000000
 m= 4 , i= 4 , k= 1 , L(i,m)= 0.666666666666667
 m= 4 , i= 4 , k= 2 , L(i,m)= 1.28571428571429
 m= 4 , i= 4 , k= 3 , L(i,m)= 0.125000000000000
 L( 4 , 4 (= 0.125000000000000
 computed U=
 1.00000000000000 1.66666666666667 2.00000000000000 2.66666666666667
 0.000000000000000 1.00000000000000 1.28571428571429 1.85714285714286
 0.000000000000000 0.000000000000000 1.00000000000000 0.812500000000000
 0.000000000000000 0.000000000000000 0.000000000000000 1.00000000000000
 computed L=
 3.00000000000000 0.000000000000000 0.000000000000000 0.000000000000000
 1.00000000000000 2.33333333333333 0.000000000000000 0.000000000000000
 3.00000000000000 -1.00000000000000 2.28571428571429 0.000000000000000
 2.00000000000000 -0.333333333333333 1.42857142857143 0.125000000000000
 list errors in L*U-A = 0
 error A( 2 , 4 )-L*U= -8.88178419700125E-16 A(i,j)= 7.00000000000000 t= 7.00000000000000
 b= 1.50000000000000 4.50000000000000 6.50000000000000 0.500000000000000
 bb= 0.500000000000000 1.71428571428571 2.93750000000000 -33.0000000000000
 x= -12.2500000000000 24.7500000000000 29.7500000000000 -33.0000000000000
 list errors in A*x-b = 0
 error Ax-b( 2 )= -4.97379915032070E-14 b(i)= 4.50000000000000 t= 4.50000000000005
 error Ax-b( 4 )= -1.42108547152020E-14 b(i)= 0.500000000000000 t= 0.500000000000014