1 |
302 |
jeremybenn |
! { dg-do run }
|
2 |
|
|
! Take the pack intrinsic through its paces, with all types that are
|
3 |
|
|
! normally accessible.
|
4 |
|
|
program main
|
5 |
|
|
implicit none
|
6 |
|
|
integer :: i
|
7 |
|
|
real(kind=4), dimension(3,3) :: r4
|
8 |
|
|
real(kind=4), dimension(9) :: vr4
|
9 |
|
|
real(kind=4), dimension(9) :: rr4
|
10 |
|
|
real(kind=8), dimension(3,3) :: r8
|
11 |
|
|
real(kind=8), dimension(9) :: vr8
|
12 |
|
|
real(kind=8), dimension(9) :: rr8
|
13 |
|
|
complex(kind=4), dimension(3,3) :: c4
|
14 |
|
|
complex(kind=4), dimension(9) :: vc4
|
15 |
|
|
complex(kind=4), dimension(9) :: rc4
|
16 |
|
|
complex(kind=8), dimension(3,3) :: c8
|
17 |
|
|
complex(kind=8), dimension(9) :: vc8
|
18 |
|
|
complex(kind=8), dimension(9) :: rc8
|
19 |
|
|
integer(kind=1), dimension(3,3) :: i1
|
20 |
|
|
integer(kind=1), dimension(9) :: vi1
|
21 |
|
|
integer(kind=1), dimension(9) :: ri1
|
22 |
|
|
integer(kind=2), dimension(3,3) :: i2
|
23 |
|
|
integer(kind=2), dimension(9) :: vi2
|
24 |
|
|
integer(kind=2), dimension(9) :: ri2
|
25 |
|
|
integer(kind=4), dimension(3,3) :: i4
|
26 |
|
|
integer(kind=4), dimension(9) :: vi4
|
27 |
|
|
integer(kind=4), dimension(9) :: ri4
|
28 |
|
|
integer(kind=8), dimension(3,3) :: i8
|
29 |
|
|
integer(kind=8), dimension(9) :: vi8
|
30 |
|
|
integer(kind=8), dimension(9) :: ri8
|
31 |
|
|
|
32 |
|
|
type i1_t
|
33 |
|
|
integer(kind=1) :: v
|
34 |
|
|
end type i1_t
|
35 |
|
|
type(i1_t), dimension(3,3) :: d_i1
|
36 |
|
|
type(i1_t), dimension(9) :: d_vi1
|
37 |
|
|
type(i1_t), dimension(9) :: d_ri1
|
38 |
|
|
|
39 |
|
|
type i4_t
|
40 |
|
|
integer(kind=4) :: v
|
41 |
|
|
end type i4_t
|
42 |
|
|
type(i4_t), dimension(3,3) :: d_i4
|
43 |
|
|
type(i4_t), dimension(9) :: d_vi4
|
44 |
|
|
type(i4_t), dimension(9) :: d_ri4
|
45 |
|
|
|
46 |
|
|
d_vi1%v = (/(i+10,i=1,9)/)
|
47 |
|
|
d_i1%v = reshape((/1_1, -1_1, 2_1, -2_1, 3_1, -3_1, 4_1, &
|
48 |
|
|
& -4_1, 5_1/), shape(i1))
|
49 |
|
|
d_ri1 = pack(d_i1,d_i1%v>0,d_vi1)
|
50 |
|
|
if (any(d_ri1%v /= (/1_1, 2_1, 3_1, 4_1, 5_1, 16_1, 17_1, 18_1, 19_1/))) &
|
51 |
|
|
& call abort
|
52 |
|
|
|
53 |
|
|
d_vi4%v = (/(i+10,i=1,9)/)
|
54 |
|
|
d_i4%v = reshape((/1_4, -1_4, 2_4, -2_4, 3_4, -3_4, 4_4, &
|
55 |
|
|
& -4_4, 5_4/), shape(d_i4))
|
56 |
|
|
d_ri4 = pack(d_i4,d_i4%v>0,d_vi4)
|
57 |
|
|
if (any(d_ri4%v /= (/1_4, 2_4, 3_4, 4_4, 5_4, 16_4, 17_4, 18_4, 19_4/))) &
|
58 |
|
|
& call abort
|
59 |
|
|
|
60 |
|
|
vr4 = (/(i+10,i=1,9)/)
|
61 |
|
|
r4 = reshape((/1.0_4, -3.0_4, 2.1_4, -4.21_4, 1.2_4, 0.98_4, -1.2_4, &
|
62 |
|
|
& -7.1_4, -9.9_4, 0.3_4 /), shape(r4))
|
63 |
|
|
rr4 = pack(r4,r4>0,vr4)
|
64 |
|
|
if (any(rr4 /= (/ 1.0_4, 2.1_4, 1.2_4, 0.98_4, 15._4, 16._4, 17._4, &
|
65 |
|
|
& 18._4, 19._4 /))) call abort
|
66 |
|
|
|
67 |
|
|
vr8 = (/(i+10,i=1,9)/)
|
68 |
|
|
r8 = reshape((/1.0_8, -3.0_8, 2.1_8, -4.21_8, 1.2_8, 0.98_8, -1.2_8, &
|
69 |
|
|
& -7.1_8, -9.9_8, 0.3_8 /), shape(r8))
|
70 |
|
|
rr8 = pack(r8,r8>0,vr8)
|
71 |
|
|
if (any(rr8 /= (/ 1.0_8, 2.1_8, 1.2_8, 0.98_8, 15._8, 16._8, 17._8, &
|
72 |
|
|
& 18._8, 19._8 /))) call abort
|
73 |
|
|
|
74 |
|
|
vc4 = (/(i+10,i=1,9)/)
|
75 |
|
|
c4 = reshape((/1.0_4, -3.0_4, 2.1_4, -4.21_4, 1.2_4, 0.98_4, -1.2_4, &
|
76 |
|
|
& -7.1_4, -9.9_4, 0.3_4 /), shape(c4))
|
77 |
|
|
rc4 = pack(c4,real(c4)>0,vc4)
|
78 |
|
|
if (any(real(rc4) /= (/ 1.0_4, 2.1_4, 1.2_4, 0.98_4, 15._4, 16._4, 17._4, &
|
79 |
|
|
& 18._4, 19._4 /))) call abort
|
80 |
|
|
if (any(aimag(rc4) /= 0)) call abort
|
81 |
|
|
|
82 |
|
|
vc8 = (/(i+10,i=1,9)/)
|
83 |
|
|
c8 = reshape((/1.0_4, -3.0_4, 2.1_4, -4.21_4, 1.2_4, 0.98_4, -1.2_4, &
|
84 |
|
|
& -7.1_4, -9.9_4, 0.3_4 /), shape(c8))
|
85 |
|
|
rc8 = pack(c8,real(c8)>0,vc8)
|
86 |
|
|
if (any(real(rc8) /= (/ 1.0_4, 2.1_4, 1.2_4, 0.98_4, 15._4, 16._4, 17._4, &
|
87 |
|
|
& 18._4, 19._4 /))) call abort
|
88 |
|
|
if (any(aimag(rc8) /= 0)) call abort
|
89 |
|
|
|
90 |
|
|
vi1 = (/(i+10,i=1,9)/)
|
91 |
|
|
i1 = reshape((/1_1, -1_1, 2_1, -2_1, 3_1, -3_1, 4_1, -4_1, 5_1/), shape(i1))
|
92 |
|
|
ri1 = pack(i1,i1>0,vi1)
|
93 |
|
|
if (any(ri1 /= (/1_1, 2_1, 3_1, 4_1, 5_1, 16_1, 17_1, 18_1, 19_1/))) &
|
94 |
|
|
& call abort
|
95 |
|
|
|
96 |
|
|
vi2 = (/(i+10,i=1,9)/)
|
97 |
|
|
i2 = reshape((/1_2, -1_2, 2_2, -2_2, 3_2, -3_2, 4_2, -4_2, 5_2/), shape(i2))
|
98 |
|
|
ri2 = pack(i2,i2>0,vi2)
|
99 |
|
|
if (any(ri2 /= (/1_2, 2_2, 3_2, 4_2, 5_2, 16_2, 17_2, 18_2, 19_2/))) &
|
100 |
|
|
& call abort
|
101 |
|
|
|
102 |
|
|
vi4 = (/(i+10,i=1,9)/)
|
103 |
|
|
i4 = reshape((/1_4, -1_4, 2_4, -2_4, 3_4, -3_4, 4_4, -4_4, 5_4/), shape(i4))
|
104 |
|
|
ri4 = pack(i4,i4>0,vi4)
|
105 |
|
|
if (any(ri4 /= (/1_4, 2_4, 3_4, 4_4, 5_4, 16_4, 17_4, 18_4, 19_4/))) &
|
106 |
|
|
& call abort
|
107 |
|
|
|
108 |
|
|
vi8 = (/(i+10,i=1,9)/)
|
109 |
|
|
i8 = reshape((/1_8, -1_8, 2_8, -2_8, 3_8, -3_8, 4_8, -4_8, 5_8/), shape(i8))
|
110 |
|
|
ri8 = pack(i8,i8>0,vi8)
|
111 |
|
|
if (any(ri8 /= (/1_8, 2_8, 3_8, 4_8, 5_8, 16_8, 17_8, 18_8, 19_8/))) &
|
112 |
|
|
& call abort
|
113 |
|
|
|
114 |
|
|
|
115 |
|
|
end program main
|