Avatar ·

¿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) } ```