Mercurial > hg > CbC > CbC_gcc

comparison gcc/testsuite/gfortran.dg/intrinsic_unpack_2.f90 @ 111:04ced10e8804

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
gcc 7
author kono
date Fri, 27 Oct 2017 22:46:09 +0900
parents
children 84e7813d76e9
comparison
equal deleted inserted replaced
68:561a7518be6b 111:04ced10e8804
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 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
44 end program