Skip to content

[ITensors] Define Sy for tJ site type #1642

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 5 commits into from
Apr 7, 2025
Merged

[ITensors] Define Sy for tJ site type #1642

Show file tree
Hide file tree
Changes from 4 commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
f43eca0
Update tj.jl
andrewkhardy Apr 4, 2025
0fe9c96
Update src/lib/SiteTypes/src/sitetypes/tj.jl
andrewkhardy Apr 7, 2025
5b6ca7b
Update test_phys_site_types.jl
andrewkhardy Apr 7, 2025
8d3db2e
Merge branch 'main' into patch-1
mtfishman Apr 7, 2025
5fceda7
Update src/lib/SiteTypes/src/sitetypes/tj.jl
mtfishman Apr 7, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
8 changes: 7 additions & 1 deletion src/lib/SiteTypes/src/sitetypes/tj.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -183,7 +183,13 @@ function op!(Op::ITensor, ::OpName"Sx", ::SiteType"tJ", s::Index)
Op[s' => 2, s => 3] = 0.5
return Op[s' => 3, s => 2] = 0.5
end

function op!(Op::ITensor, ::OpName"Sy", ::SiteType"tJ", s::Index)
Op[s' => 2, s => 3] = -0.5im
return Op[s' => 3, s => 2] = 0.5im
end
function op!(Op::ITensor, ::OpName"Sʸ", st::SiteType"tJ", s::Index)
return op!(Op, OpName("Sy"), st, s)
end
function op!(Op::ITensor, ::OpName"Sˣ", st::SiteType"tJ", s::Index)
return op!(Op, OpName("Sx"), st, s)
end
Expand Down
4 changes: 4 additions & 0 deletions test/base/test_phys_site_types.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -742,6 +742,10 @@ using ITensors.SiteTypes: has_fermion_string, op, siteind, siteinds, state
@test Sx ≈ [0.0 0.0 0; 0 0 0.5; 0 0.5 0]
Sx = Array(op(s, "Sˣ"), s', s)
@test Sx ≈ [0.0 0.0 0; 0 0 0.5; 0 0.5 0]
Sy = Array(op(s, "Sy"), s', s)
@test Sy ≈ [0.0 0.0 0; 0 0 -0.5im; 0 0.5im 0]
Sy = Array(op(s, "Sʸ"), s', s)
@test Sy ≈ [0.0 0.0 0; 0 0 -0.5im; 0 0.5im 0]
Sp = Array(op(s, "Splus"), s', s)
@test Sp ≈ [0.0 0.0 0; 0 0 1.0; 0 0 0]
Sp = Array(op(s, "Sp"), s', s)
Expand Down
Loading