Fix tmp array type for DCT and DST operator

Author: nHackel

ext/LinearOperatorFFTWExt/DCTOp.jl

@@ -34,7 +34,7 @@ returns a `DCTOpImpl <: AbstractLinearOperator` which performs a DCT on a given
 """
 function LinearOperatorCollection.DCTOp(T::Type; shape::Tuple, dcttype=2)
 
-  tmp=Array{Complex{real(T)}}(undef, shape) 
+  tmp=Array{T}(undef, shape) 
   if dcttype == 2
     plan = plan_dct!(tmp)
     iplan = plan_idct!(tmp)

ext/LinearOperatorFFTWExt/DSTOp.jl

@@ -32,7 +32,7 @@ returns a `LinearOperator` which performs a DST on a given input array.
 * `shape::Tuple`  - size of the array to transform
 """
 function LinearOperatorCollection.DSTOp(T::Type; shape::Tuple)
-  tmp=Array{Complex{real(T)}}(undef, shape)
+  tmp=Array{T}(undef, shape)
 
   plan = FFTW.plan_r2r!(tmp,FFTW.RODFT10)
   iplan = FFTW.plan_r2r!(tmp,FFTW.RODFT01)