gcc/testsuite/gfortran.dg/spread_shape_1.f90

   1 ! { dg-do compile }
   2 ! Tests the fix for PR29060 in which the shape of the result
   3 ! of SPREAD was not available to the scalarizer.
   4 !
   5 ! Contributed by Paul Thomas  <pault@gcc.gnu.org>
   6   real,dimension(:, :),pointer :: ptr
   7   real,dimension(2, 2) :: u
   8
   9   u = reshape((/0.25, 0.5, 0.75, 1.00/),(/2,2/))
  10
  11   allocate (ptr(2,2))
  12
  13 ! Original PR
  14   ptr(:, :) = u + spread ((/1.0, 2.0/), 2, size(u, 2))
  15   if (any (ptr .ne. &
  16         reshape ((/1.25, 2.50, 1.75, 3.00/), (/2, 2/)))) call abort ()
  17
  18 ! Check that the fix works correctly with the source shape after ncopies
  19   ptr(:, :) = u + spread ((/2.0, 3.0/), 1, size (u, 1))
  20   if (any (ptr .ne. &
  21         reshape ((/2.25, 2.50, 3.75, 4.00/), (/2,2/)))) call abort ()
  22 end