Что то не так с экспонентой в GAS?

@Examnes · Регистрация: 04.11.2017

Студворк — интернет-сервис помощи студентам

Я написал функцию вычисления exp(sin(x)) в виде ассемблерной вставки на C++. Я использовал компилятор G++, поэтому ассемблерные вставки написаны на GAS. Кажется никаких различий между MASM и GAS кроме синтаксиса в данном случае нет:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
#include <iostream>
#include <cmath>
using namespace std;
 
int main()
{
    float x, res, t;
    cout << "enter x: ";
    for (float x = -10; x < 100; x += 0.01)
    {
        __asm__(R"(
        finit
 
        fld %1 
        fsin
 
        fldl2e
        fmul
        
        fld     %%st
        frndint
        fsub    %%st(1), %%st
        fxch    %%st(1)
        f2xm1
        fld1
        fadd
        fscale
        fstp    %%st(1)
        fstp %4
        )"
            : "=m"(res)
            : "m"(x), "m"(four), "m"(five), "m"(t)
            :);
        cout << x << " " << t << endl;
    }
}

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

Как видно моя функция по чему - то не доходит до максимума, но в какой то момент резко опускается до минимума.
Помогите найти ошибку.

ФедосеевПавел · 20.11.2020, 19:13

Ещё раз сравните с исходным и проверьте двухоперандные инструкции - источник-приёмник в синтаксисах Intel и AT&T располагаются в разной последовательности.

А ещё лучше, приводить минимальную тестовую программу и значение, при котором видны различия в вычислениях.
Мне лень проводить за вас эту работу, поэтому - обойдётесь советом по проверке порядка операндов.

Добавлено через 35 минут
У меня в командной строке включён синтаксис Intel
Замена float на double увеличила точность

Windows Batch file
1
2
mingw32-g++.exe -Wall -fexceptions -g -masm=intel -m32  -c main.cpp -o obj\Debug\main.o
mingw32-g++.exe  -o bin\Debug\main.exe obj\Debug\main.o

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
#include <iostream>
#include <cmath>
using namespace std;
 
int main()
{
    double x, res, t, four, five;
    cout << "enter x: ";
    for (double x = -2; x < 2; x += 0.01)
    {
        asm(
        "finit \n\t"
 
        "fld %1 \n\t"
        "fsin \n\t"
 
        "fldl2e \n\t"
        "fmul \n\t"
 
        "fld     %%st \n\t"
        "frndint \n\t"
        "fsub    %%st(1), %%st \n\t"
        "fxch    %%st(1) \n\t"
        "f2xm1 \n\t"
        "fld1 \n\t"
        "fadd \n\t"
        "fscale \n\t"
        "fstp    %%st(1) \n\t"
        "fstp %4 \n\t"
 
            : "=m"(res)
            : "m"(x), "m"(four), "m"(five), "m"(t)
            :);
        cout << x << " " << exp(sin(x)) << " " << t - exp(sin(x)) << endl;
    }
}

Результат прогона (фрагмент)

Code
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
1.79 2.65401 1.96674e-016
1.8 2.64811 -2.79724e-017
1.81 2.64198 1.05384e-016
1.82 2.6356 -2.60209e-018
1.83 2.62898 1.41597e-016
1.84 2.62212 1.39428e-016
1.85 2.61503 1.73472e-017
1.86 2.60771 -2.77556e-016
1.87 2.60016 -6.09322e-017
1.88 2.59238 -1.91036e-016
1.89 2.58438 2.14672e-017
1.9 2.57616 5.72459e-017
1.91 2.56772 5.48606e-017
1.92 2.55907 -2.13805e-016
1.93 2.55021 -1.31839e-016
1.94 2.54115 1.98192e-016
1.95 2.53187 -7.84962e-017
1.96 2.5224 6.39679e-017
1.97 2.51273 -1.37911e-016
1.98 2.50287 -1.67401e-016
1.99 2.49282 -9.75782e-018

@Examnes · 20.11.2020, 21:04 **[ТС]**

Проблема действительно была в напрвлении. Спасибо, вот измененный код программы:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
#include <iostream>
#include <cmath>
using namespace std;
 
int main()
{
    float res;
    for (float x = 0; x < 100; x += 0.1)
    {
        float t;
        __asm__(R"(
        finit
        fld %1
        fsin
        fldl2e
        fmul
        fld     %%st(0)
        frndint
        fsubr   %%st(0), %%st(1)
        fld1
        fscale
        fxch 
        fstp %%st(0)
        fxch
        f2xm1
        fld1
        faddp
        fmul
        fstp %2
        )"
        : "=m"(res)
        : "m"(x), "m"(t)
        :);
        cout << exp(sin(x)) << " " << t << endl;
    }
}

Новые блоги и статьи Все статьи Все блоги /
Восстановить юзерскрипты Greasemonkey из бэкапа браузера damix 15.01.2026 Если восстановить из бэкапа профиль Firefox после переустановки винды, то список юзерскриптов в Greasemonkey будет пустым. Но восстановить их можно так. Для этого понадобится консольная утилита. . .	Изучаю kubernetes lagorue 13.01.2026 А пригодятся-ли мне знания kubernetes в России?	Сукцессия микоризы: основная теория в виде двух уравнений. anaschu 11.01.2026 https:/ / rutube. ru/ video/ 7a537f578d808e67a3c6fd818a44a5c4/	WordPad для Windows 11 Jel 10.01.2026 WordPad для Windows 11 — это приложение, которое восстанавливает классический текстовый редактор WordPad в операционной системе Windows 11. После того как Microsoft исключила WordPad из. . .
Classic Notepad for Windows 11 Jel 10.01.2026 Old Classic Notepad for Windows 11 Приложение для Windows 11, позволяющее пользователям вернуть классическую версию текстового редактора «Блокнот» из Windows 10. Программа предоставляет более. . .	Почему дизайн решает? Neotwalker 09.01.2026 В современном мире, где конкуренция за внимание потребителя достигла пика, дизайн становится мощным инструментом для успеха бренда. Это не просто красивый внешний вид продукта или сайта — это. . .	Модель микоризы: классовый агентный подход 3 anaschu 06.01.2026 aa0a7f55b50dd51c5ec569d2d10c54f6/ O1rJuneU_ls https:/ / vkvideo. ru/ video-115721503_456239114	Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ФедосеевПавел 06.01.2026 Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .