raise error if get_sigma8_0 when not calcualted

cmbant · Nov 25, 2020 · 7a85a78 · 7a85a78
1 parent 103443b
commit 7a85a78
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/camb/results.py b/camb/results.py
@@ -845,6 +845,9 @@ def get_sigma8_0(self):
 
         :return: :math:`\sigma_8` today
         """
+        nz = self.Params.Transfer.PK_num_redshifts
+        if not nz or self.Params.Transfer.PK_redshifts[nz - 1] > 1e-5:
+            raise CAMBError("sigma8 requested at z=0, but P(z=0) not calcaulted")
         return self.get_sigma8()[-1]
 
     def get_fsigma8(self):