2014-07-12 Paul Thomas <pault@gcc.gnu.org>
[official-gcc.git] / gcc / testsuite / gfortran.dg / implicit_pure_3.f90
blobd9d7734dab3525cba26267c65ed2f8efd0c9c972
1 ! { dg-do compile }
2 ! { dg-options "-O2 -fdump-tree-optimized" }
4 ! PR fortran/54556
6 ! Contributed by Joost VandeVondele
8 MODULE parallel_rng_types
10 IMPLICIT NONE
12 ! Global parameters in this module
13 INTEGER, PARAMETER :: dp=8
15 TYPE rng_stream_type
16 PRIVATE
17 CHARACTER(LEN=40) :: name
18 INTEGER :: distribution_type
19 REAL(KIND=dp), DIMENSION(3,2) :: bg,cg,ig
20 LOGICAL :: antithetic,extended_precision
21 REAL(KIND=dp) :: buffer
22 LOGICAL :: buffer_filled
23 END TYPE rng_stream_type
25 REAL(KIND=dp), DIMENSION(3,3) :: a1p0,a1p76,a1p127,&
26 a2p0,a2p76,a2p127,&
27 inv_a1,inv_a2
29 INTEGER, PARAMETER :: GAUSSIAN = 1,&
30 UNIFORM = 2
32 REAL(KIND=dp), PARAMETER :: norm = 2.328306549295727688e-10_dp,&
33 m1 = 4294967087.0_dp,&
34 m2 = 4294944443.0_dp,&
35 a12 = 1403580.0_dp,&
36 a13n = 810728.0_dp,&
37 a21 = 527612.0_dp,&
38 a23n = 1370589.0_dp,&
39 two17 = 131072.0_dp,& ! 2**17
40 two53 = 9007199254740992.0_dp,& ! 2**53
41 fact = 5.9604644775390625e-8_dp ! 1/2**24
44 CONTAINS
46 FUNCTION rn32(rng_stream) RESULT(u)
48 TYPE(rng_stream_type), POINTER :: rng_stream
49 REAL(KIND=dp) :: u
51 INTEGER :: k
52 REAL(KIND=dp) :: p1, p2
54 ! -------------------------------------------------------------------------
55 ! Component 1
57 p1 = a12*rng_stream%cg(2,1) - a13n*rng_stream%cg(1,1)
58 k = INT(p1/m1)
59 p1 = p1 - k*m1
60 IF (p1 < 0.0_dp) p1 = p1 + m1
61 rng_stream%cg(1,1) = rng_stream%cg(2,1)
62 rng_stream%cg(2,1) = rng_stream%cg(3,1)
63 rng_stream%cg(3,1) = p1
65 ! Component 2
67 p2 = a21*rng_stream%cg(3,2) - a23n*rng_stream%cg(1,2)
68 k = INT(p2/m2)
69 p2 = p2 - k*m2
70 IF (p2 < 0.0_dp) p2 = p2 + m2
71 rng_stream%cg(1,2) = rng_stream%cg(2,2)
72 rng_stream%cg(2,2) = rng_stream%cg(3,2)
73 rng_stream%cg(3,2) = p2
75 ! Combination
77 IF (p1 > p2) THEN
78 u = (p1 - p2)*norm
79 ELSE
80 u = (p1 - p2 + m1)*norm
81 END IF
83 IF (rng_stream%antithetic) u = 1.0_dp - u
85 END FUNCTION rn32
87 ! *****************************************************************************
88 FUNCTION rn53(rng_stream) RESULT(u)
90 TYPE(rng_stream_type), POINTER :: rng_stream
91 REAL(KIND=dp) :: u
93 u = rn32(rng_stream)
95 IF (rng_stream%antithetic) THEN
96 u = u + (rn32(rng_stream) - 1.0_dp)*fact
97 IF (u < 0.0_dp) u = u + 1.0_dp
98 ELSE
99 u = u + rn32(rng_stream)*fact
100 IF (u >= 1.0_dp) u = u - 1.0_dp
101 END IF
103 END FUNCTION rn53
105 END MODULE
107 ! { dg-final { scan-module-absence "parallel_rng_types" "IMPLICIT_PURE" } }
108 ! { dg-final { scan-tree-dump-times "rn32 \\(rng_stream" 3 "optimized" } }
109 ! { dg-final { cleanup-tree-dump "optimized" } }