docs: def _process_sample (self, ap1, ap2, ap3, triple, tflags):
          """We have computed one independent phase closure triple in one timeslot.
  
          """
          # Frequency-resolved:
          np.divide (triple, np.abs (triple), triple)
          phase = np.angle (triple)
  
          self.ap_spec_stats_by_ddid[self.cur_ddid].accum (ap1, phase, tflags + 0.)
          self.ap_spec_stats_by_ddid[self.cur_ddid].accum (ap2, phase, tflags + 0.)
          self.ap_spec_stats_by_ddid[self.cur_ddid].accum (ap3, phase, tflags + 0.)
  
          # Frequency-averaged:
          triple = np.dot (triple, tflags) / tflags.sum ()
          phase = np.angle (triple)
  
          self.global_stats_by_time.accum (self.cur_time, phase)
  
          self.ap_stats_by_ddid[self.cur_ddid].accum (ap1, phase)
          self.ap_stats_by_ddid[self.cur_ddid].accum (ap2, phase)
          self.ap_stats_by_ddid[self.cur_ddid].accum (ap3, phase)
          self.bp_stats_by_ddid[self.cur_ddid].accum ((ap1, ap2), phase)
          self.bp_stats_by_ddid[self.cur_ddid].accum ((ap1, ap3), phase)
          self.bp_stats_by_ddid[self.cur_ddid].accum ((ap2, ap3), phase)
  
          self.ap_time_stats_by_ddid[self.cur_ddid].accum (self.cur_time, ap1, phase)
          self.ap_time_stats_by_ddid[self.cur_ddid].accum (self.cur_time, ap2, phase)
          self.ap_time_stats_by_ddid[self.cur_ddid].accum (self.cur_time, ap3, phase)