def METHOD_NAME(self, n, rs=np.random):
    # TODO (not used yet)
    indices = rs.choice(self.linear_indices, size=n, replace=True, p=self.pdf)
    i = self.pxi[indices]
    j = self.pyi[indices]
    k = self.pzi[indices]
    return i, j, k