Tue May 20 22:31:15 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.191786 0.634101i 0.796857 0.2807i -0.621464 -0.115131i 1 : -0.625562 0.0232906i -0.433846 0.841585i 0.0853358 0.063406i 2 : -0.157328 -0.142581i 0.112121 0.466053i 0.189949 -0.0169383i 3 : 0.453755 -0.682581i 0.602077 -0.0375217i 0.506716 0.0958542i 4 : -0.988998 0.0227887i 0.199689 -0.127136i -0.705247 0.663769i 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.703488, -0.101436 ) 1 ( 0.886855, 0.295222 ) 2 ( -0.621403, -0.175209 ) 3 ( -0.127679, -0.377972 ) 4 ( 0.671591, -0.079729 ) 5 ( 0.631027, -0.627670 ) 6 ( -0.211201, 0.764469 ) 7 ( -0.384555, 0.862947 ) 8 ( 0.110678, 0.227886 ) 9 ( 0.250638, 0.787305 ) 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.0550011 -0.653537i -0.935301 -0.301224i 0.451225 -0.633868i 2 : 0.406875 0.418721i 0.68544 -0.339374i 0.651431 -0.716728i 3 : -0.68719 0.173971i 0.687163 -0.0864407i 0.776512 -0.248541i 4 : -0.535405 -0.745923i -0.576141 0.669377i 0.293045 0.0378282i 5 : -0.868205 -0.186121i -0.128381 0.0181366i 0.546387 0.585552i 6 : -0.181372 0.805624i -0.2898 -0.175048i 0.0451374 -0.704524i 7 : 0.649453 -0.440995i 0.586943 -0.773831i 0.688086 0.276878i 8 : 0.341398 -0.32277i 0.749815 -0.411236i 0.125491 -0.413465i 9 : -0.610948 0.329214i 0.0950146 -0.454288i 0.857333 0.00584544i The Fourier coefficients Y: Col: 1 2 3 Row 1 : 0.317991 0.801094i -4.54288 -3.43815i -1.13948 -1.4017i 2 : 2.09197 -2.92046i 1.58399 -3.08815i 3.77393 2.00504i 3 : 0.695708 7.28079i -1.44788 -1.39596i 2.92585 -2.07796i 4 : 0.105076 -1.476i 1.74696 -3.08488i -2.62002 -0.778302i 5 : -2.49983 3.55523i 1.62467 1.8488i 2.70853 -4.66959i 6 : -0.133301 4.48887i 3.71089 -1.54123i -3.12462 0.236011i 7 : 2.38728 -3.89083i 4.3882 1.17459i -1.60415 1.9195i 8 : 7.72603 -1.10033i -0.0402889 1.2096i 0.771804 -3.37151i 9 : -5.25714 0.120763i 0.00413727 -3.28428i -0.646928 -4.19499i The recovered data: Col: 1 2 3 Row 1 : 0.0550011 -0.653537i -0.935301 -0.301224i 0.451225 -0.633868i 2 : 0.406875 0.418721i 0.68544 -0.339374i 0.651431 -0.716728i 3 : -0.68719 0.173971i 0.687163 -0.0864407i 0.776512 -0.248541i 4 : -0.535405 -0.745923i -0.576141 0.669377i 0.293045 0.0378282i 5 : -0.868205 -0.186121i -0.128381 0.0181366i 0.546387 0.585552i 6 : -0.181372 0.805624i -0.2898 -0.175048i 0.0451374 -0.704524i 7 : 0.649453 -0.440995i 0.586943 -0.773831i 0.688086 0.276878i 8 : 0.341398 -0.32277i 0.749815 -0.411236i 0.125491 -0.413465i 9 : -0.610948 0.329214i 0.0950146 -0.454288i 0.857333 0.00584544i 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.162065 0.155607 1 -0.527798 0.428846 2 -0.193969 0.189625 3 0.723838 0.324935 4 -0.00439686 -0.991974 5 0.0726502 0.111433 6 -0.341765 -0.0673825 7 0.869814 0.107129 ...... .............. .............. 35 -0.234283 -0.412902 The Fourier coefficients: 0 -2.39185 -3.49588 1 1.67554 4.06839 2 2.22704 0.971547 3 -3.32863 0.623228 4 -3.27249 1.69992 5 4.48129 2.16391 6 -3.8453 -1.53963 7 -2.65696 -0.441187 ...... .............. .............. 35 1.4473 4.63075 The retrieved data: 0 0.162065 0.155607 1 -0.527798 0.428846 2 -0.193969 0.189625 3 0.723838 0.324935 4 -0.00439686 -0.991974 5 0.0726502 0.111433 6 -0.341765 -0.0673825 7 0.869814 0.107129 ...... .............. .............. 35 -0.234283 -0.412902 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: 3.66 2: 3.06489 3: 3.94013 4: 2.32216 5: 3.96903 6: 3.69133 7: 0.853217 8: 3.06545 9: 1.50731 10: 4.17373 The cosine coefficients: 1: -9.77939 2: -5.04841 3: -21.0621 4: 6.42003 5: 24.7636 6: -17.8466 7: -42.0073 8: 0.946816 9: 20.7204 10: 54.2881 The retrieved data: 1: 3.66 2: 3.06489 3: 3.94013 4: 2.32216 5: 3.96903 6: 3.69133 7: 0.853217 8: 3.06545 9: 1.50731 10: 4.17373 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 1.04785 1 3.98493 2 1.75816 3 4.79185 4 2.21475 5 2.28512 6 0.593425 7 3.23896 ...... .............. 34 4.4543 A (cosine) coefficients: 0 2.78276 0 -0.230716 1 0.00461761 2 0.482358 3 0.00647015 4 -0.221277 5 -0.256881 6 -0.166788 7 0.401603 8 0.644847 9 -0.333127 10 -0.737994 11 -0.0401493 12 -0.213685 13 -0.182407 14 -0.526414 15 0.222842 16 -0.588207 B (sine) coefficients: 0 0.214415 1 0.344395 2 -0.476001 3 0.0846331 4 0.533832 5 -0.212286 6 0.261741 7 -0.362398 8 -0.560517 9 -0.297185 10 -0.259043 11 0.374947 12 0.025364 13 0.290934 14 0.28419 15 0.348885 16 0.487895 The retrieved data: 0 1.04785 1 3.98493 2 1.75816 3 4.79185 4 2.21475 5 2.28512 6 0.593425 7 3.23896 ...... .............. 34 4.4543 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 1.6643 1 0.903724 2 3.63495 3 0.884005 4 1.97993 5 2.58308 6 4.0457 7 2.61798 ...... .............. 35 4.74059 A (cosine) coefficients: 0 2.57521 0 -0.0255637 1 -0.243729 2 -0.0218152 3 0.345047 4 0.148311 5 0.101767 6 -0.210109 7 0.301526 8 -0.506898 9 -0.240477 10 -0.253124 11 -0.279965 12 -0.186048 13 -0.522554 14 -0.121214 15 0.526757 16 0.526831 17 -0.249652 B (sine) coefficients: 0 0.341862 1 -0.422546 2 0.355475 3 -0.547273 4 -0.403696 5 0.234824 6 -0.349523 7 -0.197436 8 0.0374834 9 -0.393918 10 -0.297592 11 -0.638676 12 0.273944 13 -0.279744 14 0.499642 15 -0.209939 16 -0.318567 17 1.80105e-15 The retrieved data: 0 1.6643 1 0.903724 2 3.63495 3 0.884005 4 1.97993 5 2.58308 6 4.0457 7 2.61798 ...... .............. 35 4.74059 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: 0.684271 2: 3.28179 3: 1.50717 4: 3.01881 5: 2.65721 6: 3.18091 7: 0.994704 8: 2.96222 9: 0.32286 10: 4.32716 The Hartley coefficients: 1: 0.649028 2: -0.849203 3: -1.77548 4: -0.781096 5: -0.607236 6: -2.55932 7: 0.515406 8: -0.990421 9: 2.43759 10: -0.941873 The retrieved data: 1: 0.684271 2: 3.28179 3: 1.50717 4: 3.01881 5: 2.65721 6: 3.18091 7: 0.994704 8: 2.96222 9: 0.32286 10: 4.32716 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: 3.05112 2: 4.12284 3: 1.30559 4: 0.0716279 5: 3.07494 6: 0.563032 7: 3.11697 8: 2.94388 9: 1.02977 10: 2.3813 The cosine coefficients: 1: -0.0234024 2: 0.0809312 3: -0.14988 4: -0.0104379 5: 0.0282848 6: 0.0474746 7: -0.141897 8: -0.0459661 9: 0.0181968 10: -0.098281 The retrieved data: 1: 3.05112 2: 4.12284 3: 1.30559 4: 0.0716279 5: 3.07494 6: 0.563032 7: 3.11697 8: 2.94388 9: 1.02977 10: 2.3813 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.22066 2: 3.46409 3: 3.87613 4: 1.25903 5: 4.27163 6: 0.479969 7: 2.59354 8: 1.96516 9: 1.31266 10: 2.06722 The sine coefficients: 1: 0.0165593 2: -0.486886 3: -0.00466284 4: -0.296404 5: -0.114676 6: -0.291007 7: 0.114811 8: -0.107205 9: 0.0554018 10: -0.196949 The retrieved data: 1: 1.22066 2: 3.46409 3: 3.87613 4: 1.25903 5: 4.27163 6: 0.479969 7: 2.59354 8: 1.96516 9: 1.31266 10: 2.06722 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: 3.42128 2: 0.266643 3: 0.339069 4: 2.82651 5: 4.50306 6: 2.27172 7: 1.94924 8: 3.9731 9: 4.07213 10: 1.6242 The sine coefficients: 1: 2.62763 2: 216.45 3: 28.6942 4: 147.063 5: -8.95594 6: 115.472 7: 28.2458 8: 63.4907 9: 45.4735 10: 96.6203 The retrieved data: 1: 3.42128 2: 0.266643 3: 0.339069 4: 2.82651 5: 4.50306 6: 2.27172 7: 1.94924 8: 3.9731 9: 4.07213 10: 1.6242 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: 1.47608 2: 3.75082 3: 1.56824 4: 1.85092 5: 2.36206 6: 1.71666 7: 4.71244 8: 1.29013 9: 2.96133 10: 4.21901 I S(I) D(I) 0 4.225108 -2.732639 1 3.232565 -1.418996 2 2.644767 0.429369 3 4.677064 1.608889 4 3.816157 1.629413 5 3.402298 -1.840958 6 3.242302 1.394571 7 5.591748 -1.393041 8 4.278823 -1.690952 9 4.061462 -1.057182 10 5.223344 1.523444 11 2.165939 0.173383 12 4.408188 1.226434 13 2.728054 2.061979 14 4.741349 0.515168 15 2.602335 -1.329416 16 1.702246 0.331250 17 3.249110 1.509351 The retrieved data: 1: 1.47608 2: 3.75082 3: 1.56824 4: 1.85092 5: 2.36206 6: 1.71666 7: 4.71244 8: 1.29013 9: 2.96133 10: 4.21901 sftpack_test(): Normal end of execution. Tue May 20 22:31:16 2025