And so are many other sorts, I bet. The libbsd heapsort can’t sort an array that is larger than 2 Gb. This is probably because they have a line somewhere like 

 if (2*i+1 > size)

and with a 2Gb region, 2*i will overflow and return a negative number. It is a very strange experience to discover new bugs in 30 year old code.