“…The stability of the meshless method is then usually enforced via upwind approximation of fluxes similarly to Finite volume methods, cf. [4,6,9,10,23]. This technique allows to solve the systems of nonlinear hyperbolic PDE but, on the other hand, it increases complexity and overall computational costs of the method.…”