How can we be sure about the constancy of atomic clocks as in the Hafele and Keating time dilation test?
Of any one clock $\mathfrak A$ it can be measured whether its average rates were equal to each other, or by how much they differed from each other:
In regard to
any two (distinct) indications $A_{\circ F}, A_{\circ G} \in \mathcal A$ of an identifiable material participant $A$, with assigned readings $t_{\mathfrak A}[ \, A_{\circ F} \, ], t_{\mathfrak A}[ \, A_{\circ G} \, ] \in \mathbb R$ and
any two (distinct) indications $A_{\circ J}, A_{\circ K} \in \mathcal A$ of an identifiable material participant $A$, with assigned readings $t_{\mathfrak A}[ \, A_{\circ J} \, ], t_{\mathfrak A}[ \, A_{\circ K} \, ] \in \mathbb R$
clock $\mathfrak A \equiv (\mathcal A, t_{\mathfrak A})$ had run from $A_{\circ F}$ until $A_{\circ G}$ and from $A_{\circ J}$ until $A_{\circ K}$ at equal average rates if
$$(t_{\mathfrak A}[ \, A_{\circ K} \, ] - t_{\mathfrak A}[ \, A_{\circ J} \, ]) = (t_{\mathfrak A}[ \, A_{\circ G} \, ] - t_{\mathfrak A}[ \, A_{\circ F} \, ]) \, \left( \frac{\tau^A_{[ \, \circ J, \, \circ K \, ]}}{\tau^A_{[ \, \circ F, \, \circ G \, ]}} \right), $$
where $\large \left( \frac{\tau^A_{[ \, \circ J, \, \circ K \, ]}}{\tau^A_{[ \, \circ F, \, \circ G \, ]}} \right)$, i.e. the value of the ratio between $A$'s duration from its indication $A_{\circ J}$ until its indication $A_{\circ K}$ and $A$'s duration from its indication $A_{\circ F}$ until its indication $A_{\circ G}$ can and must be measured beforehand (by the familiar methods of the theory of relativity), of course.
If all average rates of clock $\mathfrak A$ between all (relevant) pairs of its indications were equal, then $\mathfrak A$ is said to have run at constant rate.
Also of interest is the comparison between the average rates of different clocks:
In regard to
any two (distinct) indications $A_{\circ J}, A_{\circ K} \in \mathcal A$ of an identifiable material participant $A$, with assigned readings $t_{\mathfrak A}[ \, A_{\circ J} \, ], t_{\mathfrak A}[ \, A_{\circ K} \, ] \in \mathbb R$ and
any two (distinct) indications $B_{\circ P}, B_{\circ Q} \in \mathcal B$ of an identifiable material participant $B$, with assigned readings $t_{\mathfrak B}[ \, B_{\circ P} \, ], t_{\mathfrak B}[ \, B_{\circ Q} \, ] \in \mathbb R$
the two clocks $\mathfrak A \equiv (\mathcal A, t_{\mathfrak A})$ and $\mathfrak B \equiv (\mathcal B, t_{\mathfrak B})$ had run at equal average rates if
$$(t_{\mathfrak B}[ \, B_{\circ Q} \, ] - t_{\mathfrak B}[ \, B_{\circ P} \, ]) = (t_{\mathfrak A}[ \, A_{\circ K} \, ] - t_{\mathfrak A}[ \, A_{\circ J} \, ]) \, \left( \frac{\tau^B_{[ \, \circ P, \, \circ Q \, ]}}{\tau^A_{[ \, \circ J, \, \circ K \, ]}} \right), $$
where $\large \left( \frac{\tau^B_{[ \, \circ P, \, \circ Q \, ]}}{\tau^A_{[ \, \circ J, \, \circ K \, ]}} \right)$, i.e. the value of the ratio between $B$'s duration from its indication $B_{\circ P}$ until its indication $B_{\circ Q}$ and $A$'s duration from its indication $A_{\circ J}$ until its indication $A_{\circ K}$ can and must be measured beforehand (by the familiar methods of the theory of relativity).