Adding MPICH support.

Descreasing allowed m2lnL variance for exploration phase (leading to overflow during matrix inversion for large dimensions).

meson.build

@@ -372,6 +372,16 @@ mpi_fortran_dep = dependency(
   required: false,
 )
 
+# MPICH is not found by the mpi dependency, so we need to check for it manually
+if not mpi_c_dep.found()
+  mpi_c_dep = dependency('mpich', required: false)
+endif
+
+# MPICH is not found by the mpi dependency, so we need to check for it manually
+if not mpi_fortran_dep.found()
+  mpi_fortran_dep = dependency('mpichf90', required: false)
+endif
+
 if mpi_c_dep.found()
   numcosmo_conf.set('HAVE_MPI', true)
   numcosmo_build_cfg.set('have_mpi_support', '#define NUMCOSMO_HAVE_MPI 1')

numcosmo/math/ncm_stats_dist.c

@@ -850,7 +850,7 @@ _ncm_stats_dist_prepare_interp (NcmStatsDist *sd, NcmVector *m2lnp)
   {
     NcmStatsDistClass *sd_class = NCM_STATS_DIST_GET_CLASS (sd);
     NcmStatsDistEval eval       = {sd, self, sd_class, NULL, m2lnp};
-    const gdouble dbl_limit     = 6.0;
+    const gdouble dbl_limit     = 2.0;
     guint i;
 
     /*