Вычислить кусочно-линейную функцию

@Tasl · Регистрация: 14.11.2021

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

Есть формула. Написал код на сопроцессоре. До этого там тоже вычисляется код. Но он верный.
Не могу понять почему не сходятся результаты. Помогите пожалуйста. И буду вам очень благодарен. Если вы сделаете это не вставками а модулем ассмеблерным. То есть ответ двумя файлами cpp и asm . Извините что много прошу просто у меня ни то ни другое не получается. Помогите пожалуйста. Буду очень признателен

Название: help.png
Просмотров: 35

Размер: 9.3 Кб

Название: help.png
Просмотров: 35

Размер: 9.3 Кб

C++

#define _USE_MATH_DEFINES
#define _CRT_SECURE_NO_WARNINGS
 
#include <iostream>
 
using namespace std;
 
 
double calc(double x)
{
    double buffer;
    int C2 = 2;
    int C6 = 6;
    __asm
    {
        finit
 
        fld x               // ST(0) = x
        fld ST              // ST(0) = x, ST(1) = x
 
        fmul ST, ST(1)      // ST(0) = x ^ 2, ST(1) = x
        fmul ST, ST         // ST(0) = x ^ 4, ST(1) = x //сразу нельзя?
 
 
        fld1                // ST(0) = 1, ST(1) = x ^ 4, ST(2) = x
        fsub ST(1), ST      // ST(0) = 1, ST(1) = x ^ 4 - 1, ST(2) = x
 
 
        fxch ST(1)          // ST(0) = x ^ 4 - 1, ST(1) = 1, ST(2) = x
        fabs                // ST(0) = abs(x ^ 4 - 1), ST(1) = 1, ST(2) = x
        fyl2x               // ST(0) = log2(abs(x ^ 4 - 1)), ST(1) = x
 
        fld1                // ST(0) = 1, ST(1) = log_2(abs(x ^ 4 - 1)), ST(2) = x
        fld st(2)           // ST(0) = x, ST(1) = 1, ST(2) = log_2(abs(x ^ 4 - 1)), ST(3) = x
        fsin                // ST(0) = sin(x), ST(1) = 1, ST(2) = log_2(abs(x ^ 4 - 1)), ST(3) = x
        fimul C2            // ST(0) = 2sin(x), ST(1) = 1, ST(2) = log_2(abs(x ^ 4 - 1)), ST(3) = x
        fmul ST, ST(2)      // ST(0) = 2sin(x)*log_2(abs(x ^ 4 - 1)), ST(1) = 1, ST(2) = log_2(abs(x ^ 4 - 1)), ST(3) = x
        faddp ST(1), ST     // ST(0) = 2sin(x)*log_2(abs(x ^ 4 - 1)) + 1 = (*), ST(1) = log_2(abs(x ^ 4 - 1)), ST(2) = x
 
        ; Вычисляем знаменатель
        fldpi               // ST(0) = PI, ST(1) = (*), ST(2) = log_2(abs(x ^ 4 - 1)), ST(3) = x
        fimul C2            // ST(0) = 2PI, ST(1) = (*), ST(2) = log_2(abs(x ^ 4 - 1)), ST(3) = x
 
        fxch ST(3)          // ST(0) = x, ST(1) = (*), ST(2) = log_2(abs(x ^ 4 - 1)), ST(3) = 2PI
        fimul C2            // ST(0) = 2x, ST(1) = (*), ST(2) = log_2(abs(x ^ 4 - 1)), ST(3) = 2PI
        fsin                // ST(0) = sin(2x), ST(1) = (*), ST(2) = log_2(abs(x ^ 4 - 1)), ST(3) = 2PI
        fldpi               // ST(0) = PI, ST(1) = sin(2x), ST(2) = (*), ST(3) = log_2(abs(x ^ 4 - 1)), ST(4) = 2PI
        fidiv C6            // ST(0) = PI/6, ST(1) = sin(2x), ST(2) = (*), ST(3) = log_2(abs(x ^ 4 - 1)), ST(4) = 2PI
        faddp ST(1), ST     // ST(0) = PI/6 + sin(2x), ST(1) = (*), ST(2) = log_2(abs(x ^ 4 - 1)), ST(3) = 2PI
        fdivp ST(3), ST     // ST(0) = (*), ST(1) = log_2(abs(x ^ 4 - 1)), ST(2) = 2PI/(PI/6 + sin(2x)) = (**)
        fxch ST(2)          // ST(0) = (**), ST(1) = log_2(abs(x ^ 4 - 1)), ST(2) = (*)
        fdivp ST(2), ST     // ST(0) = log_2(abs(x ^ 4 - 1)), ST(1) = (*)/(**) = (***)
 
        fld1                // ST(0) = 1, ST(1) = log_2(abs(x ^ 4 - 1)), ST(2) = (*)/(**) = (***)
        fpatan              // ST(0) = atan(log_2(abs(x ^ 4 - 1))/1), ST(1) = (***)
        fmul ST, ST         // ST(0) = atan^2(log_2(abs(x ^ 4 - 1))/1), ST(1) = (***)
        fsubp st(1), st     // ST(0) = (***) - atan^2()
                            // st(0) - result
    }
}
 
