Il Linguaggio Fortran 90/95

438 Operazioni di I/O su file
Al fine di confrontare l’efficienza delle operazioni di gestione di file formattati e non formattati
si consideri il seguente programma il quale crea due file contenenti ciascuno 5000 record con
quattro valori reali in doppia precisione per linea, dei quali però uno solo è di tipo formattato
mentre l’altro è di tipo non formattato. Una volta generati in maniera random i record dei
due file, essi vengono letti all’interno di due cicli separati e viene misurato il tempo speso dal
programma per effettuare questa operazione.
PROGRAM accesso_diretto
IMPLICIT NONE
INTEGER,PARAMETER :: single=SELECTED_REAL_KIND(p=6)
INTEGER,PARAMETER :: double=SELECTED_REAL_KIND(p=14)
INTEGER,PARAMETER :: max_records = 50000 ! max numero di record
INTEGER,PARAMETER :: number_of_reads = 100000 ! numero di record letti
INTEGER :: i
INTEGER :: length_fmt = 80 ! lunghezza di ciascun record formattato
INTEGER :: length_unf ! lunghezza di ciascun record non formattato
INTEGER :: irec ! puntatore al record
REAL(KIND=single) :: value
REAL(KIND=single) :: time_fmt ! tempo per lettura formattata
REAL(KIND=single) :: time_unf ! tempo per lettura non formattata
REAL(KIND=single) :: time_begin, time_end
REAL(KIND=double),DIMENSION(4) :: values ! valori del record
! Valutazione della lunghezza di ciascun record nel file non formattato
INQUIRE(IOLENGTH=length_unf) values
WRITE(*,"(A,I2)") "Lunghezza dei record non formattati: ", &
length_unf
WRITE(*,"(A,I2)") "Lunghezza dei record non formattati: ", &
length_fmt
! Connessione dei file
OPEN(UNIT=11,FILE="C:/MYPROG/OUT/myfile.fmt",ACCESS="DIRECT", &
FORM="UNFORMATTED",STATUS="REPLACE",ACTION="READWRITE", &
RECL=length_unf)
OPEN(UNIT=12,FILE="C:/MYPROG/OUT/myfile.unf",ACCESS="DIRECT", &
FORM="FORMATTED",STATUS="REPLACE",ACTION="READWRITE", &
RECL=length_fmt)
! Generazione dei record e loro inserimento nei rispettivi file
DO i=1,max_records
CALL RANDOM_NUMBER(values)
WRITE(11,REC=i) values
WRITE(12,"(4ES20.14)",REC=i) values
END DO
! Processing dei record non formattati
CALL CPU_TIME(time_begin)