A bstract The inclusive J/ ψ elliptic ( v 2 ) and triangular ( v 3 ) flow coefficients measured at forward rapidity (2 . 5 < y < 4) and the v 2 measured at midrapidity (| y | < 0 . 9) in Pb-Pb collisions at s NN = 5 . 02 TeV using the ALICE detector at the LHC are reported. The entire Pb-Pb data sample collected during Run 2 is employed, amounting to an integrated luminosity of 750 μ b − 1 at forward rapidity and 93 μ b − 1 at midrapidity. The results are obtained using the scalar product method and are reported as a function of transverse momentum p T and collision centrality. At midrapidity, the J/ ψ v 2 is in agreement with the forward rapidity measurement. The centrality averaged results indicate a positive J/ ψ v 3 with a significance of more than 5 σ at forward rapidity in the p T range 2 < p T < 5 GeV/ c . The forward rapidity v 2 , v 3 , and v 3 /v 2 results at low and intermediate p T ( p T ≲ 8 GeV/ c ) exhibit a mass hierarchy when compared to pions and D mesons, while converging into a species-independent curve at higher p T . At low and intermediate p T , the results could be interpreted in terms of a later thermalization of charm quarks compared to light quarks, while at high p T , path-length dependent effects seem to dominate. The J/ ψ v 2 measurements are further compared to a microscopic transport model calculation. Using a simplified extension of the quark scaling approach involving both light and charm quark flow components, it is shown that the D-meson v n measurements can be described based on those for charged pions and J/ ψ flow.