double condition(double x) {
    double upLimit = 0.5;
    int plus = 4;
    __asm {
        finit
        fldz // ST(0) = 0
        fld x // ST(0) = x, ST(1) = 0
        fcom
        fstsw ax
        sahf
        jb CASE_1
 
        xor ax, ax
        fld upLimit // ST(0) = 0,5 ; ST(1) = x, ST(2) = 0
        fcom
        fstsw ax
        sahf
        ja CASE_2
 
        //2^x + 3
        fxch ST(1)// ST(0) = x ; ST(1) = 0.5, ST(2) = 0
        f2xm1 // ST(0) = 2^x - 1 ; ST(1) = 0.5, ST(2) = 0
        fld plus //ST(0) = 4 , ST(1) = 2^x - 1 ; ST(2) = 0.5, ST(3) = 0
        fadd ST(1), ST(0) // ST(0) = 4 , ST(1) = 2^x + 3 ; ST(2) = 0.5, ST(3) = 0
        fxch ST(1) // ST(0) = 2^x + 3 , ST(1) = 4 ; ST(2) = 0.5, ST(3) = 0
        jmp THE_END
 
        CASE_2:
        fld1 //ST = 1,  ST(1) = 0,5 ; ST(2) = x, ST(3) = 0
        fxch ST(2)//ST = x, ST(1) = 0, 5; ST(2) = 1, ST(3) = 0
        fxch ST(1)//ST = 0,5, ST(1) = x; ST(2) = 1, ST(3) = 0
        fxch ST(2)//ST = 1, ST(1) = x; ST(2) = 0,5, ST(3) = 0
 
        fpatan // ST = atan(x),  ST(1) = x ; ST(2) = 0,5, ST(3) = 0
        fld1;//ST(0) = 1, ST = atan(x),  ST(1) = x ; ST(2) = 0,5, ST(3) = 0
        fmul ST(0), ST(1)// ST(0) = atan(x), ST = atan(x), ST(1) = 0, 5; ST(2) = x, ST(3) = 0
        fmul ST(0), ST(1)// 
        fmul ST(0), ST(1) // ST(0) = atan(x)^3, ST = atan(x), ST(1) = 0, 5; ST(2) = x, ST(3) = 0
        jmp THE_END
 
        CASE_1:
        fcos // ST(0) = cos(x), ST(1) = 0
        fmul ST(0), ST(0)
        jmp THE_END
        THE_END:
 
 
 
    }
}
 
int main() {
 
    double x, resultC, resultA, resultD;
    cout << "Input x -> ";
 
    cin >> x;
 
    cout << endl << "The first task:" << endl;
 
    resultC = (2 * sin(x) * (log2(abs(pow(x, 4) - 1))) + 1) / ((2 * M_PI) / (sin(2 * x) + (M_PI / 6))) - pow(atan(log2(abs(pow(x, 4) - 1))), 2);
 
    cout << "(C++) result -> " << resultC << endl;
 
    resultA = calc(x);
 
    cout << "(Assembler) result -> " << resultA << endl;
 
    resultD = condition(x);
 
    cout << "task 2 -> " << resultD << endl;
 
    return 0;
}

@Tasl · 23.12.2021, 18:17 **[ТС]**

Кто нибудь, помогите пожалуйста

ФедосеевПавел · 23.12.2021, 18:45

Ветвление
Электронный учебник глава 11

