дельфи с asm-вставками

[alexdelphi]

Решил тут написать программу, выводящую квадратные корни из первых _N_ чисел на экран (с использованием XMM-регистров).
Схема работает так:
1) Вводится число _N_.
2) XMM1 инициализируется как 4.0 3.0 2.0 1.0 (4 числа типа single).
3) XMM2 инициализируется как 4.0 4.0 4.0 4.0 (то же)
4) В цикле до _N_ shr 2 извлекается корень из 4 чисел сразу, после чего заносится в массив _a_, который и выводится на экран.
Проблема: первые 4 значения выводятся нормально, после чего идут одни нули.
Вопросы:
1) Как организовать сохранение XMM0, XMM1 и XMM2 в цикле? При выполнении writeln регистры XMM почему-то обнуляются.
2) Как "вытащить" значение XMM0 в массив _а_?
Код прилагается.

program Project2;
{$APPTYPE CONSOLE}
{$R *.res}
uses
  System.SysUtils;
type
  TFloat = array [1..4] of single;
var a: TFloat; i, N, j: int32;
procedure init;
asm
  mov      rax, 4                // установить счетчик
@set_xmm1:
  shufps   xmm1, xmm1, 10010011b // очистить место для нового числа
  cvtsi2ss xmm1, rax             // xmm3 = 0.0
  dec      al                    // xmm3 заполнен?
  jnz      @set_xmm1             // нет -> назад, заполнять его
  // конец цикла
  mov      cl, 4                 // установить счетчик
  mov      al, cl                // и скопировать его значение
@set_xmm2:
  cvtsi2ss xmm2, rax             // инициализация xmm2
  shufps   xmm2, xmm2, 10010011b // то же самое,
  dec      cl                    // но для
  jnz      @set_xmm2             // xmm2
end;
procedure mysqrt;
asm
  sqrtps   xmm0, xmm1
  addps    xmm1, xmm2
  movups   N, xmm0        // как сделать?
end;
begin
  readln(N);
  init;
  for i := 1 to N div 4 do
  begin
    mysqrt;
    for j := 1 to 4 do
      if i + j - 1 <= N then
        writeln(a[j]);
  end;
  readln
end.

[f.hump]

не совсем понятно почему четверки (xmm2) задвигаются на N, вместо того, чтобы задвинуть результат (xmm0) на а.

LEA rdi, [a]
MOVUPS [rdi], xmm0

[Пепел Феникса]

эмм, Delphi XE2?

[alexdelphi]

Сообщение от f.hump

не совсем понятно почему четверки (xmm2) задвигаются на N, вместо того, чтобы задвинуть результат (xmm0) на а.

LEA rdi, [a]
MOVUPS [rdi], xmm0

точно...

Сообщение от Пепел Феникса

эмм, Delphi XE2?

Да. А может, есть здесь причины использовать 2010? У меня она тоже стоит.

P.S. По идее, mysqrt должна вообще-то возвращать результат, а не модифицировать некий массив, но из-за ограничений XE2 не могу с этим справиться.

[f.hump]

По идее, mysqrt должна вообще-то возвращать результат, а не модифицировать некий массив

что за идея?

[Пепел Феникса]

P.S. По идее, mysqrt должна вообще-то возвращать результат, а не модифицировать некий массив, но из-за ограничений XE2 не могу с этим справиться.

ну для этого она как минимум должна быть функцией.
далее, я бы посмотрел как Делфи вовращает результат нужного вам типа(ну тут int32, значит поидее он в eax вернется) ну и возвращать.

[alexdelphi]

Т.е. надо написать функцию, возвращающую нужный мне тип, и под отладчиком посмотреть, в какой регистр это значение попадает?

[Пепел Феникса]

ну да, так вы сможете правильно писать функции, и чтоб использовать их по человечески.
кстати, у вас там rax, у вас точно проект 64-битный?