Mod1 for CartesianIndex

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
1 Like