An object where the size of the local magnetization $\mathbf{M}$ is fixed ($|\mathbf{M}| =M =$ Const.) can be thought of alternatively as an object where on the boundary a surface current $\mathbf{K}=\mathbf{M}\times \mathbf{n}$ flows. In the last formula, $\mathbf{n}$ is the local normal vector of the object's boundary. Additionally, an internal current $\mathbf{J}=\nabla \times \mathbf{M}$ will flow if the magnetization $\mathbf{M}$ varies in direction. $\mathbf{B}$ is a function of these two types of current via the Biot-Savart law.
Given a volume $V$ of a material with magnetization $M$, we may bound from below the strongest obtainable (local) magnetic field strength $B_{max}(M,V)$ by considering rod configurations. The previous discussion has revealed that a rod of certain length and radius which is magnetized along its axis produces the same field as a solenoid carrying a surface current $M$ going around the cylindrical boundary. A solenoid which with a smaller length $L$ (and larger cross-sectional area $A=V/L$) will not "concentrate" the magnetic field as good as a longer and slimmer solenoid. If we let $L^2/A=L^3/V$ grow large (towards the limit of an infinitely long solenoid) the field inside becomes $B=\mu_0 M$. So we have derived that for all volumes $V$
$$B_{max}(M,V)\geq \mu_0 M.$$
I think it can be argued that this inequality is in fact an equality (for all $V$) and the maximal field strength -or should we say: "supremal" field strength?- is obtained by taking these ever longer rod-shapes. I'll edit this answer when I have that question settled.
Remark: this answer of course presupposes a system size for which a continuum mechanics or mean-field approach to this question is appropriate.