Intrinsic Modules
Recent revisions of the Fortran standard support several intrinsic modules.
One must be downloaded, while the others can be USEd and will be supplied by the compiler.
ISO_VARYING_STRING
Fortran 95 did not support a variable-length string. A standardized module was defined to support a type VARYING_STRING. This module was never incorporated into the standard, so compilers do not include it, but an implementation can be found here. This module was mostly obsoleted by the variable-string capabilities of Fortran 2003, but does have a few features still lacking in the standard. A description is here. In addition to defining the standard character intrinsics for VARYING_STRING, it contains some additional functionality, such as GET and PUT to read and write a character from or into a string, REMOVE, REPLACE, SPLIT, and some other useful procedures. This module should be USEd like a programmer-written module.
Intrinsic Modules
The Fortran 2003 standard defined several intrinsics modules. They have a special form of USE:
USE, INTRINSIC :: <module>
ISO_FORTRAN_ENV
This module contains many useful variables for system parameters, some storage parameters, and KIND parameters. It also provides two intrinsic functions, COMPILER_OPTIONS and COMPILER_VERSION, which respectively return the command-line options and the compiler version used to compile the binary.
use iso_fortran_env
print *, "This executable was compiled with ",COMPILER_VERSION()
print *, "The options used were ",COMPILER_OPTIONS()
Particularly useful members of the ISO_FORTRAN_ENV module are predefined KIND parameter for specific types.
USE ISO_FORTRAN_ENV
INTEGER(int64) :: i,j
KIND Parameter | IEEE Type |
---|---|
int8 | 8-bit integer |
int16 | 16-bit integer |
int32 | 32-bit integer |
int64 | 64-bit integer |
real32 | 32-bit floating point |
real64 | 64-bit floating point |
real128 | 128-bit real |
Not all the IEEE KINDs may be supported in the hardware, particularly real128. Not all compilers support all the above KINDs, in which case it should set it to a negative value.
IEEE Modules
The three IEEE modules provide constants and procedures relating to floating-point hardware such as exceptions and arithmetic constants.
A good reference for these modules is provided by Intel for their
compiler. The NAG compiler also has useful
documentation.
IEEE Features
This module specifies the definitions of IEEE special bit patterns such as IEEE_Datatype and IEEE_Inf.
IEEE Exceptions
An exception occurs due to an illegal operation. This may include mathematically illegal operations such as taking the square root of a negative real number, dividing by zero, and so forth. Due to the finite range of floating-point numbers, other excepts are underflow and overflow. Mathematically illegal operations result in NaN (Not a Number), whereas overflow results in Inf. Operations are defined on NaNs, with the result of any arithmetic operation on a NaN being another NaN. Therefore, NaNs can easily propagate through results and it would be useful to catch them when they first occur. The IEEE Exceptions module can help with this. For example, the IEEE_SET_HALTING_MODE,HALTING) intrinsic would be invoked as follows:
USE, INTRINSIC :: IEEE_EXCEPTIONS
TYPE(IEEE_FLAG_TYPE) :: flag
LOGICAL :: halt=.true.
CALL IEEE_SET_HALTING_MODE(flag,halt)
Flag
is a variable of TYPE(IEEE_FLAG_TYPE), also defined in the module. It can be IEEE_DIVIDE_BY_ZERO, IEEE_INEXACT, IEEE_INVALID, IEEE_OVERFLOW, or IEEE_UNDERFLOW. The HALTING
argument is LOGICAL. If HALTING is set to .true.
the program will stop on occurrence of the specified MODE. Since this is often desired for INVALID (NaN), DIVIDE_BY_ZERO, and OVERFLOW, a predefined array is available.
program ieee
use, intrinsic :: ieee_exceptions
implicit none
type(ieee_flag_type), dimension(3) :: flag
logical :: halt=.true.
real :: x,y
flag=ieee_usual
call ieee_set_halting_mode(flag,halt)
x=-1.
y=sqrt(x)
print *, y
end program
IEEE Arithmetic
The IEEE_ARITHMETIC encompasses and extends the IEEE_EXCEPTIONS module. Several useful procedures and included. IEEE_SELECTED_REAL_KIND chooses only KINDs corresponding to IEEE-supported types. Other procedures determine whether arithmetic operations conform to IEEE specifications.
For more elegant handling of errors, the IEEE_IS_NAN, IEEE_IS_FINITE, and some others can be used to test a result and handle it in some manner other than halting the execution.
program ieee
use, intrinsic :: ieee_arithmetic
implicit none
real :: x,y
x=-1.
y=sqrt(x)
if ( ieee_is_nan(y) ) then
print *, "NaN"
else
print *, y
endif
end program
ISO_C_BINDING
Mixed-language programming is common, particularly mixing C with other languages. In the past, invoking C procedures from Fortran was tedious and error-prone, due to differences in conventions such as character termination, name-mangling of subprograms, and so forth. The ISO_C_BINDING module was added to simplify this.
A good reference for the content of this module is from
gfortran. Note that this module provides C equivalents to Fortran types for a given platform; in particular, the C standard does not specify the length of an int
, only a minimum, so on some platforms the default is 16 bits and on others it is 32 bits. The variables defining the correspondence can be used as KIND parameters. A subset of the most commonly used might include
Fortran Type | Module Name | C Type |
---|---|---|
INTEGER | C_INT | int |
INTEGER(int64) | C_INT64_T | int64_t |
REAL | C_FLOAT | float |
REAL(real64) | C_DOUBLE | double |
LOGICAL | C_BOOL | _Bool |
CHARACTER | C_CHAR | char |
If all the types can be matched, a C struct can be mapped to a Fortran type with the BIND(C) attribute.
USE ISO_C_BINDING
TYPE, BIND(C) :: myType
INTEGER(C_INT) :: i, j
REAL(C_DOUBLE) :: d
CHARACTER(KIND=C_CHAR) :: c
END TYPE
The KIND= keyword is required for CHARACTER because the default argument is the LEN. This would correspond to a C struct
struct {
int i, j;
double d;
char c;
} myType;
Some facts to keep in mind are that C arrays number from 0, and C strings are terminated with a NULL character (C_NULL_CHAR in the module).
Subprograms must also declare the BIND(C) attribute to set up C bindings. Caution is required since C generally passes by value, which makes a copy, whereas Fortran effectively passes by reference, i.e. the address of the memory location that holds the variable. Arrays must also receive special treatment. Pointers require particular care. For a very simple example, suppose we have a C function
int adder(int i, int* j)
and we wish to write Fortran bindings to it. The first variable is passed by value so we must add the VALUE attribute to its declaration in Fortran. The second argument is passed by reference, as is the default in Fortran. The Fortran declaration looks like
integer(c_int) function func(i,j)
use iso_c_binding, only: c_int
integer(c_int), VALUE :: i
integer(c_int) :: j
A good general discussion of Fortran-C interoperability, from which the above example is taken, is from gfortran.