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
|
! Copyright 2003 Free Software Foundation, Inc.
! Contributed by Kejia Zhao <kejia_zh@yahoo.com.cn>
!
!This file is part of the GNU Fortran 95 runtime library (libgfor).
!
!GNU libgfor is free software; you can redistribute it and/or
!modify it under the terms of the GNU Lesser General Public
!License as published by the Free Software Foundation; either
!version 2.1 of the License, or (at your option) any later version.
!
!GNU libgfor is distributed in the hope that it will be useful,
!but WITHOUT ANY WARRANTY; without even the implied warranty of
!MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
!GNU Lesser General Public License for more details.
!
!You should have received a copy of the GNU Lesser General Public
!License along with libgfor; see the file COPYING. If not,
!write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
!Boston, MA 02111-1307, USA.
!
function selected_int_kind (r)
implicit none
integer, intent (in) :: r
integer :: selected_int_kind
integer :: i
! Integer kind_range table
integer, parameter :: c = 4
type :: int_info
integer :: kind
integer :: range
end type int_info
type (int_info), parameter :: int_infos (c) = &
(/int_info (1, range (0_1)), &
int_info (2, range (0_2)), &
int_info (4, range (0_4)), &
int_info (8, range (0_8))/)
do i = 1, c
if (r <= int_infos (i) % range) then
selected_int_kind = int_infos (i) % kind
return
end if
end do
selected_int_kind = -1
return
end function
function selected_real_kind (p, r)
implicit none
integer, optional, intent (in) :: p, r
integer :: selected_real_kind
integer :: i, p2, r2
logical :: found_p, found_r
! Real kind_precision_range table
integer, parameter :: c = 2
type :: real_info
integer :: kind
integer :: precision
integer :: range
end type real_info
type (real_info) :: real_infos (c) = &
(/real_info (4, precision (0.0_4), range (0.0_4)), &
real_info (8, precision (0.0_8), range (0.0_8))/)
selected_real_kind = 0
p2 = 0
r2 = 0
found_p = .false.
found_r = .false.
if (present (p)) p2 = p
if (present (r)) r2 = r
! Assumes each type has a greater precision and range than previous one.
do i = 1, c
if (p2 <= real_infos (i) % precision) found_p = .true.
if (r2 <= real_infos (i) % range) found_r = .true.
if (found_p .and. found_r) then
selected_real_kind = real_infos (i) % kind
return
end if
end do
if (.not. (found_p)) selected_real_kind = selected_real_kind - 1
if (.not. (found_r)) selected_real_kind = selected_real_kind - 2
return
end function
|
