From 27cad277d7a294b4e1f76ea6cac51a69f9a69d78 Mon Sep 17 00:00:00 2001
From: Niklas Reinhardt <janniklas.reinhardt@aei.mpg.de>
Date: Wed, 7 Dec 2022 12:13:58 +0100
Subject: [PATCH] Correct the implementation of the the_wrt_tcb
---
lisainstrument/instrument.py | 14 +++++++-------
1 file changed, 7 insertions(+), 7 deletions(-)
diff --git a/lisainstrument/instrument.py b/lisainstrument/instrument.py
index 42db45a..1bc6cfc 100755
--- a/lisainstrument/instrument.py
+++ b/lisainstrument/instrument.py
@@ -905,14 +905,14 @@ class Instrument:
logger.debug("Computing local timer deviations")
t = self.physics_et_withinitial
- print(self.tps_wrt_tcb['1'].shape)
- print(t.shape)
self.the_wrt_tcb_withinitial = \
self.tps_wrt_tcb \
- + self.clock_offsets \
- + self.clock_freqoffsets * (t + self.tps_wrt_tcb) \
- + self.clock_freqlindrifts * (t + self.tps_wrt_tcb)**2 / 2 \
- + self.clock_freqquaddrifts * (t + self.tps_wrt_tcb)**3 / 3 \
+ + ForEachSC(lambda sc:
+ self.clock_offsets[sc] \
+ + self.clock_freqoffsets[sc] * (t + self.tps_wrt_tcb[sc]) \
+ + self.clock_freqlindrifts[sc] / 2 * (t + self.tps_wrt_tcb[sc])**2 \
+ + self.clock_freqquaddrifts[sc] / 3 * (t + self.tps_wrt_tcb[sc])**3 \
+ ) \
+ ForEachSC(lambda sc:
cumulative_trapezoid(np.broadcast_to(
self.clock_noise_fluctuations_withinitial[sc],
@@ -920,7 +920,7 @@ class Instrument:
dx=self.physics_dt, initial=0)
) \
+ ForEachSC(lambda sc:
- self.tps_wrt_tcb * self.clock_noise_fluctuations_withinitial[sc]
+ self.tps_wrt_tcb[sc] * self.clock_noise_fluctuations_withinitial[sc]
)
self.the_wrt_tcb_local = self.the_wrt_tcb_withinitial.transformed(
--
GitLab