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 0fe95f4 commit fce50b6Copy full SHA for fce50b6
src/operators.jl
@@ -90,7 +90,7 @@ function embed(basis_l::CompositeBasis, basis_r::CompositeBasis,
90
(opsb.basis_r == basis_r.bases[idxsb]) || throw(IncompatibleBases())
91
end
92
93
- S = length(operators) > 1 ? mapreduce(eltype, promote_type, operators) : Any
+ S = length(operators) > 0 ? mapreduce(eltype, promote_type, operators) : Any
94
embed_op = tensor([i ∈ indices_sb ? ops_sb[indexin(i, indices_sb)[1]] : identityoperator(T, S, basis_l.bases[i], basis_r.bases[i]) for i=1:N]...)
95
96
# Embed all joint-subspace operators.
0 commit comments