1! { dg-do run } 2! { dg-require-effective-target fortran_large_real } 3! Program to test the UNPACK intrinsic for large real type 4program 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 call abort 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) call abort 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 call abort 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 call abort 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) call abort 43 44end program 45