cmlabs.linalg.test_thomas

cmlabs.linalg.test_thomas()[source]

Test Thomas Algorithm.

Example

>>> # Test N: Thomas Algorithm
>>> # A:  [1 2 3 4 5]
>>> # B:  [1 2 3 4 5 6]
>>> # C:  [1 2 3 4 5]
>>> # F:  [1 2 3 4 5 6]
>>> thomas(A, B, C, F)
>>> # res:  [-8.774  9.774 -4.387 -1.129  5.419 -3.516]
>>> # res_scipy:  [-8.774  9.774 -4.387 -1.129  5.419 -3.516]