TypeError: in typeassert, expected Float64, got a value of type ForwardDiff.Dual{Nothing, Float64, 12}

[nilshg]

I meant this line here:

complength = zeros(length(sz),length(sz))

which should probably be

complength = zeros(eltype(ϕ_i), length(sz),length(sz))