How to become an author
Search
Write a publication
Pull to refresh

Comments 25

Лучше уж сразу назвать серию постов «Алгоритмы древнего Вавилона»
Total votes 26: ↑25 and ↓1+24
UFO just landed and posted this here
Хорошо. Ничего критичного для жизни людей не программирует.
Total votes 4: ↑2 and ↓20
while(fabs(a-b)<EPS) тут наверное обратное условие,
преобразуем, пока разница между сторонами больше ε?
Total votes 5: ↑5 and ↓0+5
Точно! Очень извиняюсь! Спасибо! Исправил
Total votes 10: ↑5 and ↓50
Ожидал увидеть очередную статью про 0x5f3759df.
Раз уж решили искать корень перебором, почему бы обычную дихотомию не использовать? Или бином Ньютона, на худой конец.
Total votes 5: ↑5 and ↓0+5
По-моему это называется итерационный аналитический алгоритм извлечения квадратного корня на основе формулы Герона.
Total votes 7: ↑7 and ↓0+7
А можете объяснить почему стандартный sqrt,

sqrt13:
double sqrt13(double n)
{
   __asm {
     fld n
     fsqrt
   }
}


и sqrt14

double inline __declspec (naked) __fastcall sqrt14(double n)
{
	_asm fld qword ptr [esp+4]
	_asm fsqrt
	_asm ret 8
}


работают разное время? Каждый следующий раза в два быстрее предыдущего (при этом точность одинаковая).

Очень уж не хочется верить в то, что в наши дни нужно писать код как в примере 14, который судя по всему ещё не будет работать на x64. Как-то я больше доверял оптимизациям компилятора.
Total votes 1: ↑1 and ↓0+1
Ну по идее обе функции должны бы заинлайниться. Да и пролог/эпилог хороший компилятор должен сам суметь решить когда выкинуть.

Total votes 2: ↑0 and ↓2-2
анализировать __asm{} блок компилятору запрещено. именно поэтому 13й вариант имеет все необходимые стандартные заголовки и не имеет права инлайниться — мало ли что ты там в asm{} блоке хочешь делать.

кстати, 14й тоже НЕ инлайнится, несмотря на __inline. иначе функция работала бы СОВСЕМ неправильно из-за ret 8 унутре
Total votes 4: ↑4 and ↓0+4
А я ожидал увидеть алгоритмы для компьютеров — быстрые, с маленькой погрешностью, эффективные.
Total votes 9: ↑9 and ↓0+9
У автора после while фигурных скобок не хватает :)
Total votes 1: ↑1 and ↓0+1
И ещё вдогонку, в приведённой функции после того как поставить скобки после условия в while, всё разно есть ещё одна ошибка,
Она связана с неучтённой погрешностью float.

Вход 1e18: codepad.org/qK4Bs558 — выполняется
Вход 1e20: codepad.org/7Rwyo0gh — бесконечный цикл
Вход 1e22: codepad.org/Ry9iRGbZ — выполняется
Total votes 3: ↑3 and ↓0+3
У меня при компиляции на gcc все представленные входы нормально считаются. Там float какого размера используется?
Total votes 2: ↑0 and ↓2-2
4 байта.
Проявление этой ошибки зависит от инструкций и их порядка. И будет скорее всего разное при разных настройках компилятора.
Можешь попробовать повычислять с числами такого же порядка, в цикле.
Total votes 1: ↑1 and ↓0+1
Sign up to leave a comment.

Articles

Show the best of all time

Your account

Sections

Information

Services

Support
© 2006–2024, Habr