  Particle-in-Cell (PIC) methods employ particle representations of unknown
fields, but also employ continuum fields for other parts of the problem. Thus
projection between particle and continuum bases is required. Moreover, we often
need to enforce conservation constraints on this projection. We derive a
mechanism for enforcement based on weak equality, and implement it in the PETSc
libraries. Scalability is demonstrated to more than 1B particles.
