julia> (σ1x, σ1y, σ1z), (σ2x, σ2y, σ2z) = @spin(σ1, σ2)
((σ1ˣ, σ1ʸ, σ1ᶻ), (σ2ˣ, σ2ʸ, σ2ᶻ))
julia> (σ1x + σ2y) * (σ2x + σ1z)
ERROR: InexactError: Int64(0 - 1im)
Stacktrace:
[1] Real at ./complex.jl:37 [inlined]
[2] convert at ./number.jl:7 [inlined]
[3] convert at /home/blegat/.julia/packages/MultivariatePolynomials/bsuXE/src/conversion.jl:43 [inlined]
[4] push! at ./array.jl:934 [inlined]
[5] mul_to_terms!(::Array{CondensedMatterSOS.SpinTerm{Int64},1}, ::CondensedMatterSOS.SpinPolynomial{Int64}, ::CondensedMatterSOS.SpinPolynomial{Int64}) at /home/blegat/.julia/packages/MultivariatePolynomials/bsuXE/src/operators.jl:195
[6] *(::CondensedMatterSOS.SpinPolynomial{Int64}, ::CondensedMatterSOS.SpinPolynomial{Int64}) at /home/blegat/.julia/packages/MultivariatePolynomials/bsuXE/src/operators.jl:201
[7] top-level scope at REPL[10]:1
