singular matrices should not be allowed

Author: Vokturz

cosmoabc/ABC_sampler.py

@@ -415,7 +415,7 @@ def UpdateWeights(self, W, previous_particle_system,
             denominator = sum(W[i3]*multivariate_normal.pdf(
                               current_particle_system[i4][:indx], 
                               previous_particle_system[i3][:indx],
-                              cov=cov1, allow_singular=True) 
+                              cov=cov1, allow_singular=False) 
                               for i3 in range(len(W)))
 
             new_weights.append(nominator/denominator)