Indexing

Indexing into a (Symmetric)Tensor{dim, order} is performed like for an Array of dimension order.

julia> A = rand(Tensor{2, 2});

julia> A[1, 2]
0.5662374165061859

julia> B = rand(SymmetricTensor{4, 2});

julia> B[1, 2, 1, 2]
0.24683718661000897

Slicing will produce a Tensor of lower order.

julia> A = rand(Tensor{2, 2});

julia> A[:, 1]
2-element Tensor{1,2,Float64,2}:
 0.5908446386657102
 0.7667970365022592

Since Tensors are immutable there is no setindex! function defined on them. Instead, use the functionality to create tensors from functions as described here. As an example, this sets the [1,2] index on a tensor to one and the rest to zero:

julia> Tensor{2, 2}((i,j) -> i == 1 && j == 2 ? 1.0 : 0.0)
2×2 Tensor{2,2,Float64,4}:
 0.0  1.0
 0.0  0.0

For symmetric tensors, note that you should only set the lower triangular part of the tensor:

julia> SymmetricTensor{2, 2}((i,j) -> i == 1 && j == 2 ? 1.0 : 0.0)
2×2 SymmetricTensor{2,2,Float64,3}:
 0.0  0.0
 0.0  0.0

julia> SymmetricTensor{2, 2}((i,j) -> i == 2 && j == 1 ? 1.0 : 0.0)
2×2 SymmetricTensor{2,2,Float64,3}:
 0.0  1.0
 1.0  0.0