Алгоритм численного поиска корня уравнения f (x) = 0 методом хорд

@DarkBoy9 · Регистрация: 06.03.2018

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

Из картинки: -1.5052 - это ожидаемый результат
У меня вывод должен действовать в StartIters, но моя программа прерывается внизу в function. Не могу понять почему она не возвращается наверх в StartIters

Assembler
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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
%include "io.inc" 
section .data
    x0 dq -2.0 
    a dq 3.0 
    x1 dq 2.0
    fmt db '%lf' , 10,13,0
section .bss
    log resq 1
    fx0 resq 1
    fxi resq 1
    xi resq 1
section .text
global CMAIN
CMAIN:
    mov ebp, esp; for correct debugging
    mov ecx, 11
    fld qword [x1]
    fstp qword [xi]
    finit
    fld qword [x0]
    fld qword [a]
    fsub st0, st1
    fabs
    fstp qword [log]
    push dword [log+4]
    push dword [log]
    call flog
    fadd st0, st1
    fstp qword [fx0]
StartIters:
    call function
    loop StartIters
    push dword [xi + 4]
    push dword [xi]
    push fmt
    call printf
    add esp, 12
    ret
flog: ; ln x
    fld qword [esp + 4]
    fldln2
    fxch
    fyl2x
    ret
iter:
    finit
    fld qword [xi]
    fld qword [fxi]
    fld qword [fx0]
    fsubr st1
    fld qword [x0]
    fxch st3
    fst st4
    fsub st0, st3
    fmul st0, st2
    fdiv st0, st1
    fsubr st0, st4
    fstp qword [xi]
    ret
function:
    finit
    fld qword [xi]
    fld qword [a]
    fsub st0, st1
    fabs
    fstp qword [log]
    push dword [log+4]
    push dword [log]
    call flog
    fadd st0, st1 
    fstp qword [fxi]
    call iter
    ret

ФедосеевПавел · 11.10.2019, 19:58

Попробуйте выполнить по шагам в отладчике ollydbg или в x64dbg. Увидите состояние стека FPU и всех регистров.

А так - вижу сброс в исходное состояние FPU несколько раз в программе - это неправильно. Для сохранения контекста применяйте fsave+frstor, как в примере в функции F Вычислить определенный интеграл методом трапеций

Добавлено через 3 минуты
При выходе из процедуры flog не восстанавливается стек.

@DarkBoy9 · 16.10.2019, 17:16 **[ТС]**

Я примерно понял о чем вы говорите. Я смог решить проблему с помощью "условия", но приобрел больше проблем. Теперь при выводе когда условие счетчика не выполняется он идет после вывода в функцию где считается логарифм, и проходит ее ровно 11 раз, превращая тем самым мой стек в хаос, но выводя правильный ответ. Преподаватель сказал убрать jmp StartIteration. И оптимизировать код. Можете помочь с этим?

Assembler
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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
%include "io.inc" 
section .data
    x0 dq -2.0 
    a dq 3.0 
    x1 dq 2.0
    fmt db '%lf' , 10,13,0
section .bss
    log resq 1
    fx0 resq 1
    fxi resq 1
    xi resq 1
section .text
global CMAIN
CMAIN:
    mov ebp, esp; for correct debugging
    mov ecx, 11
    fld qword [x1]
    fstp qword [xi]
    finit
    fld qword [x0]
    fld qword [a]
    fsub st0, st1
    fabs
    fstp qword [log]
    push dword [log+4]
    push dword [log]
    call flog
    add esp,8
    fadd st0, st1
    fstp qword [fx0]
startiteration:
    cmp ecx, 0
    jnz AFTER
    push dword [xi + 4]
    push dword [xi]
    push fmt
    call printf
    add esp, 12
    fstp st0
    fstp st0
    fstp st0
    fstp st0
    fstp st0
    ret
AFTER:
    sub ecx, 1
    call function
flog: ; ln x
    fld qword [esp + 4]
    fldln2
    fxch
    fyl2x
    ret
iteration:
    finit
    fld qword [xi]
    fld qword [fxi]
    fld qword [fx0]
    fsubr st1
    fld qword [x0]
    fxch st3
    fst st4
    fsub st0, st3
    fmul st0, st2
    fdiv st0, st1
    fsubr st0, st4
    fstp qword [xi]
    ret
function:
    finit
    fld qword [xi]
    fld qword [a]
    fsub st0, st1
    fabs
    fstp qword [log]
    push dword [log+4]
    push dword [log]
    call flog
    add esp,8
    fadd st0, st1 
    fstp qword [fxi]
    call iteration
    jmp startiteration
    ret

Программирую в SASM

Новые блоги и статьи Все статьи Все блоги /
Чистый C и C++ в Qt Creator: ставим среду без VPN и настраиваем MinGW/CMake, без фреймворка Qt 8Observer8 05.04.2026 Среду разработки Qt Creator можно установить без фреймворка Qt. Есть отдельный репозиторий для этой среды: https:/ / github. com/ qt-creator/ qt-creator, где можно скачать установщик, на вкладке Releases:. . .	AkelPad-скрипты, структуры, и немного лирики.. testuser2 05.04.2026 Такая программа, как AkelPad существует уже давно, и также давно существуют скрипты под нее. Тем не менее, прога живет, периодически что-то не спеша дополняется, улучшается. Что меня в первую очередь. . .	Отображение реквизитов в документе по условию и контроль их заполнения Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеСпецтехники", разработанного в конфигурации КА2. Данный документ берёт данные из другого нетипового документа. . .	Фото всей Земли с борта корабля Orion миссии Artemis II kumehtar 04.04.2026 Это первое подобное фото сделанное человеком за 50 лет. Снимок называют новым вариантом легендарной фотографии «The Blue Marble» 1972 года, сделанной с борта корабля «Аполлон-17». Новое фото. . .
Вывод диалогового окна перед закрытием, если документ не проведён Maks 04.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать программный контроль на предмет проведения документа. . .	Программный контроль заполнения реквизитов табличной части документа Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: 1. Реализовать контроль заполнения реквизита. . .	wmic не является внутренней или внешней командой Maks 02.04.2026 Решение: DISM / Online / Add-Capability / CapabilityName:WMIC~~~~ Отсюда: https:/ / winitpro. ru/ index. php/ 2025/ 02/ 14/ komanda-wmic-ne-naydena/	Программная установка даты и запрет ее изменения Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: при создании документов установить период списания автоматически. . .

ФедосеевПавел Модератор 8665 / 4502 / 1670 Регистрация: 01.02.2015 Сообщений: 13,926 Записей в блоге: 13
	11.10.2019, 19:58
	Попробуйте выполнить по шагам в отладчике ollydbg или в x64dbg. Увидите состояние стека FPU и всех регистров. А так - вижу сброс в исходное состояние FPU несколько раз в программе - это неправильно. Для сохранения контекста применяйте fsave+frstor, как в примере в функции F Вычислить определенный интеграл методом трапеций Добавлено через 3 минуты При выходе из процедуры flog не восстанавливается стек. 0