45 lines
1.6 KiB
Fortran
45 lines
1.6 KiB
Fortran
! { dg-do run }
|
|
! { dg-require-effective-target fortran_large_real }
|
|
! Program to test the UNPACK intrinsic for large real type
|
|
program intrinsic_unpack
|
|
implicit none
|
|
integer,parameter :: k = selected_real_kind (precision (0.0_8) + 1)
|
|
|
|
real(kind=k), dimension(3,3) :: ark, brk
|
|
complex(kind=k), dimension(3,3) :: ack, bck
|
|
|
|
logical, dimension(3, 3) :: mask
|
|
character(len=500) line1, line2
|
|
integer i
|
|
|
|
mask = reshape ((/.false.,.true.,.false.,.true.,.false.,.false.,&
|
|
&.false.,.false.,.true./), (/3, 3/));
|
|
|
|
ark = reshape ((/1._k, 0._k, 0._k, 0._k, 1._k, 0._k, 0._k, 0._k, 1._k/), &
|
|
(/3, 3/));
|
|
brk = unpack ((/2._k, 3._k, 4._k/), mask, ark)
|
|
if (any (brk .ne. reshape ((/1._k, 2._k, 0._k, 3._k, 1._k, 0._k, &
|
|
0._k, 0._k, 4._k/), (/3, 3/)))) &
|
|
call abort
|
|
write (line1,'(9F9.5)') brk
|
|
write (line2,'(9F9.5)') unpack((/2._k, 3._k, 4._k/), mask, ark)
|
|
if (line1 .ne. line2) call abort
|
|
brk = -1._k
|
|
brk = unpack ((/2._k, 3._k, 4._k/), mask, 0._k)
|
|
if (any (brk .ne. reshape ((/0._k, 2._k, 0._k, 3._k, 0._k, 0._k, &
|
|
0._k, 0._k, 4._k/), (/3, 3/)))) &
|
|
call abort
|
|
|
|
ack = reshape ((/1._k, 0._k, 0._k, 0._k, 1._k, 0._k, 0._k, 0._k, 1._k/), &
|
|
(/3, 3/));
|
|
bck = unpack ((/(2._k, 0._k), (3._k, 0._k), (4._k, 0._k)/), mask, ack)
|
|
if (any (real(bck) .ne. reshape ((/1._k, 2._k, 0._k, 3._k, 1._k, 0._k, &
|
|
0._k, 0._k, 4._k/), (/3, 3/)))) &
|
|
call abort
|
|
write (line1,'(18F9.5)') bck
|
|
write (line2,'(18F9.5)') unpack((/(2._k, 0._k), (3._k, 0._k), (4._k,0._k)/), &
|
|
mask, ack)
|
|
if (line1 .ne. line2) call abort
|
|
|
|
end program
|