Skip to content

Calibration interpolation error: `x and y arrays must be equal in length along interpolation axis`

When @cjhaster ran the following command line:

$ summarypages --webdir /home/carl-johan.haster/public_html/highQ/S190814bv/FirstSetOfRuns_SpecialInstall_q20_PESummary --samples /home/carl-johan.haster/projects/HighQ/GW190814/GW190814_highQ_SpecialInstall_q20_mcmc/NRHybSur2dq15pseudoFourPN/posterior_samples/posterior_V1H1L1_1249852257.01-347292.hdf5 /home/carl-johan.haster/projects/HighQ/GW190814/GW190814_highQ_SpecialInstall_q20_mcmc/SEOBNRv4HM_ROMpseudoFourPN/posterior_samples/posterior_V1H1L1_1249852257.01-347292.hdf5 /home/carl-johan.haster/projects/HighQ/GW190814/GW190814_highQ_SpecialInstall_q20_mcmc/IMRPhenomHMpseudoFourPN/posterior_samples/posterior_V1H1L1_1249852257.01-347292.hdf5 --psd H1:/home/carl-johan.haster/projects/HighQ/GW190814/PSDs/glitch_median_PSD_forLI_H1.dat L1:/home/carl-johan.haster/projects/HighQ/GW190814/PSDs/glitch_median_PSD_forLI_L1.dat V1:/home/carl-johan.haster/projects/HighQ/GW190814/PSDs/glitch_median_PSD_forLI_V1.dat --config /home/carl-johan.haster/projects/HighQ/GW190814/GW190814_highQ_SpecialInstall_q20_mcmc/NRHybSur2dq15pseudoFourPN/config.ini /home/carl-johan.haster/projects/HighQ/GW190814/GW190814_highQ_SpecialInstall_q20_mcmc/SEOBNRv4HM_ROMpseudoFourPN/config.ini /home/carl-johan.haster/projects/HighQ/GW190814/GW190814_highQ_SpecialInstall_q20_mcmc/IMRPhenomHMpseudoFourPN/config.ini --approximant NRHybSur2dq15 SEOBNRv4HM_ROM IMRPhenomHM --gracedb S190814bv --labels NRHybSur2dq15 SEOBNRv4HM_ROM IMRPhenomHM --gwdata H1:/home/carl-johan.haster/projects/HighQ/GW190814/GW190814_highQ_SpecialInstall_q20_mcmc/SEOBNRv4HM_ROMpseudoFourPN/caches/H-H1_HOFT_C01_CACHE-1249851715-561.lcf L1:/home/carl-johan.haster/projects/HighQ/GW190814/GW190814_highQ_SpecialInstall_q20_mcmc/SEOBNRv4HM_ROMpseudoFourPN/caches/L-L1_HOFT_C01_CACHE-1249851715-561.lcf V1:/home/carl-johan.haster/projects/HighQ/GW190814/GW190814_highQ_SpecialInstall_q20_mcmc/SEOBNRv4HM_ROMpseudoFourPN/caches/V-V1Online_CACHE-1249851715-561.lcf --no_ligo_skymap

he found the following error:

Traceback (most recent call last):
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-20200118/bin/summarypages", line 10, in <module>
    sys.exit(main())
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-20200118/lib/python3.7/site-packages/pesummary/cli/summarypages.py", line 196, in main
    args = func["input"](opts)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-20200118/lib/python3.7/site-packages/pesummary/gw/inputs.py", line 755, in __init__
    self.calibration = self.opts.calibration
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-20200118/lib/python3.7/site-packages/pesummary/gw/inputs.py", line 377, in calibration
    calibration_data = f.calibration_data_in_results_file
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-20200118/lib/python3.7/site-packages/pesummary/gw/file/formats/lalinference.py", line 114, in calibration_data_in_results_file
    return self.grab_calibration_data(grab, self.path_to_results_file)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-20200118/lib/python3.7/site-packages/pesummary/gw/file/formats/base_read.py", line 247, in grab_calibration_data
    in np.column_stack(amp_params[key])]
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-20200118/lib/python3.7/site-packages/pesummary/gw/file/formats/base_read.py", line 246, in <listcomp>
    fill_value=0, bounds_error=False)(np.log(fs)) for samp
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-20200118/lib/python3.7/site-packages/scipy/interpolate/interpolate.py", line 433, in __init__
    _Interpolator1D.__init__(self, x, y, axis=axis)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-20200118/lib/python3.7/site-packages/scipy/interpolate/polyint.py", line 60, in __init__
    self._set_yi(yi, xi=xi, axis=axis)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/conda/envs/igwn-py37-20200118/lib/python3.7/site-packages/scipy/interpolate/polyint.py", line 125, in _set_yi
    raise ValueError("x and y arrays must be equal in length along "
ValueError: x and y arrays must be equal in length along interpolation axis