rtsIDr (1158449), страница 25
Текст из файла (страница 25)
+ crtpl, dvmadr, mappl, endpl, dopl,
+ crtshg, inssh, insshd, strtsh, recvsh, sendsh, waitsh
real bptr(1)
integer dvm
integer amref, psref, mvref, plref, shgref
integer amdim(2), disaxs(2), dispar(2)
integer lshwd(2), hshwd(2), shsign(2), axis(2), coeff(2),
+ const(2)
integer lvaddr(2), lvtype(2), iiniti(2), ilasti(2), istep(2)
integer oiniti(2), olasti(2), ostep(2)
integer FDL1, FDL2, ADL1, ADL2, finssh
C Размер каждого измерения массива
C с регулярной зависимостью по данным
parameter (k = 200)
C Заголовок массива с регулярной зависимостью по данным
integer ahdr(6)
ahdr(5) = 1
ahdr(6) = 1
C Длина прямой зависимости по данным для измерения 1
FDL1 = 4
C Длина прямой зависимости по данным для измерения 2
FDL2 = 3
C Длина обратной зависимости по данным для измерения 1
ADL1 = 2
C Длина обратной зависимости по данным для измерения 2
ADL2 = 1
C Признак использования функции inssh_ для регистрации
C массива в группе границ; при нулевом значении finssh
C используется функция insshd_
finssh = 0
C Инициализация системы поддержки
dvm = linit (1)
C Создание представления абстрактной машины
C и отображение его в процессорную подсистему
amref = getam ()
psref = getps (amref)
amdim(1) = k
amdim(2) = k
mvref = crtamv (amref, 2, amdim,0)
disaxs(1) = 1
disaxs(2) = 2
dispar(1) = 0
dispar(2) = 0
dvm = distr (mvref, psref, 2, disaxs, dispar)
C Создание и распределение массива
C с регулярной зависимостью по данным
C Ширина нижней границы по измерению 1 - длина
C прямой зависимости по данным для измерения 1
lshwd(1) = FDL1
C Ширина нижней границы по измерению 2 - длина
C прямой зависимости по данным для измерения 2
lshwd(2) = FDL2
C Ширина верхней границы по измерению 1 - длина
C обратной зависимости по данным для измерения 1
hshwd(1) = ADL1
C Ширина верхней границы по измерению 2 - длина
C обратной зависимости по данным для измерения 2
hshwd(2) = ADL2
dvm = crtda (ahdr, 1, bptr, 2, 4, amdim, 0, 0, lshwd, hshwd)
axis(1) = 1
axis(2) = 2
coeff(1) = 1
coeff(2) = 1
const(1) = 0
const(2) = 0
dvm = align (ahdr, mvref, axis, coeff, const)
C Параллельный цикл инициализации массива
C с регулярной зависимостью по данным
plref = crtpl (2)
lvaddr(1) = dvmadr (j)
lvaddr(2) = dvmadr (i)
lvtype(1) = 1
lvtype(2) = 1
iiniti(1) = 0
iiniti(2) = 0
ilasti(1) = k - 1
ilasti(2) = k - 1
istep(1) = 1
istep(2) = 1
dvm = mappl (plref, ahdr, axis, coeff, const, lvaddr, lvtype
+ iiniti, ilasti, istep, oiniti, olasti, ostep)
2 if (dopl (plref) .eq. 0) goto 3
do 1 j = oiniti(1), olasti(1), ostep(1)
do 1 i = oiniti(2), olasti(2), ostep(2)
bptr( ahdr(3) + 1 + i + ahdr(2) * j ) = 1.
1 continue
goto 2
3 dvm = endpl (plref)
C Создание группы границ для обеспечения
C обратной зависимости по данным
shgref = crtshg (0)
C Нулевая ширина нижней границы по измерению 1
lshwd(1) = 0
C Нулевая ширина нижней границы по измерению 2
lshwd(2) = 0
C Ширина верхней границы по измерению 1 - длина
C обратной зависимости по данным для измерения 1
hshwd(1) = ADL1
C Ширина верхней границы по измерению 2 - длина
C обратной зависимости по данным для измерения 2
hshwd(2) = ADL2
dvm = inssh (shgref, ahdr, lshwd, hshwd, 0)
C Инициализация теневых граней для обеспечения
C обратной зависимости по данным
dvm = strtsh (shgref)
dvm = waitsh (shgref)
C Создание группы границ для обеспечения
C прямой зависимости по данным
shgref = crtshg (0)
C Ширина нижней границы по измерению 1 - длина
C прямой зависимости по данным для измерения 1
lshwd(1) = FDL1
C Ширина нижней границы по измерению 2 - длина
C прямой зависимости по данным для измерения 2
lshwd(2) = FDL2
C Нулевая ширина верхней границы по измерению 1
hshwd(1) = 0
C Нулевая ширина верхней границы по измерению 2
hshwd(2) = 0
if (finssh .eq. 0) goto 4
C Для регистрации массива в группе границ используем inssh_
dvm = inssh (shgref, ahdr, lshwd, hshwd, 0)
goto 5
C Для регистрации массива в группе границ используем insshd_
C По измерениям 1 и 2 используется только нижняя граница
4 shsign(1) = 3
shsign(2) = 3
dvm = insshd (shgref, ahdr, lshwd, hshwd, 1, shsign)
goto 5
C Прием импортируемых элементов для
С обеспечения прямой зависимости по данным
5 dvm = recvsh (shgref)
dvm = waitsh (shgref)
C ПАРАЛЛЕЛЬНЫЙ ЦИКЛ С ПРЯМОЙ И ОБРАТНОЙ
C ЗАВИСИМОСТЯМИ ПО ДАННЫМ МЕЖДУ ВИТКАМИ
plref = crtpl (2)
iiniti(1) = 1
iiniti(2) = 1
ilasti(1) = k - 2
ilasti(2) = k - 2
dvm = mappl (plref, ahdr, axis, coeff, const, lvaddr, lvtype,
+ iiniti, ilasti, istep, oiniti, olasti, ostep)
6 if (dopl (plref) .eq. 0) goto 7
do 8 j = oiniti(1), olasti(1), ostep(1)
do 8 i = oiniti(2), olasti(2), ostep(2)
bptr(ahdr(3) + 1 + i + ahdr(2) * j) =
+ ( bptr(ahdr(3) + 1 + (i - FDL2) + ahdr(2) * j) +
+ bptr(ahdr(3) + 1 + i + ahdr(2) * (j - FDL1)) +
+ bptr(ahdr(3) + 1 + (i + ADL2) + ahdr(2) * j) +
+ bptr(ahdr(3) + 1 + i + ahdr(2) * (j + ADL1)) ) / 4.
8 continue
goto 6
7 dvm = endpl (plref)
C Передача экспортируемых элементов для
C обеспечения прямой зависимости по данным
dvm = sendsh (shgref)
dvm = waitsh (shgref)
C Завершение работы с системой поддержки
dvm = lexit (0)
end
ПРОГРАММА НА ЯЗЫКЕ C
#include "dvmlib.h"
int main(int argc, char *argv[])
{
long InitPar = 0, UserRes = 0;
long StaticMV = 0, StaticDA = 0, ReDistrDA = 0,
StaticShG = 0, FullShd = 0, MaxShdCount = 1;
long ShdSignArray[2];
long i, j, Rank2 = 2, TypeSize = sizeof(float);
AMRef amref;
PSRef psref;
AMViewRef mvref;
LoopRef plref;
ShadowGroupRef fshgref, ashgref;
AddrType lvaddr[2];
long lvtype[2] = {0,0};
long amdim[2], disaxs[2], dispar[2];
long lshwdth[2], hshwdth[2], axis[2], coeff[2], cnst[2];
long iiniti[2], ilasti[2], istep[2];
long oiniti[2], olasti[2], ostep[2];
long ExtHdrSign = 0;
DVMFILE *OutFile;
long k = 200; /* размер каждого измерения массива
с регулярной зависимостью по данным */
long FDL1 = 4; /* длина прямой зависимости
по данным для измерения 1 */
long FDL2 = 3; /* длина прямой зависимости
по данным для измерения 2 */
long ADL1 = 2; /* длина обратной зависимости
по данным для измерения 1 */
long ADL2 = 1; /* длина обратной зависимости
по данным для измерения 2 */
int f_inssh = 0; /* признак использования функции inssh_ для
регистрации массива в группе границ; при
нулевом значении f_inssh используется
функция insshd_ */
long Ahdr[3]; /* заголовок массива с регулярной
зависимостью по данным */
rtl_init(InitPar, argc, argv); /* инициализация системы поддержки */
/* Создание представления абстрактной машины
и отображение его в процессорную подсистему */
amref = getam_ ();
psref = getps_ (&amref);
amdim[0] = k;
amdim[1] = k;
mvref = crtamv_ (&amref, &Rank2, amdim, &StaticMV);
disaxs[0] = 1;
disaxs[1] = 2;
dispar[0] = 0;
dispar[1] = 0;
distr_ (&mvref, &psref, &Rank2, disaxs, dispar);
/* Создание и распределение массива
с регулярной зависимостью по данным */
lshwdth[0] = FDL1; /* ширина нижней границы по измерению 1 - длина
прямой зависимости по данным для измерения 1 */
lshwdth[1] = FDL2; /* ширина нижней границы по измерению 2 - длина
прямой зависимости по данным для измерения 2 */
hshwdth[0] = ADL1; /* ширина верхней границы по измерению 1 - длина
обратной зависимости по данным для
измерения 1 */
hshwdth[1] = ADL2; /* ширина верхней границы по измерению 2 - длина
обратной зависимости по данным для
измерения 2 */
crtda_ (Ahdr, &ExtHdrSign, NULL, &Rank2, &TypeSize, amdim,
&StaticDA, &ReDistrDA, lshwdth, hshwdth);
axis[0] = 1;
axis[1] = 2;
coeff[0] = 1;
coeff[1] = 1;
cnst[0] = 0;
cnst[1] = 0;
align_ (Ahdr,&mvref,axis,coeff,cnst);
/* Параллельный цикл инициализации массива
с регулярной зависимостью по данным */
plref = crtpl_ (&Rank2);
lvaddr[0] = (AddrType)&j;
lvaddr[1] = (AddrType)&i;
iiniti[0] = 0;
iiniti[1] = 0;
ilasti[0] = k - 1;
ilasti[1] = k - 1;
istep[0] = 1;
istep[1] = 1;
mappl_ (&plref, (PatternRef *)Ahdr, axis, coeff, cnst, lvaddr,
lvtype, iiniti, ilasti, istep, oiniti, olasti, ostep);
while (dopl_ (&plref))
for ( j = oiniti[0]; j <= olasti[0]; j += ostep[0] )
for ( i = oiniti[1]; i <= olasti[1]; i += ostep[1] )
DAElm2(Ahdr, float, j, i) = 1.f;
endpl_ (&plref);
/* Создание группы границ для обеспечения
обратной зависимости по данным */
ashgref = crtshg_ (&StaticShG);
lshwdth[0] = 0; /* нулевая ширина нижней границы по измерению 1 */
lshwdth[1] = 0; /* нулевая ширина нижней границы по измерению 2 */
hshwdth[0] = ADL1; /* ширина верхней границы по измерению 1 - длина
обратной зависимости по данным для
измерения 1 */
hshwdth[1] = ADL2; /* ширина верхней границы по измерению 2 - длина
обратной зависимости по данным для
измерения 2 */
inssh_ (&ashgref, Ahdr, lshwdth, hshwdth, &FullShd);
/* Инициализация теневых граней для обеспечения
обратной зависимости по данным */
strtsh_ (&ashgref);
waitsh_ (&ashgref);
/* Создание группы границ для обеспечения
прямой зависимости по данным */
fshgref = crtshg_ (&StaticShG);
lshwdth[0] = FDL1; /* ширина нижней границы по измерению 1 - длина
прямой зависимости по данным для измерения 1 */
lshwdth[1] = FDL2; /* ширина нижней границы по измерению 2 - длина
прямой зависимости по данным для измерения 2 */
hshwdth[0] = 0; /* нулевая ширина верхней границы по измерению 1 */
hshwdth[1] = 0; /* нулевая ширина верхней границы по измерению 2 */
if(f_inssh) /* с какой функцией работаем, с inssh_ или с insshd_ */
{
/* Для регистрации массива в группе границ используем inssh_ */
inssh_ (&fshgref, Ahdr, lshwdth, hshwdth, &FullShd);
}
else
{
/* Для регистрации массива в группе границ используем insshd_ */
ShdSignArray[0] = 3; /* по измерению 1 используется
только нижняя граница */
ShdSignArray[1] = 3; /* по измерению 2 используется
только нижняя граница */
insshd_ (&fshgref, Ahdr, lshwdth, hshwdth, &MaxShdCount,
ShdSignArray);
}
/* Прием импортируемых элементов для обеспечения
прямой зависимости по данным */
recvsh_ (&fshgref);
waitsh_ (&fshgref);
/* ПАРАЛЛЕЛЬНЫЙ ЦИКЛ С ПРЯМОЙ И ОБРАТНОЙ
ЗАВИСИМОСТЯМИ ПО ДАННЫМ МЕЖДУ ВИТКАМИ */
plref = crtpl_ (&Rank2);
iiniti[0] = 1;
iiniti[1] = 1;
ilasti[0] = k - 2;
ilasti[1] = k - 2;
mappl_ (&plref, (PatternRef *)Ahdr, axis, coeff, cnst, lvaddr,
lvtype, iiniti, ilasti, istep, oiniti, olasti, ostep);
while (dopl_ (&plref))
for ( j = oiniti[0]; j <= olasti[0]; j += ostep[0] )
for ( i = oiniti[1]; i <= olasti[1]; i += ostep[1] )
DAElm2(Ahdr, float, j, i) =
( DAElm2(Ahdr, float, j, i-FDL2) +
DAElm2(Ahdr, float, j-FDL1, i) +
DAElm2(Ahdr, float, j, i+ADL2) +
DAElm2(Ahdr, float, j+ADL1, i) ) / 4.f;
endpl_ (&plref);
/* Передача экспортируемых элементов для обеспечения
прямой зависимости по данным */
sendsh_ (&fshgref);
waitsh_ (&fshgref);
/* Запись результирующего массива в файл
и завершение работы с системой поддержки */
OutFile = dvm_fopen("across.dat", "w+b");
if(OutFile)
{ dvm_dfwrite(Ahdr, 0, OutFile);
dvm_fclose(OutFile);
}
lexit_ (&UserRes);
return (int)UserRes;
}
20Литература
-
Н.А. Коновалов, В.А. Крюков, А.А. Погребцов, Ю.Л. Сазанов. C-DVM - язык разработки мобильных параллельных программ. Препринт ИПМ им. М.В. Келдыша РАН N 86, Москва, 1997 г.
-
N.A. Konovalov, V.A. Krukov, S.N. Mihailov, A.A. Pogrebtsov. Fortran-DVM language for portable parallel programs development. Proceeding of Software For Multiprocessors & Supercomputers: Theory, Practice, Experience (SMS-TPE 94), Inst. for System Programming RAS, Moscow, Sept. 1994.
-
В.А. Крюков, Р.В. Удовиченко. Средства отладки и мониторинга DVM-программ. Динамический отладчик. Отчет ИПМ РАН. Москва. 1997 г.
-
В.Е. Денисов, В.Н. Ильяков, Н.В. Ковалёва, В.А. Крюков. Отладка эффективности DVM-программ. Препринт ИПМ им. М.В. Келдыша РАН N 74, Москва, 1998 г.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
