Wed Oct 8 08:54:30 2025 sftpack_test(): python version: 3.10.12 numpy version: 1.26.4 Test sftpack(). c8mat_print_test c8mat_print prints an C8MAT. Here is a C8MAT: Col: 0 1 2 Row 0 : 10 1i 10 2i 10 3i 1 : 20 1i 20 2i 20 3i 2 : 30 1i 30 2i 30 3i 3 : 40 1i 40 2i 40 3i c8mat_print_some_test c8mat_print_some prints some of an C8MAT. Here is a C8MAT: Col: 3 4 5 Row 0 : 10 4i 10 5i 10 6i 1 : 20 4i 20 5i 20 6i 2 : 30 4i 30 5i 30 6i c8mat_uniform_01_test(): c8mat_uniform_01() computes a random C8MAT. 0 <= X <= 1 Random C8MAT: Col: 0 1 2 Row 0 : -0.386555 -0.637793i -0.347852 -0.631166i -0.454157 -0.0334784i 1 : 0.700289 -0.031857i 0.326829 0.773447i -0.485704 0.206265i 2 : -0.26855 -0.620421i -0.0348108 0.86501i -0.70587 0.237378i 3 : -0.591598 -0.672059i -0.00964125 -0.0249862i -0.539766 0.717979i 4 : 0.153472 0.89557i 0.537593 0.793277i -0.206373 0.444106i c8vec_indicator_test c8vec_indicator returns the indicator vector. The indicator vector: 0 1 -1 1 2 -2 2 3 -3 3 4 -4 4 5 -5 5 6 -6 6 7 -7 7 8 -8 8 9 -9 9 10 -10 c8vec_print_test c8vec_print prints an C8VEC. Here is a C8VEC: 0 1 2 1 3 4 2 5 6 3 7 8 c8vec_print_part_test c8vec_print_part prints part of a C8VEC. Part of the C8VEC: 0 1 -1 1 2 -2 2 3 -3 3 4 -4 4 5 -5 5 6 -6 6 7 -7 7 8 -8 ...... .............. .............. 99 100 -100 c8vec_uniform_01_test(): c8vec_uniform_01() computes pseudorandom complex values in the unit circle. 0 ( 0.136973, -0.178829 ) 1 ( 0.293454, -0.316582 ) 2 ( 0.834213, -0.373733 ) 3 ( -0.346910, 0.095523 ) 4 ( -0.310882, 0.692329 ) 5 ( 0.587201, -0.588787 ) 6 ( -0.366370, -0.437352 ) 7 ( -0.353469, 0.904121 ) 8 ( 0.878496, -0.327999 ) 9 ( 0.726129, 0.501914 ) i4_modp_test i4_modp factors a number into a multiple M and a positive remainder R. Number Divisor Multiple Remainder 107 50 2 7 107 -50 -2 7 -107 50 -3 43 -107 -50 3 43 Repeat using Python % Operator: 107 50 2 7 107 -50 -3 -43 -107 50 -3 43 -107 -50 2 -7 i4_wrap_test i4_wrap forces an integer to lie within given limits. ILO = 4 IHI = 8 I i4_wrap(I) -10 5 -9 6 -8 7 -7 8 -6 4 -5 5 -4 6 -3 7 -2 8 -1 4 0 5 1 6 2 7 3 8 4 4 5 5 6 6 7 7 8 8 9 4 10 5 11 6 12 7 13 8 14 4 15 5 16 6 17 7 18 8 19 4 20 5 r8vec_indicator1_test r8vec_indicator1 returns the 1-based indicator matrix. The 1-based indicator vector: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 r8vec_print_part_test r8vec_print_part prints part of an R8VEC. Lines 10:20 of the vector: 10: 11 11: 12 12: 13 13: 14 14: 15 15: 16 16: 17 17: 18 18: 19 19: 20 20: 21 c8mat_sft_test(): c8mat_sftf() computes the forward slow Fourier transform. c8mat_sftb() computes the backward slow Fourier transform. The data has dimension N1 = 10 by N2 = 4 The data X: Col: 1 2 3 Row 1 : 0.775894 0.213419i 0.357486 -0.34114i 0.918842 0.14402i 2 : -0.0237095 -0.451043i -0.303785 -0.410219i -0.312392 -0.711898i 3 : 0.589403 0.449646i -0.844714 0.0610534i 0.697684 0.654898i 4 : -0.508942 0.74883i -0.557157 0.0951573i -0.670043 0.37036i 5 : 0.415274 -0.464224i -0.560794 -0.304408i 0.262125 -0.227905i 6 : -0.0616147 0.57069i 0.274895 -0.821403i -0.227484 -0.332523i 7 : -0.384485 -0.597814i 0.403467 0.401047i 0.673372 -0.421197i 8 : -0.837989 -0.271514i -0.441035 0.123252i -0.107665 -0.901299i 9 : 0.164583 -0.157488i 0.0213087 -0.517409i 0.610969 0.274774i The Fourier coefficients Y: Col: 1 2 3 Row 1 : -1.83675 2.36399i 4.69333 -0.670218i 1.65719 1.33022i 2 : -1.45734 -1.82851i 0.144974 -3.40532i -0.717257 -4.61338i 3 : 1.57748 -0.861366i -0.346446 0.29254i 0.0616348 -1.11723i 4 : 0.246452 -0.520281i 1.4472 3.00382i -0.36798 0.00488623i 5 : -2.055 2.91331i 9.5048 1.87475i 5.05195 -1.65357i 6 : 0.906454 -1.53288i -2.39463 5.5271i -0.168485 -4.71549i 7 : -5.00723 0.112442i 1.40135 -1.71695i -3.03087 -0.990343i 8 : 1.13925 1.02703i 2.50635 -6.81946i 0.543594 -4.86243i 9 : 1.2232 -3.25576i -2.47181 2.97805i -3.47402 -2.72283i The recovered data: Col: 1 2 3 Row 1 : 0.775894 0.213419i 0.357486 -0.34114i 0.918842 0.14402i 2 : -0.0237095 -0.451043i -0.303785 -0.410219i -0.312392 -0.711898i 3 : 0.589403 0.449646i -0.844714 0.0610534i 0.697684 0.654898i 4 : -0.508942 0.74883i -0.557157 0.0951573i -0.670043 0.37036i 5 : 0.415274 -0.464224i -0.560794 -0.304408i 0.262125 -0.227905i 6 : -0.0616147 0.57069i 0.274895 -0.821403i -0.227484 -0.332523i 7 : -0.384485 -0.597814i 0.403467 0.401047i 0.673372 -0.421197i 8 : -0.837989 -0.271514i -0.441035 0.123252i -0.107665 -0.901299i 9 : 0.164583 -0.157488i 0.0213087 -0.517409i 0.610969 0.274774i c8vec_sft_test(): c8vec_sftf() computes the forward slow Fourier transform. c8vec_sftb() computes the backward slow Fourier transform. The number of data values, N = 36 The original data: 0 0.270075 0.280756 1 -0.00402792 0.039328 2 0.684406 0.597717 3 0.568348 0.751398 4 -0.114234 -0.585726 5 0.23403 0.134112 6 0.0762511 -0.32815 7 0.257063 0.449637 ...... .............. .............. 35 0.181242 -0.824326 The Fourier coefficients: 0 -4.07224 2.5227 1 3.26439 1.03899 2 2.51682 3.0043 3 -5.50364 0.859005 4 -2.79023 0.439126 5 -1.18244 1.02387 6 -0.914263 3.1314 7 -2.36718 1.08938 ...... .............. .............. 35 2.64742 -8.57544 The retrieved data: 0 0.270075 0.280756 1 -0.00402792 0.039328 2 0.684406 0.597717 3 0.568348 0.751398 4 -0.114234 -0.585726 5 0.23403 0.134112 6 0.0762511 -0.32815 7 0.257063 0.449637 ...... .............. .............. 35 0.181242 -0.824326 r8vec_sct_test(): r8vec_sct() does a forward or backward slow cosine transform. The number of data items is N = 256 The original data: 1: 1.52068 2: 3.75836 3: 3.5234 4: 0.735291 5: 2.59669 6: 4.29002 7: 2.63728 8: 0.606818 9: 3.52155 10: 2.44288 The cosine coefficients: 1: -36.6266 2: -2.14584 3: 21.9088 4: -30.2144 5: -7.39147 6: -29.2406 7: -4.02694 8: 24.4009 9: 37.2679 10: -3.79253 The retrieved data: 1: 1.52068 2: 3.75836 3: 3.5234 4: 0.735291 5: 2.59669 6: 4.29002 7: 2.63728 8: 0.606818 9: 3.52155 10: 2.44288 r8vec_sft_test(): r8vec_sft_n_test(): r8vec_sftf() computes the forward slow Fourier transform. r8vec_sftb() computes the backward slow Fourier transform. The number of data values, N = 35 The original data: 0 2.56225 1 0.0584862 2 2.66773 3 4.64085 4 1.85875 5 0.124641 6 1.43458 7 4.72354 ...... .............. 34 4.59226 A (cosine) coefficients: 0 2.44558 0 -0.156845 1 -0.0772948 2 -0.170103 3 0.533515 4 0.290374 5 -0.00700285 6 -0.761764 7 -0.0248241 8 0.256225 9 -0.0340988 10 -0.0270635 11 0.735754 12 -0.682267 13 0.346817 14 0.0497302 15 0.323934 16 -0.478415 B (sine) coefficients: 0 -0.0918886 1 0.150441 2 0.0578628 3 0.054358 4 0.179665 5 -0.300854 6 -0.166511 7 -0.878367 8 -0.65409 9 -0.0456513 10 -0.515882 11 0.262582 12 -0.153421 13 -0.334649 14 0.301802 15 -0.276847 16 -0.134756 The retrieved data: 0 2.56225 1 0.0584862 2 2.66773 3 4.64085 4 1.85875 5 0.124641 6 1.43458 7 4.72354 ...... .............. 34 4.59226 r8vec_sft_n_test(): r8vec_sftf() computes the forward slow Fourier transform. r8vec_sftb() computes the backward slow Fourier transform. The number of data values, N = 36 The original data: 0 4.14324 1 1.74496 2 2.44112 3 1.93853 4 3.72661 5 4.24533 6 0.51306 7 2.31783 ...... .............. 35 4.513 A (cosine) coefficients: 0 2.88426 0 0.346135 1 0.301986 2 0.0303274 3 -0.480357 4 0.602326 5 -0.553162 6 0.184591 7 0.551964 8 0.417552 9 -0.355302 10 0.0129612 11 -0.222748 12 -0.629324 13 0.528041 14 0.107685 15 0.604058 16 -0.218284 17 0.0305317 B (sine) coefficients: 0 -0.527944 1 -0.533816 2 0.270875 3 -0.304076 4 -0.171185 5 -0.204426 6 -0.512936 7 0.186653 8 -0.0703657 9 -0.317583 10 0.0688556 11 -0.135492 12 -0.0956744 13 0.374511 14 -0.460119 15 -0.0433828 16 0.402103 17 1.86694e-15 The retrieved data: 0 4.14324 1 1.74496 2 2.44112 3 1.93853 4 3.72661 5 4.24533 6 0.51306 7 2.31783 ...... .............. 35 4.513 r8vec_sht_test(): r8vec_sht() does a forward or backward slow Hartley transform. The number of data items is N = 17 The original data: 1: 3.49875 2: 0.576791 3: 4.61585 4: 0.0506913 5: 2.45443 6: 1.66166 7: 3.44323 8: 0.525687 9: 1.1149 10: 4.88229 The Hartley coefficients: 1: -0.0955496 2: 0.0170337 3: -0.451871 4: -1.14928 5: 1.07159 6: -0.851803 7: 1.65826 8: -1.49258 9: -3.28873 10: -0.976174 The retrieved data: 1: 3.49875 2: 0.576791 3: 4.61585 4: 0.0506913 5: 2.45443 6: 1.66166 7: 3.44323 8: 0.525687 9: 1.1149 10: 4.88229 r8vec_sqct_test(): r8vec_sqctf() does a forward slow quarter wave cosine transform r8vec_sqctb() does a backward slow quarter wave cosine transform. The number of data items is N = 256 The original data: 1: 4.65427 2: 3.93133 3: 2.16219 4: 4.96266 5: 4.9149 6: 3.86402 7: 4.40447 8: 4.59981 9: 3.85177 10: 3.47594 The cosine coefficients: 1: 0.079361 2: -0.0163581 3: 0.116578 4: 0.0461197 5: 0.112029 6: 0.0622258 7: -0.0512876 8: 0.0975716 9: -0.0428539 10: 0.0913004 The retrieved data: 1: 4.65427 2: 3.93133 3: 2.16219 4: 4.96266 5: 4.9149 6: 3.86402 7: 4.40447 8: 4.59981 9: 3.85177 10: 3.47594 r8vec_sqst_test(): r8vec_sqstf() does a forward slow quarter wave sine transform r8vec_sqstb() does a backward slow quarter wave sine transform. The number of data items is N = 256 The original data: 1: 1.02123 2: 1.72684 3: 3.07074 4: 1.68084 5: 3.83283 6: 4.55365 7: 4.07262 8: 0.243234 9: 0.59087 10: 3.18894 The sine coefficients: 1: 0.0905035 2: -0.535147 3: 0.00952461 4: -0.299281 5: 0.0303768 6: -0.222623 7: 0.0870694 8: -0.245521 9: 0.0310625 10: -0.197448 The retrieved data: 1: 1.02123 2: 1.72684 3: 3.07074 4: 1.68084 5: 3.83283 6: 4.55365 7: 4.07262 8: 0.243234 9: 0.59087 10: 3.18894 r8vec_sst_test(): r8vec_sst() does a forward or backward slow sine transform. The number of data items is N = 256 The original data: 1: 1.94739 2: 3.52113 3: 4.35607 4: 3.90033 5: 3.09963 6: 0.00796204 7: 1.60971 8: 3.61855 9: 0.503568 10: 4.78813 The sine coefficients: 1: 25.5655 2: 240.781 3: -27.5947 4: 118.869 5: 0.308711 6: 95.5489 7: -1.86886 8: 137.458 9: -29.4924 10: 61.0337 The retrieved data: 1: 1.94739 2: 3.52113 3: 4.35607 4: 3.90033 5: 3.09963 6: 0.00796204 7: 1.60971 8: 3.61855 9: 0.503568 10: 4.78813 r8vec_swt_test(): r8vec_swtf() computes the forward slow wavelet transform. r8vec_swtb() computes the backward slow wavelet transform. The number of data values, N = 36 The original data: 1: 4.8833 2: 1.14238 3: 4.06831 4: 3.04701 5: 4.0847 6: 1.19921 7: 0.272566 8: 0.230725 9: 4.26204 10: 0.541233 I S(I) D(I) 0 4.120387 2.867365 1 4.109310 2.185739 2 5.122034 0.483603 3 0.307347 0.438355 4 3.711952 2.081683 5 1.424560 0.853588 6 5.521588 -2.045811 7 3.622924 -0.241685 8 5.254378 0.625640 9 2.197982 2.245906 10 3.367265 -1.989085 11 3.626968 0.421407 12 3.345094 0.705294 13 2.667087 1.339738 14 3.192360 -0.260013 15 3.995688 -2.184707 16 4.244194 1.476369 17 3.386328 -1.401602 The retrieved data: 1: 4.8833 2: 1.14238 3: 4.06831 4: 3.04701 5: 4.0847 6: 1.19921 7: 0.272566 8: 0.230725 9: 4.26204 10: 0.541233 sftpack_test(): Normal end of execution. Wed Oct 8 08:54:31 2025