### Regular versus recursive functions in Fortran

This is a short post that contrasts recursive versus regular Fortran functions and gives an example factorial $n!$ recursive function.

### Basics of a Fortran function

In Fortran, unlike many other languages, a user-defined function can only calculate a single value such as an integer, an array, or a character string. Also, the name of the function should also be the name of the value that it calculates by assigning it somewhere inside the function (see example below). Lastly, a return statement is not needed in a Fortran function unless the end of the function is always reached every time it is invoked.

function add(num1, num2) ! a regular function that adds two numbers implicit none intent(in) :: num1, num2 real :: num1, num2, add !the function name "add" must be assigned to the return value add = num1 + num2 end function add

### What makes a recursive Fortran function different?

Recursive functions were added to Fortran 90 and later versions. The main syntactical difference with normal Fortran functions is the intrinsic "recursive" declaration and "result()" assignment that both go in the function declaration line. For example, using the factorial as an example:

recursive function factorial(n) result(fact) implicit none intent(in) :: n integer (kind=8):: fact, n

if (n == 0) then fact = 1 else fact = n * factorial(n-1) endif

end function factorial

**recursive**". Also, the returned value unlike in regular/non-recursive Fortran functions, is declared on the first line as well using the intrinsic "

**result()**". In this case the actual value that is returned is "

**fact**" when the function "

**factorial**" is called. Whereas in non-recursive Fortran functions the name of the function is also the name of the returned value. To use a recursive function in a program however one must still call on the function by its name but in this case the return value is distinguished from the function value. Here is a a program using the above function:

program fact implicit none integer :: n, factorial

print *, "Enter a positive integer: " read (*,*) n print *, 'factorial = ', factorial(n)

end program fact

- Guide to Fortran 2008 Programming
- Numerical Recipes in Fortran 90: Volume 2, Volume 2 of Fortran Numerical Recipes: The Art of Parallel Scientific Computing

## Comments

## Post a Comment