The Finite Element Method is a powerful tool for approximating the solutions of a large class of PDE's. Traditional FEM requires the function spaces to be linear. We discuss an extension of FEM, Geodesic Finite Elements, which allows the functions being approximated to have range in a non-linear Riemannian Manifold. We realize the set of pseudo Riemannian metrics as a symmetric space in order to pull back the matrix exponential onto it, thus endowing it with a Riemannian metric. Once done, we are able to discretize a class of pseudo-Riemannian manifolds. This has possible applications in non-linear hyperbolic PDE's.