libgomp/testsuite/libgomp.fortran/allocatable10.f90

   1 ! { dg-do run }
   2
   3   integer, allocatable :: a, b(:), c(:,:)
   4   integer :: i
   5 !$omp declare reduction (foo : integer : omp_out = omp_out + omp_in) &
   6 !$omp & initializer (omp_priv = 0)
   7   if (allocated (a) .or. allocated (b) .or. allocated (c)) call abort
   8   allocate (a, b(6:9), c(3, 8:9))
   9   a = 0
  10   b = 0
  11   c = 0
  12   if (.not.allocated (a)) call abort
  13   if (.not.allocated (b) .or. size (b) /= 4) call abort
  14   if (lbound (b, 1) /= 6 .or. ubound (b, 1) /= 9) call abort
  15   if (.not.allocated (c) .or. size (c) /= 6) call abort
  16   if (size (c, 1) /= 3 .or. size (c, 2) /= 2) call abort
  17   if (lbound (c, 1) /= 1 .or. ubound (c, 1) /= 3) call abort
  18   if (lbound (c, 2) /= 8 .or. ubound (c, 2) /= 9) call abort
  19 !$omp parallel do reduction (+:a, b, c)
  20   do i = 1, 10
  21     if (.not.allocated (a)) call abort
  22     if (.not.allocated (b) .or. size (b) /= 4) call abort
  23     if (lbound (b, 1) /= 6 .or. ubound (b, 1) /= 9) call abort
  24     if (.not.allocated (c) .or. size (c) /= 6) call abort
  25     if (size (c, 1) /= 3 .or. size (c, 2) /= 2) call abort
  26     if (lbound (c, 1) /= 1 .or. ubound (c, 1) /= 3) call abort
  27     if (lbound (c, 2) /= 8 .or. ubound (c, 2) /= 9) call abort
  28     a = a + i
  29     b = b + 2 * i
  30     c = c + 3 * i
  31   end do
  32   if (.not.allocated (a)) call abort
  33   if (.not.allocated (b) .or. size (b) /= 4) call abort
  34   if (lbound (b, 1) /= 6 .or. ubound (b, 1) /= 9) call abort
  35   if (.not.allocated (c) .or. size (c) /= 6) call abort
  36   if (size (c, 1) /= 3 .or. size (c, 2) /= 2) call abort
  37   if (lbound (c, 1) /= 1 .or. ubound (c, 1) /= 3) call abort
  38   if (lbound (c, 2) /= 8 .or. ubound (c, 2) /= 9) call abort
  39   if (a /= 55 .or. any (b /= 110) .or. any (c /= 165)) call abort
  40   a = 0
  41   b = 0
  42   c = 0
  43 !$omp parallel do reduction (foo : a, b, c)
  44   do i = 1, 10
  45     if (.not.allocated (a)) call abort
  46     if (.not.allocated (b) .or. size (b) /= 4) call abort
  47     if (lbound (b, 1) /= 6 .or. ubound (b, 1) /= 9) call abort
  48     if (.not.allocated (c) .or. size (c) /= 6) call abort
  49     if (size (c, 1) /= 3 .or. size (c, 2) /= 2) call abort
  50     if (lbound (c, 1) /= 1 .or. ubound (c, 1) /= 3) call abort
  51     if (lbound (c, 2) /= 8 .or. ubound (c, 2) /= 9) call abort
  52     a = a + i
  53     b = b + 2 * i
  54     c = c + 3 * i
  55   end do
  56   if (.not.allocated (a)) call abort
  57   if (.not.allocated (b) .or. size (b) /= 4) call abort
  58   if (lbound (b, 1) /= 6 .or. ubound (b, 1) /= 9) call abort
  59   if (.not.allocated (c) .or. size (c) /= 6) call abort
  60   if (size (c, 1) /= 3 .or. size (c, 2) /= 2) call abort
  61   if (lbound (c, 1) /= 1 .or. ubound (c, 1) /= 3) call abort
  62   if (lbound (c, 2) /= 8 .or. ubound (c, 2) /= 9) call abort
  63   if (a /= 55 .or. any (b /= 110) .or. any (c /= 165)) call abort
  64   a = 0
  65   b = 0
  66   c = 0
  67 !$omp simd reduction (+:a, b, c)
  68   do i = 1, 10
  69     if (.not.allocated (a)) call abort
  70     if (.not.allocated (b) .or. size (b) /= 4) call abort
  71     if (lbound (b, 1) /= 6 .or. ubound (b, 1) /= 9) call abort
  72     if (.not.allocated (c) .or. size (c) /= 6) call abort
  73     if (size (c, 1) /= 3 .or. size (c, 2) /= 2) call abort
  74     if (lbound (c, 1) /= 1 .or. ubound (c, 1) /= 3) call abort
  75     if (lbound (c, 2) /= 8 .or. ubound (c, 2) /= 9) call abort
  76     a = a + i
  77     b = b + 2 * i
  78     c = c + 3 * i
  79   end do
  80   if (.not.allocated (a)) call abort
  81   if (.not.allocated (b) .or. size (b) /= 4) call abort
  82   if (lbound (b, 1) /= 6 .or. ubound (b, 1) /= 9) call abort
  83   if (.not.allocated (c) .or. size (c) /= 6) call abort
  84   if (size (c, 1) /= 3 .or. size (c, 2) /= 2) call abort
  85   if (lbound (c, 1) /= 1 .or. ubound (c, 1) /= 3) call abort
  86   if (lbound (c, 2) /= 8 .or. ubound (c, 2) /= 9) call abort
  87   if (a /= 55 .or. any (b /= 110) .or. any (c /= 165)) call abort
  88   a = 0
  89   b = 0
  90   c = 0
  91 !$omp simd reduction (foo : a, b, c)
  92   do i = 1, 10
  93     if (.not.allocated (a)) call abort
  94     if (.not.allocated (b) .or. size (b) /= 4) call abort
  95     if (lbound (b, 1) /= 6 .or. ubound (b, 1) /= 9) call abort
  96     if (.not.allocated (c) .or. size (c) /= 6) call abort
  97     if (size (c, 1) /= 3 .or. size (c, 2) /= 2) call abort
  98     if (lbound (c, 1) /= 1 .or. ubound (c, 1) /= 3) call abort
  99     if (lbound (c, 2) /= 8 .or. ubound (c, 2) /= 9) call abort
 100     a = a + i
 101     b = b + 2 * i
 102     c = c + 3 * i
 103   end do
 104   if (.not.allocated (a)) call abort
 105   if (.not.allocated (b) .or. size (b) /= 4) call abort
 106   if (lbound (b, 1) /= 6 .or. ubound (b, 1) /= 9) call abort
 107   if (.not.allocated (c) .or. size (c) /= 6) call abort
 108   if (size (c, 1) /= 3 .or. size (c, 2) /= 2) call abort
 109   if (lbound (c, 1) /= 1 .or. ubound (c, 1) /= 3) call abort
 110   if (lbound (c, 2) /= 8 .or. ubound (c, 2) /= 9) call abort
 111   if (a /= 55 .or. any (b /= 110) .or. any (c /= 165)) call abort
 112 end