6 October 2025 6:18:18.069 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.120237 -0.290405 2: 0.660655 0.387224E-01 3: 0.373576 -0.192654 4: 0.127746 -0.939114 5: -0.101363 0.294780E-01 6: 0.432574 0.136442 7: -0.565756 0.203338 8: -0.362821 -0.724086E-01 ........ .............. .............. 4096: -0.314248 -0.255194 The FFT coefficients: 1: 0.245496E-02 0.752691E-02 2: 0.539668E-03 0.556899E-02 3: -0.164676E-02 0.122886E-01 4: -0.392745E-02 -0.996706E-02 5: 0.321448E-02 -0.440404E-02 6: 0.561489E-03 0.175387E-02 7: -0.112102E-01 0.105776E-02 8: -0.305711E-02 -0.118830E-03 ........ .............. .............. 4096: 0.439979E-02 0.375810E-02 The retrieved data: 1: 0.120237 -0.290405 2: 0.660655 0.387224E-01 3: 0.373576 -0.192654 4: 0.127746 -0.939113 5: -0.101363 0.294780E-01 6: 0.432574 0.136442 7: -0.565756 0.203338 8: -0.362821 -0.724086E-01 ........ .............. .............. 4096: -0.314248 -0.255194 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.194 -0.957E-01-0.142 0.478 0.284 0.464 -0.615 -0.306E-01 2 0.764 -0.588 0.443 0.319 0.720 -0.159 -0.577 0.615 3 0.902 -0.112 0.310 -0.697 -0.291 -0.587 -0.564 0.784 4 0.710E-01 0.232 0.856 0.403 -0.128 0.602 -0.602E-01-0.650 5 0.964E-02 0.198 0.618 -0.310 0.102 0.565 0.807 0.305 Col: 5 Row --- 1 0.457 -0.524 2 0.300 0.316 3 0.164 0.697 4 -0.741 -0.330 5 -0.381 0.424 Part of the FFT coefficients: Col: 1 2 3 4 Row --- 1 -0.743E-02 0.178E-01 0.146E-01 0.777E-02 0.252E-01-0.203E-02-0.111E-01 0.114E-01 2 -0.556E-02-0.579E-02 0.646E-02 0.464E-02-0.286E-02-0.157E-01 0.138E-01-0.650E-02 3 0.109E-02-0.830E-02 0.141E-01-0.349E-02-0.290E-03 0.387E-02 0.424E-02-0.970E-02 4 0.149E-01-0.253E-02 0.893E-02 0.456E-03-0.786E-02-0.402E-02-0.155E-01-0.108E-01 5 0.136E-01 0.561E-02 0.161E-01 0.838E-02 0.228E-02-0.379E-02-0.350E-02-0.713E-02 Col: 5 Row --- 1 -0.159E-01 0.987E-02 2 -0.194E-02 0.869E-02 3 -0.174E-02-0.308E-01 4 -0.555E-02 0.107E-02 5 0.274E-02 0.183E-02 Part of the retrieved data: Col: 1 2 3 4 Row --- 1 0.194 -0.957E-01-0.142 0.478 0.284 0.464 -0.615 -0.306E-01 2 0.764 -0.588 0.443 0.319 0.720 -0.159 -0.577 0.615 3 0.902 -0.112 0.310 -0.697 -0.291 -0.587 -0.564 0.784 4 0.710E-01 0.232 0.856 0.403 -0.128 0.602 -0.602E-01-0.650 5 0.964E-02 0.198 0.618 -0.310 0.102 0.565 0.807 0.305 Col: 5 Row --- 1 0.457 -0.524 2 0.300 0.316 3 0.164 0.697 4 -0.741 -0.330 5 -0.381 0.424 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.275 0.749 0.928 0.261 0.166 0.826 -0.725 -0.676 2 0.236 -0.614 0.227 0.887 0.358 -0.256 0.158 -0.837 3 -0.199 0.910E-01-0.729 0.556 0.545E-01 0.576 0.560 -0.548 4 0.661 -0.269 0.471 -0.616 -0.806 -0.283 0.402 -0.883 5 0.472 -0.546 -0.914 0.180 -0.138 0.607 -0.720 -0.380 Col: 5 Row --- 1 -0.164E-01 0.392 2 -0.418 0.400 3 -0.574 -0.791 4 0.510 -0.714 5 0.589 0.331 Part of the FFT coefficients: Col: 1 2 3 4 Row --- 1 0.414E-01-0.886E-01-0.106 0.433E-01 0.120 0.174 -0.241E-01-0.315E-01 2 0.153 0.105 0.138E-01 0.125 0.144E-01 0.917E-01-0.173 -0.747E-01 3 0.472E-01 0.401E-02 0.308E-01 0.780E-01 0.111E-01 0.440E-01-0.380E-01 0.889E-02 4 -0.820E-01 0.250E-01 0.133E-01 0.386E-01-0.850E-01-0.679E-02-0.861E-01-0.383E-01 5 0.153E-01-0.358E-01 0.248E-01 0.852E-02 0.252E-01 0.801E-01-0.203 -0.200E-01 Col: 5 Row --- 1 -0.474E-01-0.635E-01 2 -0.769E-01-0.177E-01 3 -0.307E-01 0.149 4 -0.107 0.473E-01 5 -0.117 -0.800E-01 Part of the retrieved data: Col: 1 2 3 4 Row --- 1 0.275 0.749 0.928 0.261 0.166 0.826 -0.725 -0.676 2 0.236 -0.614 0.227 0.887 0.358 -0.256 0.158 -0.837 3 -0.199 0.910E-01-0.729 0.556 0.545E-01 0.576 0.560 -0.548 4 0.661 -0.269 0.471 -0.616 -0.806 -0.283 0.402 -0.883 5 0.472 -0.546 -0.914 0.180 -0.138 0.607 -0.720 -0.380 Col: 5 Row --- 1 -0.164E-01 0.392 2 -0.418 0.400 3 -0.574 -0.791 4 0.510 -0.714 5 0.589 0.331 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.237550 2: 0.143220 3: 0.557595 4: 0.389008 5: 0.614324 6: 0.882307 7: 0.719050 8: 0.209187 ........ .............. 4096: 0.598063 The FFT coefficients: 1: 0.630870 2: -0.209008 3: 0.127423 4: -0.938022E-01 5: 0.565914E-01 6: -0.516153E-01 7: 0.451972E-01 8: -0.398003E-01 ........ .............. 4096: -0.977099E-03 The retrieved data: 1: 0.237550 2: 0.143220 3: 0.557595 4: 0.389008 5: 0.614324 6: 0.882307 7: 0.719050 8: 0.209187 ........ .............. 4096: 0.598063 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.393214E-01 0.868386 0.665147 0.816272 0.725566 2 0.565648 0.565542 0.824494 0.268706 0.383186 3 0.822565 0.976242 0.880242 0.398226 0.983264 4 0.316492E-01 0.626385 0.646968 0.373778 0.613818 5 0.522669 0.294560E-01 0.951458 0.424692 0.664426 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.683634 0.633683 0.665772 0.662227 0.619501 2 -0.287773 -0.220318 -0.190355 -0.147278 -0.138604 3 0.946720E-01 0.253703 0.206820 0.586150E-01 0.223379 4 -0.126428 -0.282037E-01 0.107696E-02 -0.139732 -0.193715E-01 5 0.894790E-01 0.661535E-01 0.120048 0.264252E-01 0.752912E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.393213E-01 0.868386 0.665147 0.816272 0.725566 2 0.565648 0.565542 0.824494 0.268706 0.383186 3 0.822565 0.976242 0.880242 0.398226 0.983264 4 0.316492E-01 0.626385 0.646969 0.373778 0.613818 5 0.522669 0.294560E-01 0.951458 0.424692 0.664426 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.806917 2: 0.647127 3: 0.492528 4: 0.232713 5: 0.715556 6: 0.366751 7: 0.493511 8: 0.548274 ........ .............. 4096: 0.599735E-01 The FFT coefficients: 1: 0.491246 2: 0.643513E-02 3: -0.104461E-01 4: -0.503332E-05 5: -0.107483E-02 6: 0.490507E-03 7: 0.358180E-02 8: 0.593174E-02 ........ .............. 4096: 0.595394E-02 The retrieved data: 1: 0.807229 2: 0.647092 3: 0.492560 4: 0.232691 5: 0.715531 6: 0.366755 7: 0.493512 8: 0.548272 ........ .............. 4096: 0.596665E-01 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.410271 0.309388 0.960547 0.959253E-01 0.454824 2 0.925642E-02 0.366944 0.575383 0.749190 0.285373 3 0.312434 0.286483 0.391445 0.690858 0.320870E-02 4 0.382324E-01 0.437629 0.716160 0.138847 0.815756 5 0.484219 0.455396 0.574479 0.716698E-02 0.246331 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.501462 0.445887 0.472604 0.472824 0.424484 2 -0.466600E-01 -0.538361E-01 0.730490E-01 -0.421746E-01 0.662763E-01 3 -0.103352 -0.757943E-02 0.535738E-01 -0.963373E-01 0.333589E-02 4 -0.213829E-02 -0.403055E-02 -0.554944E-01 0.155219E-01 -0.302751E-01 5 -0.170117 0.103010E-01 0.810304E-01 0.799423E-01 -0.810844E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.410271 0.309387 0.960547 0.959253E-01 0.454823 2 0.925605E-02 0.366945 0.575383 0.749190 0.285373 3 0.312435 0.286484 0.391445 0.690858 0.320867E-02 4 0.382324E-01 0.437629 0.716160 0.138847 0.815756 5 0.484219 0.455396 0.574479 0.716710E-02 0.246331 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.624636 2: 0.641433 3: 0.509647 4: 0.554701 5: 0.802144 6: 0.909928 7: 0.597354 8: 0.373537 ........ .............. 4096: 0.877883 The FFT coefficients: 1: 0.630724 2: -0.192255 3: 0.134711 4: -0.946780E-01 5: 0.709119E-01 6: -0.529193E-01 7: 0.524691E-01 8: -0.534126E-01 ........ .............. 4096: -0.405538E-02 The retrieved data: 1: 0.624636 2: 0.641433 3: 0.509647 4: 0.554701 5: 0.802144 6: 0.909928 7: 0.597354 8: 0.373537 ........ .............. 4096: 0.877883 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.189133 2: 0.402065 3: 0.423615 4: 0.351888E-02 5: 0.311672 6: 0.903986 7: 0.115359 8: 0.771221 ........ .............. 4096: 0.173484 The FFT coefficients: 1: 0.489936 2: 0.453254E-03 3: -0.116206E-02 4: -0.770209E-02 5: 0.482293E-02 6: 0.650081E-02 7: -0.352294E-02 8: -0.510383E-02 ........ .............. 4096: -0.530363E-02 The retrieved data: 1: 0.189133 2: 0.402065 3: 0.423615 4: 0.351888E-02 5: 0.311672 6: 0.903986 7: 0.115359 8: 0.771221 ........ .............. 4096: 0.173484 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.609305 2: 0.646674 3: 0.171295 4: 0.551438 5: 0.228903 6: 0.603435 7: 0.870343 8: 0.107605 ........ .............. 4096: 0.714984 The FFT coefficients: 1: 0.508483 2: 0.732533E-02 3: 0.143162E-02 4: 0.407829E-02 5: -0.615000E-02 6: 0.451236E-02 7: 0.389005E-02 8: 0.822787E-02 ........ .............. 4096: -0.307990E-03 The retrieved data: 1: 0.609305 2: 0.646674 3: 0.171295 4: 0.551438 5: 0.228903 6: 0.603435 7: 0.870343 8: 0.107605 ........ .............. 4096: 0.714984 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.163409 2: 0.524809 3: 0.214802 4: 0.953394 5: 0.640465 6: 0.268951 7: 0.192580 8: 0.799575 ........ .............. 4096: 0.757844 The FFT coefficients: 1: 0.639072 2: 0.509480E-02 3: 0.211259 4: -0.486685E-02 5: 0.125297 6: -0.655049E-02 7: 0.953359E-01 8: 0.358287E-02 ........ .............. 4096: 0.498609E-02 The retrieved data: 1: 0.163409 2: 0.524809 3: 0.214802 4: 0.953394 5: 0.640465 6: 0.268951 7: 0.192580 8: 0.799575 ........ .............. 4096: 0.757844 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.520866 2: 0.348540E-01 3: 0.949735 4: 0.128079 5: 0.470331 6: 0.190532 7: 0.714678 8: 0.561631 ........ .............. 4096: 0.228626 The FFT coefficients: 1: 0.493132 2: 0.553013E-02 3: -0.156647E-02 4: -0.198830E-03 5: 0.446312E-03 6: -0.127641E-02 7: 0.616680E-02 8: -0.884440E-02 ........ .............. 4096: 0.705342E-02 The retrieved data: 1: 0.520866 2: 0.348541E-01 3: 0.949735 4: 0.128079 5: 0.470331 6: 0.190532 7: 0.714678 8: 0.561631 ........ .............. 4096: 0.228626 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.541386 0.713802 0.950973 0.489065 0.846321 2 0.258205 0.167356 0.686881 0.132059 0.990121 3 0.284867 0.455791 0.555786 0.547662 0.327457 4 0.883304 0.154033 0.205622 0.822163E-01 0.860805 5 0.587640 0.919138 0.299400 0.396361 0.822996 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.494752 0.475789E-02 0.355534E-02 -0.161373E-02 -0.512108E-03 2 0. 0.331335E-02 -0.364380E-03 0.653697E-03 0.420087E-02 3 0.916433E-02 -0.134144E-01 0.214079E-01 0.928547E-02 -0.489991E-02 4 -0.104229E-02 0.603566E-02 -0.946554E-02 0.179535E-02 0.804844E-02 5 -0.292646E-02 -0.645194E-02 -0.986523E-02 -0.483085E-02 0.442636E-02 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.541386 0.713803 0.950973 0.489065 0.846321 2 0.258205 0.167356 0.686881 0.132059 0.990121 3 0.284867 0.455791 0.555786 0.547662 0.327457 4 0.883304 0.154033 0.205622 0.822163E-01 0.860805 5 0.587640 0.919138 0.299400 0.396361 0.822996 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.836581E-02 0.551298 0.834863 0.591952 0.653038 2 0.176378 0.818127 0.340675 0.355827 0.674821 3 0.355904 0.781764 0.445068 0.612677E-01 0.718127 4 0.679852 0.873624 0.345279 0.406505 0.339380 5 0.306986 0.927871 0.612592 0.374988 0.515732E-01 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.479104 0.534577 0.495900 0.424230 0.512319 2 0.400600E-01 0.331582E-01 -0.146346E-01 0.109743E-01 -0.797300E-01 3 0.645329E-01 0.593474E-01 -0.202222E-02 0.114618E-01 -0.235311E-01 4 -0.635498E-01 0.921115E-01 -0.768738E-01 0.791486E-01 -0.624550E-01 5 -0.113089 0.726904E-01 0.367127E-01 0.819146E-01 0.207836E-02 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.836580E-02 0.551298 0.834863 0.591952 0.653038 2 0.176378 0.818127 0.340675 0.355827 0.674821 3 0.355904 0.781764 0.445068 0.612678E-01 0.718127 4 0.679852 0.873624 0.345279 0.406505 0.339380 5 0.306986 0.927871 0.612592 0.374988 0.515732E-01 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.829859 2: 0.171239 3: 0.634179 4: 0.831197 5: 0.747984 6: 0.533622 7: 0.737663 8: 0.450365 ........ .............. 4096: 0.406581 The FFT coefficients: 1: 0.635879 2: 0.213375 3: 0.125391 4: 0.939981E-01 5: 0.788726E-01 6: 0.575894E-01 7: 0.429460E-01 8: 0.405307E-01 ........ .............. 4096: 0.767714E-02 The retrieved data: 1: 0.829859 2: 0.171239 3: 0.634178 4: 0.831197 5: 0.747984 6: 0.533622 7: 0.737663 8: 0.450365 ........ .............. 4096: 0.406581 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.450938 0.468528 0.827691 0.167577 0.768906 2 0.818380 0.912278 0.660934 0.566673 0.614638 3 0.440114 0.266598 0.556957 0.298481 0.891284 4 0.761677 0.688599 0.672068 0.540091 0.963287 5 0.267355 0.539677 0.924223 0.108492 0.167028 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.691485 0.481317 0.645612 0.611812 0.704458 2 0.333544 0.204742 0.262240 0.272250 0.265073 3 0.157571 0.263306 0.211849 0.173379 0.142006 4 0.225171E-01 0.137816 0.237935 0.363612E-02 0.176057 5 0.206934E-01 -0.829097E-01 0.201848 0.365910E-02 0.110226 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.450938 0.468528 0.827691 0.167577 0.768905 2 0.818380 0.912278 0.660934 0.566673 0.614638 3 0.440114 0.266598 0.556957 0.298481 0.891284 4 0.761677 0.688599 0.672068 0.540091 0.963287 5 0.267355 0.539677 0.924223 0.108492 0.167028 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.716716 2: 0.475364 3: 0.478723 4: 0.625806 5: 0.964485 6: 0.571848 7: 0.977000 8: 0.338040 ........ .............. 4096: 0.335715 The FFT coefficients: 1: 0.636029 2: 0.122078E-01 3: 0.215155 4: 0.526709E-02 5: 0.126676 6: -0.987417E-02 7: 0.933538E-01 8: 0.512637E-02 ........ .............. 4096: -0.211065E-04 The retrieved data: 1: 0.719210 2: 0.475473 3: 0.479592 4: 0.625784 5: 0.965031 6: 0.571852 7: 0.977395 8: 0.338074 ........ .............. 4096: 0.338144 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.912344 0.274916 0.971929E-01 0.222695E-01 0.503263 2 0.827798E-01 0.240565 0.247218 0.947860 0.864879 3 0.627797 0.878860 0.181765 0.656451 0.692060 4 0.581811E-01 0.611894 0.600151E-01 0.450975 0.680722 5 0.104846 0.455361 0.669805 0.262008 0.662910 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.566524 0.745715 0.624624 0.669795 0.584084 2 -0.620284E-01 -0.393946E-01 -0.491305E-01 -0.536360E-01 0.565979E-01 3 0.918730E-01 0.221493 0.609288E-01 0.146547 0.254520 4 -0.517155E-01 0.502222E-01 -0.925566E-02 -0.239402E-01 0.625759E-01 5 0.150472 0.222008 0.438953E-01 0.244171 0.105180 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.912343 0.274914 0.971923E-01 0.222677E-01 0.503263 2 0.827799E-01 0.240565 0.247218 0.947861 0.864879 3 0.627797 0.878860 0.181765 0.656451 0.692060 4 0.581816E-01 0.611894 0.600155E-01 0.450976 0.680722 5 0.104846 0.455361 0.669805 0.262008 0.662910 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.178977 -0.496662 2: -0.416236E-01 0.151921E-01 3: 0.198244 -0.203619 4: 0.672439E-02 0.175640 5: 0.704760 0.543468 6: -0.413422E-01 -0.119290 7: -0.189293 0.397423E-01 8: -0.844410E-02 0.120502 ........ .............. .............. 4096: -0.189231E-01 0.617130 The FFT coefficients: 1: -0.103565E-02 -0.249341E-02 2: 0.116177E-02 0.103153E-01 3: -0.571267E-02 -0.681815E-02 4: -0.595155E-02 0.149888E-02 5: 0.736413E-02 -0.546891E-02 6: 0.196456E-02 -0.186275E-02 7: 0.133663E-01 -0.652920E-02 8: 0.295702E-03 0.715988E-03 ........ .............. .............. 4096: -0.120936E-02 -0.238860E-02 The retrieved data: 1: -0.178977 -0.496662 2: -0.416236E-01 0.151921E-01 3: 0.198244 -0.203619 4: 0.672439E-02 0.175640 5: 0.704760 0.543468 6: -0.413422E-01 -0.119290 7: -0.189293 0.397423E-01 8: -0.844410E-02 0.120502 ........ .............. .............. 4096: -0.189231E-01 0.617130 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.171 0.601 0.599 0.329E-01-0.750 -0.255 0.480 -0.497 2 0.542 -0.285 0.351 -0.336E-01 0.310 0.683 -0.947 -0.249 3 0.591 0.153E-01 0.790 -0.363 -0.808 0.226E-01 0.610 0.697 4 -0.919 0.110E-01 0.469 0.750 0.356 0.713 0.849 0.345 5 0.715 -0.238 -0.748 0.290 -0.430 -0.626 0.595 0.657E-01 Col: 5 Row --- 1 -0.463 0.992E-01 2 -0.741 0.349 3 -0.518 0.233E-01 4 -0.440 0.385 5 0.112 -0.931 Part of the FFT coefficients: Col: 1 2 3 4 Row --- 1 -0.165E-01 0.489E-02-0.147E-01 0.106E-01-0.192E-01-0.159E-01-0.871E-02 0.245E-02 2 0.187E-01 0.456E-02-0.432E-02 0.187E-01 0.364E-02-0.884E-02 0.148E-01 0.635E-02 3 0.717E-02 0.420E-02 0.147E-01-0.146E-01-0.355E-01 0.772E-03-0.203E-02 0.231E-02 4 0.381E-02-0.247E-03 0.851E-02 0.647E-02-0.138E-01 0.252E-02 0.163E-01-0.162E-01 5 -0.517E-03 0.798E-03-0.150E-01 0.724E-02 0.175E-01 0.117E-01 0.120E-01 0.836E-03 Col: 5 Row --- 1 0.875E-02 0.199E-01 2 -0.175E-01-0.902E-02 3 -0.654E-02-0.118E-01 4 -0.522E-02-0.210E-01 5 0.991E-02 0.589E-02 Part of the retrieved data: Col: 1 2 3 4 Row --- 1 0.171 0.601 0.599 0.329E-01-0.750 -0.255 0.480 -0.497 2 0.542 -0.285 0.351 -0.336E-01 0.310 0.683 -0.947 -0.249 3 0.591 0.153E-01 0.790 -0.363 -0.808 0.226E-01 0.610 0.697 4 -0.919 0.110E-01 0.469 0.750 0.356 0.713 0.849 0.345 5 0.715 -0.238 -0.748 0.290 -0.430 -0.626 0.595 0.657E-01 Col: 5 Row --- 1 -0.463 0.992E-01 2 -0.741 0.349 3 -0.518 0.233E-01 4 -0.440 0.385 5 0.112 -0.931 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.624 -0.856E-01 0.694 0.313 0.346 0.814 0.416E-01 0.481 2 -0.417 -0.354 0.107 -0.242 0.859 0.291 -0.922 0.153 3 0.452E-01-0.595 0.865 -0.935E-01 0.954 -0.231 0.328 0.289 4 0.268E-01 0.147 0.768 -0.176 0.297E-01 0.488 0.538 -0.521 5 0.271E-02 0.507 -0.460 0.747 -0.274 -0.298 -0.173 -0.571E-01 Col: 5 Row --- 1 0.139E-01-0.388 2 -0.418 0.594 3 0.905 -0.225 4 -0.318 0.966E-01 5 -0.435 0.334 Part of the FFT coefficients: Col: 1 2 3 4 Row --- 1 0.178 -0.174 0.517E-01-0.926E-01-0.779E-01 0.104 0.592E-01 0.132E-01 2 0.106 0.747E-01 0.823E-01 0.101E-01 0.673E-01-0.103 -0.404E-01-0.833E-01 3 -0.879E-02-0.616E-01 0.225 0.928E-03-0.925E-02-0.571E-01-0.192 -0.436E-01 4 -0.784E-02 0.739E-01-0.300E-01-0.525E-01 0.112 -0.240 0.161E-01-0.316E-02 5 -0.131 -0.559E-01-0.529E-01 0.733E-01-0.179E-01 0.208 0.362E-02 0.118E-01 Col: 5 Row --- 1 -0.518E-01 0.103 2 0.129 -0.107 3 0.513E-01 0.861E-01 4 -0.134E-01 0.540E-01 5 -0.123E-01-0.342E-01 Part of the retrieved data: Col: 1 2 3 4 Row --- 1 0.624 -0.856E-01 0.694 0.313 0.346 0.814 0.416E-01 0.481 2 -0.417 -0.354 0.107 -0.242 0.859 0.291 -0.922 0.153 3 0.452E-01-0.595 0.865 -0.935E-01 0.954 -0.231 0.328 0.289 4 0.268E-01 0.147 0.768 -0.176 0.297E-01 0.488 0.538 -0.521 5 0.271E-02 0.507 -0.460 0.747 -0.274 -0.298 -0.173 -0.571E-01 Col: 5 Row --- 1 0.139E-01-0.388 2 -0.418 0.594 3 0.905 -0.225 4 -0.318 0.966E-01 5 -0.435 0.334 FFTPACK5_TEST Normal end of execution. 6 October 2025 6:18:18.085 PM