This relates to issue #6 - there is a remaining step which works around that bug, but it should not be required. Specifically the line
fill!(dA, zero(eltype(dA)))
in oneparticlebasis.jl should not be required since the 1p basis doesn't add into the storage dA but writes into it. So somewhere we are reading lots of zeros. This suggests also that there are some operations that we shouldn't be performing.
This relates to issue #6 - there is a remaining step which works around that bug, but it should not be required. Specifically the line
in
oneparticlebasis.jl
should not be required since the 1p basis doesn't add into the storagedA
but writes into it. So somewhere we are reading lots of zeros. This suggests also that there are some operations that we shouldn't be performing.