code: def _update(self):
          
          initial_time = time.time()
  
          self._updateHiddenStateTrajectories()
          self._updateEmissionProbabilities()
          self._updateTransitionMatrix()
  
          final_time = time.time()
          elapsed_time = final_time - initial_time
          logger().info("BHMM update iteration took %.3f s" % elapsed_time)
docs: Update the current model using one round of Gibbs sampling.