This paper presents the use of Powell-Sabin splines in the context of isogeometric analysis for the numerical solution of advectiondiffusion-reaction equations. Powell-Sabin splines are piecewise quadratic C functions defined on a given triangulation with a particular macro-structure. We discuss the Galerkin discretization based on a normalized Powell-Sabin B-spline basis. We focus on the accur...