В чем смысл ассемблерной конструкции (x86)?

x0 -= (x2 / x1) * x1 не ?

Декомпилируя очередную библиотеку

Но зачем???

наткнулся на странную конструкцию

Библиотека, случаем, не времён fdiv_bug?

fdiv_bug

movss

Хм...

выравнивание x2 до строго целого количества x1 в x2

https://www.felixcloutier.com/x86/DIVSS.html SSE1 не работает с целыми числами.

Я на первый коммент смотрел :) Обычно этот финт используется для целочисленной арифметики, чтобы выровнять чего-нибудь по страницам, например, или строго равномерно поделить картинку на части.

Тем не менее - суть та же самая, приведение x2 к такому значению, чтобы в x2 гарантировано было представимое в «single-precision floating-point» количество x1.

https://www.hex-rays.com/products/decompiler/v12_vs_v11.shtml запустить, он некоторые стандартные паттерны опознает

Оно и так представимо посредством присутствия в xmm0

На самом деле арифметика не целочисленная, следовательно xmm0/xmm1*xmm1==xmm0 В чем тогда подвох?

дайте файл посмотрим что скажет hexrays

Похоже на попытку оценить погрешность вычислений, так сказать, минимальную дельту при работе с числами сравнимыми с dword ptr [ebx-303BE0h]

На самом деле арифметика не целочисленная, следовательно xmm0/xmm1*xmm1==xmm0 В чем тогда подвох?

Если рассматривать совсем не целочисленный single precision float, например, то если xmm0 = 1e-23 а xmm1 = 1e23, то (xmm0/xmm1)*xmm1==0 а не xmm0. В других случаях может быть внезапно inf а не xmm0. Можешь проверить.

Но на самом деле я всё понял.

Типичный детект fdiv_bug:

4195835.0 - 3145727.0*(4195835.0/3145727.0) = 0    (Correct value)
4195835.0 - 3145727.0*(4195835.0/3145727.0) = 256  (Flawed Pentium)

В общем, это проверка на наличие FDIV_BUG скомпилённая с SSE :) в [ebx-303BE0h] должно лежать 3145727.0, а в xmm0 где-то ранее загружаться 4195835.0

В общем, это проверка на наличие FDIV_BUG скомпилённая с SSE :) в [ebx-303BE0h] должно лежать 3145727.0, а в xmm0 где-то ранее загружаться 4195835.0

Это чертовски похоже на истину, но, к сожалению, навряд ли. :( Более расширенный кусок кода выглядит так:

ucomiss xmm0, dword ptr [ebx-303BE0h] ;-7.81686689859073E-16
jbe loc_562DB1
movss xmm1, dword ptr [ebx-303BE0h]   ;-7.81686689859073E-16
movss xmm2,xmm0
divss xmm2,xmm1
mulss xmm1,xmm2
subss xmm0,xmm1
movss dword ptr [eax],xmm0
jmp loc_562DF9
loc_562DB1:
movss xmm1, dword ptr [ebx-303BDCh]   ;-5.77456881880694E-16
ucomiss xmm1,xmm0
jbe loc_562DF4
movss xmm1, dword ptr [ebx-303BE0h]   ;-7.81686689859073E-16
movss xmm2,xmm0
divss xmm2,xmm1
mulss xmm1,xmm2
addss xmm1,xmm0
movss dword ptr [eax],xmm1
jmp loc_562DF9

бинарник мы конечно не получим ? давайте я ромашек принесу и будем гадать на них

К сожалению, библиотеку полностью предоставить не могу. Да и не пойму, чем она может помочь? Возможно, я напутал с константами [ebx-303BE0h] = Hex 0x4C18E14E = Real -7.81686689859073E-16 [ebx-303BDCh] = Hex 0xCC18E14E = Real -5.77456881880694E-16 Или не?

а че толку гадать на ромашках если нельзя это посмотреть декомпилятором ?

https://www.h-schmidt.net/FloatConverter/IEEE754.html

Hex 0x4C18E14E = Real +4.00766E7
Hex 0xCC18E14E = Real -4.00766E7

Не очень понятно, как -7.81686689859073E-16 может хранится в dword. ближайшее число представимое в single precision float - -7.81686684714236017235e-16

Hex 0x4C18E14E = Real +4.00766E7
Hex 0xCC18E14E = Real -4.00766E7

Ай, спасибо, мил человек, это все ставит на место! Так решились все заботы, так узнали чё к чему! Единственный вопрос - это чего ж винхекс так в заблуждение вводит своим странным интерпретатором данных?

Если он действительно на 0x4C18E14E говорит, что число отрицательное, то он что-то совсем не то делает. Знаковый бит же сброшен.

Ну логика конструкции следующая:

if (x>4.00766E7) x-=x/4.00766E7*4.00766E7; else
{
   if (x<-4.00766E7) x+=x/4.00766E7*4.00766E7;
}

логика конструкции следующая

Я не понял, в чём эта логика заключается. Чего в итоге авторы хотели добиться.

Разъясняю для опоздавших на праздник. Данная конструкция - классический пример индийского кода в вакууме. Задумывалась проверка переполнения на +-40076600 (я думаю, многие догадались даже что это за волшебная цифра) но глупые индусы не учли что арифметика не целочисленная в итоге имеем такой безумный код

и что это за волшебная цифра?

Даю подсказку - для нормальных людей эта цифра равна 40 075 696 ;)

Длина экватора Земного шара равна 40 075 696 м.

Но мне непонятно, почему плюс-минус (а например не от нуля, до плюс длины, или плюс-минус ПОЛОВИНА длины)

https://en.wikipedia.org/wiki/Mercator_projection#/media/File:Cylindrical_Pro...

180 это не 360. А длина экватора это 360

180 это не 360. А длина экватора это 360

Так дикари-с! :D