I’m trying out `Interpolations.jl`

for matrix-valued data, for this example `Linear`

runs but `Cubic`

fails.

```
using LinearAlgebra
using Interpolations
f(t) = [cos(t) -sin(t); sin(t)*im cos(t)]
ts = 0:.01:2*pi
fs = [f(ti) for ti in ts]
itp = interpolate(fs, BSpline(Cubic(Line(OnGrid())))) # fails
itp = interpolate(fs, BSpline(Linear())) #works
ifx = scale(itp, ts)
```

The error message is `MethodError: no method matching zero(::Type{Array{Complex{Float64},2}})`

The documentation of `Interpolations.jl`

seems a bit sparse. How can one solve this problem? Thanks a lot!