Zonal operations

template<typename Count, typename Policies, typename Element, typename Zone, Rank rank>
PartitionedArray<Count, rank> zonal_diversity(Policies const &policies, PartitionedArray<Element, rank> const &array, PartitionedArray<Zone, rank> const &zones)

Return an array with per cell the number of unique element values found in the cells within the corresponding zone.