Нашел интересную статью, о том как в Quake 3 происходит вычисление обратного квадратного корня.
Вот что такое промышленное шаманство)
№ | Пользователь | Рейтинг |
---|---|---|
1 | tourist | 3993 |
2 | jiangly | 3743 |
3 | orzdevinwang | 3707 |
4 | Radewoosh | 3627 |
5 | jqdai0815 | 3620 |
6 | Benq | 3564 |
7 | Kevin114514 | 3443 |
8 | ksun48 | 3434 |
9 | Rewinding | 3397 |
10 | Um_nik | 3396 |
Страны | Города | Организации | Всё → |
№ | Пользователь | Вклад |
---|---|---|
1 | cry | 167 |
2 | Um_nik | 163 |
3 | maomao90 | 162 |
3 | atcoder_official | 162 |
5 | adamant | 159 |
6 | -is-this-fft- | 158 |
7 | awoo | 157 |
8 | TheScrasse | 154 |
9 | Dominater069 | 153 |
9 | nor | 153 |
Нашел интересную статью, о том как в Quake 3 происходит вычисление обратного квадратного корня.
Вот что такое промышленное шаманство)
Название |
---|
Это метод Ньютона нахождения корней уравнения. Довольно таки известная штука
Может и известная, но я не слышал ни разу. А вот шаманство с указателеми — реально черная магия. Не понимаю, почему минусют пост, крутой трюк на самом деле..
Вот эта строчка очень интересная
i = 0x5f3759df — ( i >> 1 );
Вот просто для всех корней берется некая константа.
Дак первое приближение берётся в виде линейной функции видимо. Не?
UPD: Наверное чуть хитрее т.к. там число из float конвертится и в старших битах будет экспонента. Т.е. приближение не по линейной функции а по экспоненциальной...
Магия же. Не просто так станут такую константу называть магической.
Вот квадратный корень из одной известной библиотечки: http://pastebin.com/hGY9YjXp.