We know that, in general, the operator norm of a matrix is larger than its spectral radius: $\|A\|_{\rm op}\ge r(A)\equiv \max_{\lambda\in\sigma(A)}\lvert \lambda\rvert$, where $\sigma(A)$ denotes the spectrum of $A$.
If $A$ is normal, these quantities are clearly the same, as $\|A\|_{\rm op}$ equals the largest singular value of $A$, which equals the largest $|\lambda|$ for $\lambda$ an eigenvalue of $A$.
I can also find examples of $\|A\|_{\rm op}=r(A)$ for $A$ not normal, by building it as direct sum of a non-normal matrix and a matrix with a large eigenvalue. E.g. $$\begin{pmatrix}2 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 1 & 1\end{pmatrix}.$$ This is not diagonalisable, but satisfies the equality.
Are there such examples with simple matrices? In other words, given a simple $A$ (i.e. a matrix composed of a single Jordan block), can we have $\|A\|_{\rm op}=r(A)$?