  Since proposed in [X. Zhang and C.-W. Shu, J. Comput. Phys., 229: 3091--3120,
2010], the Zhang--Shu framework has attracted extensive attention and motivated
many bound-preserving (BP) high-order discontinuous Galerkin and finite volume
schemes for various hyperbolic equations. A key ingredient in the framework is
the decomposition of the cell averages of the numerical solution into a convex
combination of the solution values at certain quadrature points, which helps to
rewrite high-order schemes as convex combinations of formally first-order
schemes. The classic convex decomposition originally proposed by Zhang and Shu
has been widely used over the past decade. It was verified, only for the 1D
quadratic and cubic polynomial spaces, that the classic decomposition is
optimal in the sense of achieving the mildest BP CFL condition. Yet, it
remained unclear whether the classic decomposition is optimal in multiple
dimensions. In this paper, we find that the classic multidimensional
decomposition based on the tensor product of Gauss--Lobatto and Gauss
quadratures is generally not optimal, and we discover a novel alternative
decomposition for the 2D and 3D polynomial spaces of total degree up to 2 and
3, respectively, on Cartesian meshes. Our new decomposition allows a larger BP
time step-size than the classic one, and moreover, it is rigorously proved to
be optimal to attain the mildest BP CFL condition, yet requires much fewer
nodes. The discovery of such an optimal convex decomposition is highly
nontrivial yet meaningful, as it may lead to an improvement of high-order BP
schemes for a large class of hyperbolic or convection-dominated equations, at
the cost of only a slight and local modification to the implementation code.
Several numerical examples are provided to further validate the advantages of
using our optimal decomposition over the classic one in terms of efficiency.
