From a8cbd455abea5f1978eacb19f3c21e247494d1a6 Mon Sep 17 00:00:00 2001 From: Daan Huybrechs Date: Wed, 15 Mar 2023 14:25:22 +0100 Subject: [PATCH] dct plan --- src/bases/fourier/fouriertransforms.jl | 2 +- src/operator/solvers.jl | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/bases/fourier/fouriertransforms.jl b/src/bases/fourier/fouriertransforms.jl index 4e33c49..c58ecc4 100644 --- a/src/bases/fourier/fouriertransforms.jl +++ b/src/bases/fourier/fouriertransforms.jl @@ -176,7 +176,7 @@ end # adjoint and inv. const DCTII_PLAN{T} = FFTW.DCTPlan{T,5,true} where {T} const INV_DCTII_PLAN{T} = FFTW.DCTPlan{T,4,true} where {T} -const DCTI_PLAN{T} = FFTW.r2rFFTWPlan{T,(3,),true,1} where {T} +const DCTI_PLAN{T} = FFTW.r2rFFTWPlan{T,K,true,1} where {T,K} const DCTPLANS = Union{DCTII_PLAN,INV_DCTII_PLAN,DCTI_PLAN} diff --git a/src/operator/solvers.jl b/src/operator/solvers.jl index db01a1e..b90a524 100644 --- a/src/operator/solvers.jl +++ b/src/operator/solvers.jl @@ -98,7 +98,7 @@ end # coef_dest[:] = ldiv!(solver, coef_dest) # end -adjoint(op::GenericSolverOperator) = @warn("not implemented") +adjoint(op::GenericSolverOperator) = error("adjoint not implemented for $(op)") if VERSION ≥ v"1.7-" qr_factorization(op::DictionaryOperator) = qr(matrix(op), ColumnNorm())