Numerically obtaining the inverse of a function is a common task for many scientific problems, often solved using a Newton iteration method. Here we describe an alternative scheme, based on switching variables followed by spline interpolation, which can be applied to monotonic functions under very general conditions. To optimize the algorithm, we designed a specific ultra-fast spline routine. We also derive analytically the theoretical errors of the method and test it on examples that are of interest in physics. In particular, we compute the real branch of Lambert's $W(y)$ function, which is defined as the inverse of $x \exp(x)$, and we solve Kepler's equation. In all cases, our predictions for the theoretical errors are in excellent agreement with our numerical results, and are smaller than what could be expected from the general error analysis of spline interpolation by many orders of magnitude, namely by an astonishing $3\times 10^{-22}$ factor for the computation of $W$ in the range $W(y)\in [0,10]$, and by a factor $2\times 10^{-4}$ for Kepler's problem. In our tests, this scheme is much faster than Newton-Raphson method, by a factor in the range $10^{-4}$ to $10^{-3}$ for the execution time in the examples, when the values of the inverse function over an entire interval or for a large number of points are requested. For Kepler's equation and tolerance $10^{-6}$ rad, the algorithm outperforms Newton's method for all values of the number of points $N\ge 2$.

We show that the dark energy effects can be modeled by using an \emph{Ising perfect fluid} with network interactions, whose low redshift equation of state, i.e. $\omega_0$, becomes $\omega_0=-1$ as in the $\Lambda$CDM model. In our picture, dark energy is characterized by a barotropic fluid on a lattice in the equilibrium configuration. Thus, mimicking the spin interaction by replacing the spin variable with an occupational number, the pressure naturally becomes negative. We find that the corresponding equation of state mimics the effects of a variable dark energy term, whose limiting case reduces to the cosmological constant $\Lambda$. This permits us to avoid the introduction of a vacuum energy as dark energy source by hand, alleviating the coincidence and fine tuning problems. We find fairly good cosmological constraints, by performing three tests with supernovae Ia, baryonic acoustic oscillation and cosmic microwave background measurements. Finally, we perform the AIC and BIC selection criteria, showing that our model is statistically favored with respect to the CPL parametrization.