concept-collection / numbl-open-quantum-systems
numbl-open-quantum-systems / getSci.m
11 lines · 243 BBlameHistoryRaw
1function Sci = getSci(sc,i,Nspins)
2Is = eye(2);
3Op_total = cell(1,Nspins);
4for site = 1:Nspins
5 Op_total{site} = Is+double(eq(i,site))*(sc-Is);
6end
7Sci = Op_total{1};
8for site = 2:Nspins
9 Sci = kron(Sci,Op_total{site});
10end
11end