Error raised when run `pygwb_pe`
When I run this commend pygwb_pe --path_to_file H1L1_1361652466-1361653426.pickle
. The output message is normal.
12:07 bilby INFO : Running for label 'hlv', output will be saved to './PE_Output'
12:08 bilby INFO : Analysis priors:
12:08 bilby INFO : omega_ref=LogUniform(minimum=1e-11, maximum=1e-08, name='$\\Omega_{ref}$', latex_label='$\\Omega_{ref}$', unit=None, boundary=None)
12:08 bilby INFO : alpha=Gaussian(mu=0.6666666666666666, sigma=1.5, name=None, latex_label='$\\alpha$', unit=None, boundary=None)
12:08 bilby INFO : Analysis likelihood class: <class 'pygwb.pe.PowerLawModel'>
12:08 bilby INFO : Analysis likelihood noise evidence: (55356.15206190402-331.9310295920271j)
12:08 bilby INFO : Single likelihood evaluation took 7.257e-03 s
12:08 bilby INFO : Checkpoint every check_point_delta_t = 600s
12:08 bilby WARNING : Using cached result
12:08 bilby INFO : Summary of results:
nsamples: 4047
ln_noise_evidence: 40376.692
ln_evidence: 40376.694 +/- 0.015
ln_bayes_factor: 0.002 +/- 0.015
[0.05, 0.95]
However, I got an error when I type pygwb_pe --path_to_file H1L1_1361652466-1361653426.pickle --output_dir ./PE_Output_test
12:09 bilby INFO : Running for label 'hlv', output will be saved to './PE_Output_test'
12:09 bilby INFO : Analysis priors:
12:09 bilby INFO : omega_ref=LogUniform(minimum=1e-11, maximum=1e-08, name='$\\Omega_{ref}$', latex_label='$\\Omega_{ref}$', unit=None, boundary=None)
12:09 bilby INFO : alpha=Gaussian(mu=0.6666666666666666, sigma=1.5, name=None, latex_label='$\\alpha$', unit=None, boundary=None)
12:09 bilby INFO : Analysis likelihood class: <class 'pygwb.pe.PowerLawModel'>
12:09 bilby INFO : Analysis likelihood noise evidence: (55356.15206190402-331.9310295920271j)
12:09 bilby INFO : Single likelihood evaluation took 7.103e-03 s
12:09 bilby INFO : Using sampler Dynesty with kwargs {'nlive': 1000, 'bound': 'multi', 'sample': 'rwalk', 'periodic': None, 'reflective': None, 'update_interval': None, 'first_update': None, 'npdim': None, 'rstate': None, 'queue_size': 1, 'pool': None, 'use_pool': None, 'live_points': None, 'logl_args': None, 'logl_kwargs': None, 'ptform_args': None, 'ptform_kwargs': None, 'gradient': None, 'grad_args': None, 'grad_kwargs': None, 'compute_jac': False, 'enlarge': None, 'bootstrap': None, 'walks': 10, 'facc': 0.2, 'slices': None, 'fmove': 0.9, 'max_move': 100, 'update_func': None, 'ncdim': None, 'blob': False, 'save_history': False, 'history_filename': None, 'maxiter': None, 'maxcall': None, 'dlogz': 0.1, 'logl_max': inf, 'n_effective': None, 'add_live': True, 'print_progress': True, 'print_func': <bound method Dynesty._print_func of <bilby.core.sampler.dynesty.Dynesty object at 0x7f08b8023520>>, 'save_bounds': False, 'checkpoint_file': None, 'checkpoint_every': 60, 'resume': False}
12:09 bilby INFO : Checkpoint every check_point_delta_t = 600s
12:09 bilby INFO : Using dynesty version 2.1.0
12:09 bilby INFO : Using the bilby-implemented rwalk sample method with ACT estimated walks
12:09 bilby INFO : Resume file ./PE_Output_test/hlv_resume.pickle does not exist.
12:09 bilby INFO : Generating initial points from the prior
1it [00:00, ?it/s]Traceback (most recent call last):
File "/home/hsiang-yu.huang/.conda/envs/pygwb/bin/pygwb_pe", line 190, in <module>
main()
File "/home/hsiang-yu.huang/.conda/envs/pygwb/bin/pygwb_pe", line 176, in main
hlv=bilby.run_sampler(likelihood=model,priors=priors,sampler='dynesty', npoints=1000, walks=10, maxmcmc=10000,
File "/home/hsiang-yu.huang/.conda/envs/pygwb/lib/python3.9/site-packages/bilby/core/sampler/__init__.py", line 234, in run_sampler
result = sampler.run_sampler()
File "/home/hsiang-yu.huang/.conda/envs/pygwb/lib/python3.9/site-packages/bilby/core/sampler/base_sampler.py", line 97, in wrapped
output = method(self, *args, **kwargs)
File "/home/hsiang-yu.huang/.conda/envs/pygwb/lib/python3.9/site-packages/bilby/core/sampler/dynesty.py", line 417, in run_sampler
out = self._run_external_sampler_with_checkpointing()
File "/home/hsiang-yu.huang/.conda/envs/pygwb/lib/python3.9/site-packages/bilby/core/sampler/dynesty.py", line 493, in _run_external_sampler_with_checkpointing
self.sampler.run_nested(**sampler_kwargs)
File "/home/hsiang-yu.huang/.conda/envs/pygwb/lib/python3.9/site-packages/dynesty/sampler.py", line 1016, in run_nested
for it, results in enumerate(
File "/home/hsiang-yu.huang/.conda/envs/pygwb/lib/python3.9/site-packages/dynesty/sampler.py", line 744, in sample
delta_logz = np.logaddexp(0,
TypeError: ufunc 'logaddexp' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''
1it [00:00, ?it/s]
The situation using .npz
file is quite normal, but the program shows warning
home/hsiang-yu.huang/.conda/envs/pygwb/lib/python3.9/site-packages/pygwb/baseline.py:257: UserWarning: Neither baseline nor interferometer duration is set.
warnings.warn("Neither baseline nor interferometer duration is set.")
/home/hsiang-yu.huang/.conda/envs/pygwb/lib/python3.9/site-packages/pygwb/baseline.py:454: UserWarning: Neither baseline nor interferometer sampling_frequency is set.
warnings.warn(