Sun May 24 17:47:53 2026

c83_test():
  numpy version:  1.26.4
  python version: 3.10.12
  Test c83().

c83_mv_test():
  c83_mv() multiplies a complex triadigonal matrix times a vector.
  The tridiagonal matrix A:
array([[0.        +0.j        , 0.18422485-0.797476j  ,
        0.16823333-0.7282516j ],
       [0.17921122-0.77577289j, 0.18537338-0.80244779j,
        0.21000522-0.90907454j],
       [0.19278974-0.83455184j, 0.10196643-0.44139419j,
        0.21113511-0.91396565j],
       [0.20339823-0.88047408j, 0.21411431-0.92686205j,
        0.        +0.j        ]])
  The random vector x:
array([ 0.57158354-0.01333733j,  0.50613297-0.00831643j,
       -0.63288276-0.34149783j, -0.14221992+0.47742661j])
  The product A*x:
array([ 0.17891674-0.85241586j, -0.23341439-0.36385115j,
        0.27489683+0.00499905j, -0.03368377+0.74037303j])

c83_indicator_test():
  c83_indicator() computes the indicator complex tridiagonal matrix.
  The matrix A:
array([[0.+0.j, 1.+1.j, 1.+2.j],
       [2.+1.j, 2.+2.j, 2.+3.j],
       [3.+2.j, 3.+3.j, 3.+4.j],
       [4.+3.j, 4.+4.j, 0.+0.j]])

c83_random_test():
  c83_random() computes a random complex tridiagonal matrix.
  The matrix A:
array([[0.        +0.j        , 0.84360425+0.11734004j,
        0.49283151+0.06854976j],
       [0.73120939+0.10170663j, 0.89607807+0.12463882j,
        0.27984953+0.03892531j],
       [0.96850765+0.13471332j, 0.19920644+0.02770836j,
        0.59782088+0.08315312j],
       [0.7036197 +0.09786907j, 0.85037629+0.11828199j,
        0.        +0.j        ]])

c8vec_uniform_01_test():
  c8vec_uniform_01() computes pseudorandom complex values
  in the unit circle.

       0  ( -0.221013, -0.904765 )
       1  ( 0.269543, -0.416644 )
       2  ( -0.561408, -0.646793 )
       3  ( -0.248255, 0.563299 )
       4  ( -0.438734, 0.760952 )
       5  ( -0.440082, 0.604671 )
       6  ( -0.320073, 0.560923 )
       7  ( -0.889572, 0.386665 )
       8  ( -0.020693, -0.982460 )
       9  ( 0.387219, -0.260490 )

cyclic_reduction_test():
  c83_cr_fa() factors a complex tridiagonal matrix;
  c83_cr_sl() solves a factored system.

  Matrix order N =  10
  Desired solution x:
array([ 1. +10.j,  2. +20.j,  3. +30.j,  4. +40.j,  5. +50.j,  6. +60.j,
        7. +70.j,  8. +80.j,  9. +90.j, 10.+100.j])
  Right hand side b:
array([  20.  -2.j,   40.  -4.j,   60.  -6.j,   80.  -8.j,  100. -10.j,
        120. -12.j,  140. -14.j,  160. -16.j,  180. -18.j, -889.+200.j])
  Solution x:
array([ 1. +10.j,  2. +20.j,  3. +30.j,  4. +40.j,  5. +50.j,  6. +60.j,
        7. +70.j,  8. +80.j,  9. +90.j, 10.+100.j])

c83_test():
  Normal end of execution.
Sun May 24 17:47:53 2026