Solution (source code)

= Solution

For a bounded $C^1$ <classical solution> and a smooth convex entropy $\eta$, let its <entropy flux for a scalar conservation law> satisfy $\psi'=f'\eta'$. The <chain rule> gives
$$
\partial_t\eta(u)+\partial_x\psi(u)=\eta'(u)\bigl(u_t+f'(u)u_x\bigr)=0.
$$
Integration against a nonnegative <test function>, including its initial boundary contribution, gives equality in the entropy inequality. Approximate a convex piecewise $C^1$ entropy by smooth convex functions on the bounded state range. Their values and flux antiderivatives converge uniformly, so the integrated equality passes to the limit. Therefore \b[every bounded classical solution is an entropy solution]. The bounded qualifier is required by the definition's $L^\infty$ membership; an arbitrary unbounded $C^1$ solution would not satisfy that membership condition.

Conversely, the affine entropies $\eta_+(z)=z$ and $\eta_-(z)=-z$ are both convex, with fluxes $f$ and $-f$. Their inequalities force the same linear weak integral to be both nonnegative and nonpositive. It is therefore zero for every nonnegative test, hence for every test by subtraction of nonnegative compactly supported tests. Thus $\boxed{\text{entropy solution}\Longrightarrow\text{weak solution}.}$ No assumption that the flux itself is convex is used.