Алгоритм численного поиска корня уравнения 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

Новые блоги и статьи Все статьи Все блоги /
Модель микоризы: классовый агентный подход 3 anaschu 06.01.2026 aa0a7f55b50dd51c5ec569d2d10c54f6/ O1rJuneU_ls https:/ / vkvideo. ru/ video-115721503_456239114	Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ФедосеевПавел 06.01.2026 Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .	Модель микоризы: классовый агентный подход 2 anaschu 06.01.2026 репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .	Расчёт токов в цепи постоянного тока igorrr37 05.01.2026 / * Дана цепь постоянного тока с сопротивлениями и напряжениями. Надо найти токи в ветвях. Программа составляет систему уравнений по 1 и 2 законам Кирхгофа и решает её. Последовательность действий:. . .
Новый CodeBlocs. Версия 25.03 palva 04.01.2026 Оказывается, недавно вышла новая версия CodeBlocks за номером 25. 03. Когда-то давно я возился с только что вышедшей тогда версией 20. 03. С тех пор я давно снёс всё с компьютера и забыл. Теперь. . .	Модель микоризы: классовый агентный подход anaschu 02.01.2026 Раньше это было два гриба и бактерия. Теперь три гриба, растение. И на уровне агентов добавится между грибами или бактериями взаимодействий. До того я пробовал подход через многомерные массивы,. . .	Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Programma_Boinc 28.12.2025 Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / ********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / ********/ gallery/ bKBkQFf Пост отсюда. . .	Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Нашел на реддите интересную статью под названием Anyone know where to get a free Desktop or Laptop? Ниже её машинный перевод. После долгих разбирательств я наконец-то вернула себе. . .

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