Skip to content

WIP: implementation of DiagonalTensorMap #174

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
Jutho merged 14 commits into from
Nov 18, 2024
Merged

WIP: implementation of DiagonalTensorMap #174

Changes from 1 commit
Commits
Show all changes
14 commits
Select commit Hold shift + click to select a range
41e3cbb
start implementation of DiagonalTensorMap
Jutho Nov 10, 2024
1c723ea
further progress on DiagonalTensorMap [skip ci]
Jutho Nov 12, 2024
435a2a3
more fixes and first tests
Jutho Nov 15, 2024
375d194
a bit more functionality and tests
Jutho Nov 16, 2024
3e78086
resolve conflict in Project.toml
Jutho Nov 16, 2024
17daa96
finegrain homespace permutation and selection
Jutho Nov 16, 2024
e4fdde6
fix and add more tests
Jutho Nov 17, 2024
448ca50
implement suggestions and more tests
Jutho Nov 17, 2024
e89e5c3
replace compose with mul! and compose_dest
Jutho Nov 17, 2024
486a534
Also rename testfile
lkdvos Nov 18, 2024
b9d02fa
Extend `mul!` to include `alfa` and `beta`
lkdvos Nov 18, 2024
fd9cb9f
Add TensorOperations specializations
lkdvos Nov 18, 2024
5b1e01f
Add show for Diagonal
lkdvos Nov 18, 2024
0111bc9
changelog, docs and final fixes
Jutho Nov 18, 2024
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
17 changes: 14 additions & 3 deletions src/tensors/diagonal.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -198,10 +198,21 @@ function Base.zero(d::DiagonalTensorMap)
return DiagonalTensorMap(zero.(d.data), d.domain)
end

function compose(d1::DiagonalTensorMap, d2::DiagonalTensorMap)
d1.domain == d2.domain || throw(SpaceMismatch())
return DiagonalTensorMap(d1.data .* d2.data, d1.domain)
function compose_dest(d1::DiagonalTensorMap, d2::DiagonalTensorMap)
A = promote_type(storagetype(d1), storagetype(d2))
S = spacetype(d1)
T = scalartype(A)
return DiagonalTensorMap{T,S,A}(undef, d1.domain)
end

function LinearAlgebra.mul!(dC::DiagonalTensorMap,
dA::DiagonalTensorMap,
dB::DiagonalTensorMap)
dC.domain == dA.domain == dB.domain || throw(SpaceMismatch())
dC.data .= dA.data .* dB.data
return dC
end

Base.inv(d::DiagonalTensorMap) = DiagonalTensorMap(inv.(d.data), d.domain)
function Base.:\(d1::DiagonalTensorMap, d2::DiagonalTensorMap)
d1.domain == d2.domain || throw(SpaceMismatch())
Expand Down