gcc/testsuite/gfortran.dg/intrinsic_unpack_2.f90

   1 ! { dg-do run }
   2 ! { dg-require-effective-target fortran_large_real }
   3 ! Program to test the UNPACK intrinsic for large real type
   4 program intrinsic_unpack
   5    implicit none
   6    integer,parameter :: k = selected_real_kind (precision (0.0_8) + 1)
   7
   8    real(kind=k), dimension(3,3) :: ark, brk
   9    complex(kind=k), dimension(3,3) :: ack, bck
  10
  11    logical, dimension(3, 3) :: mask
  12    character(len=500) line1, line2
  13    integer i
  14
  15    mask = reshape ((/.false.,.true.,.false.,.true.,.false.,.false.,&
  16                     &.false.,.false.,.true./), (/3, 3/));
  17
  18    ark = reshape ((/1._k, 0._k, 0._k, 0._k, 1._k, 0._k, 0._k, 0._k, 1._k/), &
  19          (/3, 3/));
  20    brk = unpack ((/2._k, 3._k, 4._k/), mask, ark)
  21    if (any (brk .ne. reshape ((/1._k, 2._k, 0._k, 3._k, 1._k, 0._k, &
  22                                0._k, 0._k, 4._k/), (/3, 3/)))) &
  23       STOP 1
  24    write (line1,'(9F9.5)') brk
  25    write (line2,'(9F9.5)') unpack((/2._k, 3._k, 4._k/), mask, ark)
  26    if (line1 .ne. line2) STOP 2
  27    brk = -1._k
  28    brk = unpack ((/2._k, 3._k, 4._k/), mask, 0._k)
  29    if (any (brk .ne. reshape ((/0._k, 2._k, 0._k, 3._k, 0._k, 0._k, &
  30       0._k, 0._k, 4._k/), (/3, 3/)))) &
  31       STOP 3
  32
  33    ack = reshape ((/1._k, 0._k, 0._k, 0._k, 1._k, 0._k, 0._k, 0._k, 1._k/), &
  34         (/3, 3/));
  35    bck = unpack ((/(2._k, 0._k), (3._k, 0._k), (4._k,   0._k)/), mask, ack)
  36    if (any (real(bck) .ne. reshape ((/1._k, 2._k, 0._k, 3._k, 1._k, 0._k, &
  37         0._k, 0._k, 4._k/), (/3, 3/)))) &
  38         STOP 4
  39    write (line1,'(18F9.5)') bck
  40    write (line2,'(18F9.5)') unpack((/(2._k, 0._k), (3._k, 0._k), (4._k,0._k)/), &
  41         mask, ack)
  42    if (line1 .ne. line2) STOP 5
  43
  44 end program