numpy 배열로 배열 채우기

Question

Kroneckerproduct에서 불특정 (한정) 행렬을 결합하고 싶습니다. 이 작업을 수행하기 위해 행렬에 행렬을 저장하려고하지만이를 수행하는 방법을 모르겠습니다. 현재 내가 가지고있는 것 :

for i in range(LNew-2): 
    for j in range(LNew-2): 
     Bulk = np.empty(shape=(LNew-1,LNew-1)) 
     if i == j: 
      Bulk[i,j]=H2 
     else: 
      Bulk[i,j]=idm 

여기서 H2와 idm은 모두 Kronecker 제품에서 결합하려는 행렬입니다. 그러나 벌크가 ndarray 객체이기 때문에 내부에 arraylike 객체를 받아들이지 않을 것이라고 생각합니다.

편집 :

이 나는이 아이디어를 사용하고자하는 기능입니다. 나는 그것을 양자 스핀 체인을위한 해밀 토니안 행렬을 만들기 위해 사용하고 있습니다. 따라서 H2는 2 개의 입자 체인에 대한 해밀턴이다. H2는 4x4 행렬이고 idm은 2x2 항등 행렬이다.

이제 세 입자 체인 np.kron이다 (H2, IDM) + np.kron (IDM, H2)

네 입자 np.kron (np.kron위한

(H2, idm) + np.kron (idm, np.kron (H2, idm)) + np.kron (idm, np.kron (idm, H2)) 등등.

def ExpandHN(LNew): 
idm = np.identity(2) 
H2 = GetH(2,'N') 
HNew = H2 
for i in range(LNew-2): 
    for j in range(LNew-2): 
     Bulk = np.empty(shape=(LNew-1,LNew-1)) 
     if i == j: 
      Bulk[i,j]=H2 
     else: 
      Bulk[i,j]=idm 
i = 0 
for i in range(LNew-2): 
    for j in range(LNew-3): 
     HNew += np.kron(Bulk[i,j],Bulk[i,j+1]) #something like this 

return HNew 

볼 수 있듯이 두 번째 for 루프 집합은 효과가 없습니다.

누군가가 완전히 다르지만 작동하는 해결책을 가지고 있다면 나는 그 점에도 만족할 것이라고 말했습니다.

Answer 1

제가 올바르게 이해한다면 질문은 numpy로 배열의 배열을 만드는 방법에 관한 것입니다.

이

Bulk = dict() 
for i in range(LNew-2): 
    for j in range(LNew-2): 
     if i == j: 
      Bulk[(i,j)]=H2 
     else: 
      Bulk[(i,j)]=idm 

키와 튜플의 사용은 당신이 행렬의 배열과 같은 인덱스를 유지 관리 할 수 ​​있습니다 : 나는 표준 파이썬 객체 dict를 사용하는 것이 좋습니다 것입니다. 두 루프 모두 외부에서 벌크를 정의해야합니다 (모든 경우).

HTH