Помощь - Поиск - Пользователи - Календарь
Полная версия: ассемблерная вставка в Delphi
Форум «Всё о Паскале» > Современный Паскаль и другие языки > Ассемблер
marwell
24.05.2012 16:17
доброго дня
необходимо вычислить интеграл методом прямоугольников (ассемблерная вставка в Delphi или в паскаль)
сам метод мне понятен
но, на каждом шаге мне нужно будет вычислять значение функции в определенной точке. хотел было использовать макросы
но не тут то было, ни паскаль ни Delphi макросы не принимают
может я неправильно описываю эти макросы?
var n,a,b: integer;
    res, h, x: real;
begin
readln(n);
a:=0;
b:=1;
asm
 fild b
 fisub a
 fidiv n
 call funct
 funct proc
...

[Error] Project1.dpr(18): Undeclared identifier: 'funct'
думал что макросы описываются след образом:
<имя> macro парам1, парам2... парам N
...
...набор команд
...
endm
IUnknown
24.05.2012 16:43
MACRO встроенным ассемблером не поддерживается. Поддерживаются инструкции, а не препроцессор. Используй функции.
marwell
24.05.2012 18:07
var n,a,b: integer;
    res, h, x: real;
begin
readln(n);// ввод количества частичных отрезков
a:=0;  // a и b - границы
b:=1;
x:=0;  // частичная сумма
res:=0;
asm
 fild b    // загрузка b на вершину стека
 fisub a   // b-a
 fidiv n   // :n
 fst h     // значение из вершины в h
 mov ecx,n
 @m:  fld x
      fadd h
      fst x
      fptan
      fdivr st(1),st
      fadd res
      fst res
      loop @m;
 end
writeln(res:5:4);
readln;
end.

почему происходит деление на ноль? я же прибавляю перед делением h
fld x 
fadd h
marwell
24.05.2012 18:28
аа, дело в 
fdivr st(1),st

кажется там на каждом шаге номер меняется?
на первом шаге надо st(5) делить на st(6)? или я неправильно понял структуру стека?
IUnknown
24.05.2012 18:29
А кто тебе сказал, что происходит деление на 0? smile.gif

Посмотри внимательно на скриншот, и скажи, что будет, когда я нажму на F7 для выполнения следующего шага (операции fdivr), и почему?
Нажмите для просмотра прикрепленного файла

Вот ответ: (Показать/Скрыть)
marwell
24.05.2012 18:33
сорри, успел отредактировать свое предыдущее сообщение до того как увидел ваш ответ
а насчет прикрепленной картинки, как то не очень понятно wacko.gif

Добавлено через 6 мин.
на каком то шаге тангенс равен 0?
IUnknown
24.05.2012 18:41
Посмотри на содержимое стека FPU. Сколько там осталось свободных регистров? Нисколько. Вот поэтому и вылет: переполнение стека FPU. Чистить стек нужно.
marwell
24.05.2012 19:00
спасибо!
var n,a,b: integer;
    res, h, x: real;
begin
readln(n);// ââîä êîëè÷åñòâà ÷àñòè÷íûõ îòðåçêîâ
a:=0;  // a è b - ãðàíèöû
b:=1;
x:=0;  // ÷àñòè÷íàÿ ñóììà
res:=0;
asm
 fild b    // çàãðóçêà b íà âåðøèíó ñòåêà
 fisub a   // b-a
 fidiv n   // :n
 fstp h     // çíà÷åíèå èç âåðøèíû â h
 mov ecx,n
 @m:  fld x
      fadd h
      fst x
      fptan
      fdivr st(1),st
      fadd res
      fstp res
      finit
      loop @m;
 end
writeln(res:5:4);
readln;
end.

осталось понять почему считает неправильно. res почему то всегда равен введенному n
IUnknown
24.05.2012 19:27
Функция какая у тебя, тангенс или котенгенс? Если тангенс - то вот так:

   asm
      fild b
      fisub a
      fidiv n
      fst h

      mov ecx, n
      fld x
   @m:
      fadd st(1)
      fld st(0) // Дублируем значение Х
      fptan
      fdivp st(1),st
      fadd res
      fstp res
      loop @m;
   end
   writeln(res:5:4);

