DVMH-Indirect (1158460), страница 4

Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 4)

! одномерное кольцо

Neighbors(i) = 2

Neighbor(I, 1) = modulo(I, nnode) + 1

Neighbor(I, 2) = modulo(I – 2, nnode) + 1

Enddo

! Новая конструкция:

!DVM$ SHADOW fi(i) with fi(neighbor(I, 1:neighbors(i)))

!DVM$ SHADOW psi(i) with psi(neighbor(I, 1:neighbors(i)))

!DVM$ REALIGN fi(i) with neighbors(i)

!DVM$ REALIGN psi(i) with neighbors(i)

! Functions

!DVM$ PARALLEL(i) on fi(i), shadow_compute

Do I = 1, nnode

Fi(i) = 1

Psi(i) = 2

enddo

! Calculate

Do it = 1, itcount

! Jacobi iteration

!DVM$ PARALLEL(i) on psi(i), shadow_renew(fi)

Do I = 1, nnode

Psi(i) = 0

Do j = 1, neighbors(i)

Psi(i) = psi(i) + fi(neighbor(I,j)) / neighbors(i)

enddo

Enddo

!DVM$ PARALLEL(i) on fi(i), shadow_renew(psi)

Do I = 1, nnode

fi(i) = 0

Do j = 1, neighbors(i)

fi(i) = fi(i) + psi(neighbor(I,j)) / neighbors(i)

enddo

Enddo

enddo

end

Расширение языков C-DVMH и Fortran-DVMH. Текущее состояние

2. indirect-rule :: indirect ( var-name )

3. derived-rule ::= derived ( derived-elem-list with derived-templ )

4. derived-elem ::= int-range-expr

5. int-range-expr ::= прозвольное целочисленное выражение + в индексных выражениях допустимы диапазоны, использование align-dummy переменных.

6. derived-templ ::= var-name [ derived-templ-axis-spec ]...

7. derived-templ-axis-spec ::= [ ] | [ @ align-dummy [ + shadow-name ]... ] | [ int-expr ]

8. Рис. 1. Формула БНФ для новых правил распределения

10. Примеры:

11. #define N 100

12. #pragma dvm array distribute[block]

13. int map[N];

14. #pragma dvm array

15. int *A;

16. A = (int *)malloc(sizeof(int) * N);

17. #pragma dvm redistribute(A[indirect(map)])

19. #pragma dvm array

20. int *C;

21. C = (int *)malloc(sizeof(int) * 3 * N);

22. #pragma dvm redistribute(C[derived(3*i, 3*i+1, 3*i+2 with A[@i])])

23. #pragma dvm redistribute(C[derived(3*i, 3*i+1, 3*i+2 with A[i])])

24. indirect-format	is INDIRECT (array-name)
    derived-format	is DERIVED (derived-elem-list) WITH derived-spec
    derived-elem	is int-range-expr
    int-range-expr	is align-dummy-use
    derived-spec	is derived-target [(derived-subscript)]
    derived-target	is array-name
    derived-subscript	is int-expr or @align-dummy[ + shadow-name ]...

26. Пример:

27. INTEGER,PARAMETER :: N=100

28. INTEGER map(N)

29. !DVM$ DISTRIBUTE map(BLOCK)

30. INTEGER,ALLOCATABLE:: A(:),C(:)

31. ALLOCATE(A(N))

32. !DVM$ REDISTRIBUTE A(INDIRECT(map))

33. ALLOCATE(C(3*N))

34. !DVM$ REDISTRIBUTE C (DERIVED(3*I,3*i+1,3*i+2) WITH A(@i))

35. !DVM$ REDISTRIBUTE C (DERIVED(3*I,3*i+1,3*i+2) WITH A(i))

37. shadow-add ::= shadow_add ( templ-name [ shadow-axis ]... = shadow-name ) [ include_to ( var-name-list ) ]

38. shadow-axis ::= [ ] | [ derived-elem-list with derived-templ ]

39. Рис. 2. Формула БНФ для задания поэлементных теневых граней

40. Пример:

41. #pragma dvm template[SZ] distribute[block]

42. void *B;

43. #pragma dvm template[SZ] distribute[derived(i-1,i,i+1 with B[@i])]

44. void *A;

45. #pragma dvm array align([i] with A[i])

46. int data[SZ];

47. int left = 2;

48. int right = 2;

49. #pragma dvm shadow_add(A[i-left with A[@i]] = "left" " 1") include_to(data)

50. #pragma dvm shadow_add(A[i+right with A[@i]] = "right") include_to(data)

51. #pragma dvm parallel([i] on A[i]) shadow_renew(data[NAM, STR(left 1)])

52. for (int i = 0; i < 0; i++) {}

53. shadow-add-directive	is SHADOW_ADD (template-name(shadow-axis-list)=shadow-name) [INCLUDE_TO (dist-array-list)]
    shadow-axis	is (derived-elem-list) WITH derived-spec

56. localize-spec ::= localize ( ref-var-name => target-var-name [ axis-specifier ]... )

57. axis-specifier ::= [ ] | [ : ]

58. Рис. 3. Формула БНФ для директивы локализации значений индексного массива

59. Пример:

60. #pragma dvm array distribute[]

61. int ref[N];

62. #pragma dvm localize(ref => A[:])

63. #pragma dvm parallel([i] on ref[i])

64. for (int i = 0; i < N; i++) {

65. A[ref[i]] = 4;

66. }

67. localize-directive	is LOCALIZE (ref-array-name =>target-array-name(target-subscript-list))
    target-subscript	is * or :

68. Пример:

69. INTEGER ref(N)

70. !DVM$ DISTRIBUTE ref(*)

71. !DVM$ LOCALIZE (ref=>A(:))

72. !DVM$ PARALLEL (i) ON ref(i)

73. DO i=1,N

74. A(ref(i)) = 4

75. ENDDO

Характеристики

Список файлов учебной работы