Allow BasicDimArrays like DimPoints, DimIndices, etc in DimStack

[asinghvi17]

this allows:

DimStack((; data = da, points = DimPoints(da)))`

very useful if you are trying to calculate e.g slope with stencils, since you then have the coordinates available to you.

[rafaqz]

Hmmm. It doesn't really fit the model of what a DimStack is.

There is a rule across DD of no duplication of dimensions. That keeps everything very simple and sane as it can be

We may need a ProtoDimindices <: AbstractArray etc that has gets turned into DimIndices, like an Array becomes a DimArray when combined with dimensions. (Yes this happens totally at compile time for every getindex)

Then if you set new dimensions, DimPoints etc layers would also be updated, etc etc

[rafaqz]

Or Proto{DimIndices}(args) I think could work for everything lazy.

[asinghvi17]

I can see where you're going with that, maybe that makes sense? But it's very confusing to tell the user that DimPoints is not a DimArray. We don't necessarily need to promote this workflow but it should either work or error informatively...

[rafaqz]

I mean, it literally cant work the way DimStack works. DimPoints is not a DimArray.

An AbstractDimArray wraps another array, it even has a type parameter A for the type of that other array - thats the fundamental split with AbstractDimArray and AbstractBasicDimArray.

You can of course just collect DimIndices etc and use them as an array. But its a very different thing. The values are no longer actually connected to the dimensions.

[rafaqz]

This approach isn't really workable, so closing