Commit 475ae7db authored by Harald Pfeiffer's avatar Harald Pfeiffer
Browse files

Merge branch 'LVCNR_odd_m_modes_bug' into 'master'

Fix odd-m modes bug in LVCNR

See merge request !697
parents 8226ca44 4f71bdeb
Pipeline #53796 passed with stages
in 73 minutes and 26 seconds
......@@ -251,9 +251,9 @@ UNUSED static UINT4 XLALSimInspiralNRWaveformGetSpinsFromHDF5FilePointer(
(file, "position2y-vs-time", ref_time);
pos2z = XLALSimInspiralNRWaveformGetInterpValueFromGroupAtPoint
(file, "position2z-vs-time", ref_time);
n_hat_x = pos2x - pos1x;
n_hat_y = pos2y - pos1y;
n_hat_z = pos2z - pos1z;
n_hat_x = pos1x - pos2x;
n_hat_y = pos1y - pos2y;
n_hat_z = pos1z - pos2z;
n_norm = sqrt(n_hat_x*n_hat_x + n_hat_y*n_hat_y + n_hat_z*n_hat_z);
n_hat_x = n_hat_x / n_norm;
n_hat_y = n_hat_y / n_norm;
......@@ -404,9 +404,9 @@ UNUSED static UINT4 XLALSimInspiralNRWaveformGetRotationAnglesFromH5File(
(filepointer, "position2y-vs-time", ref_time);
pos2z = XLALSimInspiralNRWaveformGetInterpValueFromGroupAtPoint
(filepointer, "position2z-vs-time", ref_time);
n_hat_x = pos2x - pos1x;
n_hat_y = pos2y - pos1y;
n_hat_z = pos2z - pos1z;
n_hat_x = pos1x - pos2x;
n_hat_y = pos1y - pos2y;
n_hat_z = pos1z - pos2z;
}
else
{
......
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