Approximation of solutions of Hamilton-Jacobi equations on the Heisenberg group

We propose and analyze numerical schemes for viscosity solutions of time-dependent Hamilton-Jacobi equations on the Heisenberg group. The main idea is to construct a grid compatible with the noncommutative group geometry. Under suitable assumptions on the data, the Hamiltonian and the parameters for the discrete first order scheme, we prove that the error between the viscosity solution computed at the grid nodes and the solution of the discrete problem behaves like √h where is h is the mesh step. Such an estimate is similar to those available in the Euclidean geometrical setting. The theoretical results are tested numerically on some examples for which semi-analytical formulas for the computation of geodesics are known. Other simulations are presented, for both steady and unsteady problems.