<p>Cosmography is the part of cosmology that proceeds by making minimal dynamic assumptions. That is, one does not assume the Friedmann equations (Einstein equations) unless and until absolutely necessary. On the other hand, cosmodynamics is the part of cosmology that relates the geometry to the density and pressure using the Friedmann equations. In both frameworks, we consider the amount of information and the nature of the constraints we can obtain from the Hubble flow in a FLRW universe. Indeed, the cosmological parameters contained in the Hubble relation between distance and redshift provide information on the behaviour of the universe (expansion, acceleration etc...). In the first framework, it is possible to concentrate more directly on the observational situation in a model-independent manner. We perform a number of inter-related cosmographic fits to supernova datasets, and pay particular attention to the extent to which the choice of distance scale and manner of representing the redshift scale affect the cosmological parameters. In the second framework, we use the class of w-parameter models which has become increasingly popular in the last decade. We explore the extent to which a constraint on the w-parameter leads to useful and non-trivial constraints on the Hubble flow in terms of cosmological parameters H(z), density p(z), density parameter O(z), distance scales d(z), and lookback time T(z). On another front, Numerical Relativity has experienced many breakthroughs since 2005, with full inspiral-merger-ringdown simulations now possible. One of the main goals is to provide very accurate templates of gravitational waves for ground-based and space-based interferometers. We explore the potential of a very recent and accurate numerical method, the Spectral Element Method (SEM), for Numerical Relativity, by treating a singular Schwarszchild black hole evolution as a test case. Spectral elements combine the theory of spectral and pseudo-spectral methods for high order polynomials and the variational formulation of finite elements and the associated geometric flexibility. We use the BSSN formulation of the Einstein equations with the method of the moving punctures. After applying the variational formulation to the BSSN system, we present several possible weak forms of this system and its spectral element discretization in space. We use a Runge-Kutta fourth order time discretization. The accuracy of high order methods can deteriorate in the presence of discontinuities or sharp gradients. We show that we can treat the element that contains the puncture with a filtering method to avoid artificial and spurious oscillations. These might form and propagate into the domain coming from discontinuous initial data from the BSSN system.</p>