11 April 2023 2:28:54.591 PM r83_np_test(): FORTRAN90 version: Test r83_np(). R83_NP_DET_TEST R83_NP_DET computes the determinant of an R83 matrix that was factored by R83_NP_FA. Matrix order N = 10 The factored R83 matrix: Col: 1 2 3 4 5 Row --- 1 2.00000 -1.00000 2 -0.500000 1.50000 -1.00000 3 -0.666667 1.33333 -1.00000 4 -0.750000 1.25000 -1.00000 5 -0.800000 1.20000 6 -0.833333 Col: 6 7 8 9 10 Row --- 5 -1.00000 6 1.16667 -1.00000 7 -0.857143 1.14286 -1.00000 8 -0.875000 1.12500 -1.00000 9 -0.888889 1.11111 -1.00000 10 -0.900000 1.10000 R83_NP_DET computes determinant = 11.0000 Exact determinant = 11.0000 R83_NP_FA_TEST R83_NP_FA factors a tridiagonal matrix with no pivoting, after which, R83_NP_SL can solve linear systems. Matrix order N = 10 The tridiagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.829509 2 0.956318 0.561695 0.661187E-01 3 0.415307 0.257578 0.438290E-01 4 0.109957 0.633966 0.449539 5 0.617272E-01 0.401306 6 0.754673 Col: 6 7 8 9 10 Row --- 5 0.797287 6 0.183837E-02 0.350752 7 0.897504 0.945448E-01 0.859097 8 0.136169E-01 0.840847 0.751236E-02 9 0.123104 0.260303 0.113664 10 0.912484 0.351629 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 Solution to transposed system: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 R83_NP_FS_TEST R83_NP_FS factors and solves a tridiagonal linear system. Matrix order N = 10 What is this? Col: 1 2 3 4 5 Row --- 1 0.218418 0.829509 2 0.956318 0.561695 0.661187E-01 3 0.415307 0.257578 0.438290E-01 4 0.109957 0.633966 0.449539 5 0.617272E-01 0.401306 6 0.754673 Col: 6 7 8 9 10 Row --- 5 0.797287 6 0.183837E-02 0.350752 7 0.897504 0.945448E-01 0.859097 8 0.136169E-01 0.840847 0.751236E-02 9 0.123104 0.260303 0.113664 10 0.912484 0.351629 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 R83_NP_FSS_TEST R83_NP_FSS factors a tridiagonal linear system without pivoting, and solves multiple linear systems. Matrix order N = 10 Solutions: Col: 1 2 Row --- 1 0.508301E-15 1.00000 2 1.00000 2.00000 3 2.00000 3.00000 4 3.00000 4.00000 5 4.00000 5.00000 6 5.00000 6.00000 7 6.00000 7.00000 8 7.00000 8.00000 9 8.00000 9.00000 10 9.00000 10.0000 R83_NP_ML_TEST R83_NP_ML computes A*x or A'*x where A has been factored by R83_FA. Matrix order N = 10 A*x and PLU*x: 1 1.87744 1.87744 2 2.27806 2.27806 3 1.77866 1.77866 4 5.11343 5.11343 5 7.03716 7.03716 6 6.23966 6.23966 7 12.9196 12.9196 8 6.88971 6.88971 9 4.46420 4.46420 10 11.7286 11.7286 A'*x and (PLU)'*x 1 1.35715 1.35715 2 4.39904 4.39904 3 6.03519 6.03519 4 5.59777 5.59777 5 5.81525 5.81525 6 6.74571 6.74571 7 7.50541 7.50541 8 12.9291 12.9291 9 14.2607 14.2607 10 9.79448 9.79448 R83_NP_SL_TEST R83_NP_SL solves a linear system that has been factored by R83_NP_FA. Matrix order N = 10 The tridiagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.829509 2 0.956318 0.561695 0.661187E-01 3 0.415307 0.257578 0.438290E-01 4 0.109957 0.633966 0.449539 5 0.617272E-01 0.401306 6 0.754673 Col: 6 7 8 9 10 Row --- 5 0.797287 6 0.183837E-02 0.350752 7 0.897504 0.945448E-01 0.859097 8 0.136169E-01 0.840847 0.751236E-02 9 0.123104 0.260303 0.113664 10 0.912484 0.351629 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 Solution to transposed system: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 r83_np_test(): Normal end of execution. 11 April 2023 2:28:54.591 PM