Résumé : Larger multi-qubit quantum gates allow shallower, more efficient quantum circuits, which could decrease the prohibitive effect of noise on algorithms for noisy intermediate-scale quantum (NISQ) devices and fault-tolerant error correction schemes. Such multi-qubit gates can potentially be generated by time-independent Hamiltonians comprising only physical (one- and two-local) interaction terms. Here, we present an algorithm that finds the strengths of the Hamiltonian terms by using the direction of the geodesic to the target quantum gate on the Riemannian manifold of $\mathrm{SU}(2^n)$ for $n$ qubits. Differential programming is used to determine how the Hamiltonian terms should be updated in order to follow the geodesic to the target unitary as closely as possible. We numerically compare our geodesic algorithm to gradient descent methods and show that it finds solutions with considerably fewer steps for standard multi-qubit gates such as Toffoli and Fredkin. The geodesic algorithm is then used to find previously unavailable multi-qubit gates implementing high fidelity parity checks, which could be used in a wide array of quantum codes and increase the clock speed of fault-tolerant quantum computers.