1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
|
program intrinsics_172
use iso_fortran_env, only: int32, int64
implicit none
logical :: mask(3, 4)
logical :: mask_(6, 9)
logical :: mask_3(3, 4, 3)
logical :: mask_4(4, 5, 3, 1)
integer(kind=int32) :: res_4
integer(kind=int64) :: res_8
integer, dimension(2,3) :: a, b
logical, dimension(2,3) :: mask2
integer, parameter :: c1 = count([.true., .false., .true., .false.])
integer(8), parameter :: c2 = count([.true., .false., .true., .false., .true., .false.], kind = 8)
mask = reshape([ .true., .false., .true., .false., &
.true., .false., .true., .false., &
.true., .false., .true., .true.], [3, 4])
mask_ = .false.
mask_(1, 1) = .true.
mask_(1, 2) = .true.
mask_(5, 1) = .true.
mask_(5, 2) = .true.
mask_(6, 7) = .true.
mask_(4, 8) = .true.
mask_3 = .false.
mask_3(1, 1, :) = .true.
mask_3(1, 2, :) = .true.
mask_3(3, 1, :) = .true.
mask_3(2, 2, :) = .true.
mask_3(2, 4, :) = .true.
mask_3(3, 3, :) = .true.
mask_4 = .true.
mask_4(1, 1, 1, 1) = .false.
mask_4(1, 1, 1, :) = .false.
mask_4(1, 1, 1, :) = .false.
mask_4(1, 2, 1, :) = .false.
mask_4(1, 2, 3, :) = .false.
print *, count(mask)
if (count(mask) /= 7) error stop
print *, count(mask_)
if (count(mask_) /= 6) error stop
print *, count(mask_3)
if (count(mask_3) /= 18) error stop
print *, count(mask, 1)
if (sum(count(mask, 1)) /= 7) error stop
print *, count(mask, 2)
if (sum(count(mask, 2)) /= 7) error stop
print *, sum(count(mask_, 1))
if (sum(count(mask_, 1)) /= 6) error stop
print *, shape(count(mask_, 1))
print *, sum(count(mask_, 2))
if (sum(count(mask_, 2)) /= 6) error stop
print *, shape(count(mask_, 2))
print *, sum(count(mask_3, 1))
if (sum(count(mask_3, 1)) /= 18) error stop
print *, shape(count(mask_3, 1))
print *, sum(count(mask_3, 2))
if (sum(count(mask_3, 2)) /= 18) error stop
print *, shape(count(mask_3, 2))
print *, sum(count(mask_3, 3))
if (sum(count(mask_3, 3)) /= 18) error stop
print *, sum(count(mask_4, 1))
if (sum(count(mask_4, 1)) /= 57) error stop
print *, shape(count(mask_4, 1))
res_4 = count(mask_4, kind=4)
if (res_4 /= 57) error stop
res_8 = count(mask_4, kind=8)
if (res_8 /= 57) error stop
a = reshape( [ 1, 2, 3, 4, 5, 6 ], [ 2, 3 ])
b = reshape( [ 0, 7, 3, 4, 5, 8 ], [ 2, 3 ])
print '(3i3)', a(1,:)
print '(3i3)', a(2,:)
print *
print '(3i3)', b(1,:)
print '(3i3)', b(2,:)
print *
mask2 = a.ne.b
print '(3l3)', mask2(1,:)
print '(3l3)', mask2(2,:)
print *
print '(3i3)', count(mask2)
if (count(mask2) /= 3) error stop
print *
print '(3i3)', count(mask2, 1)
print *
print '(3i3)', count(mask2, 2)
print*, c1
if (c1 /= 2) error stop
print*, c2
if (c2 /= 3) error stop
print*, kind(c2)
if (kind(c2) /= 8) error stop
end program
|
