Resolve "Bug when computing THE with respect to TCB"

Closes #115 (closed)

• Prevent Numpy's element-wise operations to compute correctly the_wrt_tcb_withinitial

lisainstrument/instrument.py

@@ -916,9 +916,9 @@ class Instrument:
         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 \
+            + self.clock_freqoffsets * (self.tps_wrt_tcb + t) \
+            + self.clock_freqlindrifts * (self.tps_wrt_tcb + t)**2 / 2 \
+            + self.clock_freqquaddrifts * (self.tps_wrt_tcb + t)**3  / 3 \
             + self.tps_wrt_tcb * self.clock_noise_fluctuations_withinitial \
             + self.integrated_clock_noise_fluctuations_withinitial