Global existence of classical solutions to the mixed initial–boundary value problem for quasilinear hyperbolic systems of diagonal form with large BV data