Raiz quadrada inversa rápida incomum de John Carmack (Quake III)

John Carmack tem uma função especial no código-fonte do Quake III que calcula a raiz quadrada inversa de um float, 4x mais rápido que o normal (float)(1.0/sqrt(x)), incluindo uma 0x5f3759dfconstante estranha . Veja o código abaixo. Alguém pode explicar linha por linha o que exatamente está...