2 November 2013 11:37:00.370 AM LAPLACIAN_PRB FORTRAN77 version Test the LAPLACIAN library. TEST01 A full-storage matrix is returned by: L1DD: Dirichlet/Dirichlet BC; L1DN: Dirichlet/Neumann BC; L1ND: Neumann/Dirichlet BC; L1NN: Neumann/Neumann BC; L1PP: Periodic BC; Using spacing H = 1.00000 L1DD: Col 1 2 3 4 5 Row 1: 2.00000 -1.00000 0.00000 0.00000 0.00000 2: -1.00000 2.00000 -1.00000 0.00000 0.00000 3: 0.00000 -1.00000 2.00000 -1.00000 0.00000 4: 0.00000 0.00000 -1.00000 2.00000 -1.00000 5: 0.00000 0.00000 0.00000 -1.00000 2.00000 L1DN: Col 1 2 3 4 5 Row 1: 2.00000 -1.00000 0.00000 0.00000 0.00000 2: -1.00000 2.00000 -1.00000 0.00000 0.00000 3: 0.00000 -1.00000 2.00000 -1.00000 0.00000 4: 0.00000 0.00000 -1.00000 2.00000 -1.00000 5: 0.00000 0.00000 0.00000 -1.00000 1.00000 L1ND: Col 1 2 3 4 5 Row 1: 1.00000 -1.00000 0.00000 0.00000 0.00000 2: -1.00000 2.00000 -1.00000 0.00000 0.00000 3: 0.00000 -1.00000 2.00000 -1.00000 0.00000 4: 0.00000 0.00000 -1.00000 2.00000 -1.00000 5: 0.00000 0.00000 0.00000 -1.00000 2.00000 L1NN: Col 1 2 3 4 5 Row 1: 1.00000 -1.00000 0.00000 0.00000 0.00000 2: -1.00000 2.00000 -1.00000 0.00000 0.00000 3: 0.00000 -1.00000 2.00000 -1.00000 0.00000 4: 0.00000 0.00000 -1.00000 2.00000 -1.00000 5: 0.00000 0.00000 0.00000 -1.00000 1.00000 L1PP: Col 1 2 3 4 5 Row 1: 2.00000 -1.00000 0.00000 0.00000 -1.00000 2: -1.00000 2.00000 -1.00000 0.00000 0.00000 3: 0.00000 -1.00000 2.00000 -1.00000 0.00000 4: 0.00000 0.00000 -1.00000 2.00000 -1.00000 5: -1.00000 0.00000 0.00000 -1.00000 2.00000 Using spacing H = 0.166667 L1DD: Col 1 2 3 4 5 Row 1: 72.0000 -36.0000 0.00000 0.00000 0.00000 2: -36.0000 72.0000 -36.0000 0.00000 0.00000 3: 0.00000 -36.0000 72.0000 -36.0000 0.00000 4: 0.00000 0.00000 -36.0000 72.0000 -36.0000 5: 0.00000 0.00000 0.00000 -36.0000 72.0000 L1DN: Col 1 2 3 4 5 Row 1: 72.0000 -36.0000 0.00000 0.00000 0.00000 2: -36.0000 72.0000 -36.0000 0.00000 0.00000 3: 0.00000 -36.0000 72.0000 -36.0000 0.00000 4: 0.00000 0.00000 -36.0000 72.0000 -36.0000 5: 0.00000 0.00000 0.00000 -36.0000 36.0000 L1ND: Col 1 2 3 4 5 Row 1: 36.0000 -36.0000 0.00000 0.00000 0.00000 2: -36.0000 72.0000 -36.0000 0.00000 0.00000 3: 0.00000 -36.0000 72.0000 -36.0000 0.00000 4: 0.00000 0.00000 -36.0000 72.0000 -36.0000 5: 0.00000 0.00000 0.00000 -36.0000 72.0000 L1NN: Col 1 2 3 4 5 Row 1: 36.0000 -36.0000 0.00000 0.00000 0.00000 2: -36.0000 72.0000 -36.0000 0.00000 0.00000 3: 0.00000 -36.0000 72.0000 -36.0000 0.00000 4: 0.00000 0.00000 -36.0000 72.0000 -36.0000 5: 0.00000 0.00000 0.00000 -36.0000 36.0000 L1PP: Col 1 2 3 4 5 Row 1: 72.0000 -36.0000 0.00000 0.00000 -36.0000 2: -36.0000 72.0000 -36.0000 0.00000 0.00000 3: 0.00000 -36.0000 72.0000 -36.0000 0.00000 4: 0.00000 0.00000 -36.0000 72.0000 -36.0000 5: -36.0000 0.00000 0.00000 -36.0000 72.0000 TEST02 The Laplacian L is applied to data U by: L1DD_APPLY for Dirichlet/Dirichlet BC; L1DN_APPLY for Dirichlet/Neumann BC; L1ND_APPLY for Neumann/Dirichlet BC; L1NN_APPLY for Neumann/Neumann BC; L1PP_APPLY for Periodic BC; Using spacing H = 0.100000 Vector U: 1: 0.90000000E-01 2: 0.16000000 3: 0.21000000 4: 0.24000000 5: 0.25000000 6: 0.24000000 7: 0.21000000 8: 0.16000000 9: 0.90000000E-01 L1DD(U): 1: 2.0000000 2: 2.0000000 3: 2.0000000 4: 2.0000000 5: 2.0000000 6: 2.0000000 7: 2.0000000 8: 2.0000000 9: 2.0000000 L1DN(U): 1: 2.0000000 2: 2.0000000 3: 2.0000000 4: 2.0000000 5: 2.0000000 6: 2.0000000 7: 2.0000000 8: 2.0000000 9: -7.0000000 L1ND(U): 1: -7.0000000 2: 2.0000000 3: 2.0000000 4: 2.0000000 5: 2.0000000 6: 2.0000000 7: 2.0000000 8: 2.0000000 9: 2.0000000 L1NN(U): 1: -7.0000000 2: 2.0000000 3: 2.0000000 4: 2.0000000 5: 2.0000000 6: 2.0000000 7: 2.0000000 8: 2.0000000 9: -7.0000000 L1PP(U): 1: -7.0000000 2: 2.0000000 3: 2.0000000 4: 2.0000000 5: 2.0000000 6: 2.0000000 7: 2.0000000 8: 2.0000000 9: -7.0000000 TEST03 Compute eigen information for the Laplacian: L1DD_EIGEN for Dirichlet/Dirichlet BC; L1DN_EIGEN for Dirichlet/Neumann BC; L1ND_EIGEN for Neumann/Dirichlet BC; L1NN_EIGEN for Neumann/Neumann BC; L1PP_EIGEN for Periodic BC; Using spacing H = 1.00000 L1DD Eigenvalues: 1: 0.26794919 2: 1.0000000 3: 2.0000000 4: 3.0000000 5: 3.7320508 L1DD Eigenvectors: Col 1 2 3 4 5 Row 1: 0.288675 0.500000 0.577350 0.500000 0.288675 2: 0.500000 0.500000 0.707050E-16 -0.500000 -0.500000 3: 0.577350 0.707050E-16 -0.577350 -0.141410E-15 0.577350 4: 0.500000 -0.500000 -0.141410E-15 0.500000 -0.500000 5: 0.288675 -0.500000 0.577350 -0.500000 0.288675 L1DD eigenerror = 0.107099E-14 L1DN Eigenvalues: 1: 0.81014053E-01 2: 0.69027853 3: 1.7153703 4: 2.8308300 5: 3.6825071 L1DN Eigenvectors: Col 1 2 3 4 5 Row 1: 0.169891 0.455734 0.596885 0.548529 0.326019 2: 0.326019 0.596885 0.169891 -0.455734 -0.548529 3: 0.455734 0.326019 -0.548529 -0.169891 0.596885 4: 0.548529 -0.169891 -0.326019 0.596885 -0.455734 5: 0.596885 -0.548529 0.455734 -0.326019 0.169891 L1DN eigenerror = 0.147635E-14 L1ND Eigenvalues: 1: 0.81014053E-01 2: 0.69027853 3: 1.7153703 4: 2.8308300 5: 3.6825071 L1ND Eigenvectors: Col 1 2 3 4 5 Row 1: 0.596885 0.548529 0.455734 0.326019 0.169891 2: 0.548529 0.169891 -0.326019 -0.596885 -0.455734 3: 0.455734 -0.326019 -0.548529 0.169891 0.596885 4: 0.326019 -0.596885 0.169891 0.455734 -0.548529 5: 0.169891 -0.455734 0.596885 -0.548529 0.326019 L1ND eigenerror = 0.143710E-14 L1NN Eigenvalues: 1: 0.0000000 2: 0.38196601 3: 1.3819660 4: 2.6180340 5: 3.6180340 L1NN Eigenvectors: Col 1 2 3 4 5 Row 1: 2.23607 0.601501 0.511667 0.371748 0.195440 2: 2.23607 0.371748 -0.195440 -0.601501 -0.511667 3: 2.23607 0.387267E-16 -0.632456 -0.116180E-15 0.632456 4: 2.23607 -0.371748 -0.195440 0.601501 -0.511667 5: 2.23607 -0.601501 0.511667 -0.371748 0.195440 L1NN eigenerror = 0.860588E-15 L1PP Eigenvalues: 1: 0.0000000 2: 1.3819660 3: 1.3819660 4: 3.6180340 5: 3.6180340 L1PP Eigenvectors: Col 1 2 3 4 5 Row 1: 0.447214 0.371748 0.511667 0.601501 0.195440 2: 0.447214 0.601501 -0.195440 -0.371748 -0.511667 3: 0.447214 0.774535E-16 -0.632456 -0.154907E-15 0.632456 4: 0.447214 -0.601501 -0.195440 0.371748 -0.511667 5: 0.447214 -0.371748 0.511667 -0.601501 0.195440 L1PP eigenerror = 0.806193E-15 Using spacing H = 0.166667 L1DD Eigenvalues: 1: 9.6461709 2: 36.000000 3: 72.000000 4: 108.00000 5: 134.35383 L1DD Eigenvectors: Col 1 2 3 4 5 Row 1: 0.288675 0.500000 0.577350 0.500000 0.288675 2: 0.500000 0.500000 0.707050E-16 -0.500000 -0.500000 3: 0.577350 0.707050E-16 -0.577350 -0.141410E-15 0.577350 4: 0.500000 -0.500000 -0.141410E-15 0.500000 -0.500000 5: 0.288675 -0.500000 0.577350 -0.500000 0.288675 L1DD eigenerror = 0.498107E-13 L1DN Eigenvalues: 1: 2.9165059 2: 24.850027 3: 61.753332 4: 101.90988 5: 132.57025 L1DN Eigenvectors: Col 1 2 3 4 5 Row 1: 0.169891 0.455734 0.596885 0.548529 0.326019 2: 0.326019 0.596885 0.169891 -0.455734 -0.548529 3: 0.455734 0.326019 -0.548529 -0.169891 0.596885 4: 0.548529 -0.169891 -0.326019 0.596885 -0.455734 5: 0.596885 -0.548529 0.455734 -0.326019 0.169891 L1DN eigenerror = 0.445168E-13 L1ND Eigenvalues: 1: 2.9165059 2: 24.850027 3: 61.753332 4: 101.90988 5: 132.57025 L1ND Eigenvectors: Col 1 2 3 4 5 Row 1: 0.596885 0.548529 0.455734 0.326019 0.169891 2: 0.548529 0.169891 -0.326019 -0.596885 -0.455734 3: 0.455734 -0.326019 -0.548529 0.169891 0.596885 4: 0.326019 -0.596885 0.169891 0.455734 -0.548529 5: 0.169891 -0.455734 0.596885 -0.548529 0.326019 L1ND eigenerror = 0.601291E-13 L1NN Eigenvalues: 1: 0.0000000 2: 13.750776 3: 49.750776 4: 94.249224 5: 130.24922 L1NN Eigenvectors: Col 1 2 3 4 5 Row 1: 2.23607 0.601501 0.511667 0.371748 0.195440 2: 2.23607 0.371748 -0.195440 -0.601501 -0.511667 3: 2.23607 0.387267E-16 -0.632456 -0.116180E-15 0.632456 4: 2.23607 -0.371748 -0.195440 0.601501 -0.511667 5: 2.23607 -0.601501 0.511667 -0.371748 0.195440 L1NN eigenerror = 0.317450E-13 L1PP Eigenvalues: 1: 0.0000000 2: 49.750776 3: 49.750776 4: 130.24922 5: 130.24922 L1PP Eigenvectors: Col 1 2 3 4 5 Row 1: 0.447214 0.371748 0.511667 0.601501 0.195440 2: 0.447214 0.601501 -0.195440 -0.371748 -0.511667 3: 0.447214 0.774535E-16 -0.632456 -0.154907E-15 0.632456 4: 0.447214 -0.601501 -0.195440 0.371748 -0.511667 5: 0.447214 -0.371748 0.511667 -0.601501 0.195440 L1PP eigenerror = 0.417302E-13 TEST04 The inverse of a full-storage matrix is returned by: L1DD_INVERSE: Dirichlet/Dirichlet BC; L1DN_INVERSE: Dirichlet/Neumann BC; L1ND_INVERSE: Neumann/Dirichlet BC; Using spacing H = 1.00000 L1DD: Col 1 2 3 4 5 Row 1: 2.00000 -1.00000 0.00000 0.00000 0.00000 2: -1.00000 2.00000 -1.00000 0.00000 0.00000 3: 0.00000 -1.00000 2.00000 -1.00000 0.00000 4: 0.00000 0.00000 -1.00000 2.00000 -1.00000 5: 0.00000 0.00000 0.00000 -1.00000 2.00000 L1DD_INVERSE: Col 1 2 3 4 5 Row 1: 0.833333 0.666667 0.500000 0.333333 0.166667 2: 0.666667 1.33333 1.00000 0.666667 0.333333 3: 0.500000 1.00000 1.50000 1.00000 0.500000 4: 0.333333 0.666667 1.00000 1.33333 0.666667 5: 0.166667 0.333333 0.500000 0.666667 0.833333 L1DD inverse error = 0.686635E-15 L1DN: Col 1 2 3 4 5 Row 1: 2.00000 -1.00000 0.00000 0.00000 0.00000 2: -1.00000 2.00000 -1.00000 0.00000 0.00000 3: 0.00000 -1.00000 2.00000 -1.00000 0.00000 4: 0.00000 0.00000 -1.00000 2.00000 -1.00000 5: 0.00000 0.00000 0.00000 -1.00000 1.00000 L1DN_INVERSE: Col 1 2 3 4 5 Row 1: 1.00000 1.00000 1.00000 1.00000 1.00000 2: 1.00000 2.00000 2.00000 2.00000 2.00000 3: 1.00000 2.00000 3.00000 3.00000 3.00000 4: 1.00000 2.00000 3.00000 4.00000 4.00000 5: 1.00000 2.00000 3.00000 4.00000 5.00000 L1DN inverse error = 0.00000 L1ND: Col 1 2 3 4 5 Row 1: 1.00000 -1.00000 0.00000 0.00000 0.00000 2: -1.00000 2.00000 -1.00000 0.00000 0.00000 3: 0.00000 -1.00000 2.00000 -1.00000 0.00000 4: 0.00000 0.00000 -1.00000 2.00000 -1.00000 5: 0.00000 0.00000 0.00000 -1.00000 2.00000 L1ND_INVERSE: Col 1 2 3 4 5 Row 1: 5.00000 4.00000 3.00000 2.00000 1.00000 2: 4.00000 4.00000 3.00000 2.00000 1.00000 3: 3.00000 3.00000 3.00000 2.00000 1.00000 4: 2.00000 2.00000 2.00000 2.00000 1.00000 5: 1.00000 1.00000 1.00000 1.00000 1.00000 L1ND inverse error = 0.00000 Using spacing H = 0.166667 L1DD: Col 1 2 3 4 5 Row 1: 72.0000 -36.0000 0.00000 0.00000 0.00000 2: -36.0000 72.0000 -36.0000 0.00000 0.00000 3: 0.00000 -36.0000 72.0000 -36.0000 0.00000 4: 0.00000 0.00000 -36.0000 72.0000 -36.0000 5: 0.00000 0.00000 0.00000 -36.0000 72.0000 L1DD_INVERSE: Col 1 2 3 4 5 Row 1: 0.231481E-01 0.185185E-01 0.138889E-01 0.925926E-02 0.462963E-02 2: 0.185185E-01 0.370370E-01 0.277778E-01 0.185185E-01 0.925926E-02 3: 0.138889E-01 0.277778E-01 0.416667E-01 0.277778E-01 0.138889E-01 4: 0.925926E-02 0.185185E-01 0.277778E-01 0.370370E-01 0.185185E-01 5: 0.462963E-02 0.925926E-02 0.138889E-01 0.185185E-01 0.231481E-01 L1DD inverse error = 0.116573E-14 L1DN: Col 1 2 3 4 5 Row 1: 72.0000 -36.0000 0.00000 0.00000 0.00000 2: -36.0000 72.0000 -36.0000 0.00000 0.00000 3: 0.00000 -36.0000 72.0000 -36.0000 0.00000 4: 0.00000 0.00000 -36.0000 72.0000 -36.0000 5: 0.00000 0.00000 0.00000 -36.0000 36.0000 L1DN_INVERSE: Col 1 2 3 4 5 Row 1: 0.277778E-01 0.277778E-01 0.277778E-01 0.277778E-01 0.277778E-01 2: 0.277778E-01 0.555556E-01 0.555556E-01 0.555556E-01 0.555556E-01 3: 0.277778E-01 0.555556E-01 0.833333E-01 0.833333E-01 0.833333E-01 4: 0.277778E-01 0.555556E-01 0.833333E-01 0.111111 0.111111 5: 0.277778E-01 0.555556E-01 0.833333E-01 0.111111 0.138889 L1DN inverse error = 0.251215E-14 L1ND: Col 1 2 3 4 5 Row 1: 36.0000 -36.0000 0.00000 0.00000 0.00000 2: -36.0000 72.0000 -36.0000 0.00000 0.00000 3: 0.00000 -36.0000 72.0000 -36.0000 0.00000 4: 0.00000 0.00000 -36.0000 72.0000 -36.0000 5: 0.00000 0.00000 0.00000 -36.0000 72.0000 L1ND_INVERSE: Col 1 2 3 4 5 Row 1: 0.138889 0.111111 0.833333E-01 0.555556E-01 0.277778E-01 2: 0.111111 0.111111 0.833333E-01 0.555556E-01 0.277778E-01 3: 0.833333E-01 0.833333E-01 0.833333E-01 0.555556E-01 0.277778E-01 4: 0.555556E-01 0.555556E-01 0.555556E-01 0.555556E-01 0.277778E-01 5: 0.277778E-01 0.277778E-01 0.277778E-01 0.277778E-01 0.277778E-01 L1ND inverse error = 0.251215E-14 TEST05 Compute upper Cholesky factors for the Laplacian: L1DD_CHOLESKY for Dirichlet/Dirichlet BC; L1DN_CHOLESKY for Dirichlet/Neumann BC; L1ND_CHOLESKY for Neumann/Dirichlet BC; L1NN_CHOLESKY for Neumann/Neumann BC; L1PP_CHOLESKY for Periodic BC; Using spacing H = 1.00000 L1DD Cholesky factor: Col 1 2 3 4 5 Row 1: 1.41421 -0.707107 0.00000 0.00000 0.00000 2: 0.00000 1.22474 -0.816497 0.00000 0.00000 3: 0.00000 0.00000 1.15470 -0.866025 0.00000 4: 0.00000 0.00000 0.00000 1.11803 -0.894427 5: 0.00000 0.00000 0.00000 0.00000 1.09545 L1DD Cholesky error = 0.888178E-15 L1DN Cholesky factor: Col 1 2 3 4 5 Row 1: 1.41421 -0.707107 0.00000 0.00000 0.00000 2: 0.00000 1.22474 -0.816497 0.00000 0.00000 3: 0.00000 0.00000 1.15470 -0.866025 0.00000 4: 0.00000 0.00000 0.00000 1.11803 -0.894427 5: 0.00000 0.00000 0.00000 0.00000 0.447214 L1DN Cholesky error = 0.777156E-15 L1ND Cholesky factor: Col 1 2 3 4 5 Row 1: 1.00000 -1.00000 0.00000 0.00000 0.00000 2: 0.00000 1.00000 -1.00000 0.00000 0.00000 3: 0.00000 0.00000 1.00000 -1.00000 0.00000 4: 0.00000 0.00000 0.00000 1.00000 -1.00000 5: 0.00000 0.00000 0.00000 0.00000 1.00000 L1ND Cholesky error = 0.00000 L1NN Cholesky factor: Col 1 2 3 4 5 Row 1: 1.00000 -1.00000 0.00000 0.00000 0.00000 2: 0.00000 1.00000 -1.00000 0.00000 0.00000 3: 0.00000 0.00000 1.00000 -1.00000 0.00000 4: 0.00000 0.00000 0.00000 1.00000 -1.00000 5: 0.00000 0.00000 0.00000 0.00000 0.00000 L1NN Cholesky error = 0.00000 L1PP Cholesky factor: Col 1 2 3 4 5 Row 1: 1.41421 -0.707107 0.00000 0.00000 -0.707107 2: 0.00000 1.22474 -0.816497 0.00000 -0.408248 3: 0.00000 0.00000 1.15470 -0.866025 -0.288675 4: 0.00000 0.00000 0.00000 1.11803 -1.11803 5: 0.00000 0.00000 0.00000 0.00000 0.00000 L1PP Cholesky error = 0.111300E-14 Using spacing H = 0.166667 L1DD Cholesky factor: Col 1 2 3 4 5 Row 1: 8.48528 -4.24264 0.00000 0.00000 0.00000 2: 0.00000 7.34847 -4.89898 0.00000 0.00000 3: 0.00000 0.00000 6.92820 -5.19615 0.00000 4: 0.00000 0.00000 0.00000 6.70820 -5.36656 5: 0.00000 0.00000 0.00000 0.00000 6.57267 L1DD Cholesky error = 0.375984E-13 L1DN Cholesky factor: Col 1 2 3 4 5 Row 1: 8.48528 -4.24264 0.00000 0.00000 0.00000 2: 0.00000 7.34847 -4.89898 0.00000 0.00000 3: 0.00000 0.00000 6.92820 -5.19615 0.00000 4: 0.00000 0.00000 0.00000 6.70820 -5.36656 5: 0.00000 0.00000 0.00000 0.00000 2.68328 L1DN Cholesky error = 0.375984E-13 L1ND Cholesky factor: Col 1 2 3 4 5 Row 1: 6.00000 -6.00000 0.00000 0.00000 0.00000 2: 0.00000 6.00000 -6.00000 0.00000 0.00000 3: 0.00000 0.00000 6.00000 -6.00000 0.00000 4: 0.00000 0.00000 0.00000 6.00000 -6.00000 5: 0.00000 0.00000 0.00000 0.00000 6.00000 L1ND Cholesky error = 0.00000 L1NN Cholesky factor: Col 1 2 3 4 5 Row 1: 6.00000 -6.00000 0.00000 0.00000 0.00000 2: 0.00000 6.00000 -6.00000 0.00000 0.00000 3: 0.00000 0.00000 6.00000 -6.00000 0.00000 4: 0.00000 0.00000 0.00000 6.00000 -6.00000 5: 0.00000 0.00000 0.00000 0.00000 0.00000 L1NN Cholesky error = 0.00000 L1PP Cholesky factor: Col 1 2 3 4 5 Row 1: 8.48528 -4.24264 0.00000 0.00000 -4.24264 2: 0.00000 7.34847 -4.89898 0.00000 -2.44949 3: 0.00000 0.00000 6.92820 -5.19615 -1.73205 4: 0.00000 0.00000 0.00000 6.70820 -6.70820 5: 0.00000 0.00000 0.00000 0.00000 0.00000 L1PP Cholesky error = 0.414314E-13 TEST06 Compute LU factors for the Laplacian: L1DD_LU for Dirichlet/Dirichlet BC; L1DN_LU for Dirichlet/Neumann BC; L1ND_LU for Neumann/Dirichlet BC; L1NN_LU for Neumann/Neumann BC; L1PP_LU for Periodic BC; Using spacing H = 1.00000 L1DD L factor: Col 1 2 3 4 5 Row 1: 1.00000 0.00000 0.00000 0.00000 0.00000 2: -0.500000 1.00000 0.00000 0.00000 0.00000 3: 0.00000 -0.666667 1.00000 0.00000 0.00000 4: 0.00000 0.00000 -0.750000 1.00000 0.00000 5: 0.00000 0.00000 0.00000 -0.800000 1.00000 L1DD U factor: Col 1 2 3 4 5 Row 1: 2.00000 -1.00000 0.00000 0.00000 0.00000 2: 0.00000 1.50000 -1.00000 0.00000 0.00000 3: 0.00000 0.00000 1.33333 -1.00000 0.00000 4: 0.00000 0.00000 0.00000 1.25000 -1.00000 5: 0.00000 0.00000 0.00000 0.00000 1.20000 L1DD LU error = 0.00000 L1DN L factor: Col 1 2 3 4 5 Row 1: 1.00000 0.00000 0.00000 0.00000 0.00000 2: -0.500000 1.00000 0.00000 0.00000 0.00000 3: 0.00000 -0.666667 1.00000 0.00000 0.00000 4: 0.00000 0.00000 -0.750000 1.00000 0.00000 5: 0.00000 0.00000 0.00000 -0.800000 1.00000 L1DN U factor: Col 1 2 3 4 5 Row 1: 2.00000 -1.00000 0.00000 0.00000 0.00000 2: 0.00000 1.50000 -1.00000 0.00000 0.00000 3: 0.00000 0.00000 1.33333 -1.00000 0.00000 4: 0.00000 0.00000 0.00000 1.25000 -1.00000 5: 0.00000 0.00000 0.00000 0.00000 0.200000 L1DN LU error = 0.00000 L1ND L factor: Col 1 2 3 4 5 Row 1: 1.00000 0.00000 0.00000 0.00000 0.00000 2: -1.00000 1.00000 0.00000 0.00000 0.00000 3: 0.00000 -1.00000 1.00000 0.00000 0.00000 4: 0.00000 0.00000 -1.00000 1.00000 0.00000 5: 0.00000 0.00000 0.00000 -1.00000 1.00000 L1ND U factor: Col 1 2 3 4 5 Row 1: 1.00000 -1.00000 0.00000 0.00000 0.00000 2: 0.00000 1.00000 -1.00000 0.00000 0.00000 3: 0.00000 0.00000 1.00000 -1.00000 0.00000 4: 0.00000 0.00000 0.00000 1.00000 -1.00000 5: 0.00000 0.00000 0.00000 0.00000 1.00000 L1ND LU error = 0.00000 L1NN L factor: Col 1 2 3 4 5 Row 1: 1.00000 0.00000 0.00000 0.00000 0.00000 2: -1.00000 1.00000 0.00000 0.00000 0.00000 3: 0.00000 -1.00000 1.00000 0.00000 0.00000 4: 0.00000 0.00000 -1.00000 1.00000 0.00000 5: 0.00000 0.00000 0.00000 -1.00000 1.00000 L1NN U factor: Col 1 2 3 4 5 Row 1: 1.00000 -1.00000 0.00000 0.00000 0.00000 2: 0.00000 1.00000 -1.00000 0.00000 0.00000 3: 0.00000 0.00000 1.00000 -1.00000 0.00000 4: 0.00000 0.00000 0.00000 1.00000 -1.00000 5: 0.00000 0.00000 0.00000 0.00000 0.00000 L1NN LU error = 0.00000 L1PP L factor: Col 1 2 3 4 5 Row 1: 1.00000 0.00000 0.00000 0.00000 0.00000 2: -0.500000 1.00000 0.00000 0.00000 0.00000 3: 0.00000 -0.666667 1.00000 0.00000 0.00000 4: 0.00000 0.00000 -0.750000 1.00000 0.00000 5: -0.500000 -0.333333 -0.250000 -1.00000 1.00000 L1PP U factor: Col 1 2 3 4 5 Row 1: 2.00000 -1.00000 0.00000 0.00000 -1.00000 2: 0.00000 1.50000 -1.00000 0.00000 -0.500000 3: 0.00000 0.00000 1.33333 -1.00000 -0.333333 4: 0.00000 0.00000 0.00000 1.25000 -1.25000 5: 0.00000 0.00000 0.00000 0.00000 0.00000 L1PP LU error = 0.00000 Using spacing H = 0.166667 L1DD L factor: Col 1 2 3 4 5 Row 1: 6.00000 0.00000 0.00000 0.00000 0.00000 2: -3.00000 6.00000 0.00000 0.00000 0.00000 3: 0.00000 -4.00000 6.00000 0.00000 0.00000 4: 0.00000 0.00000 -4.50000 6.00000 0.00000 5: 0.00000 0.00000 0.00000 -4.80000 6.00000 L1DD U factor: Col 1 2 3 4 5 Row 1: 12.0000 -6.00000 0.00000 0.00000 0.00000 2: 0.00000 9.00000 -6.00000 0.00000 0.00000 3: 0.00000 0.00000 8.00000 -6.00000 0.00000 4: 0.00000 0.00000 0.00000 7.50000 -6.00000 5: 0.00000 0.00000 0.00000 0.00000 7.20000 L1DD LU error = 0.710543E-14 L1DN L factor: Col 1 2 3 4 5 Row 1: 6.00000 0.00000 0.00000 0.00000 0.00000 2: -3.00000 6.00000 0.00000 0.00000 0.00000 3: 0.00000 -4.00000 6.00000 0.00000 0.00000 4: 0.00000 0.00000 -4.50000 6.00000 0.00000 5: 0.00000 0.00000 0.00000 -4.80000 6.00000 L1DN U factor: Col 1 2 3 4 5 Row 1: 12.0000 -6.00000 0.00000 0.00000 0.00000 2: 0.00000 9.00000 -6.00000 0.00000 0.00000 3: 0.00000 0.00000 8.00000 -6.00000 0.00000 4: 0.00000 0.00000 0.00000 7.50000 -6.00000 5: 0.00000 0.00000 0.00000 0.00000 1.20000 L1DN LU error = 0.100486E-13 L1ND L factor: Col 1 2 3 4 5 Row 1: 6.00000 0.00000 0.00000 0.00000 0.00000 2: -6.00000 6.00000 0.00000 0.00000 0.00000 3: 0.00000 -6.00000 6.00000 0.00000 0.00000 4: 0.00000 0.00000 -6.00000 6.00000 0.00000 5: 0.00000 0.00000 0.00000 -6.00000 6.00000 L1ND U factor: Col 1 2 3 4 5 Row 1: 6.00000 -6.00000 0.00000 0.00000 0.00000 2: 0.00000 6.00000 -6.00000 0.00000 0.00000 3: 0.00000 0.00000 6.00000 -6.00000 0.00000 4: 0.00000 0.00000 0.00000 6.00000 -6.00000 5: 0.00000 0.00000 0.00000 0.00000 6.00000 L1ND LU error = 0.00000 L1NN L factor: Col 1 2 3 4 5 Row 1: 6.00000 0.00000 0.00000 0.00000 0.00000 2: -6.00000 6.00000 0.00000 0.00000 0.00000 3: 0.00000 -6.00000 6.00000 0.00000 0.00000 4: 0.00000 0.00000 -6.00000 6.00000 0.00000 5: 0.00000 0.00000 0.00000 -6.00000 6.00000 L1NN U factor: Col 1 2 3 4 5 Row 1: 6.00000 -6.00000 0.00000 0.00000 0.00000 2: 0.00000 6.00000 -6.00000 0.00000 0.00000 3: 0.00000 0.00000 6.00000 -6.00000 0.00000 4: 0.00000 0.00000 0.00000 6.00000 -6.00000 5: 0.00000 0.00000 0.00000 0.00000 0.00000 L1NN LU error = 0.00000 L1PP L factor: Col 1 2 3 4 5 Row 1: 6.00000 0.00000 0.00000 0.00000 0.00000 2: -3.00000 6.00000 0.00000 0.00000 0.00000 3: 0.00000 -4.00000 6.00000 0.00000 0.00000 4: 0.00000 0.00000 -4.50000 6.00000 0.00000 5: -3.00000 -2.00000 -1.50000 -6.00000 6.00000 L1PP U factor: Col 1 2 3 4 5 Row 1: 12.0000 -6.00000 0.00000 0.00000 -6.00000 2: 0.00000 9.00000 -6.00000 0.00000 -3.00000 3: 0.00000 0.00000 8.00000 -6.00000 -2.00000 4: 0.00000 0.00000 0.00000 7.50000 -7.50000 5: 0.00000 0.00000 0.00000 0.00000 0.00000 L1PP LU error = 0.00000 LAPLACIAN_PRB Normal end of execution. 2 November 2013 11:37:00.382 AM