Wrong distributive scaling properties

[apkille]

julia> @op A
A

julia> (1+im)*A
1 + 1imA

[apkille]

Could be that printing functionalities are misleading here:

julia> QuantumSymbolics.arguments((1+im)*A)
2-element Vector{Any}:
 1 + 1im
    A

Will address in upcoming PR.

[apkille]

Also, we currently have

julia> @bra b; @ket k;

julia> (b*k)*A
ERROR: TypeError: non-boolean (SymbolicUtils.BasicSymbolic{Bool}) used in boolean context
Stacktrace:
 [1] *(c::SBraKet, x::SOperator)
   @ QuantumSymbolics ~/Documents/Julia Packages/QuantumJulia/QuantumSymbolics.jl/src/QSymbolicsBase/basic_ops_homogeneous.jl:33
 [2] top-level scope
   @ REPL[16]:1

[apkille]

Could be that printing functionalities are misleading here:

julia> QuantumSymbolics.arguments((1+im)*A)
2-element Vector{Any}:
 1 + 1im
    A

Will address in upcoming PR.

Another printing issue is when you multiply a tensor product with an operator, say A*(B ⊗ C), which currently gives AB⊗C, rather than the preferred A(B ⊗ C).