code: def query(self, query_samples):
  
          
  
          self.sampled_topics = np.zeros((self.samples, self.N),
                                         dtype=np.int)
  
          for s in range(self.samples):            
  
              self.sampled_topics[s, :] = \
                  samplers_lda.sampler_query(self.docid, self.tokens,
                                             self.topic_seed,
                                             np.ascontiguousarray(
                                                 self.tt[:, :, s],
                                                 dtype=np.float),
                                             self.N, self.K, self.D,
                                             self.alpha, query_samples)
  
              print("Sample %d queried" % s)
  
          self.dt = np.zeros((self.D, self.K, self.samples))
  
          for s in range(self.samples):            
              self.dt[:, :, s] = \
                 samplers_lda.dt_comp(self.docid, self.sampled_topics[s, :],
                                      self.N, self.K, self.D, self.alpha)
docs: Query docs with query_samples number of Gibbs
          sampling iterations.