3 ! Test the fix for PR41113 and PR41117, in which unnecessary calls
4 ! to internal_pack and internal_unpack were being generated.
6 ! Contributed by Joost VandeVondele <jv244@cam.ac.uk>
10 REAL :: data(10) = [(i
, i
= 1, 10)]
13 SUBROUTINE S1(data, i
, chksum
)
14 REAL, DIMENSION(*) :: data
16 real :: subsum
, chksum
19 subsum
= subsum
+ data(j
)
21 if (abs(subsum
- chksum
) > 1e-6) STOP 1
29 real :: data1(10) = [(i
, i
= 1, 10)]
30 REAL :: data(-4:5,-4:5) = reshape ([(real(i
), i
= 1, 100)], [10,10])
33 CALL S1(d
%data, 10, sum (d
%data))
34 CALL S1(data1
, 10, sum (data1
))
38 CALL S1(data(:,i
), 10, sum (data(:,i
)))
41 ! With the fix for PR41113/7 this is the only time that _internal_pack
42 ! was called. The final part of the fix for PR43072 put paid to it too.
44 CALL S1(data(-2:,i
), 8, sum (data(-2:,i
)))
47 CALL S1(data(:,i
:i
+1), 20, sum (reshape (data(:,i
:i
+1), [20])))
50 CALL S1(data(2,i
), 1, data(2,i
))