Provides a c++ code to do a fast MCMC search with 4 params quickly. Performs trilinear interpolation saves the chains and check for convergence. However, use has to provide the4d grid and grid of models.