code: def get_matrix(self, indices):
          
          new = numpy.empty(self.samples1.shape)
          for idx in range(len(indices)):
              if indices[idx]:
                  new[idx] = self.samples1[idx]
              else:
                  new[idx] = self.samples2[idx]
  
          if self.poly:
              new = self.poly(*new)
          return new
docs: Retrieve Saltelli matrix.