add SEOBNRv5_ROM time length estimator (gwastro#4451)

pycbc/pnutils.py

@@ -524,6 +524,9 @@ def _get_imr_duration(m1, m2, s1z, s2z, f_low, approximant="SEOBNRv4"):
         # NB for no clear reason this function has f_low as first argument
         time_length = lalsimulation.SimIMRSEOBNRv4ROMTimeOfFrequency(
                            f_low, m1 * lal.MSUN_SI, m2 * lal.MSUN_SI, s1z, s2z)
+    elif approximant == 'SEOBNRv5_ROM':
+        time_length = lalsimulation.SimIMRSEOBNRv5ROMTimeOfFrequency(
+                           f_low, m1 * lal.MSUN_SI, m2 * lal.MSUN_SI, s1z, s2z)
     elif approximant == 'SPAtmplt' or approximant == 'TaylorF2':
         chi = lalsimulation.SimInspiralTaylorF2ReducedSpinComputeChi(
             m1, m2, s1z, s2z

pycbc/waveform/waveform.py

@@ -1030,6 +1030,11 @@ def seobnrv4_length_in_time(**kwds):
     """
     return get_imr_length("SEOBNRv4", **kwds)
 
+def seobnrv5_length_in_time(**kwds):
+    """Stub for holding the calculation of SEOBNRv5_ROM waveform duration.
+    """
+    return get_imr_length("SEOBNRv5_ROM", **kwds)
+
 def imrphenomd_length_in_time(**kwds):
     """Stub for holding the calculation of IMRPhenomD waveform duration.
     """
@@ -1091,6 +1096,7 @@ def get_hm_length_in_time(lor_approx, maxm_default, **kwargs):
 _filter_time_lengths["SEOBNRv4HM_ROM"] = seobnrv4hm_length_in_time
 _filter_time_lengths["SEOBNRv4"] = seobnrv4_length_in_time
 _filter_time_lengths["SEOBNRv4P"] = seobnrv4_length_in_time
+_filter_time_lengths["SEOBNRv5_ROM"] = seobnrv5_length_in_time
 _filter_time_lengths["IMRPhenomC"] = imrphenomd_length_in_time
 _filter_time_lengths["IMRPhenomD"] = imrphenomd_length_in_time
 _filter_time_lengths["IMRPhenomPv2"] = imrphenomd_length_in_time