Il Linguaggio Fortran 90/95

5.15 Procedure ricorsive 283
RECURSIVE SUBROUTINE quick_sort(a)
IMPLICIT NONE
INTEGER, DIMENSION(:), INTENT(inout) :: a
END SUBROUTINE quick_sort
END INTERFACE
PRINT "(’Lista di partenza:’,/,(I2,1X),/)", a
CALL quick_sort(a)
PRINT "(’Lista ordinata:’,/,(I2,1X))", a
END PROGRAM prova_quick_sort
L’output del programmino di test fornisce questo risultato:
Lista di partenza:
11 6 2 3 13 0 24 -2 9 4
Lista ordinata:
-2 0 2 3 4 6 9 11 13 24
a riprova della bontà della routine.
In tutti gli esempi sviluppati in questo paragrafo una procedura ricorsiva era tale perché
invocava direttamente sé stessa (si parla in tal caso di ricorsione diretta). Tuttavia, una procedura
può essere invocata anche per ricorsione indiretta, cioè può chiamare sé stessa attraverso
l’invocazione di un’altra procedura che la richiami. Per illustrare questo concetto può essere
utile riferirsi ad un esempio concreto. Si supponga si voler eseguire l’integrazione di una funzione
bidimensionale f(x, y) su di un dominio rettangolare avendo, però, a disposizione soltanto
una procedura di integrazione monodimensionale del tipo seguente:
RECURSIVE FUNCTION integrate(f,bounds)
! Scopo: Integrare f da bounds(1) a bounds(2)
IMPLICIT NONE
REAL :: integrate
INTERFACE
FUNCTION f(x)
REAL :: f
REAL, INTENT(IN) :: x
END FUNCTION f
END INTERFACE
REAL, DIMENSION(2),INTENT(IN) :: bounds
...
END FUNCTION integrate
Ebbene, in questo caso si potrebbe scrivere una funzione di modulo che riceva il valore di x
come argomento ed il valore di y dal modulo stesso per host association:
MODULE func