(логика очень проста: в стеке держим только приращение в ST1, и текущее значение X в ST0, все остальное сразу после вычисления убираем). Если функция - КОтангенс - поменяй fdivp на fdivrp
marwell
24.05.2012 20:00
Цитата(IUnknown @ 24.05.2012 15:27) *

Функция какая у тебя, тангенс или котенгенс? Если тангенс - то вот так:

   asm
      fild b
      fisub a
      fidiv n
      fst h

      mov ecx, n
      fld x
   @m:
      fadd st(1)
      fld st(0) // Дублируем значение Х
      fptan
      fdivp st(1),st
      fadd res
      fstp res
      loop @m;
   end
   writeln(res:5:4);

(логика очень проста: в стеке держим только приращение в ST1, и текущее значение X в ST0, все остальное сразу после вычисления убираем). Если функция - КОтангенс - поменяй fdivp на fdivrp

спасибо! да, тангенс нужен. И логика понятна.
Только Delphi ругается на эту строчку
fadd st(1)

говорит [Error] Project1.dpr(25): Invalid combination of opcode and operands
IUnknown
24.05.2012 20:29
FPC не ругается, в Дельфи сейчас посмотрю.

Update
Да, замени ту строку вот на эту:
      fadd st, st(1)
marwell
24.05.2012 20:49
Цитата(IUnknown @ 24.05.2012 16:29) *

FPC не ругается, в Дельфи сейчас посмотрю.

Update
Да, замени ту строку вот на эту:
      fadd st, st(1)


сделал. Запускается без ошибок. Вот только ответ... Ответ должен быть приблизительно 0.261168 да? А прога при n=100 выдает 62,3434. Чем больше n, тем дальше ответ растет wacko.gif
IUnknown
24.05.2012 21:41
Не верю smile.gif Ты просто не доделал метод прямоугольников. Вот код, вычисляющий значение и на Ассемблере, и на Паскале:

program Project1;

{$APPTYPE CONSOLE}
uses math;

var n, a, b: integer;
    res, h, x: real;

    i : Integer;
begin
   // readln(n);
   n := 15000;

   a := 0;
   b := 1;
   x := 0;
   res := 0;
   asm
      fild b
      fisub a
      fidiv n
      fst h

      mov ecx, n
      fld x
   @m:
      fadd st, st(1)
      fld st(0)
      fptan
      fdivp st(1), st
      fadd res
      fstp res
      loop @m;

      fld res
      fmul st, st(2)
      fstp res
   end
   writeln(res:5:4);

   // Test
   h := (b - a) / n;
   x := 0;
   res := 0;
   for i := 1 to n do
   begin
      x := x + h;
      res := res + tan(x);
   end;
   res := res * h;
   writeln(res:5:4);

   readln;
end.
Чем выше N - тем выше точность. А так значения абсолютно одинаковые.
marwell
24.05.2012 21:57
Цитата(IUnknown @ 24.05.2012 17:41) *

Не верю smile.gif Ты просто не доделал метод прямоугольников. Вот код, вычисляющий значение и на Ассемблере, и на Паскале:

program Project1;

{$APPTYPE CONSOLE}
uses math;

var n, a, b: integer;
    res, h, x: real;

    i : Integer;
begin
   // readln(n);
   n := 15000;

   a := 0;
   b := 1;
   x := 0;
   res := 0;
   asm
      fild b
      fisub a
      fidiv n
      fst h

      mov ecx, n
      fld x
   @m:
      fadd st, st(1)
      fld st(0)
      fptan
      fdivp st(1), st
      fadd res
      fstp res
      loop @m;

      fld res
      fmul st, st(2)
      fstp res
   end
   writeln(res:5:4);

   // Test
   h := (b - a) / n;
   x := 0;
   res := 0;
   for i := 1 to n do
   begin
      x := x + h;
      res := res + tan(x);
   end;
   res := res * h;
   writeln(res:5:4);

   readln;
end.
Чем выше N - тем выше точность. А так значения абсолютно одинаковые.

большое и уже не первое спасибо smile.gif
Это текстовая версия — только основной контент. Для просмотра полной версии этой страницы, пожалуйста, нажмите сюда.