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
|
! { dg-do run }
! PR48615 Invalid UP/DOWN rounding with E and ES descriptors
! Test case provided by Thomas Henlich.
program pr48615
call checkfmt("(RU,F17.0)", 2.5, " 3.")
call checkfmt("(RU,-1P,F17.1)", 2.5, " 0.3")
call checkfmt("(RU,E17.1)", 2.5, " 0.3E+01")
call checkfmt("(RU,1P,E17.0)", 2.5, " 3.E+00")
call checkfmt("(RU,ES17.0)", 2.5, " 3.E+00")
call checkfmt("(RU,EN17.0)", 2.5, " 3.E+00")
call checkfmt("(RU,F2.0)", 2.0, "2.")
call checkfmt("(RU,F6.4)", 2.0, "2.0000")
call checkfmt("(RU,1P,E6.0E2)", 2.0, "2.E+00")
call checkfmt("(RU,1P,E7.1E2)", 2.5, "2.5E+00")
call checkfmt("(RU,1P,E10.4E2)", 2.5, "2.5000E+00")
call checkfmt("(RU,1P,G6.0E2)", 2.0, "2.E+00")
call checkfmt("(RU,1P,G10.4E2)", 2.3456e5, "2.3456E+05")
call checkfmt("(RU,F2.0)", 0.09, "1.") ! 0.
call checkfmt("(RD,F3.0)", -0.09, "-1.") ! -0.
call checkfmt("(RU,F2.0)", 2.0, "2.") ! 3.
call checkfmt("(RD,F3.0)", -2.0, "-2.") ! -3.
call checkfmt("(RU,F6.4)", 2.0, "2.0000") ! 2.0001
call checkfmt("(RD,F7.4)", -2.0, "-2.0000") ! -2.0001
call checkfmt("(RU,1P,E6.0E2)", 2.0, "2.E+00") ! 3.E+00
call checkfmt("(RD,1P,E7.0E2)", -2.0, "-2.E+00") ! -3.E+00
call checkfmt("(RU,1P,E7.1E2)", 2.5, "2.5E+00") ! 2.6E+00
call checkfmt("(RD,1P,E8.1E2)", -2.5, "-2.5E+00") ! -2.6E+00
call checkfmt("(RU,1P,E10.4E2)", 2.5, "2.5000E+00") ! 2.5001E+00
call checkfmt("(RD,1P,E11.4E2)", -2.5, "-2.5000E+00") ! -2.5001E+00
call checkfmt("(RU,1P,G6.0E2)", 2.0, "2.E+00") ! 3.E+00
call checkfmt("(RD,1P,G7.0E2)", -2.0, "-2.E+00") ! -3.E+00
call checkfmt("(RU,1P,G10.4E2)", 2.3456e5, "2.3456E+05") ! 2.3457E+05
call checkfmt("(RD,1P,G11.4E2)", -2.3456e5, "-2.3456E+05") ! -2.3457E+05
call checkfmt("(RD,F17.0)", 2.5, " 2.")
call checkfmt("(RD,-1P,F17.1)", 2.5, " 0.2")
call checkfmt("(RD,E17.1)", 2.5, " 0.2E+01")
call checkfmt("(RD,1P,E17.0)", 2.5, " 2.E+00")
call checkfmt("(RD,ES17.0)", 2.5, " 2.E+00")
call checkfmt("(RD,EN17.0)", 2.5, " 2.E+00")
call checkfmt("(RC,F17.0)", 2.5, " 3.")
call checkfmt("(RC,-1P,F17.1)", 2.5, " 0.3")
call checkfmt("(RC,E17.1)", 2.5, " 0.3E+01")
call checkfmt("(RC,1P,E17.0)", 2.5, " 3.E+00")
call checkfmt("(RC,ES17.0)", 2.5, " 3.E+00")
call checkfmt("(RC,EN17.0)", 2.5, " 3.E+00")
call checkfmt("(RN,F17.0)", 2.5, " 2.")
call checkfmt("(RN,-1P,F17.1)", 2.5, " 0.2")
call checkfmt("(RN,E17.1)", 2.5, " 0.2E+01")
call checkfmt("(RN,1P,E17.0)", 2.5, " 2.E+00")
call checkfmt("(RN,ES17.0)", 2.5, " 2.E+00")
call checkfmt("(RN,EN17.0)", 2.5, " 2.E+00")
call checkfmt("(RZ,F17.0)", 2.5, " 2.")
call checkfmt("(RZ,-1P,F17.1)", 2.5, " 0.2")
call checkfmt("(RZ,E17.1)", 2.5, " 0.2E+01")
call checkfmt("(RZ,1P,E17.0)", 2.5, " 2.E+00")
call checkfmt("(RZ,ES17.0)", 2.5, " 2.E+00")
call checkfmt("(RZ,EN17.0)", 2.5, " 2.E+00")
call checkfmt("(RZ,F17.0)", -2.5, " -2.")
call checkfmt("(RZ,-1P,F17.1)", -2.5, " -0.2")
call checkfmt("(RZ,E17.1)", -2.5, " -0.2E+01")
call checkfmt("(RZ,1P,E17.0)", -2.5, " -2.E+00")
call checkfmt("(RZ,ES17.0)", -2.5, " -2.E+00")
call checkfmt("(RZ,EN17.0)", -2.5, " -2.E+00")
call checkfmt("(RN,F17.0)", -2.5, " -2.")
call checkfmt("(RN,-1P,F17.1)", -2.5, " -0.2")
call checkfmt("(RN,E17.1)", -2.5, " -0.2E+01")
call checkfmt("(RN,1P,E17.0)", -2.5, " -2.E+00")
call checkfmt("(RN,ES17.0)", -2.5, " -2.E+00")
call checkfmt("(RN,EN17.0)", -2.5, " -2.E+00")
call checkfmt("(RC,F17.0)", -2.5, " -3.")
call checkfmt("(RC,-1P,F17.1)", -2.5, " -0.3")
call checkfmt("(RC,E17.1)", -2.5, " -0.3E+01")
call checkfmt("(RC,1P,E17.0)", -2.5, " -3.E+00")
call checkfmt("(RC,ES17.0)", -2.5, " -3.E+00")
call checkfmt("(RC,EN17.0)", -2.5, " -3.E+00")
call checkfmt("(RU,E17.1)", nearest(2.0, 1.0), " 0.3E+01")
call checkfmt("(RD,E17.1)", nearest(3.0, -1.0), " 0.2E+01")
contains
subroutine checkfmt(fmt, x, cmp)
character(len=*), intent(in) :: fmt
real, intent(in) :: x
character(len=*), intent(in) :: cmp
character(len=20) :: s
write(s, fmt) x
if (s /= cmp) call abort
!if (s /= cmp) print "(a,1x,a,' expected: ',1x)", fmt, s, cmp
end subroutine
end program
|
