add a stage1_density_discr to QBXLayerPotentialSource

At the moment this just forwards to the density_discr as before,
but it allows marking parts of the code that need stage1 explicitly.

The next step would be to differentiate between the user-provided
density_discr and the stage1 discr and only solve densities on the
density_discr, not on stage1