classy-szfast 0.0.25.post16__py3-none-any.whl → 0.0.25.post18__py3-none-any.whl

classy_szfast/classy_sz.py

@@ -24,6 +24,19 @@ class classy_sz(classy):
                 self.log, "Could not find CLASS_SZ. Check error message above.")
         
         from classy_sz import Class, CosmoSevereError, CosmoComputationError
+
+        # import classy_sz
+        # import os
+        # path_to_data = os.getenv("PATH_TO_CLASS_SZ_DATA")
+        # print(path_to_data)
+
+        # print(classy_sz.__file__)
+
+        # import classy_szfast
+        # print(classy_szfast.__file__)
+        # import sys; sys.exit()
+        # print(sys.path)
+
         
         global CosmoComputationError, CosmoSevereError
 
@@ -255,6 +268,7 @@ class classy_sz(classy):
                         args_names=[],
                         args=[])
 
+        ## added lensing fwith limber integral
         if "cl_lens_lens" in requirements:
                 # make sure cobaya still runs as it does for standard classy
                 requirements.pop("cl_lens_lens")
@@ -446,10 +460,32 @@ class classy_sz(classy):
                 else:
 
                     start_time = time.time()
+                    # self.classy = Class()
+                    # print('pars in classy',self.classy.pars)
                     self.classy.compute_class_szfast(likelihood_mode=True)
                     end_time = time.time()
                     # self.log.info("Execution time of class_szfast: {:.5f} seconds".format(end_time - start_time))
-                    # print('pars in classy',self.classy.pars)
+                    
+                    # some debug printouts
+                    # cl_kk_hm = self.classy.cl_kk
+                    # ell = np.asarray(cl_kk_hm()['ell'])
+                    # fac = ell*(ell+1.)/2./np.pi
+                    # cl_kk_hf = np.asarray(cl_kk_hm()['hf'])/fac
+                    # print('cl_kk_hf',cl_kk_hf[0])
+                    # print('ell',ell[0])
+                    # print('pp',self.classy.lensed_cl()['pp'][2])
+                    # print('angular_distance',self.classy.angular_distance(0.1))
+                    # print('Omega_Lambda',self.classy.Omega_Lambda())
+                    # print('Hubble',self.classy.Hubble(0))
+                    # print('Omega_m',self.classy.Omega_m())
+                    # print('Neff',self.classy.Neff())
+                    # print('m_ncdm_tot',self.classy.get_current_derived_parameters(['m_ncdm_tot']))
+                    # print('Neff',self.classy.get_current_derived_parameters(['Neff']))
+                    # print('Omega_m',self.classy.get_current_derived_parameters(['Omega_m']))
+                    # print('h',self.classy.get_current_derived_parameters(['h']))
+                    # print('m_ncdm_in_eV',self.classy.get_current_derived_parameters(['m_ncdm_in_eV']))
+                    # print("\n\n-----------------------------------\n\n")
+                    # import sys; sys.exit()
             else:
 
                 self.classy.compute()
@@ -546,6 +582,21 @@ class classy_sz(classy):
         state["derived_extra"] = deepcopy(d_extra)
 
 
+    def set(self, params_values_dict):
+        # If no output requested, remove arguments that produce an error
+        # (e.g. complaints if halofit requested but no Cl's computed.) ?????
+        # Needed for facilitating post-processing
+        if not self.extra_args["output"]:
+            for k in ["non_linear"]:
+                self.extra_args.pop(k, None)
+        # Prepare parameters to be passed: this-iteration + extra
+        args = {self.translate_param(p): v for p, v in params_values_dict.items()}
+        args.update(self.extra_args)
+        # Generate and save
+        # print("Setting parameters: %r", args)
+        self.classy.set(**args)
+
+
 
     @classmethod
     def is_installed(cls, **kwargs):

{classy_szfast-0.0.25.post16.dist-info → classy_szfast-0.0.25.post18.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: classy_szfast
-Version: 0.0.25.post16
+Version: 0.0.25.post18
 Summary: The accelerator of the class_sz code from https://github.com/CLASS-SZ
 Maintainer-email: Boris Bolliet <bb667@cam.ac.uk>
 License: MIT

{classy_szfast-0.0.25.post16.dist-info → classy_szfast-0.0.25.post18.dist-info}/RECORD

@@ -1,5 +1,5 @@
 classy_szfast/__init__.py,sha256=E2thrL0Z9oXFfdzwcsu-xbOytudLFTlRlPqVFGlPPPg,279
-classy_szfast/classy_sz.py,sha256=uUQtksZ5GZJTi01bowLaM1EUEXaQad8cMiNK_iXnFpI,22839
+classy_szfast/classy_sz.py,sha256=oy2MzoWsbgzzrr52ocn1knA4Y1u9UlOuslNfyufLSNs,25268
 classy_szfast/classy_szfast.py,sha256=PFb4Mqy3ACTRwEDAgN3m_-2batQhIZ1z2prT8zSGu5U,44556
 classy_szfast/config.py,sha256=v6DGcBHmfn5JtuO48dKyXCh-Dmn0uwOF_izvVOJFnqw,279
 classy_szfast/cosmopower.py,sha256=ooYK2BDOZSo3XtGHfPtjXHxr5UW-yVngLPkb5gpvTx8,2351
@@ -14,7 +14,7 @@ classy_szfast/custom_bias/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJW
 classy_szfast/custom_bias/custom_bias.py,sha256=aR2t5RTIwv7P0m2bsEU0Eq6BTkj4pG10AebH6QpG4qM,486
 classy_szfast/custom_profiles/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 classy_szfast/custom_profiles/custom_profiles.py,sha256=4LZwb2XoqwCyWNmW2s24Z7AJdmgVdaRG7yYaBYe-d9Q,1188
-classy_szfast-0.0.25.post16.dist-info/METADATA,sha256=HK53ITw__TpeiD5IBctJ1PwOUM8HPiCNbVWQvZs40D0,572
-classy_szfast-0.0.25.post16.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
-classy_szfast-0.0.25.post16.dist-info/top_level.txt,sha256=hRgqpilUck4lx2KkaWI2y9aCDKqF6pFfGHfNaoPFxv0,14
-classy_szfast-0.0.25.post16.dist-info/RECORD,,
+classy_szfast-0.0.25.post18.dist-info/METADATA,sha256=fZBmdgjo1wb0mZb5_3f-l2BBJYkOP2uiIvZWVSrEnBE,572
+classy_szfast-0.0.25.post18.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
+classy_szfast-0.0.25.post18.dist-info/top_level.txt,sha256=hRgqpilUck4lx2KkaWI2y9aCDKqF6pFfGHfNaoPFxv0,14
+classy_szfast-0.0.25.post18.dist-info/RECORD,,