separate mul! and ctmul! function calls out of op constructor

MagneticResonanceImaging · Jul 11, 2023 · 58b2698 · 58b2698
1 parent 397e69f
commit 58b2698
Showing 1 changed file with 5 additions and 2 deletions.
diff --git a/MRIOperators/src/FieldmapNFFTOp.jl b/MRIOperators/src/FieldmapNFFTOp.jl
@@ -134,10 +134,13 @@ function Base.copy(S::FieldmapNFFTOp{T,Nothing,Function,D}) where {T,D}
  D_ = length(shape)
  circTraj = S.circTraj
 
+ mul! = (res,x) -> produ!(res,x,x_tmp,shape,plans,idx,cparam,circTraj,d,p)
+ ctmul! = (res,y) -> ctprodu!(res,y,y_tmp,shape,plans,idx,cparam,circTraj,d,p)
+
  return FieldmapNFFTOp{T,Nothing,Function,D_}(S.nrow, S.ncol, false, false
- , (res,x) -> produ!(res,x,x_tmp,shape,plans,idx,cparam,circTraj,d,p)
+ , mul!
  , nothing
- , (res,y) -> ctprodu!(res,y,y_tmp,shape,plans,idx,cparam,circTraj,d,p), 0, 0, 0, false, false, false, Complex{T}[], Complex{T}[]
+ , ctmul!, 0, 0, 0, false, false, false, Complex{T}[], Complex{T}[]
  , plans, idx, circTraj, shape, cparam)
 end