Обращение к предыдущему/следующему элементу массива по его адресу.

[InternetStranger]

Прошу помочь разобраться.

Предыстория вопроса:
Значит, есть у меня массив двумерный (of real), площадью > 1e6.
В рекурентной формуле интенсивно пользуюсь его элементами.
Поскольку читаю его последовательно, то можно оптимизировать обращение к массиву. (Определить указатель, и читать как из файла, смещая указатель)
Простейший пример реализации этого алгоритма:

	var
   u : array[1..10] of real;

Type
    TRealPointer = ^real;

procedure TForm1.Button1Click(Sender: TObject);
var
   i          : integer;
   s          : string;
   my_pointer : TRealPointer;   // указатель на участок памяти длиной 8 байт
begin
     // Тестовые данные в массив
     for i:=1 to 10 do u[i]:= 0.1*i;

     // типы двойной точности "real" имеют в памяти размер 8 байт.
     // массив - это последовательное расположение в памяти этих блоков - элементов- в данном случае по 8 байт
     my_pointer := @u; // или так @u[1]; // запишем адрес начала массива
                       // доступ к значению в памяти, на которое сейчас указывает указатель - через "разыменование" указателя my_pointer^ 

     // выведет первый элемент массива 
     label2.Caption := FloatToStr( my_pointer^ );

     // а теперь распечатаем весь массив:
     // чтобы получить адрес след. элемента массива u[i+1], скажем из адреса элемента u[i], нужно задать смещение на длину блока,
     // в данном случае на 8 байт:  " @u[i] + $08 " (след. элемент)   ; " @u[i] - $08 " (пред. элемент)
     s:= '';
     for i:=1 to 10 do begin
         s:= s + FloatToStr( my_pointer^ ) + '  |  ';
         inc(my_pointer,$01); // -> сдвигаем указатель вправо к след. элементу  (на 8 байт) . вот эта строка преобразуется всего в одну ассемблерную строчку, потому нужно ЯВНО указывать смещение: либо число, либо константу (НЕ типизированную)
     end;
     // выведет весь массив
     label4.Caption := s ;
end;

Видно, что указатель my_pointer в пределах i-й итерации указывает на область памяти, содержащую значение u[i]

Собственно вопрос:
Внутри цикла, на i-й итерации, необходимо обратиться к трем элементам u[i-1], u[i], u[i+1].
Не могу сообразить как указать явно указать смещение вперед/назад относительно my_pointer.
Охота что-то вроде my_pointer-$08, my_pointer, my_pointer+$08 чтобы эти конструкции можно было подставить в одно выражение
(Требуется извлечь максимальную эффективность! )
Можно, конечно, определить дополнительную переменную и получится что-то вроде

     ...
     for i:=2 to 9 do begin
	 prev_ptr:= my_pointer; 
	 next_ptr:= my_pointer;
	 inc(next_ptr); dec(prev_ptr);
         s:= s + FloatToStr( my_pointer^ ) + '  |  Сумма окрестностных элементов '+ FloatToStr( prev_ptr^ + next_ptr^ );
         inc(my_pointer,$01); // -> сдвигаем указатель вправо к след. элементу  (на 8 байт) . вот эта строка преобразуется всего в одну ассемблерную строчку, потому нужно ЯВНО указывать смещение: либо число, либо константу (НЕ типизированную)
     end;
     ...
end;

Но тут целых две операции присваивания, да еще инкрейз/декрейз. Что-то не нравится мне. проще (и может быть даже быстрее) обращаться к элементу массива
Работа с указателями, конечно, у Дельфина довольно специфична

[maxionans]

Возможно, так подойдет:

  S := S + FloatToStr( TRealPointer(Cardinal(my_pointer) - 8)^ ) +
              FloatToStr( my_pointer^ ) +
              FloatToStr( TRealPointer(Cardinal(my_pointer) + 8)^ );

[InternetStranger]

Спасибо )) Отличный способ. Все как положено - жесткое приведение типов. Не мог сообразить, к какому же типу нужно привести указатели ))

Изучаю соответствующий ассемблерный код

Кста, а почему Cardinal-то ? Это из каких-то соображений, или просто так?
Он вроде переменной длины?

[maxionans]

Кста, а почему Cardinal-то ? Это из каких-то соображений, или просто так?

Да, из соображений того, что Cardinal - это четырехбайтовый целочисленный тип без знака, каким и должен быть указатель в IA32.

Он вроде переменной длины?

Что значит "переменной длины"?

[InternetStranger]

Да это я глупость какую-то прочитал на сайте одном. не так понял.

Тип Cardinal это тип целых чисел, чей размер не гарантирован, базовый беззнаковый целый тип в Delphi, и на настоящий момент он имеет тот же размер, что и LongWord: 32 бита.

И щас не понимаю, что это значит.

[maxionans]

Это они видимо говорят о том, что, когда выйдет delphi, способная генерировать код для 64-битной архитектуры, в ней Cardinal может стать 8-байтовым целочисленным типов без знака.

[InternetStranger]

Протестировал на скорость выполнения.
К сожалению, компилятор не правильно понимает, чего от него хотят (может я просто не так спрашивал)).

В любом случае, следующий способ хоть и не очень удобен, но явно эффективнее:

 
...
my_pointer:= @u[i];
prev_ptr:= my_pointer; 
next_ptr:= my_pointer;
inc(next_ptr); dec(prev_ptr);

for i:=2 to 9 do begin
         s:= s + FloatToStr( my_pointer^ ) + '  |  Сумма окрестностных элементов '+ FloatToStr( prev_ptr^ + next_ptr^ );
         inc(my_pointer); inc(next_ptr); inc(prev_ptr);
end;
...

Т.е. хоть и приходится поддерживать 3 указателя (не шибко эстетично), но за то явно указано, по какому адресу обратиться )

PS: За предложенный алгоритм тоже спасибо.
maxionans.reputacia.inc(+1);

[maxionans]

А можно поинтересоваться отчего такая щепетильность к оптимизации выполнения? У вас там гигабайты обрабатываемых данных или же вы пишете для калькулятора?

[InternetStranger]

из "калькуляторов" уже давно вырос ))
вот собственно моя задача Статья

Под скромненькой фразой во введении "С помощью численных методов рассмотрена..." кроется многолетняя борьба за каждый такт, за каждый квант процессорного времени )) И все во имя отечественной науки.

вот это тоже ради нее

[maxionans]

нет, я про калькулятор говорил не в том смысле, что вы создаете программу-калькулятор, а в том, что компьютер, для которого вы пишите, по производительности сравним с калькулятором.