Инструкции сравнения берите те, которые сразу изменяют флаги CPU, и выполняйте условный переход, как при сравнении беззнаковых чисел - т.е. один из jb/jbe/je/jae/ja

Для сравнения можете загружать константы из переменных или fldz (загрузка 0) или последовательность для 0,5

Assembler

fld1
fld1
fadd st(0), st(0)
fdivp st(1), st(0)

@Left-guy · 26.12.2021, 12:56

ФедосеевПавел, Павел, каким образлм брать com, если я работаю с сопроцессором, а там fcom только

ФедосеевПавел · 26.12.2021, 13:31

О каком com речь?

Новые блоги и статьи
Какой локальный веб-сервер выбрать InfoMaster 19.01.2025 В современной веб-разработке локальные веб-серверы играют ключевую роль, предоставляя разработчикам надежную среду для создания, тестирования и отладки веб-приложений без необходимости использования. . .	Почему планшеты и iPad уже не так популярны, как раньше InfoMaster 19.01.2025 Эра революционных инноваций История планшетных компьютеров началась задолго до того, как эти устройства стали привычными спутниками нашей повседневной жизни. В начале 1990-х годов появились первые. . .	Как самому прошить BIOS ноутбука InfoMaster 19.01.2025 BIOS (Basic Input/ Output System) представляет собой важнейший компонент любого компьютера или ноутбука, который обеспечивает базовое взаимодействие между аппаратным и программным обеспечением. . .	Какой Linux выбрать для домашнего компьютера InfoMaster 19.01.2025 Современные реалии выбора операционной системы В современном мире выбор операционной системы для домашнего компьютера становится все более важным решением, которое может существенно повлиять на. . .	Как объединить два словаря одним выражением в Python InfoMaster 19.01.2025 В мире программирования на Python работа со словарями является неотъемлемой частью разработки. Словари представляют собой мощный инструмент для хранения и обработки данных в формате "ключ-значение". . . .	Как без исключения проверить существование файла в Python InfoMaster 19.01.2025 При разработке программного обеспечения на Python часто возникает необходимость проверить существование файла перед выполнением операций с ним. Это критически важная задача, которая помогает избежать. . .
Как определить, содержит ли строка подстроку в JavaScript InfoMaster 19.01.2025 При разработке веб-приложений часто возникает необходимость выполнять различные операции со строками, среди которых особое место занимает поиск подстрок. JavaScript предоставляет несколько встроенных. . .	Что такое метаклассы в Python InfoMaster 19.01.2025 Метаклассы в Python представляют собой один из самых мощных и одновременно сложных механизмов языка, позволяющий программистам контролировать процесс создания классов. По своей сути, метакласс. . .	Как удалить свойство из объекта JavaScript InfoMaster 19.01.2025 В современной веб-разработке объекты JavaScript играют фундаментальную роль в организации и структурировании данных. Они представляют собой контейнеры, которые хранят связанные данные и. . .	Какая разница между String и string в C# InfoMaster 19.01.2025 В языке программирования C# существует интересная особенность: для работы со строками можно использовать как String, так и string. Эта двойственность часто вызывает вопросы у разработчиков, особенно. . .	Как в Git откатить репозиторий к предыдущему коммиту InfoMaster 19.01.2025 В современной разработке программного обеспечения система контроля версий Git стала неотъемлемой частью рабочего процесса, предоставляя разработчикам мощные инструменты для управления изменениями в. . .	Как работают замыкания (closure) в JavaScript InfoMaster 19.01.2025 В мире современной веб-разработки замыкания (closures) представляют собой один из фундаментальных концептов языка JavaScript, который часто вызывает затруднения у начинающих разработчиков, но при. . .

@Tasl 0 / 0 / 0 Регистрация: 14.11.2021 Сообщений: 20
	23.12.2021, 18:17 [ТС]	2
	Кто нибудь, помогите пожалуйста 0

@Left-guy 0 / 0 / 0 Регистрация: 26.12.2021 Сообщений: 2
	26.12.2021, 12:56	4
	ФедосеевПавел, Павел, каким образлм брать com, если я работаю с сопроцессором, а там fcom только 0

ФедосеевПавел Модератор 8561 / 4409 / 1655 Регистрация: 01.02.2015 Сообщений: 13,699 Записей в блоге: 9
	26.12.2021, 13:31	5
	О каком `com` речь? 0

Опции темы