Dark energy or apparent acceleration due to a relativistic cosmological model more complex than the Friedmann-Lemaitre-Robertson-Walker model?

We use the Szekeres inhomogeneous relativistic models in order to fit supernovacombined data sets. We show that with a choice of the spatial curvature function that is guided by current observations, the models fit the supernovadata almost as well as the Lambda-CDM modelwithout requiring a dark energy component. The Szekeres models were originally derived as an exact solution to Einstein's equations with a general metricthat has no symmetries, and they are regarded as good candidates to model the true lumpy universe that we observe. The null geodesics in these models are not radial. The best fit modelfound is also consistent with the requirement of spatial flatness at cosmic microwave backgroundscales. The first results presented here seem to encourage further investigations of apparent acceleration using various inhomogeneous models, and other constraints from cosmic microwave backgroundand large structure need to be explored next.