We present a comparative analysis of the continuous Galerkin (CG) and discontinuous Galerkin (DG) finite element (FE) methods, focusing on their efficiency in solving advection–diffusion problems in the Summation By Parts–Simultaneous Approximation Term (SBP–SAT) framework. Our findings provide insights into the practical advantages and limitations of each approach, guiding future applications in high–order numerical simulations of real world applications.