¿Se puede implementar la raíz cuadrada inversa rápida de Quake en Rust?
¡Hola a todos! Estoy aprendiendo Rust y me topé con la legendaria función InvSqrt() del juego Quake III Arena, que calcula 1/sqrt(x) increíblemente rápido mediante trucos de bits. En C/C++ se ve así:
```c
float Q_rsqrt(float number) {
long i;
float x2, y;
const float threehalfs = 1.5F;
x2 = number * 0.5F;
y = number;
i = * ( long * ) &y; // type punning malvado
i = 0x5f3759df - ( i >> 1 ); // constante mágica
y = * ( float * ) &i;
y = y * ( threehalfs - ( x2 * y * y ) ); // una iteración de Newton
return y;
}
```
El problema es que en Rust, ese casteo directo de punteros (`*(long*)&y`) se considera comportamiento indefinido (UB) y está prohibido. Intenté hacerlo con `transmute`:
```rust
fn inv_sqrt(x: f32) -> f32 {
let x2 = x * 0.5;
let mut y = x;
let mut i: i32 = unsafe { std::mem::transmute(y) };
i = 0x5f3759df - (i >> 1);
y = unsafe { std::mem::transmute(i) };
y * (1.5 - x2 * y * y)
}
```