great work on parallelizing it! I know how difficult it can be.
qp = pf.h.quad_proj(0, "Density")
I'll try to test it out later today. One comment - for add_projection we have this:
but above it's in the other order. Do we want to make these more consistent?