I am using Julia version 1.1. I am working a lot with matrices that can be constructed from smaller matrices, e.g., the Pauli matrices. It is not clear to me how to efficiently construct big matrices by using a set of smaller matrices in Julia, i.e., directly write the smaller matrix into a certain index position.
Julias kron
is not satisfactory since I would need to generate several "big matrices" to get my final result. E.g. I would like to create something like this (this is only a very small example)
sy = [[0 -im]; [im 0]]
M = [[0 sy adjoint(sy)];
[adjoint(sy) 0 sy];
[sy adjoint(sy) 0]]
It would be possible to do so by doing two kronecker products adding the two results. However this would be a huge waste, especially if the matrices got bigger.
I also already tried to work with the package BlockArrays.jl
but realised that it does not fullfill my need.
In the end I want to be able to address "matrix blocks" of my big matrix so that I can directly assign the construction matrices to the right position, for the above example this would look like the following (I did not use a loop here to make my point clearer):
M[1, 2] = sy
M[1, 3] = adjoint(sy)
M[2, 1] = adjoint(sy)
M[2, 3] = sy
M[3, 1] = sy
M[3, 2] = adjoint(sy)
I realise that this means reducing my original big array indices to something like array "block indices".
I thought about doing this with views, where I create a matrix of SubArrays
that I can then address with the matrix block index notation e.g.
S0 = view(M, 1:2, 1:2)
S1 = view(M, 1:2, 2:4)
S2 = view(M, 1:2, 4:6)
...
Viewmatrix = [[S0 S1 S2]; [S3 S4 S5]; [S6 S7 S8]]
Viewmatrix[1, 2] .= sy
Viewmatrix[1, 3] .= adjoint(sy)
...
Now it is unclear to me how one would actually go about this and write such a view matrix in general or if this is even a feasable way to address the problem. If there is a better way to approach this problem I would like to know it.
BlockArray.jl
always only divide the matrix into 2x2 blocks of arbitrary size. I need an arbitrary amount of blocks depending on the size of my submatrix and the size of the "big" matrix.