Solve linear system of equation with gauss elimination
This code is taken from [Asinari code MIXLBM.f90]
(http://staff.polito.it/pietro.asinari/rome08/index.html)
This routine serves as a wrapper to call other routines responsible for
dynamic load balancing in Musubi. The dynamic load balancing algorithm
follows the SPARTA algorithm in an iterative way to perform balancing for
Multilevel simulations as well as complex boundary conditions.