Skip to content

GitLab

Commit 73f6a8ef authored by Marc Arene's avatar Marc Arene
Browse files

Describing why I left some commented lines

parent c2078762
Hide whitespace changes
Inline Side-by-side
Showing with 4 additions and 0 deletions
Library/likelihood_gradient.py
View file @ 73f6a8ef
......@@ -374,6 +374,7 @@ class GWTransientLikelihoodGradient(LikelihoodGradient):
ifo.exp_two_pi_f_dt_at_point = ifo.get_fd_time_translation(self.likelihood.parameters)
ifo.template_at_point = ifo.get_detector_response(waveform_polarizations, self.likelihood.parameters, exp_two_pi_f_dt=ifo.exp_two_pi_f_dt_at_point)
ifo.dh = {}
# # [DEBUG]: old way of computing inner products which does not use the strain_data.window_factor
# ifo.s_inner_h_at_point = nwip(ifo.fd_strain, ifo.template_at_point, ifo.psd_array, ifo.strain_data.duration)
# ifo.h_inner_h_at_point = np.real(nwip(ifo.template_at_point, ifo.template_at_point, ifo.psd_array, ifo.strain_data.duration))
ifo.s_inner_h_at_point = ifo.inner_product(ifo.template_at_point)
......@@ -504,6 +505,7 @@ class GWTransientLikelihoodGradient(LikelihoodGradient):
dh = (template_plus - template_minus) / (2 * offset)
ifo.dh[traj_param_key] = dh
# # [DEBUG]: old way of computing inner products which does not use the strain_data.window_factor
# ifo.s_inner_h_plus = nwip(ifo.fd_strain, template_plus, ifo.psd_array, ifo.strain_data.duration)
# ifo.s_inner_h_minus = nwip(ifo.fd_strain, template_minus, ifo.psd_array, ifo.strain_data.duration)
# ifo.h_inner_dh = nwip(ifo.template_at_point, dh, ifo.psd_array, ifo.strain_data.duration)
......@@ -555,6 +557,7 @@ class GWTransientLikelihoodGradient(LikelihoodGradient):
dh = (template_plus - ifo.template_at_point) / offset
ifo.dh[traj_param_key] = dh
# # [DEBUG]: old way of computing inner products which does not use the strain_data.window_factor
# ifo.s_inner_h_plus = nwip(ifo.fd_strain, template_plus, ifo.psd_array, ifo.strain_data.duration)
# ifo.h_inner_dh = nwip(ifo.template_at_point, dh, ifo.psd_array, ifo.strain_data.duration)
ifo.s_inner_h_plus = ifo.inner_product(template_plus)
......@@ -602,6 +605,7 @@ class GWTransientLikelihoodGradient(LikelihoodGradient):
dh = (ifo.template_at_point - template_minus) / offset
ifo.dh[traj_param_key] = dh
# # [DEBUG]: old way of computing inner products which does not use the strain_data.window_factor
# ifo.s_inner_h_minus = nwip(ifo.fd_strain, template_minus, ifo.psd_array, ifo.strain_data.duration)
# ifo.h_inner_dh = nwip(ifo.template_at_point, dh, ifo.psd_array, ifo.strain_data.duration)
ifo.s_inner_h_minus = ifo.inner_product(template_minus)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment