Skip to content
Snippets Groups Projects
New look: Off

Set likelihood to -inf if waveform_generator fails

    • Closed Issue created by Sylvia Biscoveanu @sylvia.biscoveanu

    Various XLAL issues can occur, like if the sampler steps outside the prior on its MCMC walk:

    SMPLR 0 init evolve:  35%|ââââ      | 35/100 [14:37<29:02, 26.81s/it] XLAL Error - XLALSimIMRSpinAlignedEOBModes: Mass ratio larger than 100!           [109/869]
SEOBNRv2, SEOBNRv2_opt, SEOBNRv4, and SEOBNRv4_opt are only available for mass ratios up to 100.                                                        [108/869]
XLAL Error - XLALSimIMRSpinAlignedEOBModes (LALSimIMRSpinAlignedEOB.c:993): Invalid argument                                                            [107/869]
XLAL Error - XLALSimIMRSpinAlignedEOBWaveformAll (LALSimIMRSpinAlignedEOB.c:3217): Internal function call failed: Invalid argument                      [106/869]
XLAL Error - XLALSimIMRSpinAlignedEOBWaveform (LALSimIMRSpinAlignedEOB.c:635): Internal function call failed: Invalid argument                          [105/869]
XLAL Error - XLALSimInspiralChooseTDWaveform (LALSimInspiral.c:968): Internal function call failed: Invalid argument                                    [104/869]
XLAL Error - XLALSimInspiralTDFromTD (LALSimInspiral.c:1894): Internal function call failed: Invalid argument                                           [103/869]
XLAL Error - XLALSimInspiralTD (LALSimInspiral.c:2133): Internal function call failed: Invalid argument                                                 [102/869]
XLAL Error - XLALSimInspiralFD (LALSimInspiral.c:2346): Internal function call failed: Invalid argument

    One solution is to catch the exit code returned by the LALSim error and include an exception in the waveform generator such that if it fails, the likelihood is set to -inf to indicate that this region of parameter space isn't allowed, and so that sampling can continue. LALSim returns specific exit codes for different errors: https://lscsoft.docs.ligo.org/lalsuite/lal/group___x_l_a_l_error__h.html#ggacd0e128cb3ee1cc61f98d9e1add4d267a785106db6dbfb84fcfc4c99ae74b306d

    Designs

    Child items ...

    • Activity

    • Sylvia Biscoveanu changed milestone to %Future

      changed milestone to %Future

    • Sylvia Biscoveanu added 1 deleted label

      added 1 deleted label

    • Colm Talbot
      Colm Talbot @colm.talbot ·
      Maintainer

      We could either:

      • build this into the waveform generator, and return None if an exception is raised, this should then return a -np.inf ln likelihood.
      • do this in the likelihood when calling the waveform generator.

      I'm not sure which is the better solution.

      Regardless of the method, I agree that we should definitely pipe the full error message from lalsimulation to some kind of output for debugging.

      Edited by Colm Talbot
    • Gregory Ashton added 1 deleted label

      added 1 deleted label

    • Gregory Ashton removed 1 deleted label

      removed 1 deleted label

    • Gregory Ashton added 1 deleted label

      added 1 deleted label

    • Serguei Ossokine
      Serguei Ossokine @serguei.ossokine ·
      Developer

      This would very useful, as there are several scenarios where a waveform model may refuse to generate a waveform for legitimate reasons, e.g. for SEOBNR* models if the waveform were to start too close to merger, the waveform generator would return an XLAL_EDOM. lalinference currently catches that and sets the likelihood to -inf as far as I know.

    • Sylvia Biscoveanu mentioned in merge request !653 (merged)

      mentioned in merge request !653 (merged)

    • Moritz Huebner changed milestone to %0.6.4

      changed milestone to %0.6.4

    • Colm Talbot
      Colm Talbot @colm.talbot ·
      Maintainer

      Closing after !653 (merged)

    • Colm Talbot removed 1 deleted label

      removed 1 deleted label

    • Colm Talbot closed

      closed

    • Roberto Cotesta
      Roberto Cotesta @roberto.cotesta ·
      Reporter

      Hi @sylvia.biscoveanu @colm.talbot I realised that the errors are not well handled if they appear at the beginning of the PE run. Is there anything different that is happening at the beginning of the run?

    • Colm Talbot
      Colm Talbot @colm.talbot ·
      Maintainer

      @roberto.cotesta could you provide a minimal example of the failure?

      The only thing I could imagine at the beginning of the run is that there are parts of the prior space where the waveform is failing. Is there a reason not to exclude forbidden regions in the prior?

    • Serguei Ossokine
      Serguei Ossokine @serguei.ossokine ·
      Developer

      @colm.talbot I think the issue is that sometimes the "forbidden" region is almost impossible to define in terms of the priors we use, since it uses some condition which results in a highly non trivial shape in parameter space

    • Roberto Cotesta
      Roberto Cotesta @roberto.cotesta ·
      Reporter

      Hi @colm.talbot this is the error I get when using IMRPhenomPv3HM. It's not easy to remove these regions because they depend on the spin directions and in the code to some high PN term.

      11:22 bilby INFO    : sampling_keys=['mass_1', 'mass_2', 'dec', 'ra', 'cos_theta_jn', 'psi', 'phase', 'a_1', 'a_2', 'tilt_1', 'tilt_2', 'phi_12', 'phi_jl', 'recalib_H1_amplitude_0', 'recalib_H1_amplitude_1', 'recalib_H1_amplitude_2
', 'recalib_H1_amplitude_3', 'recalib_H1_amplitude_4', 'recalib_H1_amplitude_5', 'recalib_H1_amplitude_6', 'recalib_H1_amplitude_7', 'recalib_H1_amplitude_8', 'recalib_H1_amplitude_9', 'recalib_H1_phase_0', 'recalib_H1_phase_1', 'r
ecalib_H1_phase_2', 'recalib_H1_phase_3', 'recalib_H1_phase_4', 'recalib_H1_phase_5', 'recalib_H1_phase_6', 'recalib_H1_phase_7', 'recalib_H1_phase_8', 'recalib_H1_phase_9', 'recalib_L1_amplitude_0', 'recalib_L1_amplitude_1', 'reca
lib_L1_amplitude_2', 'recalib_L1_amplitude_3', 'recalib_L1_amplitude_4', 'recalib_L1_amplitude_5', 'recalib_L1_amplitude_6', 'recalib_L1_amplitude_7', 'recalib_L1_amplitude_8', 'recalib_L1_amplitude_9', 'recalib_L1_phase_0', 'recal
ib_L1_phase_1', 'recalib_L1_phase_2', 'recalib_L1_phase_3', 'recalib_L1_phase_4', 'recalib_L1_phase_5', 'recalib_L1_phase_6', 'recalib_L1_phase_7', 'recalib_L1_phase_8', 'recalib_L1_phase_9', 'recalib_V1_amplitude_0', 'recalib_V1_a
mplitude_1', 'recalib_V1_amplitude_2', 'recalib_V1_amplitude_3', 'recalib_V1_amplitude_4', 'recalib_V1_amplitude_5', 'recalib_V1_amplitude_6', 'recalib_V1_amplitude_7', 'recalib_V1_amplitude_8', 'recalib_V1_amplitude_9', 'recalib_V
1_phase_0', 'recalib_V1_phase_1', 'recalib_V1_phase_2', 'recalib_V1_phase_3', 'recalib_V1_phase_4', 'recalib_V1_phase_5', 'recalib_V1_phase_6', 'recalib_V1_phase_7', 'recalib_V1_phase_8', 'recalib_V1_phase_9', 'time_jitter']
11:22 bilby INFO    : Periodic keys: ['ra', 'psi', 'phase', 'phi_12', 'phi_jl', 'time_jitter']
11:22 bilby INFO    : Reflective keys: ['recalib_H1_amplitude_0', 'recalib_H1_amplitude_1', 'recalib_H1_amplitude_2', 'recalib_H1_amplitude_3', 'recalib_H1_amplitude_4', 'recalib_H1_amplitude_5', 'recalib_H1_amplitude_6', 'recalib_
H1_amplitude_7', 'recalib_H1_amplitude_8', 'recalib_H1_amplitude_9', 'recalib_H1_phase_0', 'recalib_H1_phase_1', 'recalib_H1_phase_2', 'recalib_H1_phase_3', 'recalib_H1_phase_4', 'recalib_H1_phase_5', 'recalib_H1_phase_6', 'recalib
_H1_phase_7', 'recalib_H1_phase_8', 'recalib_H1_phase_9', 'recalib_L1_amplitude_0', 'recalib_L1_amplitude_1', 'recalib_L1_amplitude_2', 'recalib_L1_amplitude_3', 'recalib_L1_amplitude_4', 'recalib_L1_amplitude_5', 'recalib_L1_ampli
tude_6', 'recalib_L1_amplitude_7', 'recalib_L1_amplitude_8', 'recalib_L1_amplitude_9', 'recalib_L1_phase_0', 'recalib_L1_phase_1', 'recalib_L1_phase_2', 'recalib_L1_phase_3', 'recalib_L1_phase_4', 'recalib_L1_phase_5', 'recalib_L1_
phase_6', 'recalib_L1_phase_7', 'recalib_L1_phase_8', 'recalib_L1_phase_9', 'recalib_V1_amplitude_0', 'recalib_V1_amplitude_1', 'recalib_V1_amplitude_2', 'recalib_V1_amplitude_3', 'recalib_V1_amplitude_4', 'recalib_V1_amplitude_5',
 'recalib_V1_amplitude_6', 'recalib_V1_amplitude_7', 'recalib_V1_amplitude_8', 'recalib_V1_amplitude_9', 'recalib_V1_phase_0', 'recalib_V1_phase_1', 'recalib_V1_phase_2', 'recalib_V1_phase_3', 'recalib_V1_phase_4', 'recalib_V1_phas
e_5', 'recalib_V1_phase_6', 'recalib_V1_phase_7', 'recalib_V1_phase_8', 'recalib_V1_phase_9']
11:22 bilby INFO    : Using priors:
11:22 bilby INFO    : mass_1: Uniform(minimum=1.0, maximum=60.0, name='mass_1', latex_label='$m_{1}$', unit=None, boundary=None)
11:22 bilby INFO    : mass_2: Uniform(minimum=1.0, maximum=60.0, name='mass_2', latex_label='$m_{2}$', unit=None, boundary=None)
11:22 bilby INFO    : mass_ratio: Constraint(minimum=0.02, maximum=1, name='mass_ratio', latex_label='$q$', unit=None)
11:22 bilby INFO    : dec: Cosine(name='dec', latex_label='$\\mathrm{DEC}$', unit=None, minimum=-1.5707963267948966, maximum=1.5707963267948966, boundary=None)
11:22 bilby INFO    : ra: Uniform(minimum=0, maximum=6.283185307179586, name='ra', latex_label='$\\mathrm{RA}$', unit=None, boundary='periodic')
11:22 bilby INFO    : cos_theta_jn: Uniform(minimum=-1, maximum=1, name='cos_theta_jn', latex_label='$\\cos\\theta_{JN}$', unit=None, boundary=None)
11:22 bilby INFO    : psi: Uniform(minimum=0, maximum=3.141592653589793, name='psi', latex_label='$\\psi$', unit=None, boundary='periodic')
11:22 bilby INFO    : phase: Uniform(minimum=0, maximum=6.283185307179586, name='phase', latex_label='$\\phi$', unit=None, boundary='periodic')
11:22 bilby INFO    : a_1: Uniform(minimum=0, maximum=0.99, name='a_1', latex_label='$a_1$', unit=None, boundary=None)
11:22 bilby INFO    : a_2: Uniform(minimum=0, maximum=0.99, name='a_2', latex_label='$a_2$', unit=None, boundary=None)
11:22 bilby INFO    : tilt_1: Sine(name='tilt_1', latex_label='$\\theta_1$', unit=None, minimum=0, maximum=3.141592653589793, boundary=None)
11:22 bilby INFO    : tilt_2: Sine(name='tilt_2', latex_label='$\\theta_2$', unit=None, minimum=0, maximum=3.141592653589793, boundary=None)
11:22 bilby INFO    : phi_12: Uniform(minimum=0, maximum=6.283185307179586, name='phi_12', latex_label='$\\Delta\\phi$', unit=None, boundary='periodic')
11:22 bilby INFO    : phi_jl: Uniform(minimum=0, maximum=6.283185307179586, name='phi_jl', latex_label='$\\phi_{JL}$', unit=None, boundary='periodic')
11:22 bilby INFO    : luminosity_distance: 100.0
11:22 bilby INFO    : geocent_time: 
11:22 bilby INFO    : recalib_H1_amplitude_0: Gaussian(mu=-0.00948324596377989, sigma=0.031058514685322173, name='recalib_H1_amplitude_0', latex_label='$A^H1_0$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_amplitude_1: Gaussian(mu=-0.007187518147561477, sigma=0.01719412533843237, name='recalib_H1_amplitude_1', latex_label='$A^H1_1$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_amplitude_2: Gaussian(mu=0.006901770564401236, sigma=0.018818067610333087, name='recalib_H1_amplitude_2', latex_label='$A^H1_2$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_amplitude_3: Gaussian(mu=0.005939568154794071, sigma=0.012133856273908183, name='recalib_H1_amplitude_3', latex_label='$A^H1_3$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_amplitude_4: Gaussian(mu=0.0006369035121663344, sigma=0.009613940831560006, name='recalib_H1_amplitude_4', latex_label='$A^H1_4$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_amplitude_5: Gaussian(mu=0.004085633411212843, sigma=0.012802830537223612, name='recalib_H1_amplitude_5', latex_label='$A^H1_5$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_amplitude_6: Gaussian(mu=0.0064092629565881375, sigma=0.007904933472510041, name='recalib_H1_amplitude_6', latex_label='$A^H1_6$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_amplitude_7: Gaussian(mu=0.004884799174715116, sigma=0.0069481071637376305, name='recalib_H1_amplitude_7', latex_label='$A^H1_7$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_amplitude_8: Gaussian(mu=0.009394240376282932, sigma=0.008659290300670265, name='recalib_H1_amplitude_8', latex_label='$A^H1_8$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_amplitude_9: Gaussian(mu=0.008156793527082942, sigma=0.009337956219428126, name='recalib_H1_amplitude_9', latex_label='$A^H1_9$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_phase_0: Gaussian(mu=-0.01761507245582772, sigma=0.026220668789912444, name='recalib_H1_phase_0', latex_label='$\\phi^H1_0$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_phase_1: Gaussian(mu=0.002322530308909263, sigma=0.0146614040128832, name='recalib_H1_phase_1', latex_label='$\\phi^H1_1$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_phase_2: Gaussian(mu=0.004005566416632068, sigma=0.01345338471458025, name='recalib_H1_phase_2', latex_label='$\\phi^H1_2$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_phase_3: Gaussian(mu=-0.0036850571718003203, sigma=0.01590503154906307, name='recalib_H1_phase_3', latex_label='$\\phi^H1_3$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_phase_4: Gaussian(mu=0.00034550454900019736, sigma=0.014296787711954561, name='recalib_H1_phase_4', latex_label='$\\phi^H1_4$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_phase_5: Gaussian(mu=0.004761551834495761, sigma=0.0056336088951495795, name='recalib_H1_phase_5', latex_label='$\\phi^H1_5$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_phase_6: Gaussian(mu=0.0027124660483071436, sigma=0.00451876081632824, name='recalib_H1_phase_6', latex_label='$\\phi^H1_6$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_phase_7: Gaussian(mu=0.004777433487351768, sigma=0.004735669948299628, name='recalib_H1_phase_7', latex_label='$\\phi^H1_7$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_phase_8: Gaussian(mu=0.006368308764655624, sigma=0.0069754016970722706, name='recalib_H1_phase_8', latex_label='$\\phi^H1_8$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_phase_9: Gaussian(mu=0.003882043952897555, sigma=0.007855524479463492, name='recalib_H1_phase_9', latex_label='$\\phi^H1_9$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_H1_frequency_0: DeltaFunction(peak=19.999999999999996, name='recalib_H1_frequency_0', latex_label='$f^H1_0$', unit=None)
11:22 bilby INFO    : recalib_H1_frequency_1: DeltaFunction(peak=30.970547307245077, name='recalib_H1_frequency_1', latex_label='$f^H1_1$', unit=None)
11:22 bilby INFO    : recalib_H1_frequency_2: DeltaFunction(peak=47.95874002551529, name='recalib_H1_frequency_2', latex_label='$f^H1_2$', unit=None)
11:22 bilby INFO    : recalib_H1_frequency_3: DeltaFunction(peak=74.2654213378045, name='recalib_H1_frequency_3', latex_label='$f^H1_3$', unit=None)
11:22 bilby INFO    : recalib_H1_frequency_4: DeltaFunction(peak=115.00203724174807, name='recalib_H1_frequency_4', latex_label='$f^H1_4$', unit=None)
11:22 bilby INFO    : recalib_H1_frequency_5: DeltaFunction(peak=178.08380174125588, name='recalib_H1_frequency_5', latex_label='$f^H1_5$', unit=None)
11:22 bilby INFO    : recalib_H1_frequency_6: DeltaFunction(peak=275.76764032408107, name='recalib_H1_frequency_6', latex_label='$f^H1_6$', unit=None)
11:22 bilby INFO    : recalib_H1_frequency_7: DeltaFunction(peak=427.0337375232151, name='recalib_H1_frequency_7', latex_label='$f^H1_7$', unit=None)
11:22 bilby INFO    : recalib_H1_frequency_8: DeltaFunction(peak=661.2734284876209, name='recalib_H1_frequency_8', latex_label='$f^H1_8$', unit=None)
11:22 bilby INFO    : recalib_H1_frequency_9: DeltaFunction(peak=1024.0, name='recalib_H1_frequency_9', latex_label='$f^H1_9$', unit=None)
11:22 bilby INFO    : recalib_L1_amplitude_0: Gaussian(mu=0.009031710548076489, sigma=0.02302311541461535, name='recalib_L1_amplitude_0', latex_label='$A^L1_0$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_amplitude_1: Gaussian(mu=0.007271849076313515, sigma=0.04649882239229492, name='recalib_L1_amplitude_1', latex_label='$A^L1_1$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_amplitude_2: Gaussian(mu=-0.0011990308836519752, sigma=0.02879659101816254, name='recalib_L1_amplitude_2', latex_label='$A^L1_2$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_amplitude_3: Gaussian(mu=0.0017450923568419186, sigma=0.009704522378558627, name='recalib_L1_amplitude_3', latex_label='$A^L1_3$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_amplitude_4: Gaussian(mu=0.0044650289911179475, sigma=0.014251440343065344, name='recalib_L1_amplitude_4', latex_label='$A^L1_4$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_amplitude_5: Gaussian(mu=0.00873676782288403, sigma=0.0164648127270254, name='recalib_L1_amplitude_5', latex_label='$A^L1_5$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_amplitude_6: Gaussian(mu=0.008048922318671797, sigma=0.009725357490388588, name='recalib_L1_amplitude_6', latex_label='$A^L1_6$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_amplitude_7: Gaussian(mu=0.007351057817572623, sigma=0.00872998772773446, name='recalib_L1_amplitude_7', latex_label='$A^L1_7$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_amplitude_8: Gaussian(mu=0.005851811277763345, sigma=0.008572093142218103, name='recalib_L1_amplitude_8', latex_label='$A^L1_8$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_amplitude_9: Gaussian(mu=0.003109979729194971, sigma=0.009111455218091404, name='recalib_L1_amplitude_9', latex_label='$A^L1_9$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_phase_0: Gaussian(mu=-0.0013733005644884517, sigma=0.025778104464724477, name='recalib_L1_phase_0', latex_label='$\\phi^L1_0$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_phase_1: Gaussian(mu=-0.011747589262647283, sigma=0.045567252891636364, name='recalib_L1_phase_1', latex_label='$\\phi^L1_1$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_phase_2: Gaussian(mu=-0.0038841428191166818, sigma=0.029071803827067505, name='recalib_L1_phase_2', latex_label='$\\phi^L1_2$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_phase_3: Gaussian(mu=0.0008707616361380174, sigma=0.02548854674214132, name='recalib_L1_phase_3', latex_label='$\\phi^L1_3$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_phase_4: Gaussian(mu=0.003397299776513891, sigma=0.01948102741400912, name='recalib_L1_phase_4', latex_label='$\\phi^L1_4$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_phase_5: Gaussian(mu=0.003381357273364288, sigma=0.006632734457708485, name='recalib_L1_phase_5', latex_label='$\\phi^L1_5$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_phase_6: Gaussian(mu=0.004953494339707422, sigma=0.0065748086559958285, name='recalib_L1_phase_6', latex_label='$\\phi^L1_6$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_phase_7: Gaussian(mu=0.008334779089412624, sigma=0.006570491212245042, name='recalib_L1_phase_7', latex_label='$\\phi^L1_7$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_phase_8: Gaussian(mu=0.01362871029754724, sigma=0.006820099800675018, name='recalib_L1_phase_8', latex_label='$\\phi^L1_8$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_phase_9: Gaussian(mu=0.02079306113951578, sigma=0.008028982231097917, name='recalib_L1_phase_9', latex_label='$\\phi^L1_9$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_L1_frequency_0: DeltaFunction(peak=30.000000000000004, name='recalib_L1_frequency_0', latex_label='$f^L1_0$', unit=None)
11:22 bilby INFO    : recalib_L1_frequency_1: DeltaFunction(peak=44.40935277910318, name='recalib_L1_frequency_1', latex_label='$f^L1_1$', unit=None)
11:22 bilby INFO    : recalib_L1_frequency_2: DeltaFunction(peak=65.73968714196128, name='recalib_L1_frequency_2', latex_label='$f^L1_2$', unit=None)
11:22 bilby INFO    : recalib_L1_frequency_3: DeltaFunction(peak=97.31523192917435, name='recalib_L1_frequency_3', latex_label='$f^L1_3$', unit=None)
11:22 bilby INFO    : recalib_L1_frequency_4: DeltaFunction(peak=144.05688218409838, name='recalib_L1_frequency_4', latex_label='$f^L1_4$', unit=None)
11:22 bilby INFO    : recalib_L1_frequency_5: DeltaFunction(peak=213.24909670571103, name='recalib_L1_frequency_5', latex_label='$f^L1_5$', unit=None)
11:22 bilby INFO    : recalib_L1_frequency_6: DeltaFunction(peak=315.67514551430025, name='recalib_L1_frequency_6', latex_label='$f^L1_6$', unit=None)
11:22 bilby INFO    : recalib_L1_frequency_7: DeltaFunction(peak=467.2976300246428, name='recalib_L1_frequency_7', latex_label='$f^L1_7$', unit=None)
11:22 bilby INFO    : recalib_L1_frequency_8: DeltaFunction(peak=691.7461768201063, name='recalib_L1_frequency_8', latex_label='$f^L1_8$', unit=None)
11:22 bilby INFO    : recalib_L1_frequency_9: DeltaFunction(peak=1024.0, name='recalib_L1_frequency_9', latex_label='$f^L1_9$', unit=None)
11:22 bilby INFO    : recalib_V1_amplitude_0: Gaussian(mu=0.0, sigma=0.05000000000000005, name='recalib_V1_amplitude_0', latex_label='$A^V1_0$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_amplitude_1: Gaussian(mu=0.0, sigma=0.050000000000000024, name='recalib_V1_amplitude_1', latex_label='$A^V1_1$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_amplitude_2: Gaussian(mu=0.0, sigma=0.050000000000000044, name='recalib_V1_amplitude_2', latex_label='$A^V1_2$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_amplitude_3: Gaussian(mu=0.0, sigma=0.05000000000000004, name='recalib_V1_amplitude_3', latex_label='$A^V1_3$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_amplitude_4: Gaussian(mu=0.0, sigma=0.05000000000000005, name='recalib_V1_amplitude_4', latex_label='$A^V1_4$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_amplitude_5: Gaussian(mu=0.0, sigma=0.050000000000000044, name='recalib_V1_amplitude_5', latex_label='$A^V1_5$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_amplitude_6: Gaussian(mu=0.0, sigma=0.05000000000000003, name='recalib_V1_amplitude_6', latex_label='$A^V1_6$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_amplitude_7: Gaussian(mu=0.0, sigma=0.05000000000000003, name='recalib_V1_amplitude_7', latex_label='$A^V1_7$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_amplitude_8: Gaussian(mu=0.0, sigma=0.050000000000000044, name='recalib_V1_amplitude_8', latex_label='$A^V1_8$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_amplitude_9: Gaussian(mu=0.0, sigma=0.050000000000000044, name='recalib_V1_amplitude_9', latex_label='$A^V1_9$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_phase_0: Gaussian(mu=0.0, sigma=0.034929197639720425, name='recalib_V1_phase_0', latex_label='$\\phi^V1_0$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_phase_1: Gaussian(mu=0.0, sigma=0.03496078300818482, name='recalib_V1_phase_1', latex_label='$\\phi^V1_1$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_phase_2: Gaussian(mu=0.0, sigma=0.03503640793642491, name='recalib_V1_phase_2', latex_label='$\\phi^V1_2$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_phase_3: Gaussian(mu=0.0, sigma=0.03521708937024615, name='recalib_V1_phase_3', latex_label='$\\phi^V1_3$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_phase_4: Gaussian(mu=0.0, sigma=0.03564662257025354, name='recalib_V1_phase_4', latex_label='$\\phi^V1_4$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_phase_5: Gaussian(mu=0.0, sigma=0.036656116163713476, name='recalib_V1_phase_5', latex_label='$\\phi^V1_5$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_phase_6: Gaussian(mu=0.0, sigma=0.038970427043521064, name='recalib_V1_phase_6', latex_label='$\\phi^V1_6$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_phase_7: Gaussian(mu=0.0, sigma=0.044027143041336855, name='recalib_V1_phase_7', latex_label='$\\phi^V1_7$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_phase_8: Gaussian(mu=0.0, sigma=0.05426593740433323, name='recalib_V1_phase_8', latex_label='$\\phi^V1_8$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_phase_9: Gaussian(mu=0.0, sigma=0.07319891790170277, name='recalib_V1_phase_9', latex_label='$\\phi^V1_9$', unit=None, boundary='reflective')
11:22 bilby INFO    : recalib_V1_frequency_0: DeltaFunction(peak=19.999999999999996, name='recalib_V1_frequency_0', latex_label='$f^V1_0$', unit=None)
11:22 bilby INFO    : recalib_V1_frequency_1: DeltaFunction(peak=30.970547307245077, name='recalib_V1_frequency_1', latex_label='$f^V1_1$', unit=None)
11:22 bilby INFO    : recalib_V1_frequency_2: DeltaFunction(peak=47.95874002551529, name='recalib_V1_frequency_2', latex_label='$f^V1_2$', unit=None)
11:22 bilby INFO    : recalib_V1_frequency_3: DeltaFunction(peak=74.2654213378045, name='recalib_V1_frequency_3', latex_label='$f^V1_3$', unit=None)
11:22 bilby INFO    : recalib_V1_frequency_4: DeltaFunction(peak=115.00203724174807, name='recalib_V1_frequency_4', latex_label='$f^V1_4$', unit=None)
11:22 bilby INFO    : recalib_V1_frequency_5: DeltaFunction(peak=178.08380174125588, name='recalib_V1_frequency_5', latex_label='$f^V1_5$', unit=None)
11:22 bilby INFO    : recalib_V1_frequency_6: DeltaFunction(peak=275.76764032408107, name='recalib_V1_frequency_6', latex_label='$f^V1_6$', unit=None)
11:22 bilby INFO    : recalib_V1_frequency_7: DeltaFunction(peak=427.0337375232151, name='recalib_V1_frequency_7', latex_label='$f^V1_7$', unit=None)
11:22 bilby INFO    : recalib_V1_frequency_8: DeltaFunction(peak=661.2734284876209, name='recalib_V1_frequency_8', latex_label='$f^V1_8$', unit=None)
11:22 bilby INFO    : recalib_V1_frequency_9: DeltaFunction(peak=1024.0, name='recalib_V1_frequency_9', latex_label='$f^V1_9$', unit=None)
11:22 bilby INFO    : time_jitter: Uniform(minimum=-0.00048828125, maximum=0.00048828125, name=None, latex_label=None, unit=None, boundary='periodic')
11:22 bilby INFO    : Initialize NestedSampler with {
 "bound": "multi",
 "enlarge": 1.5,
 "facc": 0.5,
 "first_update": {
  "min_eff": 10,
  "min_ncall": 3000
 },
 "maxmcmc": 5000,
 "nact": 5,
 "nlive": 1500,
 "sample": "rwalk",
 "vol_check": 8,
 "vol_dec": 0.5,
 "walks": 100
}
XLAL Error - InitializeSystem (../../../lalsimulation/lib/LALSimInspiralFDPrecAngles_internals.c:270): Omegaz5 = 1456.7873362282709877 which is larger than expected. Not generating a waveform here.

XLAL Error - InitializeSystem (../../../lalsimulation/lib/LALSimInspiralFDPrecAngles_internals.c:270): Input domain error
XLAL Error - XLALOrbitalAngMom3PNSpinning (../../../lalsimulation/lib/LALSimInspiralFDPrecAngles.c:198): InitializeSystem failed
XLAL Error - XLALOrbitalAngMom3PNSpinning (../../../lalsimulation/lib/LALSimInspiralFDPrecAngles.c:198): Internal function call failed: Input domain error
XLAL Error - PhenomInternal_OrbAngMom3PN (../../../lalsimulation/lib/LALSimIMRPhenomInternalUtils.c:283): XLALOrbitalAngMom3PNSpinning failed.
XLAL Error - PhenomInternal_OrbAngMom3PN (../../../lalsimulation/lib/LALSimIMRPhenomInternalUtils.c:283): Internal function call failed: Input domain error
XLAL Error - InitializeSystem (../../../lalsimulation/lib/LALSimInspiralFDPrecAngles_internals.c:270): Omegaz5 = 1456.7873362282709877 which is larger than expected. Not generating a waveform here.

XLAL Error - InitializeSystem (../../../lalsimulation/lib/LALSimInspiralFDPrecAngles_internals.c:270): Input domain error
XLAL Error - init_PhenomPv3HM_Storage (../../../lalsimulation/lib/LALSimIMRPhenomPv3HM.c:274): InitializeSystem failed
XLAL Error - init_PhenomPv3HM_Storage (../../../lalsimulation/lib/LALSimIMRPhenomPv3HM.c:274): Internal function call failed: Input domain error
XLAL Error - XLALSimIMRPhenomPv3HMGetHplusHcross (../../../lalsimulation/lib/LALSimIMRPhenomPv3HM.c:385): init_PhenomPv3HM_Storage failed
XLAL Error - XLALSimIMRPhenomPv3HMGetHplusHcross (../../../lalsimulation/lib/LALSimIMRPhenomPv3HM.c:385): Internal function call failed: Input domain error
XLAL Error - XLALSimInspiralChooseFDWaveform (../../../lalsimulation/lib/LALSimInspiral.c:1960): Internal function call failed: Input domain error
11:22 bilby WARNING : Evaluating the waveform failed with error: Internal function call failed: Input domain error
The parameters were {'mass_1': 3.2808496489702265e+31, 'mass_2': 1.919389755604094e+31, 'spin_1': (array(0.10447106), array(-0.25623486), array(-0.37294073)), 'spin_2': (array(0.34334453), array(-0.25623486), array(-0.83691068)), '
luminosity_distance': 3.085677581491367e+24, 'iota': array(1.49187458), 'phase': 3.1723843378852434, 'eccentricity': 0.0, 'start_frequency': 20}
Likelihood will be set to -inf.
/home/sossokine/envs/modern3/lib/python3.7/site-packages/dynesty/dynesty.py:373: UserWarning: A note of caution: having `nlive < ndim * (ndim + 1) // 2` may result in unconstrained bounding distributions.
  warnings.warn("A note of caution: "
Traceback (most recent call last):
  File "/home/sossokine/envs/modern3/bin/parallel_bilby_analysis", line 11, in <module>
    load_entry_point('parallel-bilby==0.1.2', 'console_scripts', 'parallel_bilby_analysis')()
  File "/home/sossokine/envs/modern3/lib/python3.7/site-packages/pkg_resources/__init__.py", line 490, in load_entry_point
    return get_distribution(dist).load_entry_point(group, name)
  File "/home/sossokine/envs/modern3/lib/python3.7/site-packages/pkg_resources/__init__.py", line 2853, in load_entry_point
    return ep.load()
  File "/home/sossokine/envs/modern3/lib/python3.7/site-packages/pkg_resources/__init__.py", line 2444, in load
    return self.resolve()
  File "/home/sossokine/envs/modern3/lib/python3.7/site-packages/pkg_resources/__init__.py", line 2450, in resolve
    module = __import__(self.module_name, fromlist=['__name__'], level=0)
  File "/home/sossokine/envs/modern3/lib/python3.7/site-packages/parallel_bilby-0.1.2-py3.7.egg/parallel_bilby/analysis.py", line 525, in <module>
    **init_sampler_kwargs,
  File "/home/sossokine/envs/modern3/lib/python3.7/site-packages/dynesty/dynesty.py", line 523, in NestedSampler
    live_points[1][i]))
ValueError: The log-likelihood (nan) of live point 575 located at u=[0.5449461  0.05996337 0.8592629  0.08377273 0.04117
 0.08467092 0.91669568 0.86717991 0.64842079 0.11786119
 0.37710353 0.95683466 0.06912789 0.95326059 0.68091307
 0.81790561 0.4073056  0.07928778 0.25911268 0.73948133
 0.43235686 0.92296221 0.6949312  0.10239935 0.32454368
 0.4398925  0.17883582 0.14853475 0.4201956  0.80304055
 0.59645636 0.15355395 0.79257409 0.87897312 0.88215923
 0.02313054 0.59382942 0.41553409 0.60373097 0.0344003
 0.46100025 0.46073657 0.66226719 0.91259949 0.05415461
 0.95149725 0.04952897 0.30672676 0.12056345 0.10227468
 0.18351813 0.1160018  0.3904681  0.71525295 0.3708249
 0.21956906 0.31540079 0.40344707 0.2567768  0.23387467
 0.24299617 0.07769591 0.65898277 0.03638956 0.04588957
 0.61372448 0.70930895 0.28808388 0.40068907 0.86197617
 0.92110922 0.07780648 0.60313486 0.94172446] v=[ 3.31518199e+01  4.53783890e+00  8.01680371e-01  5.26359585e-01
 -9.17660001e-01  2.66001525e-01  5.75976882e+00  8.58508114e-01
  6.41936586e-01  7.00875896e-01  1.32245868e+00  6.01196951e+00
  4.34343355e-01  4.26119983e-02  8.98080803e-04  2.39775177e-02
  3.09440335e-03 -1.29175852e-02 -4.18606253e-03  1.14822318e-02
  3.70100209e-03  2.17361749e-02  1.29180032e-02 -5.08628249e-02
 -4.34885439e-03  1.97085280e-03 -1.83146768e-02 -1.45623042e-02
  3.62698175e-03  6.56485405e-03  5.93381462e-03 -7.55734146e-04
  1.02873227e-02  3.59657357e-02  6.24124782e-02 -5.85907333e-02
  4.04901256e-03  1.42474381e-03  1.30672811e-02 -9.64857245e-03
  6.49626880e-03  5.00678927e-03  6.92456932e-03  3.36059794e-02
 -8.49212760e-02  4.43618993e-02 -4.11710160e-02 -6.44354053e-03
 -4.39336664e-03 -3.38794176e-03  2.40794109e-03  5.47723460e-03
  1.85602076e-02  2.84398357e-02 -1.64834691e-02 -3.86824721e-02
 -2.40299467e-02 -1.22217365e-02 -3.26657203e-02 -3.63072949e-02
 -3.48348581e-02 -7.10370929e-02  2.04844258e-02 -6.26704335e-02
 -5.89469056e-02  1.01269106e-02  1.94175422e-02 -1.99261463e-02
 -9.22135638e-03  4.24481845e-02  6.21915160e-02 -7.70566476e-02
  1.91393023e-02  4.31371545e-04] is invalid.
      Edited by Roberto Cotesta
    • Colm Talbot
      Colm Talbot @colm.talbot ·
      Maintainer

      Thanks, both. That's a reasonable reason. Is there by any chance a function that quickly evaluates whether a waveform can be generated that we could call as a constraint on the prior? Or does it require actually evaluating the waveform?

      This error looks weird. It looks like the value isn't propagated through the sampler properly. I'll dig in a little and try to reproduce the error.

    • Colm Talbot reopened

      reopened

    • Colm Talbot
      Colm Talbot @colm.talbot ·
      Maintainer

      I'm beginning to suspect that this is actually not the source of the issue. I think the issue is actually unrelated and there's some other effect causing the likelihood to return NaN.

    • Roberto Cotesta
      Roberto Cotesta @roberto.cotesta ·
      Reporter

      The developer of the waveform model are probably the best to ask. It looks like they are checking the value of this coefficient Omegaz5(here). Maybe they can provide the function to compute it without computing the entire waveform.

    • Colm Talbot
      Colm Talbot @colm.talbot ·
      Maintainer

      I'm going to reclose this issue, as I'm reasonably convinced it isn't causing this failure.

    • Colm Talbot closed

      closed

    • Michael Williams unassigned @sylvia.biscoveanu

      unassigned @sylvia.biscoveanu

    Please register or sign in to reply