17 November 2007 9:06:35.591 AM FFTPACK5_PRB FORTRAN90 version A set of tests for FFTPACK5. TEST01 For single precision complex fast Fourier transforms, 1D, CFFT1I initializes the transforms, CFFT1F does a forward transforms; CFFT1B does a backward transforms. The number of data items is N = 4096 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 9 0.158831 0.390522 10 0.503769 -0.474654 The FFT coefficients: 1 -0.779193E-03 -0.169766E-02 2 0.269171E-02 0.171117E-01 3 -0.985016E-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 9 -0.710866E-02 0.339033E-04 10 0.166492E-02 0.526741E-02 The retrieved data: 1 -0.122847 -0.187107E-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 9 0.158831 0.390522 10 0.503769 -0.474654 TEST02 For single precision complex fast Fourier transforms, 2D, CFFT2I initializes the transforms, CFFT2F does a forward transforms; CFFT2B does a backward transforms. The data is stored in an L by M array, with L = 32 M = 64 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 single precision complex fast Fourier transforms, 1D, multiple CFFTMI initializes the transforms, CFFTMF does a forward transforms; CFFTMB does a backward transforms. The number of sequences is LOT = 6 The length of each sequence is N = 32 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 single precision real fast cosine transforms, 1D, COSQ1I initializes the transforms, COSQ1F does a forward transforms; COSQ1B does a backward transforms. The number of data items is N = 4096 The original data: 1 0.15441427E-01 2 0.52405602 3 0.80930066 4 0.91671026 5 0.14975941 6 0.64746253E-02 7 0.81902778 8 0.40008637 9 0.25180870 10 0.14908369 The FFT coefficients: 1 0.62831903 2 -.21445246 3 0.12118940 4 -.86556278E-01 5 0.69241442E-01 6 -.67744911E-01 7 0.46076778E-01 8 -.56036085E-01 9 0.34167781E-01 10 -.45675509E-01 The retrieved data: 1 0.15441358E-01 2 0.52405584 3 0.80930054 4 0.91671020 5 0.14975947 6 0.64746439E-02 7 0.81902778 8 0.40008652 9 0.25180852 10 0.14908373 TEST05 For single precision real fast cosine transforms, 1D, multiple COSQMI initializes the transforms, COSQMF does a forward transforms; COSQMB does a backward transforms. The number of sequences is LOT = 6 The length of each sequence is N = 32 Part of the original data: Col 1 2 3 4 5 Row 1 0.154414E-01 0.819757 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.800813E-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.438616E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.154414E-01 0.819757 0.644887 0.353603E-01 0.907628 2 0.524056 0.647485 0.616000 0.300492 0.505506 3 0.809301 0.289558 0.103113 0.362978 0.325764E-01 4 0.916710 0.606128 0.160017E-01 0.573576 0.511787 5 0.149759 0.195894 0.938865 0.100744 0.611072 TEST06 For single precision real fast cosine transforms, 1D, COST1I initializes the transforms, COST1F does a forward transforms; COST1B does a backward transforms. The number of data items is N = 4096 The original data: 1 0.15441427E-01 2 0.52405602 3 0.80930066 4 0.91671026 5 0.14975941 6 0.64746253E-02 7 0.81902778 8 0.40008637 9 0.25180870 10 0.14908369 The FFT coefficients: 1 0.49683216 2 -.82119294E-02 3 -.15197121E-03 4 -.67906347E-02 5 0.10358667E-01 6 -.14575866E-01 7 0.10416001E-02 8 -.14773483E-01 9 -.29696987E-02 10 -.13238335E-01 The retrieved data: 1 0.15211732E-01 2 0.52423316 3 0.80931890 4 0.91670543 5 0.14977117 6 0.64959014E-02 7 0.81902009 8 0.40008989 9 0.25179392 10 0.14908217 TEST07 For single precision real fast cosine transforms, 1D, multiple COSTMI initializes the transforms, COSTMF does a forward transforms; COSTMB does a backward transforms. The number of sequences is LOT = 6 The length of each sequence is N = 32 Part of the original data: Col 1 2 3 4 5 Row 1 0.154414E-01 0.819757 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.476391 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.899186E-01 -0.775478E-01 -0.566894E-01 -0.768210E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.154415E-01 0.819757 0.644887 0.353613E-01 0.907629 2 0.524056 0.647486 0.616000 0.300491 0.505505 3 0.809301 0.289558 0.103113 0.362978 0.325767E-01 4 0.916709 0.606128 0.160022E-01 0.573576 0.511787 5 0.149760 0.195894 0.938864 0.100744 0.611072 TEST075 For double precision real fast cosine transforms, 1D, DCOSQ1I initializes the transforms, DCOSQ1F does a forward transforms; DCOSQ1B does a backward transforms. The number of data items is N = 4096 The original data: 1 0.15441427E-01 2 0.52405601 3 0.80930069 4 0.91671028 5 0.14975941 6 0.64746253E-02 7 0.81902779 8 0.40008638 9 0.25180872 10 0.14908369 The FFT coefficients: 1 0.62831907 2 -.21445245 3 0.12118939 4 -.86556279E-01 5 0.69241438E-01 6 -.67744910E-01 7 0.46076784E-01 8 -.56036087E-01 9 0.34167785E-01 10 -.45675506E-01 The retrieved data: 1 0.15441427E-01 2 0.52405601 3 0.80930069 4 0.91671028 5 0.14975941 6 0.64746253E-02 7 0.81902779 8 0.40008638 9 0.25180872 10 0.14908369 TEST08 For double precision real fast Fourier transforms, 1D, DFFT1I initializes the transforms, DFFT1F does a forward transforms; DFFT1B does a backward transforms. The number of data items is N = 4096 The original data: 1 0.15441427E-01 2 0.52405601 3 0.80930069 4 0.91671028 5 0.14975941 6 0.64746253E-02 7 0.81902779 8 0.40008638 9 0.25180872 10 0.14908369 The FFT coefficients: 1 0.49678912 2 -.23654481E-03 3 -.50048965E-03 4 0.10274917E-01 5 0.37467898E-02 6 0.93833513E-03 7 -.30616504E-02 8 -.30867563E-02 9 -.88267511E-02 10 0.20676679E-02 The retrieved data: 1 0.15441427E-01 2 0.52405601 3 0.80930069 4 0.91671028 5 0.14975941 6 0.64746253E-02 7 0.81902779 8 0.40008638 9 0.25180872 10 0.14908369 TEST09 For single precision real fast Fourier transforms, 1D, RFFT1I initializes the transforms, RFFT1F does a forward transforms; RFFT1B does a backward transforms. The number of data items is N = 4096 The original data: 1 0.15441427E-01 2 0.52405602 3 0.80930066 4 0.91671026 5 0.14975941 6 0.64746253E-02 7 0.81902778 8 0.40008637 9 0.25180870 10 0.14908369 The FFT coefficients: 1 0.49678910 2 -.23654522E-03 3 -.50048961E-03 4 0.10274917E-01 5 0.37467894E-02 6 0.93833520E-03 7 -.30616506E-02 8 -.30867567E-02 9 -.88267513E-02 10 0.20676679E-02 The retrieved data: 1 0.15441388E-01 2 0.52405608 3 0.80930066 4 0.91671026 5 0.14975947 6 0.64747632E-02 7 0.81902766 8 0.40008625 9 0.25180873 10 0.14908370 TEST10 For single precision real fast Fourier transforms, 2D, RFFT2I initializes the transforms, RFFT2F does a forward transforms; RFFT2B does a backward transforms. The L by M data is stored in an LDIM by M array, with L = 32 LDIM = 34 M = 64 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.179073E-02 -0.706918E-02 0.360346E-02 -0.493869E-02 2 0. -0.924539E-03 0.234106E-02 0.805991E-02 0.363998E-02 3 0.139190E-02 -0.857094E-02 0.182175E-01 0.411902E-02 -0.107389E-01 4 -0.430501E-02 0.202246E-01 0.202669E-02 -0.613034E-02 0.115338E-02 5 -0.912771E-02 0.286053E-02 -0.677017E-02 -0.529703E-02 0.193206E-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.826156E-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 single precision real fast Fourier transforms, 1D, multiple RFFTMI initializes the transforms, RFFTMF does a forward transforms; RFFTMB does a backward transforms. The number of sequences is LOT = 6 The length of each sequence is N = 32 Part of the original data: Col 1 2 3 4 5 Row 1 0.154414E-01 0.819757 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.975671E-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.647485 0.615999 0.300492 0.505506 3 0.809301 0.289558 0.103113 0.362978 0.325765E-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 single precision real fast sine transforms, 1D, SINQ1I initializes the transforms, SINQ1F does a forward transforms; SINQ1B does a backward transforms. The number of data items is N = 4096 The original data: 1 0.15441427E-01 2 0.52405602 3 0.80930066 4 0.91671026 5 0.14975941 6 0.64746253E-02 7 0.81902778 8 0.40008637 9 0.25180870 10 0.14908369 The FFT coefficients: 1 0.63644505 2 0.20595019 3 0.12816204 4 0.81669316E-01 5 0.84993973E-01 6 0.51287916E-01 7 0.56464050E-01 8 0.34127738E-01 9 0.37419993E-01 10 0.18238375E-01 The retrieved data: 1 0.15441403E-01 2 0.52405614 3 0.80930072 4 0.91671032 5 0.14975937 6 0.64746141E-02 7 0.81902784 8 0.40008628 9 0.25180870 10 0.14908376 TEST13 For single precision real fast sine transforms, 1D, multiple SINQMI initializes the transforms, SINQMF does a forward transforms; SINQMB does a backward transforms. The number of sequences is LOT = 6 The length of each sequence is N = 32 Part of the original data: Col 1 2 3 4 5 Row 1 0.154414E-01 0.819757 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.776554E-03 -0.136147E-02 0.497461E-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.647485 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.611072 TEST14 For single precision real fast sine transforms, 1D, SINT1I initializes the transforms, SINT1F does a forward transforms; SINT1B does a backward transforms. The number of data items is N = 4096 The original data: 1 0.15441427E-01 2 0.52405602 3 0.80930066 4 0.91671026 5 0.14975941 6 0.64746253E-02 7 0.81902778 8 0.40008637 9 0.25180870 10 0.14908369 The FFT coefficients: 1 0.63174033 2 -.50661963E-03 3 0.20499378 4 0.37456874E-02 5 0.13336287 6 -.30891683E-02 7 0.94992638E-01 8 -.88718403E-02 9 0.68625592E-01 10 -.18615942E-01 The retrieved data: 1 0.18086851E-01 2 0.52409309 3 0.81036377 4 0.91665316 5 0.15040898 6 0.64504296E-02 7 0.81948835 8 0.40001872 9 0.25209671 10 0.14910299 TEST15 For single precision real fast sine transforms, 1D, multiple SINTMI initializes the transforms, SINTMF does a forward transforms; SINTMB does a backward transforms. The number of sequences is LOT = 6 The length of each sequence is N = 32 Part of the original data: Col 1 2 3 4 5 Row 1 0.154414E-01 0.819757 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.352925E-02 -0.271148E-01 0.160491 -0.154556E-01 -0.438414E-01 3 0.232019 0.327303 0.237239 0.175305 0.269578 4 0.455154E-02 -0.826640E-01 -0.586933E-01 -0.426873E-01 0.609585E-01 5 0.991900E-01 0.496054E-01 0.679210E-01 0.702973E-01 0.464044E-01 Part of the retrieved data: Col 1 2 3 4 5 Row 1 0.154410E-01 0.819755 0.644885 0.353586E-01 0.907628 2 0.524056 0.647486 0.616000 0.300492 0.505506 3 0.809301 0.289558 0.103112 0.362978 0.325764E-01 4 0.916710 0.606128 0.160021E-01 0.573577 0.511788 5 0.149760 0.195894 0.938865 0.100744 0.611072 FFTPACK5_PRB Normal end of execution. 17 November 2007 9:06:35.675 AM