Pensieve header: Testing the commutativity of two independent Nwu’s.

SetDirectory["C:\\drorbn\\AcademicPensieve\\People\\Vo"];

NOE-0

(Subscript[R, i_,j_]^(+)):=𝔼[(Subscript[b, i])(Subscript[c, j])+(Subscript[b, i]^(-1))((E^(Subscript[b, i]))-1)(Subscript[u, i])(Subscript[w, j])]; (Subscript[R, i_,j_]^(-)):=𝔼[-(Subscript[b, i])(Subscript[c, j])+(Subscript[b, i]^(-1))((E^(-(Subscript[b, i])))-1)(Subscript[u, i])(Subscript[w, j])];

CF[ω_.𝔼[Q_]] := Simplify[ω]𝔼[Simplify[Q]];
𝔼 /: 𝔼[Q1_]𝔼[Q2_] := CF@𝔼[Q1+Q2];
ω1_.𝔼[Q1_] ≡ ω2_.𝔼[Q2_] := Simplify[ω1==ω2∧Q1==Q2];

(Subscript[N, (Subscript[u, i_])(Subscript[c, j_])->k_])[ω_.𝔼[Q_]] := CF[
ω 𝔼[(E^(-γ))β (Subscript[u, k])+γ (Subscript[c, k])+(Q/.(Subscript[c, j])|(Subscript[u, i])->0)]/.{γ->(∂_(Subscript[c, j]))Q,β->(∂_(Subscript[u, i]))Q}];
(Subscript[N, (Subscript[w, i_])(Subscript[c, j_])->k_])[ω_.𝔼[Q_]] := CF[
ω 𝔼[(E^γ)α (Subscript[w, k])+γ (Subscript[c, k])+(Q/.(Subscript[c, j])|(Subscript[w, i])->0)]/.{γ->(∂_(Subscript[c, j]))Q,α->(∂_(Subscript[w, i]))Q}];
(Subscript[N, (Subscript[w, i_])(Subscript[u, j_])->k_])[ω_.𝔼[Q_]] := CF[
ν ω 𝔼[-(Subscript[b, k])ν α β+ν β (Subscript[u, k])+ν δ (Subscript[u, k])(Subscript[w, k])+ν α (Subscript[w, k]) +(Q/.(Subscript[w, i])|(Subscript[u, j])->0)] /. ν-> (((1+(Subscript[b, k]) δ))^(-1)) /. {α->(∂_(Subscript[w, i]))Q/.(Subscript[u, j])->0, β->(∂_(Subscript[u, j]))Q/.(Subscript[w, i])->0, δ->(∂_((Subscript[w, i]),(Subscript[u, j])))Q}];

(Subscript[m, i_,j_->k_])[ω_.𝔼[Q_]] :=CF[Module[{x},
(ω 𝔼[Q] /. (Subscript[b, i|j])->(Subscript[b, k]) // (Subscript[N, (Subscript[w, i])(Subscript[c, j])->x])//(Subscript[N, (Subscript[u, i])(Subscript[c, x])->x])//(Subscript[N, (Subscript[w, x])(Subscript[u, j])->x])) /. {(Subscript[c, i])->(Subscript[c, k]),(Subscript[w, j])->(Subscript[w, k]),(Subscript[y_, x]):>(Subscript[y, k])}]]

Z0=(((Subscript[R, 1,2]^(+))//(Subscript[m, 1,2->1])) /. (Subscript[x_, 1]):>(Subscript[x, 1])+(Subscript[x, 2]))((Subscript[R, 3,5]^(+))//(Subscript[m, 3,5->3]))((Subscript[R, 4,6]^(+))//(Subscript[m, 4,6->1]))

Z0 // (Subscript[m, 1,3->1])//(Subscript[m, 2,4->2])

Z0 // (Subscript[m, 3,1->1])//(Subscript[m, 4,2->2])

R[p_,i_,j_] := 𝔼[p (Subscript[b, i])(Subscript[c, j])+(Subscript[b, i]^(-1))((E^(p (Subscript[b, i])))-1)(Subscript[u, i])(Subscript[w, j])];

R[-1,i,j]

R[1/2,1,2]R[1/2,3,4] //(Subscript[m, 1,3->1])//(Subscript[m, 2,4->2])