We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 01019af commit fe980c0Copy full SHA for fe980c0
src/factorizations/truncation.jl
@@ -91,7 +91,7 @@ for f! in (:eig_trunc!, :eigh_trunc!)
91
(D, V)::Tuple{DiagonalTensorMap,AbstractTensorMap},
92
strategy::TruncationStrategy)
93
ind = MAK.findtruncated(diagview(D), strategy)
94
- V_truncated = spacetype(D)(c => length(I) for (c, I) in ind)
+ V_truncated = truncate_space(space(D, 1), ind)
95
96
D̃ = DiagonalTensorMap{scalartype(D)}(undef, V_truncated)
97
truncate_diagonal!(D̃, D, ind)
0 commit comments