julia> A = Tensor(Dense(SparseTriangle{2}(Element(0.0))), collect(reshape(1:27, 3, 3, 3)))
Dense [:,:,1:3]
├─ [:, :, 1]: SparseTriangle{2} (0.0) [:,1:3]
│  ├─ [1, 1]: 1.0
│  ├─ [2, 1]: 4.0
│  ├─ [2, 2]: 5.0
│  ├─ [3, 1]: 7.0
│  ├─ [3, 2]: 8.0
│  └─ [3, 3]: 9.0
├─ [:, :, 2]: SparseTriangle{2} (0.0) [:,1:3]
│  ├─ [1, 1]: 10.0
│  ├─ [2, 1]: 13.0
│  ├─ [2, 2]: 14.0
│  ├─ [3, 1]: 16.0
│  ├─ [3, 2]: 17.0
│  └─ [3, 3]: 18.0
└─ [:, :, 3]: SparseTriangle{2} (0.0) [:,1:3]
   ├─ [1, 1]: 19.0
   ├─ [2, 1]: 22.0
   ├─ [2, 2]: 23.0
   ├─ [3, 1]: 25.0
   ├─ [3, 2]: 26.0
   └─ [3, 3]: 27.0

julia> C = Scalar(0)
Scalar{0, Int64}(0)
julia> @finch for k = _, j = _, i = _
        C[] += A[i, j, k]
    end
NamedTuple()

