From 72d3de6516a962ed399f21eb6ebd5e9e80eadab7 Mon Sep 17 00:00:00 2001
From: joanvallve <jvallve@iri.upc.edu>
Date: Mon, 21 Jun 2021 16:41:17 +0200
Subject: [PATCH] bug when calling computeSagnacCorrection

---
 src/tdcp.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/tdcp.cpp b/src/tdcp.cpp
index 110c0f2..5d5ff39 100644
--- a/src/tdcp.cpp
+++ b/src/tdcp.cpp
@@ -248,8 +248,8 @@ TdcpOutput Tdcp(SnapshotPtr               snapshot_r,
         // sagnac corrections
         if (tdcp_params.sagnac_correction != 0)
         {
-            double sagnac_corr_r = computeSagnacCorrection(x_r,s_r); //OMGE*(s_r.col(row)(0)*x_r(1)-s_r.col(row)(1)*x_r(0))/CLIGHT;
-            double sagnac_corr_k = computeSagnacCorrection(x_r,s_k); //OMGE*(s_k.col(row)(0)*(x_r(1)+d_0(1))-s_k.col(row)(1)*(x_r(0)+d_0(0)))/CLIGHT;
+            double sagnac_corr_r = computeSagnacCorrection(x_r,s_r.col(row)); //OMGE*(s_r.col(row)(0)*x_r(1)-s_r.col(row)(1)*x_r(0))/CLIGHT;
+            double sagnac_corr_k = computeSagnacCorrection(x_r,s_k.col(row)); //OMGE*(s_k.col(row)(0)*(x_r(1)+d_0(1))-s_k.col(row)(1)*(x_r(0)+d_0(0)))/CLIGHT;
 
             if (tdcp_params.sagnac_correction == 1)
                 drho_m(row) += -sagnac_corr_k + sagnac_corr_r;
-- 
GitLab