|
111
|
1 ! { dg-do run }
|
|
|
2 ! { dg-add-options ieee }
|
|
|
3 ! { dg-skip-if "NaN not supported" { spu-*-* } }
|
|
|
4 real :: a(3,3), b(3), nan, minf, pinf, h
|
|
|
5 logical :: l, l2
|
|
|
6 logical :: l3(3,3), l4(3,3), l5(3,3)
|
|
|
7
|
|
|
8 nan = 0.0
|
|
|
9 minf = 0.0
|
|
|
10 pinf = 0.0
|
|
|
11 nan = 0.0/nan
|
|
|
12 minf = -1.0/minf
|
|
|
13 pinf = 1.0/pinf
|
|
|
14 h = huge(h)
|
|
|
15 l = .false.
|
|
|
16 l2 = .true.
|
|
|
17 l3 = .false.
|
|
|
18 l4 = .true.
|
|
|
19 l5 = .true.
|
|
|
20 l5(1,1) = .false.
|
|
|
21 l5(1,2) = .false.
|
|
|
22 l5(2,3) = .false.
|
|
|
23 a = reshape ((/ nan, nan, nan, pinf, pinf, pinf, pinf, minf, pinf /), (/ 3, 3 /))
|
|
|
24 if (minval (a).ne.minf) call abort
|
|
|
25 if (any (minloc (a).ne.(/ 2, 3 /))) call abort
|
|
|
26 b = minval (a, dim = 1)
|
|
|
27 if (.not.isnan(b(1))) call abort
|
|
|
28 b(1) = 0.0
|
|
|
29 if (any (b.ne.(/ 0.0, pinf, minf /))) call abort
|
|
|
30 if (any (minloc (a, dim = 1).ne.(/ 1, 1, 2 /))) call abort
|
|
|
31 b = minval (a, dim = 2)
|
|
|
32 if (any (b.ne.(/ pinf, minf, pinf /))) call abort
|
|
|
33 if (any (minloc (a, dim = 2).ne.(/ 2, 3, 2 /))) call abort
|
|
|
34 if (minval (a, mask = l).ne.h) call abort
|
|
|
35 if (any (minloc (a, mask = l).ne.(/ 0, 0 /))) call abort
|
|
|
36 b = minval (a, dim = 1, mask = l)
|
|
|
37 if (any (b.ne.(/ h, h, h /))) call abort
|
|
|
38 if (any (minloc (a, dim = 1, mask = l).ne.(/ 0, 0, 0 /))) call abort
|
|
|
39 b = minval (a, dim = 2, mask = l)
|
|
|
40 if (any (b.ne.(/ h, h, h /))) call abort
|
|
|
41 if (any (minloc (a, dim = 2, mask = l).ne.(/ 0, 0, 0 /))) call abort
|
|
|
42 if (minval (a, mask = l3).ne.h) call abort
|
|
|
43 if (any (minloc (a, mask = l3).ne.(/ 0, 0 /))) call abort
|
|
|
44 b = minval (a, dim = 1, mask = l3)
|
|
|
45 if (any (b.ne.(/ h, h, h /))) call abort
|
|
|
46 if (any (minloc (a, dim = 1, mask = l3).ne.(/ 0, 0, 0 /))) call abort
|
|
|
47 b = minval (a, dim = 2, mask = l3)
|
|
|
48 if (any (b.ne.(/ h, h, h /))) call abort
|
|
|
49 if (any (minloc (a, dim = 2, mask = l3).ne.(/ 0, 0, 0 /))) call abort
|
|
|
50 if (minval (a, mask = l2).ne.minf) call abort
|
|
|
51 if (minval (a, mask = l4).ne.minf) call abort
|
|
|
52 if (any (minloc (a, mask = l2).ne.(/ 2, 3 /))) call abort
|
|
|
53 if (any (minloc (a, mask = l4).ne.(/ 2, 3 /))) call abort
|
|
|
54 b = minval (a, dim = 1, mask = l2)
|
|
|
55 if (.not.isnan(b(1))) call abort
|
|
|
56 b(1) = 0.0
|
|
|
57 if (any (b.ne.(/ 0.0, pinf, minf /))) call abort
|
|
|
58 if (any (minloc (a, dim = 1, mask = l2).ne.(/ 1, 1, 2 /))) call abort
|
|
|
59 b = minval (a, dim = 2, mask = l2)
|
|
|
60 if (any (b.ne.(/ pinf, minf, pinf /))) call abort
|
|
|
61 if (any (minloc (a, dim = 2, mask = l2).ne.(/ 2, 3, 2 /))) call abort
|
|
|
62 b = minval (a, dim = 1, mask = l4)
|
|
|
63 if (.not.isnan(b(1))) call abort
|
|
|
64 b(1) = 0.0
|
|
|
65 if (any (b.ne.(/ 0.0, pinf, minf /))) call abort
|
|
|
66 if (any (minloc (a, dim = 1, mask = l2).ne.(/ 1, 1, 2 /))) call abort
|
|
|
67 b = minval (a, dim = 2, mask = l4)
|
|
|
68 if (any (b.ne.(/ pinf, minf, pinf /))) call abort
|
|
|
69 if (any (minloc (a, dim = 2, mask = l2).ne.(/ 2, 3, 2 /))) call abort
|
|
|
70 if (minval (a, mask = l5).ne.pinf) call abort
|
|
|
71 if (any (minloc (a, mask = l5).ne.(/ 2, 2 /))) call abort
|
|
|
72 b = minval (a, dim = 1, mask = l5)
|
|
|
73 if (.not.isnan(b(1))) call abort
|
|
|
74 b(1) = 0.0
|
|
|
75 if (any (b.ne.(/ 0.0, pinf, pinf /))) call abort
|
|
|
76 if (any (minloc (a, dim = 1, mask = l5).ne.(/ 2, 2, 1 /))) call abort
|
|
|
77 b = minval (a, dim = 2, mask = l5)
|
|
|
78 if (any (b.ne.(/ pinf, pinf, pinf /))) call abort
|
|
|
79 if (any (minloc (a, dim = 2, mask = l5).ne.(/ 3, 2, 2 /))) call abort
|
|
|
80 a = nan
|
|
|
81 if (.not.isnan(minval (a))) call abort
|
|
|
82 if (minval (a, mask = l).ne.h) call abort
|
|
|
83 if (.not.isnan(minval (a, mask = l2))) call abort
|
|
|
84 if (minval (a, mask = l3).ne.h) call abort
|
|
|
85 if (.not.isnan(minval (a, mask = l4))) call abort
|
|
|
86 if (.not.isnan(minval (a, mask = l5))) call abort
|
|
|
87 if (any (minloc (a).ne.(/ 1, 1 /))) call abort
|
|
|
88 if (any (minloc (a, mask = l).ne.(/ 0, 0 /))) call abort
|
|
|
89 if (any (minloc (a, mask = l2).ne.(/ 1, 1 /))) call abort
|
|
|
90 if (any (minloc (a, mask = l3).ne.(/ 0, 0 /))) call abort
|
|
|
91 if (any (minloc (a, mask = l4).ne.(/ 1, 1 /))) call abort
|
|
|
92 if (any (minloc (a, mask = l5).ne.(/ 2, 1 /))) call abort
|
|
|
93 a = pinf
|
|
|
94 if (minval (a).ne.pinf) call abort
|
|
|
95 if (minval (a, mask = l).ne.h) call abort
|
|
|
96 if (minval (a, mask = l2).ne.pinf) call abort
|
|
|
97 if (minval (a, mask = l3).ne.h) call abort
|
|
|
98 if (minval (a, mask = l4).ne.pinf) call abort
|
|
|
99 if (minval (a, mask = l5).ne.pinf) call abort
|
|
|
100 if (any (minloc (a).ne.(/ 1, 1 /))) call abort
|
|
|
101 if (any (minloc (a, mask = l).ne.(/ 0, 0 /))) call abort
|
|
|
102 if (any (minloc (a, mask = l2).ne.(/ 1, 1 /))) call abort
|
|
|
103 if (any (minloc (a, mask = l3).ne.(/ 0, 0 /))) call abort
|
|
|
104 if (any (minloc (a, mask = l4).ne.(/ 1, 1 /))) call abort
|
|
|
105 if (any (minloc (a, mask = l5).ne.(/ 2, 1 /))) call abort
|
|
|
106 a = nan
|
|
|
107 a(1,3) = pinf
|
|
|
108 if (minval (a).ne.pinf) call abort
|
|
|
109 if (minval (a, mask = l).ne.h) call abort
|
|
|
110 if (minval (a, mask = l2).ne.pinf) call abort
|
|
|
111 if (minval (a, mask = l3).ne.h) call abort
|
|
|
112 if (minval (a, mask = l4).ne.pinf) call abort
|
|
|
113 if (minval (a, mask = l5).ne.pinf) call abort
|
|
|
114 if (any (minloc (a).ne.(/ 1, 3 /))) call abort
|
|
|
115 if (any (minloc (a, mask = l).ne.(/ 0, 0 /))) call abort
|
|
|
116 if (any (minloc (a, mask = l2).ne.(/ 1, 3 /))) call abort
|
|
|
117 if (any (minloc (a, mask = l3).ne.(/ 0, 0 /))) call abort
|
|
|
118 if (any (minloc (a, mask = l4).ne.(/ 1, 3 /))) call abort
|
|
|
119 if (any (minloc (a, mask = l5).ne.(/ 1, 3 /))) call abort
|
|
|
120 end
|
