ler 0.4.1__py3-none-any.whl → 0.4.2__py3-none-any.whl

ler/__init__.py

@@ -51,7 +51,7 @@ set_multiprocessing_start_method()
 __author__ = 'hemanta_ph <hemantaphurailatpam@gmail.com>'
 
 # The version as used in the setup.py
-__version__ = "0.4.1"
+__version__ = "0.4.2"
 
 # add __file__
 import os

ler/gw_source_population/cbc_source_parameter_distribution.py

@@ -46,7 +46,7 @@ class CBCSourceParameterDistribution(CBCSourceRedshiftDistribution):
         Dictionary of prior sampler functions and its input parameters.
         Check for available priors and corresponding input parameters by running,
         >>> from ler.gw_source_population import CBCSourceParameterDistribution
-        >>> cbc = CompactBinaryPopulation()
+        >>> cbc = CBCSourceParameterDistribution()
         >>> cbc.available_gw_prior_list_and_its_params()
         # To check the current chosen priors and its parameters, run,
         >>> print("default priors=",cbc.gw_param_samplers)
@@ -77,7 +77,7 @@ class CBCSourceParameterDistribution(CBCSourceRedshiftDistribution):
 
     Instance Attributes
     ----------
-    CompactBinaryPopulation has the following instance attributes:\n
+    CBCSourceParameterDistribution has the following instance attributes:\n
     +-------------------------------------+----------------------------------+
     | Atrributes                          | Type                             |
     +=====================================+==================================+
@@ -114,7 +114,7 @@ class CBCSourceParameterDistribution(CBCSourceRedshiftDistribution):
 
     Instance Methods
     ----------
-    CompactBinaryPopulation has the following instance methods:\n
+    CBCSourceParameterDistribution has the following instance methods:\n
     +-------------------------------------+----------------------------------+
     | Methods                             | Type                             |
     +=====================================+==================================+
@@ -1136,7 +1136,7 @@ class CBCSourceParameterDistribution(CBCSourceRedshiftDistribution):
         Examples
         ----------
         >>> from ler.gw_source_population import CBCSourceParameterDistribution
-        >>> cbc = CompactBinaryPopulation()
+        >>> cbc = CBCSourceParameterDistribution()
         >>> priors = cbc.available_gw_prior_list_and_its_params
         >>> priors.keys()  # type of priors
         dict_keys(['merger_rate_density', 'source_frame_masses', 'spin', 'geocent_time', 'ra', 'phase', 'psi', 'theta_jn'])

ler/gw_source_population/cbc_source_redshift_distribution.py

@@ -25,7 +25,7 @@ from .jit_functions import merger_rate_density_bbh_popI_II_oguri2018, star_forma
 
 class CBCSourceRedshiftDistribution(object):
     """Class to generate a population of source galaxies.
-    This class is inherited by :class:`~ler.ler.CompactBinaryPopulation` and :class:`~ler.ler.LensGalaxyParameterDistribution` class.
+    This class is inherited by :class:`~ler.ler.CBCSourceParameterDistribution` and :class:`~ler.ler.LensGalaxyParameterDistribution` class.
 
     Parameters
     ----------

ler/image_properties/image_properties.py

@@ -1,7 +1,7 @@
 # -*- coding: utf-8 -*-
 """
 This module contains the LensGalaxyPopulation class, which is used to sample lens galaxy parameters, source parameters conditioned on the source being strongly lensed, image properties, and lensed SNRs. \n
-The class inherits from the CompactBinaryPopulation class, which is used to sample source parameters. \n
+The class inherits from the CBCSourceParameterDistribution class, which is used to sample source parameters. \n
 """
 
 import warnings
@@ -341,6 +341,8 @@ class ImageProperties():
         ####################################
         lens_parameters.update(image_parameters)
         lens_parameters["n_images"] = n_images
+        lens_parameters["x_source"] = x_source
+        lens_parameters["y_source"] = y_source
 
         return lens_parameters
 

ler/lens_galaxy_population/lens_galaxy_parameter_distribution.py

@@ -2,7 +2,7 @@
 """
 This module contains the LensGalaxyPopulation class, which is used to sample lens galaxy parameters, source parameters conditioned on the source being strongly lensed. \n
 The class inherits from the ImageProperties class, which is used calculate image properties (magnification, timedelays, source position, image position, morse phase). \n
