Нашел интересную статью, о том как в Quake 3 происходит вычисление обратного квадратного корня.
Вот что такое промышленное шаманство)
№ | Пользователь | Рейтинг |
---|---|---|
1 | tourist | 3803 |
2 | jiangly | 3707 |
3 | Benq | 3627 |
4 | ecnerwala | 3584 |
5 | orzdevinwang | 3573 |
6 | Geothermal | 3569 |
6 | cnnfls_csy | 3569 |
8 | Radewoosh | 3542 |
9 | jqdai0815 | 3532 |
10 | gyh20 | 3447 |
Страны | Города | Организации | Всё → |
№ | Пользователь | Вклад |
---|---|---|
1 | awoo | 163 |
2 | maomao90 | 160 |
3 | adamant | 159 |
4 | maroonrk | 152 |
5 | -is-this-fft- | 150 |
6 | atcoder_official | 148 |
6 | SecondThread | 148 |
8 | nor | 147 |
9 | TheScrasse | 146 |
10 | Petr | 144 |
Нашел интересную статью, о том как в Quake 3 происходит вычисление обратного квадратного корня.
Вот что такое промышленное шаманство)
Название |
---|
Это метод Ньютона нахождения корней уравнения. Довольно таки известная штука
Может и известная, но я не слышал ни разу. А вот шаманство с указателеми — реально черная магия. Не понимаю, почему минусют пост, крутой трюк на самом деле..
Вот эта строчка очень интересная
i = 0x5f3759df — ( i >> 1 );
Вот просто для всех корней берется некая константа.
Дак первое приближение берётся в виде линейной функции видимо. Не?
UPD: Наверное чуть хитрее т.к. там число из float конвертится и в старших битах будет экспонента. Т.е. приближение не по линейной функции а по экспоненциальной...
Магия же. Не просто так станут такую константу называть магической.
Вот квадратный корень из одной известной библиотечки: http://pastebin.com/hGY9YjXp.