diff --git a/benchmark.jl b/benchmark.jl index 12c1b9f..5ed311d 100644 --- a/benchmark.jl +++ b/benchmark.jl @@ -15,7 +15,7 @@ end T=Float32 -function V_test(r2::T)::T +function V_test(r2) return -4*exp(-r2/4) end @@ -31,7 +31,7 @@ n_image=1 for L::T in 5.0:14.0 println("Constructing H operator...") - @time H=HOperator{T}(V_test,3,3,N,L,convert(T,μ),n_image) + @time H=HOperator{T}(V_test,3,3,N,L,convert(T,0),convert(T,μ),n_image) println("Applying H 1000 times...") if GPU_mode v=CUDA.rand(Complex{T},vectorDims(H)...)