Replace this with, using ntuple:
function mod1(CI::CartesianIndex{N}, MOD::CartesianIndex{N}) where N
vals = ntuple(i -> Base.mod1(CI[i], MOD[i]), N)
CartesianIndex(vals...)
end
Replace this with, using ntuple:
function mod1(CI::CartesianIndex{N}, MOD::CartesianIndex{N}) where N
vals = ntuple(i -> Base.mod1(CI[i], MOD[i]), N)
CartesianIndex(vals...)
end