On entry: the data must be stored in
X as if in a two-dimensional array of dimension
1:M,0:N+1; each of the
m sequences is stored in a
row of the array.
In other words, if the data values of the
pth sequence to be transformed are denoted by
xjp, for
j=0,1,…,n-1 and
p=1,2,…,m, then the first
mn elements of the array
X must contain the values
The
n+1th and
n+2th elements of each row
xnp , x n+1 p , for
p=1,2,…,m, are required as workspace. These
2m elements may contain arbitrary values as they are set to zero by the routine.
On exit: the
m quarter-wave cosine transforms stored as if in a two-dimensional array of dimension
1:M,0:N+1. Each of the
m transforms is stored in a
row of the array, overwriting the corresponding original sequence.
If the
n components of the
pth quarter-wave cosine transform are denoted by
x^kp, for
k=0,1,…,n-1 and
p=1,2,…,m, then the
mn+2 elements of the array
X contain the values