  A finite set of the Euclidean space is called an $s$-distance set provided
the number of Euclidean distances in the set is $s$. Determining the largest
possible $s$-distance set for the Euclidean space of a given dimension is
challenging. This problem was solved only when dealing with small values of $s$
and dimensions. Lison\v{e}k (1997) achieved the classification of the largest
2-distance sets for dimensions up to $7$, using computer assistance and graph
representation theory. In this study, we consider a theory analogous to these
results of Lison\v{e}k for the pseudo-Euclidean space $\mathbb{R}^{p,q}$. We
consider an $s$-indefinite-distance set in a pseudo-Euclidean space that uses
the value \[ || x-y ||=(x_1-y_1)^2 +\cdots +(x_p
-y_p)^2-(x_{p+1}-y_{p+1})^2-\cdots -(x_{p+q}-y_{p+q})^2 \] instead of the
Euclidean distance. We develop a representation theory for symmetric matrices
in the context of $s$-indefinite-distance sets, which includes or improves the
results of Euclidean $s$-distance sets with large $s$ values. Moreover, we
classify the largest possible $2$-indefinite-distance sets for small
dimensions.
