Get subdomain assemble working for point integrals
Fix assemble of scalars in parallel
 -- only include contribution from vertex with the lowest rank from shared processes