Carolina Cuesta
Carolina Cuesta
Are there any news on this thread? It'd be a useful thing to have!
That one is closer 
It seems to be coming from the ensure_physical flag  But then it doesn't work if you sample from the non physical field, I get this error message, `ValueError: lam...
Yes, I tried that too but it is still problematic  This is the code I ran, ```python import matplotlib.pyplot as plt import numpy as np from astropy.cosmology import Planck18...
This is for 500 samples  Std is large
It seems to be a matter of amplitude. When P(k) = 1/k, the variance is low,  But when P(k) = 5.e2/k (to match the amplitude of halofit's power spectrum),...