-Either the class takes in initialized CompactBinaryPopulation class as input or inherits the CompactBinaryPopulation class with default params (if no input) \n
+Either the class takes in initialized CBCSourceParameterDistribution class as input or inherits the CBCSourceParameterDistribution class with default params (if no input) \n
 """
 
 import warnings
@@ -181,8 +181,8 @@ class LensGalaxyParameterDistribution(CBCSourceParameterDistribution, ImagePrope
 
     # Attributes
     cbc_pop = None
-    """:class:`~CompactBinaryPopulation` class\n
-    This is an already initialized class that contains a function (CompactBinaryPopulation.sample_gw_parameters) that actually samples the source parameters. 
+    """:class:`~CBCSourceParameterDistribution` class\n
+    This is an already initialized class that contains a function (CBCSourceParameterDistribution.sample_gw_parameters) that actually samples the source parameters. 
     """
 
     z_min = None
@@ -313,7 +313,7 @@ class LensGalaxyParameterDistribution(CBCSourceParameterDistribution, ImagePrope
             dictionary of parameters to initialize the parent classes
         """
 
-        # initialization of CompactBinaryPopulation class
+        # initialization of CBCSourceParameterDistribution class
         # it also initializes the CBCSourceRedshiftDistribution class
         # list of relevant initialized instances,
         # 1. self.sample_source_redshift
@@ -639,7 +639,7 @@ class LensGalaxyParameterDistribution(CBCSourceParameterDistribution, ImagePrope
         def zs_function(zs_sl):
             # get zs
             # self.sample_source_redshifts from CBCSourceRedshiftDistribution class
-            zs = self.sample_zs(size)  # this function is from CompactBinaryPopulation class
+            zs = self.sample_zs(size)  # this function is from CBCSourceParameterDistribution class
             # put strong lensing condition with optical depth
             tau = self.strong_lensing_optical_depth(zs)
             tau_max = self.strong_lensing_optical_depth(np.array([z_max]))[0] # tau increases with z

ler/rates/gwrates.py

@@ -84,7 +84,7 @@ class GWRATES(CBCSourceParameterDistribution):
         
     Instance Attributes
     ----------
-    LeR class has the following attributes, \n
+    LeR class has the following attributes:\n
     +-------------------------------------+----------------------------------+
     | Atrributes                          | Type                             |
     +=====================================+==================================+
@@ -121,7 +121,7 @@ class GWRATES(CBCSourceParameterDistribution):
 
     Instance Methods
     ----------
-    LeR class has the following methods, \n
+    LeR class has the following methods:\n
     +-------------------------------------+----------------------------------+
     | Methods                             | Description                      |
     +=====================================+==================================+
@@ -506,7 +506,7 @@ class GWRATES(CBCSourceParameterDistribution):
             dictionary of parameters to initialize the parent classes
         """
 
-        # initialization of CompactBinaryPopulation class
+        # initialization of CBCSourceParameterDistribution class
         # it also initializes the CBCSourceRedshiftDistribution class
         input_params = dict(
             z_min=self.z_min,

ler/rates/ler.py

@@ -76,6 +76,23 @@ class LeR(LensGalaxyParameterDistribution):
     interpolator_directory : `str`
         directory to store the interpolators.
         default interpolator_directory = './interpolator_pickle'. This is used for storing the various interpolators related to `ler` and `gwsnr` package.
+    create_new_interpolator : `bool` or `dict`
+        default create_new_interpolator = False.
+        if True, the all interpolators (including `gwsnr`'s)will be created again.
+        if False, the interpolators will be loaded from the interpolator_directory if they exist.
+        if dict, you can specify which interpolators to create new. Complete example (change any of them to True), create_new_interpolator = create_new_interpolator = dict(
+            redshift_distribution=dict(create_new=False, resolution=1000),
+            z_to_luminosity_distance=dict(create_new=False, resolution=1000),
+            velocity_dispersion=dict(create_new=False, resolution=1000),
+            axis_ratio=dict(create_new=False, resolution=1000),
+            optical_depth=dict(create_new=False, resolution=200),
+            z_to_Dc=dict(create_new=False, resolution=1000),
+            Dc_to_z=dict(create_new=False, resolution=1000),
+            angular_diameter_distance=dict(create_new=False, resolution=1000),
+            differential_comoving_volume=dict(create_new=False, resolution=1000),
+            Dl_to_z=dict(create_new=False, resolution=1000),
+            gwsnr=False,
+        )
     ler_directory : `str`
         directory to store the parameters.
         default ler_directory = './ler_data'. This is used for storing the parameters of the simulated events.
@@ -98,7 +115,7 @@ class LeR(LensGalaxyParameterDistribution):
         
     Instance Attributes
     ----------
-    LeR class has the following attributes, \n
+    LeR class has the following attributes: \n
     +-------------------------------------+----------------------------------+
     | Atrributes                          | Type                             |
     +=====================================+==================================+
@@ -141,7 +158,7 @@ class LeR(LensGalaxyParameterDistribution):
 
     Instance Methods
     ----------
-    LeR class has the following methods, \n
+    LeR class has the following methods:\n
     +-------------------------------------+----------------------------------+
     | Methods                             | Description                      |
     +=====================================+==================================+
@@ -194,7 +211,7 @@ class LeR(LensGalaxyParameterDistribution):
     |                                     | ratio between lensed and         |
     |                                     | unlensed events.                 |
     +-------------------------------------+----------------------------------+
-    |:meth:`~rate_comparision_with_rate_calculation                          |
+    |:meth:`~rate_comparison_with_rate_calculation`                          |
     +-------------------------------------+----------------------------------+
     |                                     | Function to calculate rates for  |
     |                                     | unleesed and lensed events and   |
@@ -437,6 +454,7 @@ class LeR(LensGalaxyParameterDistribution):
         list_of_detectors=None,
         json_file_names=None,
         interpolator_directory="./interpolator_pickle",
+        create_new_interpolator=False,
         ler_directory="./ler_data",
         verbose=True,
         **kwargs,
@@ -453,6 +471,7 @@ class LeR(LensGalaxyParameterDistribution):
         if json_file_names:
             self.json_file_names.update(json_file_names)
         self.interpolator_directory = interpolator_directory
+        kwargs["create_new_interpolator"] = create_new_interpolator
         self.ler_directory = ler_directory
         # create directory if not exists
         if not os.path.exists(ler_directory):
@@ -820,6 +839,21 @@ class LeR(LensGalaxyParameterDistribution):
                 if key in input_params:
                     input_params[key] = value
         self.snr_calculator_dict = input_params
+
+        # dealing with create_new_interpolator param
+        if isinstance(input_params["create_new_interpolator"], bool):
+            pass
+        elif isinstance(input_params["create_new_interpolator"], dict):
+            # check input_params["gwsnr"] exists
+            if "gwsnr" in input_params["create_new_interpolator"]:
+                if isinstance(input_params["create_new_interpolator"]["gwsnr"], bool):
+                    input_params["create_new_interpolator"] = input_params["create_new_interpolator"]["gwsnr"]
+                else:
+                    raise ValueError("create_new_interpolator['gwsnr'] should be a boolean.")
+            else:
+                raise ValueError("create_new_interpolator should be a boolean or a dictionary with 'gwsnr' key.")
+
+        # initialization of GWSNR class
         gwsnr = GWSNR(
                     npool=input_params["npool"],
                     mtot_min=input_params["mtot_min"],
@@ -835,7 +869,7 @@ class LeR(LensGalaxyParameterDistribution):
                     psds=input_params["psds"],
                     ifos=input_params["ifos"],
                     interpolator_dir=input_params["interpolator_dir"],
-                    # create_new_interpolator=input_params["create_new_interpolator"],
+                    create_new_interpolator=input_params["create_new_interpolator"],
                     gwsnr_verbose=input_params["gwsnr_verbose"],
                     multiprocessing_verbose=input_params["multiprocessing_verbose"],
                     mtot_cut=input_params["mtot_cut"],
@@ -910,7 +944,7 @@ class LeR(LensGalaxyParameterDistribution):
             resume = False (default) or True.
             if True, the function will resume from the last batch.
         save_batch : `bool`
-            if True, the function will save the parameters in batches. if False, the function will save all the parameters at the end of sampling. save_batch=False is faster.
+            if True, the function will save the parameters in batches. if False (default), the function will save all the parameters at the end of sampling. save_batch=False is faster.
         output_jsonfile : `str`
             json file name for storing the parameters.
             default output_jsonfile = 'unlensed_params.json'. Note that this file will be stored in the self.ler_directory.
@@ -1670,7 +1704,7 @@ class LeR(LensGalaxyParameterDistribution):
 
         return snr_hit
 
-    def rate_comparision_with_rate_calculation(
+    def rate_comparison_with_rate_calculation(
         self,
         unlensed_param=None,
         snr_threshold_unlensed=8.0,
@@ -1733,7 +1767,7 @@ class LeR(LensGalaxyParameterDistribution):
         >>> ler = LeR()
         >>> ler.unlensed_cbc_statistics();
         >>> ler.lensed_cbc_statistics();
-        >>> rate_ratio, unlensed_param, lensed_param = ler.rate_comparision_with_rate_calculation()
+        >>> rate_ratio, unlensed_param, lensed_param = ler.rate_comparison_with_rate_calculation()
         """
 
         # call json_file_ler_param and add the results

ler/utils/utils.py

@@ -5,6 +5,7 @@ This module contains helper routines for other modules in the ler package.
 
 import os
 import pickle
+import h5py
 import numpy as np
 import json
 from scipy.interpolate import interp1d
@@ -52,6 +53,66 @@ class NumpyEncoder(json.JSONEncoder):
             return obj.tolist()
         return json.JSONEncoder.default(self, obj)
 
+def load_pickle(file_name):
+    """Load a pickle file.
+
+    Parameters
+    ----------
+    file_name : `str`
+        pickle file name for storing the parameters.
+
+    Returns
+    ----------
+    param : `dict`
+    """
+    with open(file_name, "rb") as handle:
+        param = pickle.load(handle)
+
+    return param
+
+def save_pickle(file_name, param):
+    """Save a dictionary as a pickle file.
+
+    Parameters
+    ----------
+    file_name : `str`
+        pickle file name for storing the parameters.
+    param : `dict`
+        dictionary to be saved as a pickle file.
+    """
+    with open(file_name, "wb") as handle:
+        pickle.dump(param, handle, protocol=pickle.HIGHEST_PROTOCOL)
+
+# hdf5
+def load_hdf5(file_name):
+    """Load a hdf5 file.
+
+    Parameters
+    ----------
+    file_name : `str`
+        hdf5 file name for storing the parameters.
+
+    Returns
+    ----------
+    param : `dict`
+    """
+
+    return h5py.File(file_name, 'r')
+
+def save_hdf5(file_name, param):
+    """Save a dictionary as a hdf5 file.
+
+    Parameters
+    ----------
+    file_name : `str`
+        hdf5 file name for storing the parameters.
+    param : `dict`
+        dictionary to be saved as a hdf5 file.
+    """
+    with h5py.File(file_name, 'w') as f:
+        for key, value in param.items():
+            f.create_dataset(key, data=value)
+
 def load_json(file_name):
     """Load a json file.
 

{ler-0.4.1.dist-info → ler-0.4.2.dist-info}/METADATA

@@ -1,7 +1,7 @@
 Metadata-Version: 2.1
 Name: ler
-Version: 0.4.1
-Summary: Gravitational waves Lensing Rates
+Version: 0.4.2
+Summary: LVK (LIGO-Virgo-KAGRA collaboration) Event (compact-binary mergers) Rate calculator and simulator
 Home-page: https://github.com/hemantaph/ler
 Author: Hemantakumar
 Author-email: hemantaphurailatpam@gmail.com
@@ -9,20 +9,20 @@ License: MIT
 Requires-Python: >=3.10
 Description-Content-Type: text/markdown
 License-File: LICENSE
-Requires-Dist: setuptools>=65.5.0
-Requires-Dist: matplotlib>=3.4.2
-Requires-Dist: pycbc>=2.0.4
-Requires-Dist: numpy>=1.18
-Requires-Dist: numba>=0.57.1
-Requires-Dist: bilby>=1.0.2
-Requires-Dist: gwsnr>=0.3.2
-Requires-Dist: scipy<1.14.0
-Requires-Dist: lenstronomy>=1.10.4
-Requires-Dist: astropy>=5.1
-Requires-Dist: tqdm>=4.64.1
-Requires-Dist: pointpats>=2.3
-Requires-Dist: shapely>=2.0.1
-Requires-Dist: gwcosmo==2.1.0
+Requires-Dist: setuptools >=65.5.0
+Requires-Dist: matplotlib >=3.4.2
+Requires-Dist: pycbc >=2.0.4
+Requires-Dist: numpy >=1.18
+Requires-Dist: numba >=0.57.1
+Requires-Dist: bilby >=1.0.2
+Requires-Dist: gwsnr >=0.3.2
+Requires-Dist: scipy <1.14.0
+Requires-Dist: lenstronomy >=1.10.4
+Requires-Dist: astropy >=5.1
+Requires-Dist: tqdm >=4.64.1
+Requires-Dist: pointpats >=2.3
+Requires-Dist: shapely >=2.0.1
+Requires-Dist: gwcosmo ==2.1.0
 
 # LeR
 [![DOI](https://zenodo.org/badge/626733473.svg)](https://zenodo.org/badge/latestdoi/626733473) [![PyPI version](https://badge.fury.io/py/ler.svg)](https://badge.fury.io/py/ler) [![DOCS](https://readthedocs.org/projects/ler/badge/?version=latest)](https://ler.readthedocs.io/en/latest/)

{ler-0.4.1.dist-info → ler-0.4.2.dist-info}/RECORD

@@ -1,25 +1,25 @@
-ler/__init__.py,sha256=aC-PEncn43TwcmsVQdUYivq0I8ylmw8wKclBoEWVQv4,2174
+ler/__init__.py,sha256=KUIsDaNuHiBdV2gljkp4QbqUCwyMxx5gdZk3r7mVrK8,2174
 ler/gw_source_population/__init__.py,sha256=HG0ve5wTpBDN2fNMxHLnoOqTz-S0jXM_DsWEJ5PEHAw,126
-ler/gw_source_population/cbc_source_parameter_distribution.py,sha256=6rKmyX6DLawgXsgbrjLHN_Mgd1j5QcaMHFHEIHVLHh0,67457
-ler/gw_source_population/cbc_source_redshift_distribution.py,sha256=o2qAM_-9SeLxxfGwqXrdVWTCeEAaXVan_OPDd4jrplg,28559
+ler/gw_source_population/cbc_source_parameter_distribution.py,sha256=A6tEvCMB7sqEQPlEBR6g7nqQMN8uPkiWiqT7LQLocfw,67485
+ler/gw_source_population/cbc_source_redshift_distribution.py,sha256=F2MBdYyXxh0kr5No778wNc4esBs7MES3ej2FrLaZD0M,28566
 ler/gw_source_population/jit_functions.py,sha256=aQV9mv3IY5b3OLiPeXmoLWJ_TbFUS9M1OgnIyIY3eX4,8668
 ler/image_properties/__init__.py,sha256=XfJFlyZuOrKODT-z9WxjR9mI8eT399YJV-jzcJKTqGo,71
-ler/image_properties/image_properties.py,sha256=A6PDoeK28C1zKLgOjLAfPqo8CFROaHK4OdbvI34j_hM,25376
+ler/image_properties/image_properties.py,sha256=CI-0_Jj4u5vf9KhZcbA9Nmc1naxEfbFYCMuRjUWi-oA,25477
 ler/image_properties/multiprocessing_routine.py,sha256=hYnQTM7PSIj3X-5YNDqMxH9UgeXHUPPdLG70h_r6sEY,18333
 ler/lens_galaxy_population/__init__.py,sha256=TXk1nwiYy0tvTpKs35aYK0-ZK63g2JLPyGG_yfxD0YU,126
 ler/lens_galaxy_population/jit_functions.py,sha256=tCTcr4FWyQXH7SQlHsUWeZBpv4jnG00DsBIljdWFs5M,8472
-ler/lens_galaxy_population/lens_galaxy_parameter_distribution.py,sha256=8dj1hoef_cjkb7tbXwR2pyUl_Tyyel_Ze3OAJ4DTUBw,48335
+ler/lens_galaxy_population/lens_galaxy_parameter_distribution.py,sha256=tRNy4DtSUPGfm5AdmQ4tbD-hHDwQ_HeOgq_mkDN5qe8,48377
 ler/lens_galaxy_population/mp.py,sha256=TPnFDEzojEqJzE3b0g39emZasHeeaeXN2q7JtMcgihk,6387
 ler/lens_galaxy_population/optical_depth.py,sha256=rZ_Inpw7ChpFdDLp3kJrCmA0PL3RxN6T_W_NTFhj_ko,42542
 ler/rates/__init__.py,sha256=N4li9NouSVjZl5HIhyuiKKRyrpUgQkBZaUeDgL1m4ic,43
-ler/rates/gwrates.py,sha256=akw6rKAkETr_ERmymVJx3APRXs0XqqFccZ-LIzAV4jM,58465
+ler/rates/gwrates.py,sha256=Ek_wqAEFV1Uffu5TtiihdLJuSaEqw2KJhNpNFUWdjSw,58470
 ler/rates/ler copy.py,sha256=BlnGlRISUwiWUhUNwp32_lvh7tHdT-d1VDhFelwKO_c,101873
-ler/rates/ler.py,sha256=LlpUJKw5KkPacz_mm5voLDLqLkmr8FD4bQ68jcQeSqY,105771
+ler/rates/ler.py,sha256=38z4EGRs33-B6AUaIVAbcWXYRBZW0l0EixyZ1ngrsH4,107884
 ler/utils/__init__.py,sha256=JWF9SKoqj1BThpV_ynfoyUeU06NQQ45DHCUGaaMSp_8,42
 ler/utils/plots.py,sha256=uq-usKRnEymtOSAPeHFOfMQW1XX76_WP2aBkT40RvLo,15664
-ler/utils/utils.py,sha256=45fdEJ6B125kCp_lXl3IFulY61hKyOX1JewSQ_LLQSU,31464
-ler-0.4.1.dist-info/LICENSE,sha256=9LeXXC3WaHBpiUGhLVgOVnz0F12olPma1RX5zgpfp8Q,1081
-ler-0.4.1.dist-info/METADATA,sha256=nmop_EQhbQPN2M6_ISxap_6LUvSA1i2BFkA0T5ob4Sk,6514
-ler-0.4.1.dist-info/WHEEL,sha256=GV9aMThwP_4oNCtvEC2ec3qUYutgWeAzklro_0m4WJQ,91
-ler-0.4.1.dist-info/top_level.txt,sha256=VWeWLF_gNMjzquGmqrLXqp2J5WegY86apTUimMTh68I,4
-ler-0.4.1.dist-info/RECORD,,
+ler/utils/utils.py,sha256=wxcHWHNdjdU9fi4RftDez4PuIeu5m7OCL-lS-L6j5b8,32798
+ler-0.4.2.dist-info/LICENSE,sha256=9LeXXC3WaHBpiUGhLVgOVnz0F12olPma1RX5zgpfp8Q,1081
+ler-0.4.2.dist-info/METADATA,sha256=FL3i3pflIZDyQ4ZpXT945zepKway86BfCgh8DUogG8g,6592
+ler-0.4.2.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
+ler-0.4.2.dist-info/top_level.txt,sha256=VWeWLF_gNMjzquGmqrLXqp2J5WegY86apTUimMTh68I,4
+ler-0.4.2.dist-info/RECORD,,

{ler-0.4.1.dist-info → ler-0.4.2.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (75.1.0)
+Generator: bdist_wheel (0.43.0)
 Root-Is-Purelib: true
 Tag: py3-none-any