getSci.m at c33ae98dafbc4f154b8372938c2493a50f4386a1 - concept-collection/numbl-open-quantum-systems

concept-collection / numbl-open-quantum-systems

11 lines · 243 BBlame History Raw

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