Abstract We present a semi-Lagrangian method for the numerical resolution of Vlasov-type equations on multi-patch meshes. Following N. Crouseilles et al. [ A parallel Vlasov solver based local cubic spline interpolation patches. Journal Computational Physics (2009)], we employ with Hermite boundary conditions between The derivative reconstruction is adapted to cope non-uniform meshes as well non-conforming situations. In conforming case, constraint number points each patch, found in previous studies, removed; however, small global system must now be solved. that representations coincide corresponding reconstruction. Alternatively, can choose not apply and derivatives approximated. influence most distant diminishes per patch increases. For uniform configurations, study explicit asymptotic behavior this has been led. validated using two-dimensional guiding-center model an O-point. All results are carried out Gyselalib++ library.