13 May 2013 10:20:36.701 AM FFTPACK5.1D_PRB FORTRAN90 version Test the FFTPACK5.1D library. TEST01 For complex double 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.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 TEST02 For complex double 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 TEST03 For complex double 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 TEST04 For real double 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.154414E-01 2: 0.524056 3: 0.809301 4: 0.916710 5: 0.149759 6: 0.647463E-02 7: 0.819028 8: 0.400086 ........ .............. 4096: 0.625828 The FFT coefficients: 1: 0.628319 2: -0.214452 3: 0.121189 4: -0.865563E-01 5: 0.692414E-01 6: -0.677449E-01 7: 0.460768E-01 8: -0.560361E-01 ........ .............. 4096: 0.290026E-02 The retrieved data: 1: 0.154414E-01 2: 0.524056 3: 0.809301 4: 0.916710 5: 0.149759 6: 0.647463E-02 7: 0.819028 8: 0.400086 ........ .............. 4096: 0.625828 TEST05 For real double 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.154414E-01 0.819756 0.644887 0.353603E-01 0.907628 2: 0.524056 0.647486 0.615999 0.300492 0.505506 3: 0.809301 0.289558 0.103113 0.362978 0.325764E-01 4: 0.916710 0.606128 0.160016E-01 0.573577 0.511787 5: 0.149759 0.195894 0.938865 0.100744 0.611073 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1: 0.554079 0.581825 0.632356 0.560119 0.566316 2: -0.157676 -0.178926 -0.772595E-01 -0.237567 -0.199914 3: 0.161507 0.222528 0.430641E-01 0.800814E-01 0.201477 4: -0.871489E-01 -0.206338 -0.195369 -0.130487 -0.850273E-01 5: 0.607968E-01 0.512716E-01 0.165265E-01 0.317574E-01 -0.438617E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1: 0.154414E-01 0.819756 0.644887 0.353603E-01 0.907628 2: 0.524056 0.647486 0.615999 0.300492 0.505506 3: 0.809301 0.289558 0.103113 0.362978 0.325764E-01 4: 0.916710 0.606128 0.160016E-01 0.573577 0.511787 5: 0.149759 0.195894 0.938865 0.100744 0.611073 TEST06 For real double 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.154414E-01 2: 0.524056 3: 0.809301 4: 0.916710 5: 0.149759 6: 0.647463E-02 7: 0.819028 8: 0.400086 ........ .............. 4096: 0.625828 The FFT coefficients: 1: 0.496832 2: -0.821193E-02 3: -0.151919E-03 4: -0.679062E-02 5: 0.103586E-01 6: -0.145759E-01 7: 0.104165E-02 8: -0.147735E-01 ........ .............. 4096: 0.225528E-02 The retrieved data: 1: 0.154414E-01 2: 0.524056 3: 0.809301 4: 0.916710 5: 0.149759 6: 0.647463E-02 7: 0.819028 8: 0.400086 ........ .............. 4096: 0.625828 TEST07 For real double 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.154414E-01 0.819756 0.644887 0.353603E-01 0.907628 2: 0.524056 0.647486 0.615999 0.300492 0.505506 3: 0.809301 0.289558 0.103113 0.362978 0.325764E-01 4: 0.916710 0.606128 0.160016E-01 0.573577 0.511787 5: 0.149759 0.195894 0.938865 0.100744 0.611073 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1: 0.433520 0.476392 0.460543 0.446714 0.458791 2: 0.111741E-01 -0.403864E-01 0.127365 -0.321794E-01 -0.366949E-01 3: 0.493879E-01 0.121119 0.123673E-01 -0.477149E-01 0.645467E-01 4: 0.216718E-01 -0.444270E-01 -0.124038 -0.345698E-01 0.476763E-01 5: -0.202301E-01 -0.899187E-01 -0.775478E-01 -0.566893E-01 -0.768211E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1: 0.154414E-01 0.819756 0.644887 0.353603E-01 0.907628 2: 0.524056 0.647486 0.615999 0.300492 0.505506 3: 0.809301 0.289558 0.103113 0.362978 0.325764E-01 4: 0.916710 0.606128 0.160016E-01 0.573577 0.511787 5: 0.149759 0.195894 0.938865 0.100744 0.611073 TEST09 For real double precision fast cosine 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.154414E-01 2: 0.524056 3: 0.809301 4: 0.916710 5: 0.149759 6: 0.647463E-02 7: 0.819028 8: 0.400086 ........ .............. 4096: 0.625828 The FFT coefficients: 1: 0.496789 2: -0.236545E-03 3: -0.500490E-03 4: 0.102749E-01 5: 0.374679E-02 6: 0.938335E-03 7: -0.306165E-02 8: -0.308676E-02 ........ .............. 4096: 0.218022E-02 The retrieved data: 1: 0.154414E-01 2: 0.524056 3: 0.809301 4: 0.916710 5: 0.149759 6: 0.647463E-02 7: 0.819028 8: 0.400086 ........ .............. 4096: 0.625828 TEST10 For real double 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 = 253 LENWRK = 4352 Part of the original data: Col 1 2 3 4 5 Row 1: 0.154414E-01 0.289558 0.938865 0.670009 0.989057 2: 0.524056 0.606128 0.506925 0.840560 0.826155E-01 3: 0.809301 0.195894 0.895098 0.289515 0.519208 4: 0.916710 0.397889 0.912458 0.874942 0.320776 5: 0.149759 0.316762 0.682392 0.146499 0.284288 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1: 0.493520 0.179072E-02 -0.924538E-03 -0.706917E-02 0.234106E-02 2: 0.695952E-03 0.843925E-02 0.958180E-03 -0.402839E-02 0.371127E-02 3: 0.215250E-02 -0.119762E-02 -0.371299E-02 0.964164E-02 -0.145895E-02 4: -0.456386E-02 -0.421337E-02 -0.522234E-03 -0.366154E-02 -0.334043E-03 5: -0.605849E-02 0.152034E-02 -0.697590E-02 -0.493230E-02 -0.269114E-02 Part of the retrieved data: Col 1 2 3 4 5 Row 1: 0.154414E-01 0.289558 0.938865 0.670009 0.989057 2: 0.524056 0.606128 0.506925 0.840560 0.826155E-01 3: 0.809301 0.195894 0.895098 0.289515 0.519208 4: 0.916710 0.397889 0.912458 0.874942 0.320776 5: 0.149759 0.316762 0.682392 0.146499 0.284288 TEST11 For real double 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.154414E-01 0.819756 0.644887 0.353603E-01 0.907628 2: 0.524056 0.647486 0.615999 0.300492 0.505506 3: 0.809301 0.289558 0.103113 0.362978 0.325764E-01 4: 0.916710 0.606128 0.160016E-01 0.573577 0.511787 5: 0.149759 0.195894 0.938865 0.100744 0.611073 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1: 0.432976 0.481375 0.458128 0.438682 0.465870 2: 0.449657E-01 0.121036 0.213464E-01 -0.639088E-01 0.698868E-01 3: -0.705120E-02 -0.386068E-01 0.162617 -0.637046E-02 -0.537158E-01 4: -0.227500E-01 -0.975672E-01 -0.898767E-01 -0.751189E-01 -0.466340E-01 5: 0.102067E-01 -0.736704E-01 -0.379095E-01 -0.266392E-01 0.746719E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1: 0.154414E-01 0.819756 0.644887 0.353603E-01 0.907628 2: 0.524056 0.647486 0.615999 0.300492 0.505506 3: 0.809301 0.289558 0.103113 0.362978 0.325764E-01 4: 0.916710 0.606128 0.160016E-01 0.573577 0.511787 5: 0.149759 0.195894 0.938865 0.100744 0.611073 TEST12 For real double 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.154414E-01 2: 0.524056 3: 0.809301 4: 0.916710 5: 0.149759 6: 0.647463E-02 7: 0.819028 8: 0.400086 ........ .............. 4096: 0.625828 The FFT coefficients: 1: 0.636445 2: 0.205950 3: 0.128162 4: 0.816693E-01 5: 0.849940E-01 6: 0.512879E-01 7: 0.564641E-01 8: 0.341277E-01 ........ .............. 4096: 0.170207E-02 The retrieved data: 1: 0.154414E-01 2: 0.524056 3: 0.809301 4: 0.916710 5: 0.149759 6: 0.647463E-02 7: 0.819028 8: 0.400086 ........ .............. 4096: 0.625828 TEST13 For real double 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.154414E-01 0.819756 0.644887 0.353603E-01 0.907628 2: 0.524056 0.647486 0.615999 0.300492 0.505506 3: 0.809301 0.289558 0.103113 0.362978 0.325764E-01 4: 0.916710 0.606128 0.160016E-01 0.573577 0.511787 5: 0.149759 0.195894 0.938865 0.100744 0.611073 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1: 0.542828 0.613997 0.541278 0.587242 0.593326 2: 0.166503 0.118491 0.311751 0.194124 0.134563 3: 0.141454 0.225297 0.233000 0.909497E-01 0.146886 4: 0.112766 0.948606E-01 0.480827E-01 0.604473E-01 0.170748 5: 0.399389E-01 -0.478320E-01 -0.776566E-03 -0.136147E-02 0.497462E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1: 0.154414E-01 0.819756 0.644887 0.353603E-01 0.907628 2: 0.524056 0.647486 0.615999 0.300492 0.505506 3: 0.809301 0.289558 0.103113 0.362978 0.325764E-01 4: 0.916710 0.606128 0.160016E-01 0.573577 0.511787 5: 0.149759 0.195894 0.938865 0.100744 0.611073 TEST14 For real double 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.154414E-01 2: 0.524056 3: 0.809301 4: 0.916710 5: 0.149759 6: 0.647463E-02 7: 0.819028 8: 0.400086 ........ .............. 4096: 0.625828 The FFT coefficients: 1: 0.631740 2: -0.506637E-03 3: 0.204992 4: 0.374571E-02 5: 0.133361 6: -0.308919E-02 7: 0.949909E-01 8: -0.887183E-02 ........ .............. 4096: 0.835070E-04 The retrieved data: 1: 0.154414E-01 2: 0.524056 3: 0.809301 4: 0.916710 5: 0.149759 6: 0.647463E-02 7: 0.819028 8: 0.400086 ........ .............. 4096: 0.625828 TEST15 For real double 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.154414E-01 0.819756 0.644887 0.353603E-01 0.907628 2: 0.524056 0.647486 0.615999 0.300492 0.505506 3: 0.809301 0.289558 0.103113 0.362978 0.325764E-01 4: 0.916710 0.606128 0.160016E-01 0.573577 0.511787 5: 0.149759 0.195894 0.938865 0.100744 0.611073 Part of the FFT coefficients: Col 1 2 3 4 5 Row 1: 0.535771 0.566881 0.587510 0.590318 0.564129 2: -0.352924E-02 -0.271147E-01 0.160491 -0.154556E-01 -0.438414E-01 3: 0.232019 0.327303 0.237239 0.175305 0.269578 4: 0.455153E-02 -0.826640E-01 -0.586933E-01 -0.426873E-01 0.609585E-01 5: 0.991899E-01 0.496055E-01 0.679210E-01 0.702973E-01 0.464044E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1: 0.154414E-01 0.819756 0.644887 0.353603E-01 0.907628 2: 0.524056 0.647486 0.615999 0.300492 0.505506 3: 0.809301 0.289558 0.103113 0.362978 0.325764E-01 4: 0.916710 0.606128 0.160016E-01 0.573577 0.511787 5: 0.149759 0.195894 0.938865 0.100744 0.611073 FFTPACK5.1D_PRB Normal end of execution. 13 May 2013 10:20:36.714 AM