Feda /dev/ - No Javascript allowed

https://www.campusmvp.es/recursos/post/la-linea-de-codigo-que-cambio-la-industria-de-los-graficos-3d.aspx

#52529 Si me permites una corrección: los aparejadores no usan la grúa. Son «arquitectos técnicos», y una de sus funciones es llevar el control material de la obra.

#52564 si no te has molestado en leerlo estás que me voy a molestar yo en explicártelo. En el artículo lo dicen

#52568 #52567 Lo único que dice el artículo es:

Esa fórmula es una manera aparentemente "mágica" de calcular la raíz cuadrada inversa aproximada, una operación común para determinar distancias entre vectores en modelado 3D en tiempo real.

Yo también puedo poner un hola mundo y decir que de forma mágica calcula si te has comido las tres pollas anuales o aún no.

O a lo mejor es que vais de desus de la vida pero tampoco sabéis explicar esa línea de código, no pasa nada.

edit: En la wiki lo explica un poco mejor:

Básicamente el numerito 0x5f3759df es una constante que representa la raiz cuadrada de 2 elevado a 127, que estaría bien que alguien explicase por qué este número y no la raíz cuadrada de mi abuela en bicicleta por ejemplo.

Por lo que dice, se pasa un número en coma flotante de 32 bits y coge la mitad de use número, entonces coge esos bits y los trata en formato integer en vez de float, nosequé historias de un bit a la derecha y entonces le resta la constante de mi abuela en bicicleta.

Pero para mí sigue siendo todo algún tipo de magia negra oscura.

edit2: Usando a chatGPT:
i = 0x5f3759df - ( i >> 1 ); // what the fuck?
i tiene un valor que es el que quieren calcular la raiz inversa cuadrada.
La constante ya lo hemos explicado.
i >> 1 lo que hace es desplazar un bit a la derecha el valor de i, como está en binario, esto equivale a dividirlo por 2

Por tanto, constante - i/2 --> un número cercano a la raiz cuadrada inversa de i.

#52569

Pero para mí sigue siendo todo algún tipo de magia negra oscura.

Por eso en su caso sí se puede decir que es una línea que ha cambiado el mundo y tú solo puedes memear con un hola mundo

#52572 El código en sí lo que hace es una única iteración del método de Newton sobre la función 1 / x ** 2 - c, con c = número para el que quieres saber la inversa de la raíz cuadrada, de hecho con la fórmula en la mano, si f(x) = 1 / x ** 2 - c

x' = x - f(x) / f'(x) es (x/2) * (3 - c * x ** 2), que es la línea

 y * ( threehalfs - ( x2 * y * y ) )

La fumada con bits es para encontrar un x inicial tal que con una única iteración ya te acerques lo bastante al resultado.

Siendo francos no se cómo cojones se le ocurrió a nadie tal número. La resta de i como integer shiftado a la derecha se me ocurre que es porque en IEE754 los números están representado como exponente y mantisa, y con eso divides el exponente por 2, lo cual parece bastante razonable.

La constante mágica ya sí que confieso que se escapa a mi entendimiento, en su día en la uni intenté entender por qué y no llegué muy lejos xD

La idea de la inversa entiendo que es porque la puedes calcular dividiendo únicamente por constantes, que se puede expresar como un producto, y así te ahorras dividir en coma flotante, que es bastante lento, y en la época probablemente te jodía la performance varios órdenes de magnitud.

#52576Seyriuu:

desplazas bit a la derecha y obtienes la mitad del número.

y sabes por qué se hace así?

O por qué en esta parte del código multiplica por 0.5 para dividir entre 2?

long i;
float x2, y;
const float threehalfs = 1.5F;

x2 = number * 0.5F;

#52580 pero mira el contexto:
1) usuario pone la línea magica que según el texto hace tal
2) yo pregunto eso como funciona
3) se me responde diciendo que ya lo pone en el artículo
4) yo hago una parodia del artículo para explicar que el artículo no explica la línea de código, explica el resultado final

Desde mi perspectiva, los soberbios habéis sido vosotros, yo desde un primer momento he dejado claro que NO sabía como funcionaba esa línea de código y posteriormente he investigado usando wikipedia y chatgpt.

Y sí me había leído el artículo, que nunca llega a explicar qué hacía esa línea en concreto, si no el resultado final.

#52582 no lo estaría, pero la duda original era entender qué hacía esa línea de código (restar una constante a la mitad de un valor, lo cual es el resultado aproximado de su raiz cuadrada inversa) tampoco se merecía la respuesta que ha tenido.

No podemos ser rudos y luego quejarnos de que nos contestan con rudeza, la verdad.

#52583 te pido disculpas si he estado borde contigo, un saludo

#52584sergioRG:

Me da a entender que tú tampoco sabías lo que hacía esa función, y que de hecho ni te acordabas del método de Newton.

lo he mirado después de leer el artículo en otra página que lo explican bien y cómo llegan al método de newton, pero no lo he puesto para que nadie me saltase al cuello la verdad

Eran estos dos liks:

https://breq.dev/2021/03/17/5F3759DF

https://betterexplained.com/articles/understanding-quakes-fast-inverse-square-root/