Skip to content
Snippets Groups Projects
Commit c29102cd authored by Wanting Niu's avatar Wanting Niu :speech_balloon:
Browse files

gstlal_inspiral_bank_splitter: fix the calculation of spin1 and spin2 using numpy.dot

parent 727a03f4
No related branches found
No related tags found
1 merge request!369Document the template_duration in svd_bank match to the actual waveform
Pipeline #482028 failed
......@@ -276,9 +276,9 @@ for outputrows in banks_subbins:
row.template_duration = chirptime.imr_time(options.f_low,
lal.MSUN_SI * row.mass1,
lal.MSUN_SI * row.mass2,
(row.spin1x**2 + row.spin1y**2 + row.spin1z**2)**.5,
(row.spin2x**2 + row.spin2y**2 + row.spin2z**2)**.5,
f_max = min(row.f_final, 2 * chirptime.ringf(lal.MSUN_SI * row.mass1 + lal.MSUN_SI * row.mass2, chirptime.overestimate_j_from_chi(max((row.spin1x**2 + row.spin1y**2 + row.spin1z**2)**.5, (row.spin2x**2 + row.spin2y**2 + row.spin2z**2)**.5))) if approximant in templates.gstlal_IMR_approximants else spawaveform.ffinal(row.mass1, row.mass2, 'bkl_isco')))
numpy.dot(row.spin1x**2 + row.spin1y**2 + row.spin1z**2, row.spin1x**2 + row.spin1y**2 + row.spin1z**2)**.5,
numpy.dot(row.spin2x**2 + row.spin2y**2 + row.spin2z**2, row.spin2x**2 + row.spin2y**2 + row.spin2z**2)**.5,
f_max = min(row.f_final, 2 * chirptime.ringf(lal.MSUN_SI * row.mass1 + lal.MSUN_SI * row.mass2, chirptime.overestimate_j_from_chi(max(numpy.dot(row.spin1x**2 + row.spin1y**2 + row.spin1z**2, row.spin1x**2 + row.spin1y**2 + row.spin1z**2)**.5, numpy.dot(row.spin2x**2 + row.spin2y**2 + row.spin2z**2, row.spin2x**2 + row.spin2y**2 + row.spin2z**2)**.5))) if approximant in templates.gstlal_IMR_approximants else spawaveform.ffinal(row.mass1, row.mass2, 'bkl_isco')))
# template_duration may be overwritten by svd_bank if the flow is increased in order to satisfy the target maximum duration
outputrows[ind] = (approximant, rows)
num_banks += len(outputrows)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment