Implement a reinterpret function

[tpadioleau]

In some situations it can be useful to reinterpret the support of a ChunkSpan:

ChunkSpan<T, DiscreteDomain<DDim0, DDim1>> chk;
DiscreteDomain<DDim2, DDim3> dom;
ChunkSpan<T, DiscreteDomain<DDim2, DDim3>> reinterpreted_chk(chk.data_handle(), dom);

This way is unsafe because the constructor cannot check the bounds.

I would like a function that takes an input a ChunkSpan and a DiscreteDomain/DiscreteElement and returns a reinterpreted ChunkSpan with a check of the bounds in case of a DiscreteDomain.