25 February 2023 6:28:44.698 AM 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.122847 -0.187108E-01 2: 0.779209 -0.449593 3: 0.386667 0.157388E-01 4: -0.732450 0.531549 5: 0.297286 0.404264 6: 0.289634 -0.752138 7: 0.171588 0.682045 8: 0.564249 0.401782 ........ .............. .............. 4096: 0.278509E-01 0.326458 The FFT coefficients: 1: -0.779193E-03 -0.169766E-02 2: 0.269171E-02 0.171117E-01 3: -0.985017E-03 -0.557189E-02 4: 0.783427E-02 0.472466E-02 5: -0.128650E-01 0.158053E-02 6: 0.131372E-02 0.200234E-01 7: 0.124140E-01 0.197644E-02 8: 0.107536E-02 -0.101186E-02 ........ .............. .............. 4096: -0.611606E-02 0.546467E-02 The retrieved data: 1: -0.122847 -0.187108E-01 2: 0.779209 -0.449593 3: 0.386667 0.157388E-01 4: -0.732450 0.531549 5: 0.297286 0.404264 6: 0.289634 -0.752138 7: 0.171588 0.682045 8: 0.564249 0.401782 ........ .............. .............. 4096: 0.278508E-01 0.326457 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.123 -0.187E-01-0.599 -0.535 -0.952 -0.329E-01 0.724 0.229 2 0.779 -0.450 0.319 0.322E-01-0.180 -0.134E-01 0.528 0.310E-01 3 0.387 0.157E-01-0.968 -0.421E-01-0.223 0.749 -0.172 0.136 4 -0.732 0.532 0.807 -0.495 0.458 -0.518E-01 0.620 0.644 5 0.297 0.404 0.802 -0.196 0.864 0.493 -0.574 0.172 Col: 5 Row --- 1 -0.568E-02-0.749 2 -0.400E-01 0.990 3 0.810 0.212 4 -0.847E-01 0.733 5 0.378 0.489 Part of the FFT coefficients: Col: 1 2 3 4 Row --- 1 -0.117E-01 0.704E-02 0.361E-02 0.952E-04-0.193E-01-0.441E-02 0.280E-01 0.525E-02 2 0.155E-01 0.140E-02-0.596E-02 0.615E-02 0.320E-02-0.133E-02 0.990E-02 0.117E-01 3 0.926E-02 0.659E-02 0.139E-01 0.586E-02 0.120E-02-0.142E-02-0.849E-02-0.104E-02 4 0.978E-02 0.587E-02-0.220E-02 0.114E-03-0.399E-03 0.994E-02-0.288E-02 0.786E-02 5 -0.116E-01 0.816E-02-0.604E-02-0.210E-01-0.244E-01 0.487E-02 0.996E-02-0.532E-02 Col: 5 Row --- 1 -0.581E-02-0.437E-02 2 -0.379E-02-0.249E-01 3 0.354E-02-0.919E-02 4 -0.268E-01 0.921E-02 5 -0.143E-02 0.837E-02 Part of the retrieved data: Col: 1 2 3 4 Row --- 1 -0.123 -0.187E-01-0.599 -0.535 -0.952 -0.329E-01 0.724 0.229 2 0.779 -0.450 0.319 0.322E-01-0.180 -0.134E-01 0.528 0.310E-01 3 0.387 0.157E-01-0.968 -0.421E-01-0.223 0.749 -0.172 0.136 4 -0.732 0.532 0.807 -0.495 0.458 -0.518E-01 0.620 0.644 5 0.297 0.404 0.802 -0.196 0.864 0.493 -0.574 0.172 Col: 5 Row --- 1 -0.568E-02-0.749 2 -0.400E-01 0.990 3 0.810 0.212 4 -0.847E-01 0.733 5 0.378 0.489 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.123 -0.187E-01-0.599 -0.535 -0.952 -0.329E-01 0.724 0.229 2 0.779 -0.450 0.319 0.322E-01-0.180 -0.134E-01 0.528 0.310E-01 3 0.387 0.157E-01-0.968 -0.421E-01-0.223 0.749 -0.172 0.136 4 -0.732 0.532 0.807 -0.495 0.458 -0.518E-01 0.620 0.644 5 0.297 0.404 0.802 -0.196 0.864 0.493 -0.574 0.172 Col: 5 Row --- 1 -0.568E-02-0.749 2 -0.400E-01 0.990 3 0.810 0.212 4 -0.847E-01 0.733 5 0.378 0.489 Part of the FFT coefficients: Col: 1 2 3 4 Row --- 1 -0.133 0.349E-01-0.103 0.582E-01-0.567E-01 0.170 0.107 -0.209E-01 2 0.273E-01-0.109 -0.513E-01-0.119 -0.712E-01 0.209E-01 0.751E-01 0.101 3 -0.785E-01-0.728E-01-0.671E-01-0.115 0.525E-01-0.768E-01 0.726E-01 0.318E-02 4 0.127E-02 0.821E-02 0.800E-02 0.736E-01-0.845E-01-0.264E-01 0.978E-01 0.114 5 -0.601E-01-0.134 -0.145 -0.203E-01-0.552E-01-0.110 -0.839E-01 0.765E-01 Col: 5 Row --- 1 0.244E-01 0.171 2 0.224 0.133 3 0.176 0.137 4 -0.146 -0.957E-01 5 -0.420E-01-0.452E-01 Part of the retrieved data: Col: 1 2 3 4 Row --- 1 -0.123 -0.187E-01-0.599 -0.535 -0.952 -0.329E-01 0.724 0.229 2 0.779 -0.450 0.319 0.322E-01-0.180 -0.134E-01 0.528 0.310E-01 3 0.387 0.157E-01-0.968 -0.421E-01-0.223 0.749 -0.172 0.136 4 -0.732 0.532 0.807 -0.495 0.458 -0.518E-01 0.620 0.644 5 0.297 0.404 0.802 -0.196 0.864 0.493 -0.574 0.172 Col: 5 Row --- 1 -0.568E-02-0.749 2 -0.400E-01 0.990 3 0.810 0.212 4 -0.847E-01 0.733 5 0.378 0.489 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.539534 2: 0.416468 3: 0.193611 4: 0.817742 5: 0.719557 6: 0.497288 7: 0.869629 8: 0.704313 ........ .............. 4096: 0.157169 The FFT coefficients: 1: 0.632680 2: -0.204740 3: 0.134645 4: -0.875416E-01 5: 0.583805E-01 6: -0.628708E-01 7: 0.584447E-01 8: -0.341073E-01 ........ .............. 4096: 0.224411E-03 The retrieved data: 1: 0.539534 2: 0.416468 3: 0.193611 4: 0.817742 5: 0.719557 6: 0.497288 7: 0.869629 8: 0.704313 ........ .............. 4096: 0.157169 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.480174 0.845840 0.845505 0.795895 0.298671 2 0.883769 0.870402 0.634498 0.139651E-01 0.612482 3 0.981831 0.424082 0.190107 0.711233 0.230011 4 0.866995 0.984562 0.539882 0.762234 0.109115 5 0.274882 0.281187 0.512555 0.519369 0.926991 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.780645 0.672595 0.720443 0.627566 0.515074 2 -0.154687 -0.184961 -0.304978 -0.244869 -0.179657 3 0.194195 0.104132 0.860862E-01 0.220580 0.788737E-01 4 -0.193197E-01 -0.995412E-03 -0.473579E-01 -0.419295E-01 -0.630928E-01 5 -0.289852E-01 0.126964 0.141807 0.386468E-01 0.874671E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.480174 0.845840 0.845505 0.795895 0.298671 2 0.883769 0.870402 0.634498 0.139651E-01 0.612482 3 0.981831 0.424082 0.190107 0.711233 0.230011 4 0.866995 0.984562 0.539882 0.762234 0.109115 5 0.274882 0.281187 0.512555 0.519369 0.926991 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.299255 2: 0.338710 3: 0.185129 4: 0.689363 5: 0.867619 6: 0.782228 7: 0.675736 8: 0.134861 ........ .............. 4096: 0.127278 The FFT coefficients: 1: 0.504828 2: -0.413500E-02 3: -0.705105E-03 4: 0.729934E-02 5: -0.735645E-02 6: -0.701813E-02 7: -0.176098E-02 8: -0.472289E-02 ........ .............. 4096: -0.232263E-02 The retrieved data: 1: 0.299216 2: 0.338769 3: 0.185142 4: 0.689370 5: 0.867658 6: 0.782210 7: 0.675743 8: 0.134854 ........ .............. 4096: 0.127310 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.641528 0.662216 0.851953 0.250036 0.535345 2 0.610206 0.304597 0.797175 0.127638E-01 0.692552E-01 3 0.679809 0.910579E-01 0.505227 0.543552 0.410298 4 0.168023 0.305452 0.482583E-01 0.331507 0.336956 5 0.905025 0.494696 0.124296 0.544585 0.831872 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.459098 0.516222 0.367013 0.461798 0.479833 2 0.572270E-01 -0.141414 -0.664907E-01 -0.563128E-01 -0.477650E-01 3 -0.518515E-01 -0.350065E-01 0.127038 -0.938872E-01 0.582555E-01 4 0.423522E-01 0.168424 -0.199881E-02 0.132102 0.278599E-01 5 0.883156E-01 -0.161942 0.573737E-01 -0.656550E-01 -0.103712 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.641528 0.662214 0.851953 0.250035 0.535345 2 0.610206 0.304597 0.797174 0.127640E-01 0.692549E-01 3 0.679809 0.910584E-01 0.505227 0.543552 0.410298 4 0.168024 0.305452 0.482590E-01 0.331508 0.336956 5 0.905024 0.494696 0.124296 0.544584 0.831872 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.783004 2: 0.429467 3: 0.555037 4: 0.202918E-01 5: 0.383178 6: 0.195411E-01 7: 0.257777E-01 8: 0.115422 ........ .............. 4096: 0.640982 The FFT coefficients: 1: 0.635444 2: -0.220632 3: 0.116371 4: -0.954503E-01 5: 0.737503E-01 6: -0.640337E-01 7: 0.390152E-01 8: -0.411695E-01 ........ .............. 4096: -0.126226E-01 The retrieved data: 1: 0.783004 2: 0.429467 3: 0.555037 4: 0.202918E-01 5: 0.383178 6: 0.195411E-01 7: 0.257777E-01 8: 0.115422 ........ .............. 4096: 0.640982 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.274335 2: 0.255898E-01 3: 0.590385 4: 0.213442 5: 0.154537 6: 0.118647 7: 0.371873 8: 0.444302 ........ .............. 4096: 0.974435 The FFT coefficients: 1: 0.500507 2: -0.301298E-02 3: -0.373750E-02 4: -0.103363E-01 5: -0.478553E-02 6: -0.374341E-02 7: -0.408529E-02 8: -0.805673E-02 ........ .............. 4096: -0.105803E-03 The retrieved data: 1: 0.274335 2: 0.255898E-01 3: 0.590385 4: 0.213442 5: 0.154537 6: 0.118647 7: 0.371873 8: 0.444302 ........ .............. 4096: 0.974435 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.203852 2: 0.440061 3: 0.111566 4: 0.216969 5: 0.371172 6: 0.395827 7: 0.468025 8: 0.905393 ........ .............. 4096: 0.587722 The FFT coefficients: 1: 0.502383 2: -0.168977E-02 3: 0.179176E-02 4: -0.961130E-02 5: 0.534654E-02 6: -0.783438E-02 7: -0.362094E-03 8: -0.216722E-02 ........ .............. 4096: 0.399631E-02 The retrieved data: 1: 0.203852 2: 0.440061 3: 0.111566 4: 0.216969 5: 0.371172 6: 0.395827 7: 0.468025 8: 0.905393 ........ .............. 4096: 0.587722 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.948544 2: 0.228583 3: 0.380337E-01 4: 0.942235 5: 0.990162 6: 0.113311 7: 0.276048 8: 0.807852 ........ .............. 4096: 0.491501 The FFT coefficients: 1: 0.635523 2: -0.375535E-03 3: 0.222153 4: -0.863166E-02 5: 0.136237 6: -0.121956E-01 7: 0.807226E-01 8: 0.109851E-01 ........ .............. 4096: -0.129016E-02 The retrieved data: 1: 0.948544 2: 0.228583 3: 0.380337E-01 4: 0.942235 5: 0.990162 6: 0.113311 7: 0.276048 8: 0.807852 ........ .............. 4096: 0.491501 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.833493 2: 0.321451 3: 0.734789 4: 0.362489 5: 0.938502E-01 6: 0.125501E-01 7: 0.693220 8: 0.647398 ........ .............. 4096: 0.890332 The FFT coefficients: 1: 0.505929 2: 0.424052E-02 3: -0.583533E-02 4: -0.170163E-02 5: -0.118607E-02 6: -0.778442E-02 7: 0.486312E-02 8: -0.329760E-02 ........ .............. 4096: -0.907481E-03 The retrieved data: 1: 0.833493 2: 0.321451 3: 0.734789 4: 0.362489 5: 0.938501E-01 6: 0.125499E-01 7: 0.693220 8: 0.647398 ........ .............. 4096: 0.890332 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.796579 0.168208 0.380958 0.717646 0.470663 2 0.799316 0.832680 0.676921 0.343789 0.342172E-01 3 0.201249 0.873256 0.697109 0.402565E-01 0.415608 4 0.584365 0.789817 0.337861 0.346513 0.111743 5 0.218751 0.543482 0.161350 0.480548 0.742245 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.486962 0.706137E-02 0.926804E-02 0.328156E-02 0.156774E-02 2 0. -0.372434E-02 -0.243291E-02 0.541377E-02 0.441244E-02 3 0.682429E-02 -0.153819E-01 0.179983E-02 0.871157E-02 -0.240382E-02 4 0.221088E-02 -0.156267E-01 0.156163E-01 -0.101578E-01 0.251903E-02 5 0.431919E-03 -0.105081E-01 0.863711E-02 0.640720E-02 0.565592E-03 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.796579 0.168208 0.380958 0.717646 0.470663 2 0.799316 0.832680 0.676921 0.343789 0.342173E-01 3 0.201249 0.873256 0.697109 0.402564E-01 0.415608 4 0.584365 0.789817 0.337861 0.346513 0.111743 5 0.218751 0.543482 0.161350 0.480548 0.742245 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.833369 0.776870 0.513855 0.578250 0.926766 2 0.500557 0.283205 0.659361 0.782724 0.561755 3 0.695131 0.679943 0.928252E-01 0.395978 0.663067 4 0.610411 0.993588 0.378326 0.912780E-01 0.928057E-01 5 0.110912 0.212234 0.220060 0.396343 0.477923 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.453140 0.497418 0.407217 0.507289 0.512047 2 -0.158100 0.596114E-01 -0.155569 -0.369332E-01 0.252753E-01 3 -0.127501E-01 0.409192E-01 0.332282E-01 0.153468 -0.488889E-01 4 0.134184E-01 0.725468E-02 0.543634E-01 -0.392605E-01 0.436864E-01 5 0.745535E-01 -0.209704E-01 -0.984666E-02 -0.929179E-01 -0.565715E-02 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.833369 0.776870 0.513855 0.578250 0.926766 2 0.500557 0.283205 0.659361 0.782724 0.561755 3 0.695131 0.679943 0.928252E-01 0.395978 0.663067 4 0.610411 0.993588 0.378326 0.912780E-01 0.928057E-01 5 0.110912 0.212234 0.220060 0.396343 0.477923 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.912553 2: 0.243193 3: 0.143173 4: 0.226523 5: 0.525968 6: 0.457689 7: 0.711849 8: 0.383520 ........ .............. 4096: 0.281544 The FFT coefficients: 1: 0.635503 2: 0.207360 3: 0.137049 4: 0.886762E-01 5: 0.595402E-01 6: 0.522006E-01 7: 0.563963E-01 8: 0.432355E-01 ........ .............. 4096: 0.156242E-02 The retrieved data: 1: 0.912553 2: 0.243193 3: 0.143173 4: 0.226523 5: 0.525968 6: 0.457689 7: 0.711849 8: 0.383520 ........ .............. 4096: 0.281544 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.576052 0.954240 0.350392 0.262384 0.968487 2 0.129080E-01 0.261012 0.226828 0.331330 0.476049 3 0.628773 0.438460 0.448746 0.380222 0.314307 4 0.920156E-01 0.633983 0.583190 0.945970E-01 0.405892 5 0.869633 0.686490 0.927990 0.367519 0.682225E-01 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.556490 0.705520 0.646630 0.583854 0.694991 2 0.143360 0.225367 0.124201 0.206582 -0.713499E-01 3 0.117603 0.235490 0.193251 0.237787E-01 0.520636E-01 4 0.185424 0.131638 0.196058 -0.463332E-01 0.138453 5 0.266792E-01 0.230225E-01 0.726005E-01 0.631765E-01 0.582444E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.576052 0.954240 0.350392 0.262384 0.968487 2 0.129080E-01 0.261012 0.226828 0.331330 0.476049 3 0.628773 0.438460 0.448746 0.380222 0.314307 4 0.920155E-01 0.633983 0.583190 0.945970E-01 0.405892 5 0.869633 0.686489 0.927990 0.367519 0.682224E-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.481496 2: 0.396902 3: 0.737643E-01 4: 0.999704 5: 0.227493 6: 0.626439 7: 0.574397 8: 0.825406 ........ .............. 4096: 0.905353 The FFT coefficients: 1: 0.630500 2: -0.637816E-02 3: 0.217164 4: -0.193237E-02 5: 0.125503 6: 0.311987E-02 7: 0.830325E-01 8: -0.124317E-01 ........ .............. 4096: -0.651245E-02 The retrieved data: 1: 0.483932 2: 0.396899 3: 0.745303E-01 4: 0.999483 5: 0.228143 6: 0.626389 7: 0.574711 8: 0.825405 ........ .............. 4096: 0.907535 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.923100 0.379382 0.162137E-02 0.519200 0.323778 2 0.474944 0.839333 0.657888 0.544317 0.503467 3 0.515133E-01 0.112484 0.422173E-01 0.537279 0.659370E-02 4 0.691902 0.138806 0.582501 0.631411 0.938553E-02 5 0.772650 0.829866 0.830888 0.650464 0.251231 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1 0.629365 0.600179 0.498184 0.757081 0.689841 2 -0.585092E-01 -0.746914E-01 -0.285941E-01 -0.334700E-01 -0.315580E-01 3 0.315045 0.650125E-01 0.272100 0.303344 0.138853 4 -0.513820E-01 -0.365363E-01 -0.611473E-02 0.308574E-01 -0.144775 5 0.192458 0.205793 0.149558 0.201013 0.974980E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.923099 0.379378 0.162048E-02 0.519200 0.323776 2 0.474945 0.839334 0.657888 0.544316 0.503468 3 0.515133E-01 0.112484 0.422172E-01 0.537279 0.659356E-02 4 0.691902 0.138807 0.582501 0.631411 0.938597E-02 5 0.772650 0.829865 0.830888 0.650464 0.251231 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.122847 -0.187108E-01 2: 0.779210 -0.449592 3: 0.386667 0.157388E-01 4: -0.732450 0.531549 5: 0.297286 0.404264 6: 0.289634 -0.752138 7: 0.171588 0.682045 8: 0.564249 0.401782 ........ .............. .............. 4096: 0.278509E-01 0.326458 The FFT coefficients: 1: -0.779161E-03 -0.169767E-02 2: 0.269171E-02 0.171117E-01 3: -0.985015E-03 -0.557189E-02 4: 0.783426E-02 0.472466E-02 5: -0.128650E-01 0.158052E-02 6: 0.131372E-02 0.200234E-01 7: 0.124140E-01 0.197644E-02 8: 0.107536E-02 -0.101186E-02 ........ .............. .............. 4096: -0.611606E-02 0.546467E-02 The retrieved data: 1: -0.122847 -0.187108E-01 2: 0.779210 -0.449592 3: 0.386667 0.157388E-01 4: -0.732450 0.531549 5: 0.297286 0.404264 6: 0.289634 -0.752138 7: 0.171588 0.682045 8: 0.564249 0.401782 ........ .............. .............. 4096: 0.278509E-01 0.326458 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.123 -0.187E-01-0.599 -0.535 -0.952 -0.329E-01 0.724 0.229 2 0.779 -0.450 0.319 0.322E-01-0.180 -0.134E-01 0.528 0.310E-01 3 0.387 0.157E-01-0.968 -0.421E-01-0.223 0.749 -0.172 0.136 4 -0.732 0.532 0.807 -0.495 0.458 -0.518E-01 0.620 0.644 5 0.297 0.404 0.802 -0.196 0.864 0.493 -0.574 0.172 Col: 5 Row --- 1 -0.568E-02-0.749 2 -0.400E-01 0.990 3 0.810 0.212 4 -0.847E-01 0.733 5 0.378 0.489 Part of the FFT coefficients: Col: 1 2 3 4 Row --- 1 -0.117E-01 0.704E-02 0.361E-02 0.952E-04-0.193E-01-0.441E-02 0.280E-01 0.525E-02 2 0.155E-01 0.140E-02-0.596E-02 0.615E-02 0.320E-02-0.133E-02 0.990E-02 0.117E-01 3 0.926E-02 0.659E-02 0.139E-01 0.586E-02 0.120E-02-0.142E-02-0.849E-02-0.104E-02 4 0.978E-02 0.587E-02-0.220E-02 0.114E-03-0.399E-03 0.994E-02-0.288E-02 0.786E-02 5 -0.116E-01 0.816E-02-0.604E-02-0.210E-01-0.244E-01 0.487E-02 0.996E-02-0.532E-02 Col: 5 Row --- 1 -0.581E-02-0.437E-02 2 -0.379E-02-0.249E-01 3 0.354E-02-0.919E-02 4 -0.268E-01 0.921E-02 5 -0.143E-02 0.837E-02 Part of the retrieved data: Col: 1 2 3 4 Row --- 1 -0.123 -0.187E-01-0.599 -0.535 -0.952 -0.329E-01 0.724 0.229 2 0.779 -0.450 0.319 0.322E-01-0.180 -0.134E-01 0.528 0.310E-01 3 0.387 0.157E-01-0.968 -0.421E-01-0.223 0.749 -0.172 0.136 4 -0.732 0.532 0.807 -0.495 0.458 -0.518E-01 0.620 0.644 5 0.297 0.404 0.802 -0.196 0.864 0.493 -0.574 0.172 Col: 5 Row --- 1 -0.568E-02-0.749 2 -0.400E-01 0.990 3 0.810 0.212 4 -0.847E-01 0.733 5 0.378 0.489 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.123 -0.187E-01-0.599 -0.535 -0.952 -0.329E-01 0.724 0.229 2 0.779 -0.450 0.319 0.322E-01-0.180 -0.134E-01 0.528 0.310E-01 3 0.387 0.157E-01-0.968 -0.421E-01-0.223 0.749 -0.172 0.136 4 -0.732 0.532 0.807 -0.495 0.458 -0.518E-01 0.620 0.644 5 0.297 0.404 0.802 -0.196 0.864 0.493 -0.574 0.172 Col: 5 Row --- 1 -0.568E-02-0.749 2 -0.400E-01 0.990 3 0.810 0.212 4 -0.847E-01 0.733 5 0.378 0.489 Part of the FFT coefficients: Col: 1 2 3 4 Row --- 1 -0.133 0.349E-01-0.103 0.582E-01-0.567E-01 0.170 0.107 -0.209E-01 2 0.273E-01-0.109 -0.513E-01-0.119 -0.712E-01 0.209E-01 0.751E-01 0.101 3 -0.785E-01-0.728E-01-0.671E-01-0.115 0.525E-01-0.768E-01 0.726E-01 0.318E-02 4 0.127E-02 0.821E-02 0.800E-02 0.736E-01-0.845E-01-0.264E-01 0.978E-01 0.114 5 -0.601E-01-0.134 -0.145 -0.203E-01-0.552E-01-0.110 -0.839E-01 0.765E-01 Col: 5 Row --- 1 0.244E-01 0.171 2 0.224 0.133 3 0.176 0.137 4 -0.146 -0.957E-01 5 -0.420E-01-0.452E-01 Part of the retrieved data: Col: 1 2 3 4 Row --- 1 -0.123 -0.187E-01-0.599 -0.535 -0.952 -0.329E-01 0.724 0.229 2 0.779 -0.450 0.319 0.322E-01-0.180 -0.134E-01 0.528 0.310E-01 3 0.387 0.157E-01-0.968 -0.421E-01-0.223 0.749 -0.172 0.136 4 -0.732 0.532 0.807 -0.495 0.458 -0.518E-01 0.620 0.644 5 0.297 0.404 0.802 -0.196 0.864 0.493 -0.574 0.172 Col: 5 Row --- 1 -0.568E-02-0.749 2 -0.400E-01 0.990 3 0.810 0.212 4 -0.847E-01 0.733 5 0.378 0.489 FFTPACK5_TEST Normal end of execution. 25 February 2023 6:28:44.718 AM