v6.0: release MagIC 6.0

Support parallel tridiagonal and pentadiagonal solvers when finite differences are employed. This feature is only available for 2nd order F.D., it then makes uses of ghost zones to ensure boundary conditions. In this configuration, the all-to-all communications are not used except for calling the output routines. This yields an overall 20-40% speed-up depending on the spatial resolution.