Let $E$ be a supersingular curve over a field of characteristic $p$ with endomorphism ring $\mathcal O_D$ which is a maximal order in a division ring $D$ over $Q$ ramified at $p$ and $\infty$.
The norm defines a (positive definite) quadratic form in $4$ variables, let us call it $Q(x,y,z,w)$. Given an integer $n$, is it known what the number of solutions to $Q(x,y,z,w) = n$ is? I would like this to depend only on $p$, $n$ and the size of the unit group in the maximal order.
For instance, if $D$ were the Hamitonian quaternions (over $\mathbb Q$), then this is classical and Wikipedia has a formula essentially in terms of $\sigma(n)$, the sum of divisors on $n$ here. They don't give the formula for the maximal order but I am really interested in the maximal order and even just the case of $n$ prime might be good enough.