I have a set of data, e.g. damped waves in a finite box of size L, so that the data is periodic in space but not time. At a particular instant in time, i.e. taking a cross section of data
at a particular time, which I'll call dataSlice
, I can create an interpolating function
intFn = Interpolation[dataSlice, InterpolationOrder -> 1, PeriodicInterpolation -> True]
and I don't run into trouble. However, if I attempt a periodic interpolation on the whole data set, it fails
intFn = Interpolation[dataSlice, InterpolationOrder -> 1, PeriodicInterpolation -> True]
(*In dimension 1 the data at the endpoints of the fundamental period are not equal*)
That makes sense, of course. For my data, it would make sense to do a semi-periodic interpolation. Does such a function exist? I would like it to function like
intFn = Interpolation[dataSlice, InterpolationOrder -> 1, PeriodicInterpolation -> {False,True,True,True}]
so the first (time) dimension would be non-periodic, and the spatial dimension is periodic.