@@ -809,13 +809,16 @@ end program demo_dzfftb
809
809
810
810
## DCT type-1 (DCT-1)
811
811
812
- ### Initialize DCT-1: ` dcosti `
812
+ ### Initialize DCT-1: ` dcosti ` or ` dct_t1i `
813
813
814
814
#### Description
815
815
816
816
Initializes the array ` wsave ` which is used in subroutine ` dcost ` .
817
817
The prime factorization of ` n ` together with a tabulation of the trigonometric functions are computed and stored in ` wsave ` .
818
818
819
+ The two procedures are completely equivalent and expect the same arguments.
820
+ It is a matter of personal preference which one you choose to use.
821
+
819
822
#### Status
820
823
821
824
Experimental
@@ -851,18 +854,21 @@ program demo_dcosti
851
854
end program demo_dcosti
852
855
```
853
856
854
- ### Compute DCT-1: ` dcost `
857
+ ### Compute DCT-1: ` dcost ` or ` dct_t1 `
855
858
856
859
#### Description
857
860
858
861
Computes the DCT-1 of the input real data.
859
862
The transform is defined below at output parameter ` x ` .
860
863
864
+ The two procedures are completely equivalent and expect the same arguments.
865
+ It is a matter of personal preference which one you choose to use.
866
+
861
867
For real input data ` x ` of length ` n ` , the DCT-1 of ` x ` is equivalent, up to a
862
868
scaling factor, to the DFT of the even extension of ` x ` with length ` 2*(n-1) ` ,
863
869
where the first and last entries of the original data are not repeated in the
864
- extension. For example, the DCT-1 of input data * abcde* (size \[ 5 \] ) is
865
- equivalent to the DFT of data * abcdedcb* (size \[ 2* 4=8\] ).
870
+ extension. For example, the DCT-1 of input data * abcde* (size \( 5 \) ) is
871
+ equivalent to the DFT of data * abcdedcb* (size \( 2* 4=8\) ).
866
872
867
873
Also, ` dcost ` is the unnormalized inverse of itself. This means that a call of
868
874
` dcost ` followed by another call of ` dcost ` will multiply the input sequence ` x `
@@ -932,7 +938,7 @@ end program demo_dcost
932
938
933
939
## DCT of types 2, 3 (DCT-2, 3), a.k.a "Quarter" cosine transforms
934
940
935
- ### Initialize DCT-2, 3: ` dcosqi `
941
+ ### Initialize DCT-2, 3: ` dcosqi ` or ` dct_t23i `
936
942
937
943
#### Description
938
944
@@ -941,6 +947,9 @@ The prime factorization of `n` together with
941
947
a tabulation of the trigonometric functions are computed and
942
948
stored in ` wsave ` .
943
949
950
+ The two procedures are completely equivalent and expect the same arguments.
951
+ It is a matter of personal preference which one you choose to use.
952
+
944
953
#### Status
945
954
946
955
Experimental
@@ -978,13 +987,16 @@ program demo_dcosqi
978
987
end program demo_dcosqi
979
988
```
980
989
981
- ### Compute DCT-3: ` dcosqf `
990
+ ### Compute DCT-3: ` dcosqf ` or ` dct_t3 `
982
991
983
992
#### Description
984
993
985
994
Computes the DCT-3 of the input real data.
986
995
The transform is defined below at output parameter ` x ` .
987
996
997
+ The two procedures are completely equivalent and expect the same arguments.
998
+ It is a matter of personal preference which one you choose to use.
999
+
988
1000
Also, ` dcosqf ` (DCT-3) is the unnormalized inverse of ` dcosqb ` (DCT-2), since a
989
1001
call of ` dcosqf ` followed by a call of ` dcosqb ` will multiply the input sequence
990
1002
` x ` by ` 4*n ` .
@@ -1049,13 +1061,16 @@ program demo_dcosqf
1049
1061
end program demo_dcosqf
1050
1062
```
1051
1063
1052
- ### Compute DCT-2: ` dcosqb `
1064
+ ### Compute DCT-2: ` dcosqb ` or ` dct_t2 `
1053
1065
1054
1066
#### Description
1055
1067
1056
1068
Computes the DCT-2 of the input real data.
1057
1069
The transform is defined below at output parameter ` x ` .
1058
1070
1071
+ The two procedures are completely equivalent and expect the same arguments.
1072
+ It is a matter of personal preference which one you choose to use.
1073
+
1059
1074
For real input data ` x ` of length ` n ` , the DCT-2 of ` x ` is equivalent, up to a
1060
1075
scaling factor, to the DFT of the even extension of ` x ` with length ` 4*n ` ,
1061
1076
where all the even-frequency entries are zero.
0 commit comments