Tue Aug 16 12:06:57 2022 r83_np_test(): Python version: 3.6.9 Test r83_np(). r83_np_det_test(): r83_np_det() computes the determinant of a tridiagonal matrix factored by r83_np_fa. Matrix order N = 10 The factored R83 matrix: Col: 0 1 2 3 4 Row --- 0: 2 -1 1: -0.5 1.5 -1 2: -0.666667 1.33333 -1 3: -0.75 1.25 -1 4: -0.8 1.2 5: -0.833333 Col: 5 6 7 8 9 Row --- 4: -1 5: 1.16667 -1 6: -0.857143 1.14286 -1 7: -0.875 1.125 -1 8: -0.888889 1.11111 -1 9: -0.9 1.1 r83_np_det() computes determinant = 10.999999999999998 Exact determinant is = 11 r83_np_fa_test(): r83_np_fa() factors a tridiagonal matrix with no pivoting Matrix order N = 10 The tridiagonal matrix: Col: 0 1 2 3 4 Row --- 0: 0.543081 0.885338 1: 0.299569 0.342346 0.568441 2: 0.292904 0.293687 0.649803 3: 0.0724743 0.563777 0.0547356 4: 0.854588 0.404391 5: 0.970195 Col: 5 6 7 8 9 Row --- 4: 0.851008 5: 0.956989 0.977807 6: 0.776726 0.957145 0.200744 7: 0.383979 0.488641 0.518432 8: 0.344702 0.849382 0.0555662 9: 0.545662 0.440529 Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 Solution to tranposed system: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 r83_np_fs_test(): r83_np_fs() factors and solves a tridiagonal linear system. Matrix order N = 10 Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 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: 0: 1.37674 1.37674 1: 2.84778 2.84778 2: 2.8591 2.8591 3: 6.8071 6.8071 4: 4.87787 4.87787 5: 8.13162 8.13162 6: 12.9133 12.9133 7: 5.86067 5.86067 8: 15.2796 15.2796 9: 3.75188 3.75188 A'*x and (PLU)'*x 0: 0.713395 0.713395 1: 2.83836 2.83836 2: 4.31231 4.31231 3: 2.69593 2.69593 4: 9.64095 9.64095 5: 8.08702 8.08702 6: 6.72158 6.72158 7: 9.50925 9.50925 8: 17.9782 17.9782 9: 7.06252 7.06252 r83_np_sl_test(): r83_np_sl() solves a linear system after the tridiagonal matrix has been factored by r83_np_fa(). Matrix order N = 10 The tridiagonal matrix: Col: 0 1 2 3 4 Row --- 0: 0.799796 0.127965 1: 0.639924 0.154573 0.602818 2: 0.849302 0.34825 0.891448 3: 0.358966 0.853392 0.794225 4: 0.281341 0.6686 5: 0.537398 Col: 5 6 7 8 9 Row --- 4: 0.321289 5: 0.897704 0.0481679 6: 0.955134 0.569751 0.965453 7: 0.285188 0.309461 0.574781 8: 0.634949 0.948259 0.952792 9: 0.894094 0.818907 Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 Solution to tranposed system: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 r83_np_test(): Normal end of execution. Tue Aug 16 12:06:57 2022