9 May 2025 9:14:05.084 PM fftpack5_test(): Fortran90 version Test fftpack5(). cfft1_test For complex single precision fast Fourier transforms, 1D, CFFT1I initializes the transform, CFFT1F does a forward transform; CFFT1B does a backward transform. The number of data items is N = 4096 LENSAV = 8208 LENWRK = 8192 The original data: 1: 0.186384 -0.107270E-01 2: 0.153105 -0.732133E-02 3: -0.400113 -0.340877 4: 0.116562 -0.132258 5: -0.105840 -0.900648E-01 6: 0.955881E-01 -0.396886 7: -0.868819 -0.192563 8: -0.416828 -0.315053 ........ .............. .............. 4096: -0.127660 -0.453328E-01 The FFT coefficients: 1: 0.227399E-02 0.749479E-02 2: 0.453210E-02 0.403717E-02 3: 0.659163E-02 0.621137E-02 4: 0.872048E-02 0.415778E-02 5: 0.378696E-02 -0.635333E-02 6: -0.104865E-01 0.344505E-02 7: -0.550223E-02 -0.746335E-03 8: -0.267224E-02 0.357275E-02 ........ .............. .............. 4096: 0.398312E-02 0.228714E-02 The retrieved data: 1: 0.186384 -0.107270E-01 2: 0.153105 -0.732128E-02 3: -0.400113 -0.340877 4: 0.116562 -0.132258 5: -0.105840 -0.900649E-01 6: 0.955881E-01 -0.396886 7: -0.868819 -0.192563 8: -0.416828 -0.315053 ........ .............. .............. 4096: -0.127660 -0.453326E-01 cfft2_test For complex single precision fast Fourier transforms, 2D, CFFT2I initializes the transform, CFFT2F does a forward transform; CFFT2B does a backward transform. The data is stored in an L by M array, with L = 32 M = 64 LENSAV = 211 LENWRK = 4096 Part of the original data: Col: 1 2 3 4 Row --- 1 0.372 -0.527 -0.716 -0.345 -0.224 0.613 0.162 -0.219 2 0.442 0.524 -0.840E-01 0.749 -0.588 0.770 -0.477 -0.370 3 0.475 0.469 -0.749 -0.469 0.449 0.893E-01-0.332E-01 0.862 4 -0.936 -0.184 0.847 -0.356 -0.953E-01-0.731 -0.318 -0.501 5 0.916 0.360 0.397 -0.823 0.348 -0.698 0.474 0.583 Col: 5 Row --- 1 0.733 -0.523 2 -0.734 -0.628 3 0.349 -0.743E-01 4 -0.463 -0.704E-01 5 -0.575 -0.400 Part of the FFT coefficients: Col: 1 2 3 4 Row --- 1 0.176E-01 0.454E-02 0.267E-02-0.130E-01 0.863E-02-0.102E-01 0.148E-02-0.694E-02 2 0.772E-02 0.829E-02-0.326E-02-0.176E-02 0.102E-01 0.488E-02 0.813E-02-0.281E-03 3 0.940E-02-0.930E-02 0.601E-02-0.228E-02-0.102E-01 0.238E-02 0.863E-02 0.127E-02 4 0.619E-03 0.104E-01 0.109E-01-0.911E-02-0.297E-03 0.219E-01-0.780E-02-0.111E-01 5 0.156E-02 0.672E-02-0.213E-01 0.417E-02-0.639E-02-0.731E-02 0.634E-02-0.941E-02 Col: 5 Row --- 1 0.147E-01 0.126E-01 2 -0.598E-02 0.593E-02 3 -0.111E-01-0.795E-02 4 -0.367E-02 0.191E-01 5 0.273E-02-0.163E-01 Part of the retrieved data: Col: 1 2 3 4 Row --- 1 0.372 -0.527 -0.716 -0.345 -0.224 0.613 0.162 -0.219 2 0.442 0.524 -0.840E-01 0.749 -0.588 0.770 -0.477 -0.370 3 0.475 0.469 -0.749 -0.469 0.449 0.893E-01-0.332E-01 0.862 4 -0.936 -0.184 0.847 -0.356 -0.953E-01-0.731 -0.318 -0.501 5 0.916 0.360 0.397 -0.823 0.348 -0.698 0.474 0.583 Col: 5 Row --- 1 0.733 -0.523 2 -0.734 -0.628 3 0.349 -0.743E-01 4 -0.463 -0.704E-01 5 -0.575 -0.400 cfftm_test For complex single precision fast Fourier transforms, 1D, multiple CFFTMI initializes the transform, CFFTMF does a forward transform; CFFTMB does a backward transform. The number of sequences is LOT = 6 The length of each sequence is N = 32 LENC = 192 LENSAV = 73 LENWRK = 384 Part of the original data: Col: 1 2 3 4 Row --- 1 0.674E-01 0.262 -0.647 0.175E-01-0.277 -0.164 0.453 -0.275 2 0.297 -0.230 -0.278 0.139 -0.134 -0.687 0.599E-01-0.754 3 0.460 0.444 -0.453 0.246 -0.606 -0.298 0.608 0.490 4 0.392 0.266 0.645 0.634E-01 0.108 0.776 -0.398 0.856 5 0.159 0.762 0.871 0.406E-01 0.315E-01 0.160 -0.238 0.596 Col: 5 Row --- 1 0.614 -0.547 2 -0.461 -0.480 3 -0.705 -0.356 4 -0.265 0.480 5 -0.832 0.765E-01 Part of the FFT coefficients: Col: 1 2 3 4 Row --- 1 0.145 -0.335E-01-0.319E-01 0.315E-01-0.125E-02 0.188E-01-0.106 0.152 2 0.113 0.348E-01 0.869E-01-0.531E-01 0.838E-01-0.768E-01-0.302E-01-0.935E-01 3 0.540E-01-0.248E-01-0.216 -0.104 0.139 -0.476E-02 0.130 -0.617E-01 4 0.854E-01 0.288E-01 0.959E-01-0.152 -0.468E-01-0.291E-01 0.103 -0.981E-01 5 0.320E-02-0.427E-01 0.227E-02-0.387E-01-0.818E-01 0.584E-01 0.243E-02-0.528E-01 Col: 5 Row --- 1 -0.951E-01-0.723E-01 2 0.322E-01 0.229 3 -0.282E-01 0.144 4 0.152 -0.746E-01 5 -0.547E-01 0.120E-01 Part of the retrieved data: Col: 1 2 3 4 Row --- 1 0.674E-01 0.262 -0.647 0.175E-01-0.277 -0.164 0.453 -0.275 2 0.297 -0.230 -0.278 0.139 -0.134 -0.687 0.599E-01-0.754 3 0.460 0.444 -0.453 0.246 -0.606 -0.298 0.608 0.490 4 0.392 0.266 0.645 0.634E-01 0.108 0.776 -0.398 0.856 5 0.159 0.762 0.871 0.406E-01 0.315E-01 0.160 -0.238 0.596 Col: 5 Row --- 1 0.614 -0.547 2 -0.461 -0.480 3 -0.705 -0.356 4 -0.265 0.480 5 -0.832 0.765E-01 cosq1_test For real single precision fast cosine transforms, 1D, COSQ1I initializes the transform, COSQ1F does a forward transform; COSQ1B does a backward transform. The number of data items is N = 4096 LENSAV = 8208 LENWRK = 4096 The original data: 1: 0.706044 2: 0.292692 3: 0.225037 4: 0.837388 5: 0.232472 6: 0.919223E-01 7: 0.566928 8: 0.337849 ........ .............. 4096: 0.969999 The FFT coefficients: 1: 0.632126 2: -0.217230 3: 0.137904 4: -0.865639E-01 5: 0.705304E-01 6: -0.595712E-01 7: 0.409178E-01 8: -0.452051E-01 ........ .............. 4096: -0.398737E-02 The retrieved data: 1: 0.706043 2: 0.292692 3: 0.225037 4: 0.837388 5: 0.232472 6: 0.919224E-01 7: 0.566928 8: 0.337849 ........ .............. 4096: 0.969999 cosqm_test(): For real single precision fast cosine transform, 1D, multiple COSQMI initializes the transform, COSQMF does a forward transform; COSQMB does a backward transform. The number of sequences is LOT = 6 The length of each sequence is N = 32 LENR = 192 LENSAV = 73 LENWRK = 192 Part of the original data: Col 1 2 3 4 5 Row 1 0.951819 0.111432 0.718531 0.956881 0.125137 2 0.160636 0.173878 0.384385 0.266248 0.412317 3 0.549095 0.582308 0.380510 0.456122 0.981501 4 0.934515 0.713881 0.301994 0.405541 0.367536 5 0.557744 0.161113E-01 0.550743 0.659004 0.567438 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.757674 0.538201 0.635623 0.573292 0.592222 2 -0.320999 -0.215495 -0.172701 -0.194795 -0.132948 3 0.275829 0.170690 0.182863 0.108360 0.171101 4 -0.175555 -0.194525 -0.107255 -0.899716E-01 -0.769766E-01 5 0.124022 0.177872E-02 -0.840820E-01 0.541914E-01 0.609913E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.951819 0.111432 0.718531 0.956881 0.125137 2 0.160636 0.173878 0.384385 0.266248 0.412317 3 0.549095 0.582308 0.380510 0.456122 0.981501 4 0.934515 0.713881 0.301994 0.405541 0.367536 5 0.557744 0.161114E-01 0.550742 0.659004 0.567438 cost1_test For real single precision fast cosine transforms, 1D, COST1I initializes the transform, COST1F does a forward transform; COST1B does a backward transform. The number of data items is N = 4096 LENSAV = 8208 LENWRK = 4095 The original data: 1: 0.626457 2: 0.559142 3: 0.985273 4: 0.511272 5: 0.494264 6: 0.572294 7: 0.918338 8: 0.319133 ........ .............. 4096: 0.381620 The FFT coefficients: 1: 0.500399 2: 0.602747E-02 3: -0.879294E-02 4: 0.142742E-02 5: -0.813149E-02 6: 0.946650E-02 7: 0.122038E-01 8: 0.881228E-03 ........ .............. 4096: -0.518478E-02 The retrieved data: 1: 0.626584 2: 0.559193 3: 0.985277 4: 0.511250 5: 0.494284 6: 0.572311 7: 0.918332 8: 0.319128 ........ .............. 4096: 0.381490 costm_test(): For real single precision fast cosine transforms, 1D, multiple COSTMI initializes the transform, COSTMF does a forward transform; COSTMB does a backward transform. The number of sequences is LOT = 6 The length of each sequence is N = 32 LENR = 192 LENSAV = 73 LENWRK = 198 Part of the original data: Col 1 2 3 4 5 Row 1 0.654684 0.771326 0.936701 0.701799 0.752087 2 0.585156 0.996604 0.703819 0.315298 0.820469 3 0.874833 0.884333 0.494847 0.724295 0.327840 4 0.112074 0.739293 0.358927E-01 0.996162 0.382762 5 0.824903 0.470468 0.348855 0.680367 0.287583 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.612138 0.572121 0.507711 0.440315 0.557659 2 0.446469E-01 0.892617E-02 -0.995454E-02 0.124544 0.110315E-01 3 -0.106717 0.100231 -0.978530E-02 0.115690 -0.206837E-01 4 -0.633240E-01 0.140863 -0.853348E-01 -0.690699E-02 0.984862E-01 5 0.839630E-01 0.958036E-01 0.565770E-01 0.191823E-01 -0.156480 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.654685 0.771325 0.936702 0.701800 0.752086 2 0.585156 0.996604 0.703819 0.315298 0.820469 3 0.874832 0.884332 0.494847 0.724294 0.327840 4 0.112075 0.739294 0.358937E-01 0.996161 0.382763 5 0.824901 0.470468 0.348854 0.680367 0.287583 dcosq1_test For real double precision fast cosine transform, 1D, DCOSQ1I initializes the transform, DCOSQ1F does a forward transform; DCOSQ1B does a backward transform. The number of data items is N = 4096 The original data: 1: 0.172295 2: 0.920697E-01 3: 0.318619 4: 0.238460 5: 0.227795 6: 0.248477 7: 0.253556 8: 0.626010 ........ .............. 4096: 0.119007 The FFT coefficients: 1: 0.639307 2: -0.214121 3: 0.126543 4: -0.926403E-01 5: 0.660180E-01 6: -0.667198E-01 7: 0.569842E-01 8: -0.448732E-01 ........ .............. 4096: 0.790149E-02 The retrieved data: 1: 0.172295 2: 0.920697E-01 3: 0.318619 4: 0.238460 5: 0.227795 6: 0.248477 7: 0.253556 8: 0.626010 ........ .............. 4096: 0.119007 dcost1_test For double precision fast cosine transforms, 1D, DCOST1I initializes the transforms, DCOST1F does a forward transforms; DCOST1B does a backward transforms. The number of data items is N = 4096 The original data: 1: 0.519663 2: 0.692729 3: 0.955357 4: 0.666035E-02 5: 0.561664 6: 0.370985 7: 0.184872E-01 8: 0.784523 ........ .............. 4096: 0.699190 The FFT coefficients: 1: 0.502141 2: 0.134623E-02 3: 0.601895E-02 4: -0.315048E-02 5: -0.544620E-03 6: 0.567853E-03 7: -0.688361E-02 8: 0.574635E-02 ........ .............. 4096: -0.254470E-02 The retrieved data: 1: 0.519663 2: 0.692729 3: 0.955357 4: 0.666035E-02 5: 0.561664 6: 0.370985 7: 0.184872E-01 8: 0.784523 ........ .............. 4096: 0.699190 dfft1_test For real double precision fast Fourier transform, 1D, DFFT1I initializes the transform, DFFT1F does a forward transform; DFFT1B does a backward transform. The number of data items is N = 4096 LENSAV = 4112 LENWRK = 4096 The original data: 1: 0.209585 2: 0.742596 3: 0.731516 4: 0.285760 5: 0.924417 6: 0.817090 7: 0.654617 8: 0.958241 ........ .............. 4096: 0.324592 The FFT coefficients: 1: 0.502937 2: -0.537305E-03 3: -0.265300E-02 4: 0.726416E-02 5: 0.473760E-02 6: -0.787573E-02 7: 0.743814E-03 8: -0.506395E-02 ........ .............. 4096: 0.105674E-01 The retrieved data: 1: 0.209585 2: 0.742596 3: 0.731516 4: 0.285760 5: 0.924417 6: 0.817090 7: 0.654617 8: 0.958241 ........ .............. 4096: 0.324592 DSINT1_TEST For double precision fast sine transforms, 1D, DSINT1I initializes the transforms, DSINT1F does a forward transforms; DSINT1B does a backward transforms. The number of data items is N = 4096 The original data: 1: 0.344096 2: 0.154080 3: 0.793470 4: 0.233125 5: 0.689708 6: 0.234382 7: 0.971090 8: 0.490974 ........ .............. 4096: 0.108788 The FFT coefficients: 1: 0.628989 2: 0.820065E-03 3: 0.211881 4: -0.384422E-02 5: 0.119869 6: 0.485903E-02 7: 0.883894E-01 8: -0.202404E-02 ........ .............. 4096: 0.246687E-02 The retrieved data: 1: 0.344096 2: 0.154080 3: 0.793470 4: 0.233125 5: 0.689708 6: 0.234382 7: 0.971090 8: 0.490974 ........ .............. 4096: 0.108788 rfft1_test For real single precision fast Fourier transforms, 1D, RFFT1I initializes the transform, RFFT1F does a forward transform; RFFT1B does a backward transform. The number of data items is N = 4096 LENSAV = 4112 LENWRK = 4096 The original data: 1: 0.266527 2: 0.533744 3: 0.767436E-01 4: 0.202738 5: 0.960959 6: 0.278685 7: 0.450325E-02 8: 0.323797 ........ .............. 4096: 0.786361 The FFT coefficients: 1: 0.504168 2: -0.532625E-02 3: 0.130908E-01 4: 0.525131E-02 5: -0.215366E-02 6: 0.350756E-02 7: 0.220849E-02 8: 0.619585E-02 ........ .............. 4096: 0.597599E-02 The retrieved data: 1: 0.266527 2: 0.533744 3: 0.767436E-01 4: 0.202738 5: 0.960959 6: 0.278685 7: 0.450313E-02 8: 0.323797 ........ .............. 4096: 0.786361 rfft2_test For real single precision fast Fourier transform, 2D, RFFT2I initializes the transform, RFFT2F does a forward transform; RFFT2B does a backward transform. The L by M data is stored in an LDIM by M array, with L = 32 LDIM = 34 M = 64 LENSAV = 207 LENWRK = 4352 Part of the original data: Col 1 2 3 4 5 Row 1 0.622865E-01 0.320458 0.896809 0.813736 0.497912 2 0.732348 0.722404 0.264486 0.225999 0.823771 3 0.844100 0.592508 0.557380 0.261820 0.478705 4 0.238016 0.156924 0.895208 0.462358 0.520321 5 0.762519 0.891760 0.495887 0.811086 0.513774 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.505362 0.992125E-02 -0.831857E-03 0.405797E-03 0.172506E-02 2 0. 0.398182E-02 -0.423837E-03 0.262383E-02 -0.432488E-02 3 0.901262E-02 -0.112430E-02 -0.969403E-02 -0.656296E-02 -0.185144E-01 4 -0.113886E-02 0.490277E-03 0.799331E-02 0.140899E-02 -0.175332E-01 5 -0.590240E-02 -0.308549E-02 -0.522697E-02 0.560926E-02 -0.169952E-02 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.622866E-01 0.320458 0.896809 0.813736 0.497911 2 0.732348 0.722404 0.264486 0.225999 0.823771 3 0.844100 0.592508 0.557380 0.261821 0.478705 4 0.238016 0.156924 0.895208 0.462358 0.520321 5 0.762519 0.891760 0.495887 0.811086 0.513774 rfftm_test() For real single precision fast Fourier transform, 1D, multiple RFFTMI initializes the transform, RFFTMF does a forward transform; RFFTMB does a backward transform. The number of sequences is LOT = 6 The length of each sequence is N = 32 LENR = 192 LENSAV = 41 LENWRK = 192 Part of the original data: Col 1 2 3 4 5 Row 1 0.408042 0.640253 0.222508 0.899309 0.497544 2 0.468089 0.102367 0.238483 0.683750 0.134382 3 0.307225 0.480574E-01 0.517002 0.228738 0.381537 4 0.600365E-01 0.210910E-01 0.193528 0.301932 0.241840 5 0.272247 0.351568 0.883346 0.533116 0.853918 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.440958 0.463954 0.446600 0.496083 0.447818 2 -0.591730E-01 0.100173E-01 -0.511691E-01 -0.319143E-01 -0.219869E-01 3 0.443544E-01 -0.204759 -0.122072 -0.785558E-02 0.157843E-01 4 0.250543E-01 -0.338897E-01 0.283327E-01 -0.104266 0.636178E-01 5 -0.178334 -0.142662 0.258904E-02 0.860088E-01 -0.522408E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.408042 0.640253 0.222507 0.899309 0.497544 2 0.468089 0.102367 0.238483 0.683750 0.134382 3 0.307225 0.480574E-01 0.517002 0.228738 0.381537 4 0.600365E-01 0.210910E-01 0.193528 0.301932 0.241840 5 0.272247 0.351568 0.883346 0.533116 0.853918 sinq1_test For real single precision fast sine transforms, 1D, SINQ1I initializes the transform, SINQ1F does a forward transform; SINQ1B does a backward transform. The number of data items is N = 4096 LENSAV = 8208 LENWRK = 4096 The original data: 1: 0.977781 2: 0.538126 3: 0.269470 4: 0.369379 5: 0.953215 6: 0.745763 7: 0.939792 8: 0.232244 ........ .............. 4096: 0.115617 The FFT coefficients: 1: 0.640417 2: 0.217110 3: 0.128037 4: 0.100200 5: 0.752530E-01 6: 0.471672E-01 7: 0.371505E-01 8: 0.437916E-01 ........ .............. 4096: -0.576210E-02 The retrieved data: 1: 0.977781 2: 0.538126 3: 0.269470 4: 0.369379 5: 0.953215 6: 0.745763 7: 0.939792 8: 0.232244 ........ .............. 4096: 0.115617 sinqm_test() For real single precision fast sine transforms, 1D, multiple SINQMI initializes the transform, SINQMF does a forward transform; SINQMB does a backward transform. The number of sequences is LOT = 6 The length of each sequence is N = 32 LENR = 192 LENSAV = 73 LENWRK = 192 Part of the original data: Col 1 2 3 4 5 Row 1 0.914544 0.824281 0.784808 0.458871 0.164998 2 0.859975 0.385673 0.909163 0.297626 0.940197 3 0.728160 0.892325 0.532980 0.463218 0.158476E-01 4 0.627823 0.782550 0.607877E-01 0.548621 0.649564 5 0.335164 0.849782 0.643413 0.374109 0.513583E-01 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.605518 0.571246 0.842222 0.730025 0.540563 2 0.256463 0.167831 0.197549 0.114784 0.292728 3 0.176980 0.159716 0.264892 0.241633 0.125879 4 0.108374 0.155128 -0.471758E-01 0.335264E-01 0.137625 5 0.114735 0.756007E-01 -0.190540E-01 0.184967E-01 0.360905E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.914544 0.824281 0.784808 0.458871 0.164998 2 0.859974 0.385673 0.909163 0.297626 0.940197 3 0.728160 0.892325 0.532980 0.463218 0.158477E-01 4 0.627823 0.782550 0.607877E-01 0.548621 0.649564 5 0.335164 0.849782 0.643413 0.374109 0.513582E-01 sint1_test For real single precision fast sine transforms, 1D, SINT1I initializes the transform, SINT1F does a forward transform; SINT1B does a backward transform. The number of data items is N = 4096 LENSAV = 6160 LENWRK = 8194 The original data: 1: 0.191825 2: 0.686901 3: 0.224347 4: 0.305209 5: 0.455097 6: 0.217538 7: 0.685046E-01 8: 0.750933 ........ .............. 4096: 0.412250 The FFT coefficients: 1: 0.634556 2: 0.292887E-02 3: 0.211986 4: -0.906883E-02 5: 0.126241 6: 0.601850E-02 7: 0.942761E-01 8: -0.289637E-02 ........ .............. 4096: -0.111413E-01 The retrieved data: 1: 0.195018 2: 0.686806 3: 0.225335 4: 0.305053 5: 0.455868 6: 0.217496 7: 0.689398E-01 8: 0.750926 ........ .............. 4096: 0.415313 sintm_test For real single precision fast sine transforms, 1D, multiple SINTMI initializes the transform, SINTMF does a forward transform; SINTMB does a backward transform. The number of sequences is LOT = 6 The length of each sequence is N = 32 LENR = 192 LENSAV = 57 LENWRK = 408 Part of the original data: Col 1 2 3 4 5 Row 1 0.517969 0.685957 0.896646 0.569269 0.902925E-01 2 0.145529 0.798701 0.244238 0.170255 0.728352 3 0.923690 0.997748 0.258038 0.141020 0.862115 4 0.249461 0.307719 0.431161 0.645979 0.248401 5 0.243665 0.489461 0.708804 0.462194 0.604454 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.577505 0.537015 0.636630 0.499258 0.545787 2 0.534118E-01 0.102753E-01 0.308733E-02 0.107174 0.281272E-01 3 0.167100 0.159030 0.197303 0.250751 0.117732 4 -0.250818E-01 0.628145E-01 -0.294144E-01 -0.811610E-01 0.124471 5 0.134519 0.198189 0.141446 0.940636E-01 0.175470 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.517967 0.685956 0.896645 0.569270 0.902921E-01 2 0.145530 0.798701 0.244238 0.170255 0.728352 3 0.923690 0.997748 0.258038 0.141021 0.862116 4 0.249461 0.307719 0.431161 0.645980 0.248401 5 0.243664 0.489461 0.708804 0.462195 0.604454 zfft1_test For complex double precision fast Fourier transforms, 1D, ZFFT1I initializes the transform, ZFFT1F does a forward transform; ZFFT1B does a backward transform. The number of data items is N = 4096 LENSAV = 8204 LENWRK = 8192 The original data: 1: 0.606234E-01 0.989959 2: -0.884357E-01 -0.194141 3: -0.309992 -0.831882 4: -0.920826 0.338469 5: 0.254542E-01 -0.253750 6: -0.650790E-01 -0.366139E-02 7: 0.746834 0.428833 8: 0.290077 -0.460828E-01 ........ .............. .............. 4096: 0.160006 -0.190727 The FFT coefficients: 1: 0.398121E-02 0.718901E-02 2: 0.923545E-03 0.295203E-02 3: 0.989183E-04 0.231283E-02 4: 0.326863E-02 -0.390890E-02 5: -0.866010E-02 -0.417581E-02 6: 0.458634E-02 -0.846394E-02 7: -0.996783E-02 -0.370458E-02 8: -0.246905E-03 0.745489E-02 ........ .............. .............. 4096: -0.443660E-02 -0.979399E-02 The retrieved data: 1: 0.606234E-01 0.989959 2: -0.884357E-01 -0.194141 3: -0.309992 -0.831882 4: -0.920826 0.338469 5: 0.254542E-01 -0.253750 6: -0.650790E-01 -0.366139E-02 7: 0.746834 0.428833 8: 0.290077 -0.460828E-01 ........ .............. .............. 4096: 0.160006 -0.190727 zfft2_test For complex double precision fast Fourier transforms, 2D, ZFFT2I initializes the transform, ZCFFT2F does a forward transform; ZFFT2B does a backward transform. The data is stored in an L by M array, with L = 32 M = 64 LENSAV = 207 LENWRK = 4096 Part of the original data: Col: 1 2 3 4 Row --- 1 -0.399 -0.830E-01-0.186E-01 0.789 -0.558 -0.309 0.275 -0.564 2 0.189 -0.962 0.172 -0.261 0.598E-01-0.337 -0.573E-01-0.158 3 0.808 -0.459 0.747 0.311 0.754 -0.614 -0.745 0.574 4 0.626 0.178 0.602 -0.428 0.441 0.592 0.198 0.880 5 0.184E-01 0.120 0.534 -0.549 0.678E-01-0.913 -0.682 0.614 Col: 5 Row --- 1 0.347 0.713 2 0.772 -0.156 3 0.266 -0.223 4 -0.452 -0.680 5 -0.890 0.264 Part of the FFT coefficients: Col: 1 2 3 4 Row --- 1 0.930E-02-0.146E-01-0.102E-01-0.738E-02 0.184E-01 0.250E-02 0.721E-02-0.120E-01 2 -0.129E-01-0.143E-01 0.597E-02-0.111E-01 0.259E-02-0.128E-01-0.160E-03-0.141E-01 3 0.697E-02 0.223E-01-0.367E-02 0.871E-02 0.105E-01 0.493E-02-0.837E-02-0.248E-02 4 0.896E-02-0.101E-01-0.148E-01 0.783E-02 0.355E-02-0.519E-04 0.166E-01-0.595E-02 5 -0.357E-02-0.971E-02 0.163E-01-0.669E-03-0.136E-01 0.273E-03 0.156E-01-0.896E-03 Col: 5 Row --- 1 0.213E-02 0.121E-01 2 -0.198E-02 0.644E-03 3 -0.898E-02 0.770E-02 4 0.276E-02-0.298E-01 5 0.522E-02-0.127E-01 Part of the retrieved data: Col: 1 2 3 4 Row --- 1 -0.399 -0.830E-01-0.186E-01 0.789 -0.558 -0.309 0.275 -0.564 2 0.189 -0.962 0.172 -0.261 0.598E-01-0.337 -0.573E-01-0.158 3 0.808 -0.459 0.747 0.311 0.754 -0.614 -0.745 0.574 4 0.626 0.178 0.602 -0.428 0.441 0.592 0.198 0.880 5 0.184E-01 0.120 0.534 -0.549 0.678E-01-0.913 -0.682 0.614 Col: 5 Row --- 1 0.347 0.713 2 0.772 -0.156 3 0.266 -0.223 4 -0.452 -0.680 5 -0.890 0.264 zfftm_test For complex double precision fast Fourier transforms, 1D, multiple ZFFTMI initializes the transform, ZFFTMF does a forward transform; ZFFTMB does a backward transform. The number of sequences is LOT = 6 The length of each sequence is N = 32 LENC = 192 LENSAV = 71 LENWRK = 384 Part of the original data: Col: 1 2 3 4 Row --- 1 -0.527E-01-0.352 0.248 0.487 0.207 -0.468E-01 0.169 -0.529 2 -0.478 -0.834 0.386 0.128 0.277 -0.110 -0.437 -0.596E-01 3 0.576 -0.535E-01 0.113 -0.444 -0.122 0.294 0.287 0.319 4 0.725 0.341 -0.478 0.256 -0.496 -0.359 0.796 0.364 5 0.155 0.141 -0.150 -0.667 -0.143E-01-0.341 -0.152 -0.344 Col: 5 Row --- 1 0.925 0.338 2 -0.822 0.396 3 -0.321E-01 0.936 4 0.719E-01-0.310 5 -0.821 0.204E-01 Part of the FFT coefficients: Col: 1 2 3 4 Row --- 1 0.479E-01-0.122 0.517E-02-0.142E-01 0.674E-01-0.824E-01-0.381E-03 0.443E-01 2 0.745E-01-0.108 0.439E-01 0.391E-01-0.365E-01 0.416E-01-0.297E-01-0.126 3 0.207 0.736E-01 0.562E-01-0.563E-01-0.972E-02 0.149 0.806E-01-0.622E-01 4 -0.745E-01-0.101 0.324E-01 0.127 0.147E-01-0.179E-01 0.127E-01 0.149 5 -0.242 0.534E-01 0.821E-01 0.152 0.247E-01 0.114 -0.182E-01-0.145 Col: 5 Row --- 1 -0.133 0.790E-01 2 -0.442E-01 0.101E-01 3 0.268E-01 0.109 4 0.698E-01-0.122 5 0.201 0.452E-01 Part of the retrieved data: Col: 1 2 3 4 Row --- 1 -0.527E-01-0.352 0.248 0.487 0.207 -0.468E-01 0.169 -0.529 2 -0.478 -0.834 0.386 0.128 0.277 -0.110 -0.437 -0.596E-01 3 0.576 -0.535E-01 0.113 -0.444 -0.122 0.294 0.287 0.319 4 0.725 0.341 -0.478 0.256 -0.496 -0.359 0.796 0.364 5 0.155 0.141 -0.150 -0.667 -0.143E-01-0.341 -0.152 -0.344 Col: 5 Row --- 1 0.925 0.338 2 -0.822 0.396 3 -0.321E-01 0.936 4 0.719E-01-0.310 5 -0.821 0.204E-01 FFTPACK5_TEST Normal end of execution. 9 May 2025 9:14:05.100 PM