15 September 2021 8:27:07.182 AM LAPLACIAN_TEST FORTRAN90 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. -1. 0. 0. 0. 2: -1. 2. -1. 0. 0. 3: 0. -1. 2. -1. 0. 4: 0. 0. -1. 2. -1. 5: 0. 0. 0. -1. 2. L1DN: Col 1 2 3 4 5 Row 1: 2. -1. 0. 0. 0. 2: -1. 2. -1. 0. 0. 3: 0. -1. 2. -1. 0. 4: 0. 0. -1. 2. -1. 5: 0. 0. 0. -1. 1. L1ND: Col 1 2 3 4 5 Row 1: 1. -1. 0. 0. 0. 2: -1. 2. -1. 0. 0. 3: 0. -1. 2. -1. 0. 4: 0. 0. -1. 2. -1. 5: 0. 0. 0. -1. 2. L1NN: Col 1 2 3 4 5 Row 1: 1. -1. 0. 0. 0. 2: -1. 2. -1. 0. 0. 3: 0. -1. 2. -1. 0. 4: 0. 0. -1. 2. -1. 5: 0. 0. 0. -1. 1. L1PP: Col 1 2 3 4 5 Row 1: 2. -1. 0. 0. -1. 2: -1. 2. -1. 0. 0. 3: 0. -1. 2. -1. 0. 4: 0. 0. -1. 2. -1. 5: -1. 0. 0. -1. 2. Using spacing H = 0.166667 L1DD: Col 1 2 3 4 5 Row 1: 72. -36. 0. 0. 0. 2: -36. 72. -36. 0. 0. 3: 0. -36. 72. -36. 0. 4: 0. 0. -36. 72. -36. 5: 0. 0. 0. -36. 72. L1DN: Col 1 2 3 4 5 Row 1: 72. -36. 0. 0. 0. 2: -36. 72. -36. 0. 0. 3: 0. -36. 72. -36. 0. 4: 0. 0. -36. 72. -36. 5: 0. 0. 0. -36. 36. L1ND: Col 1 2 3 4 5 Row 1: 36. -36. 0. 0. 0. 2: -36. 72. -36. 0. 0. 3: 0. -36. 72. -36. 0. 4: 0. 0. -36. 72. -36. 5: 0. 0. 0. -36. 72. L1NN: Col 1 2 3 4 5 Row 1: 36. -36. 0. 0. 0. 2: -36. 72. -36. 0. 0. 3: 0. -36. 72. -36. 0. 4: 0. 0. -36. 72. -36. 5: 0. 0. 0. -36. 36. L1PP: Col 1 2 3 4 5 Row 1: 72. -36. 0. 0. -36. 2: -36. 72. -36. 0. 0. 3: 0. -36. 72. -36. 0. 4: 0. 0. -36. 72. -36. 5: -36. 0. 0. -36. 72. 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.108488E-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.147139E-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.163165E-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.968414E-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.920414E-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.474347E-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.446274E-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.567516E-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.415822E-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. -1. 0. 0. 0. 2: -1. 2. -1. 0. 0. 3: 0. -1. 2. -1. 0. 4: 0. 0. -1. 2. -1. 5: 0. 0. 0. -1. 2. 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. 0.666667 0.333333 3: 0.500000 1. 1.50000 1. 0.500000 4: 0.333333 0.666667 1. 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. -1. 0. 0. 0. 2: -1. 2. -1. 0. 0. 3: 0. -1. 2. -1. 0. 4: 0. 0. -1. 2. -1. 5: 0. 0. 0. -1. 1. L1DN_INVERSE: Col 1 2 3 4 5 Row 1: 1. 1. 1. 1. 1. 2: 1. 2. 2. 2. 2. 3: 1. 2. 3. 3. 3. 4: 1. 2. 3. 4. 4. 5: 1. 2. 3. 4. 5. L1DN inverse error = 0.00000 L1ND: Col 1 2 3 4 5 Row 1: 1. -1. 0. 0. 0. 2: -1. 2. -1. 0. 0. 3: 0. -1. 2. -1. 0. 4: 0. 0. -1. 2. -1. 5: 0. 0. 0. -1. 2. L1ND_INVERSE: Col 1 2 3 4 5 Row 1: 5. 4. 3. 2. 1. 2: 4. 4. 3. 2. 1. 3: 3. 3. 3. 2. 1. 4: 2. 2. 2. 2. 1. 5: 1. 1. 1. 1. 1. L1ND inverse error = 0.00000 Using spacing H = 0.166667 L1DD: Col 1 2 3 4 5 Row 1: 72. -36. 0. 0. 0. 2: -36. 72. -36. 0. 0. 3: 0. -36. 72. -36. 0. 4: 0. 0. -36. 72. -36. 5: 0. 0. 0. -36. 72. 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.116739E-14 L1DN: Col 1 2 3 4 5 Row 1: 72. -36. 0. 0. 0. 2: -36. 72. -36. 0. 0. 3: 0. -36. 72. -36. 0. 4: 0. 0. -36. 72. -36. 5: 0. 0. 0. -36. 36. 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.236037E-14 L1ND: Col 1 2 3 4 5 Row 1: 36. -36. 0. 0. 0. 2: -36. 72. -36. 0. 0. 3: 0. -36. 72. -36. 0. 4: 0. 0. -36. 72. -36. 5: 0. 0. 0. -36. 72. 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.257035E-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. 0. 0. 2: 0. 1.22474 -0.816497 0. 0. 3: 0. 0. 1.15470 -0.866025 0. 4: 0. 0. 0. 1.11803 -0.894427 5: 0. 0. 0. 0. 1.09545 L1DD Cholesky error = 0.800593E-15 L1DN Cholesky factor: Col 1 2 3 4 5 Row 1: 1.41421 -0.707107 0. 0. 0. 2: 0. 1.22474 -0.816497 0. 0. 3: 0. 0. 1.15470 -0.866025 0. 4: 0. 0. 0. 1.11803 -0.894427 5: 0. 0. 0. 0. 0.447214 L1DN Cholesky error = 0.777156E-15 L1ND Cholesky factor: Col 1 2 3 4 5 Row 1: 1. -1. 0. 0. 0. 2: 0. 1. -1. 0. 0. 3: 0. 0. 1. -1. 0. 4: 0. 0. 0. 1. -1. 5: 0. 0. 0. 0. 1. L1ND Cholesky error = 0.00000 L1NN Cholesky factor: Col 1 2 3 4 5 Row 1: 1. -1. 0. 0. 0. 2: 0. 1. -1. 0. 0. 3: 0. 0. 1. -1. 0. 4: 0. 0. 0. 1. -1. 5: 0. 0. 0. 0. 0. L1NN Cholesky error = 0.00000 L1PP Cholesky factor: Col 1 2 3 4 5 Row 1: 1.41421 -0.707107 0. 0. -0.707107 2: 0. 1.22474 -0.816497 0. -0.408248 3: 0. 0. 1.15470 -0.866025 -0.288675 4: 0. 0. 0. 1.11803 -1.11803 5: 0. 0. 0. 0. 0. L1PP Cholesky error = 0.102056E-14 Using spacing H = 0.166667 L1DD Cholesky factor: Col 1 2 3 4 5 Row 1: 8.48528 -4.24264 0. 0. 0. 2: 0. 7.34847 -4.89898 0. 0. 3: 0. 0. 6.92820 -5.19615 0. 4: 0. 0. 0. 6.70820 -5.36656 5: 0. 0. 0. 0. 6.57267 L1DD Cholesky error = 0.375984E-13 L1DN Cholesky factor: Col 1 2 3 4 5 Row 1: 8.48528 -4.24264 0. 0. 0. 2: 0. 7.34847 -4.89898 0. 0. 3: 0. 0. 6.92820 -5.19615 0. 4: 0. 0. 0. 6.70820 -5.36656 5: 0. 0. 0. 0. 2.68328 L1DN Cholesky error = 0.375984E-13 L1ND Cholesky factor: Col 1 2 3 4 5 Row 1: 6. -6. 0. 0. 0. 2: 0. 6. -6. 0. 0. 3: 0. 0. 6. -6. 0. 4: 0. 0. 0. 6. -6. 5: 0. 0. 0. 0. 6. L1ND Cholesky error = 0.00000 L1NN Cholesky factor: Col 1 2 3 4 5 Row 1: 6. -6. 0. 0. 0. 2: 0. 6. -6. 0. 0. 3: 0. 0. 6. -6. 0. 4: 0. 0. 0. 6. -6. 5: 0. 0. 0. 0. 0. L1NN Cholesky error = 0.00000 L1PP Cholesky factor: Col 1 2 3 4 5 Row 1: 8.48528 -4.24264 0. 0. -4.24264 2: 0. 7.34847 -4.89898 0. -2.44949 3: 0. 0. 6.92820 -5.19615 -1.73205 4: 0. 0. 0. 6.70820 -6.70820 5: 0. 0. 0. 0. 0. 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. 0. 0. 0. 0. 2: -0.500000 1. 0. 0. 0. 3: 0. -0.666667 1. 0. 0. 4: 0. 0. -0.750000 1. 0. 5: 0. 0. 0. -0.800000 1. L1DD U factor: Col 1 2 3 4 5 Row 1: 2. -1. 0. 0. 0. 2: 0. 1.50000 -1. 0. 0. 3: 0. 0. 1.33333 -1. 0. 4: 0. 0. 0. 1.25000 -1. 5: 0. 0. 0. 0. 1.20000 L1DD LU error = 0.00000 L1DN L factor: Col 1 2 3 4 5 Row 1: 1. 0. 0. 0. 0. 2: -0.500000 1. 0. 0. 0. 3: 0. -0.666667 1. 0. 0. 4: 0. 0. -0.750000 1. 0. 5: 0. 0. 0. -0.800000 1. L1DN U factor: Col 1 2 3 4 5 Row 1: 2. -1. 0. 0. 0. 2: 0. 1.50000 -1. 0. 0. 3: 0. 0. 1.33333 -1. 0. 4: 0. 0. 0. 1.25000 -1. 5: 0. 0. 0. 0. 0.200000 L1DN LU error = 0.00000 L1ND L factor: Col 1 2 3 4 5 Row 1: 1. 0. 0. 0. 0. 2: -1. 1. 0. 0. 0. 3: 0. -1. 1. 0. 0. 4: 0. 0. -1. 1. 0. 5: 0. 0. 0. -1. 1. L1ND U factor: Col 1 2 3 4 5 Row 1: 1. -1. 0. 0. 0. 2: 0. 1. -1. 0. 0. 3: 0. 0. 1. -1. 0. 4: 0. 0. 0. 1. -1. 5: 0. 0. 0. 0. 1. L1ND LU error = 0.00000 L1NN L factor: Col 1 2 3 4 5 Row 1: 1. 0. 0. 0. 0. 2: -1. 1. 0. 0. 0. 3: 0. -1. 1. 0. 0. 4: 0. 0. -1. 1. 0. 5: 0. 0. 0. -1. 1. L1NN U factor: Col 1 2 3 4 5 Row 1: 1. -1. 0. 0. 0. 2: 0. 1. -1. 0. 0. 3: 0. 0. 1. -1. 0. 4: 0. 0. 0. 1. -1. 5: 0. 0. 0. 0. 0. L1NN LU error = 0.00000 L1PP L factor: Col 1 2 3 4 5 Row 1: 1. 0. 0. 0. 0. 2: -0.500000 1. 0. 0. 0. 3: 0. -0.666667 1. 0. 0. 4: 0. 0. -0.750000 1. 0. 5: -0.500000 -0.333333 -0.250000 -1. 1. L1PP U factor: Col 1 2 3 4 5 Row 1: 2. -1. 0. 0. -1. 2: 0. 1.50000 -1. 0. -0.500000 3: 0. 0. 1.33333 -1. -0.333333 4: 0. 0. 0. 1.25000 -1.25000 5: 0. 0. 0. 0. 0. L1PP LU error = 0.00000 Using spacing H = 0.166667 L1DD L factor: Col 1 2 3 4 5 Row 1: 1. 0. 0. 0. 0. 2: -0.500000 1. 0. 0. 0. 3: 0. -0.666667 1. 0. 0. 4: 0. 0. -0.750000 1. 0. 5: 0. 0. 0. -0.800000 1. L1DD U factor: Col 1 2 3 4 5 Row 1: 72. -36. 0. 0. 0. 2: 0. 54. -36. 0. 0. 3: 0. 0. 48. -36. 0. 4: 0. 0. 0. 45. -36. 5: 0. 0. 0. 0. 43.2000 L1DD LU error = 0.00000 L1DN L factor: Col 1 2 3 4 5 Row 1: 1. 0. 0. 0. 0. 2: -0.500000 1. 0. 0. 0. 3: 0. -0.666667 1. 0. 0. 4: 0. 0. -0.750000 1. 0. 5: 0. 0. 0. -0.800000 1. L1DN U factor: Col 1 2 3 4 5 Row 1: 72. -36. 0. 0. 0. 2: 0. 54. -36. 0. 0. 3: 0. 0. 48. -36. 0. 4: 0. 0. 0. 45. -36. 5: 0. 0. 0. 0. 7.20000 L1DN LU error = 0.00000 L1ND L factor: Col 1 2 3 4 5 Row 1: 1. 0. 0. 0. 0. 2: -1. 1. 0. 0. 0. 3: 0. -1. 1. 0. 0. 4: 0. 0. -1. 1. 0. 5: 0. 0. 0. -1. 1. L1ND U factor: Col 1 2 3 4 5 Row 1: 36. -36. 0. 0. 0. 2: 0. 36. -36. 0. 0. 3: 0. 0. 36. -36. 0. 4: 0. 0. 0. 36. -36. 5: 0. 0. 0. 0. 36. L1ND LU error = 0.00000 L1NN L factor: Col 1 2 3 4 5 Row 1: 1. 0. 0. 0. 0. 2: -1. 1. 0. 0. 0. 3: 0. -1. 1. 0. 0. 4: 0. 0. -1. 1. 0. 5: 0. 0. 0. -1. 1. L1NN U factor: Col 1 2 3 4 5 Row 1: 36. -36. 0. 0. 0. 2: 0. 36. -36. 0. 0. 3: 0. 0. 36. -36. 0. 4: 0. 0. 0. 36. -36. 5: 0. 0. 0. 0. 0. L1NN LU error = 0.00000 L1PP L factor: Col 1 2 3 4 5 Row 1: 1. 0. 0. 0. 0. 2: -0.500000 1. 0. 0. 0. 3: 0. -0.666667 1. 0. 0. 4: 0. 0. -0.750000 1. 0. 5: -0.500000 -0.333333 -0.250000 -1. 1. L1PP U factor: Col 1 2 3 4 5 Row 1: 72. -36. 0. 0. -36. 2: 0. 54. -36. 0. -18. 3: 0. 0. 48. -36. -12. 4: 0. 0. 0. 45. -45. 5: 0. 0. 0. 0. 0. L1PP LU error = 0.00000 LAPLACIAN_TEST Normal end of execution. 15 September 2021 8:27:07.183 AM