I am interested in generating a list of n 4-vectors (t,x,y,z) such that -t^2+x^2+y^2+z^2=0 for each vector and the sum of the n 4-vectors equals zero. All of the t,x,y,z are real. I, particular, I am interested in solutions where none of the entries are zero, thus excluding the trivial solution and the case in which thing are collinear, e.g. the four vectors (plus/minus 1,0,0,plus/minus 1). I'm interested in this question because it is relevant to the study of scattering amplitudes in the center of mass frame. Is there an efficient algorithmic implementation of this, or, barring that, an implementation that works for n~10 to 20?
Thanks, Ning