Avviso di runtime Fortran: array temporaneo

Question

Ricevo l'avviso di runtime fortran "Un array temporaneo è stato creato" durante l'esecuzione del mio codice (compilato con gfortran) e vorrei sapere se esiste un modo migliore per risolvere questo avviso.

mio codice originale è qualcosa di simile:

allocate(flx_est(lsign,3)) 
allocate(flx_err(lsign,3)) 
do i=1,lsign 
call combflx_calc(flx_est(i,:),flx_err(i,:)) 
enddo 

all'interno della subroutine definisco le variabili in questo modo:

subroutine combflx_calc(flx_est,flx_err) 
use,intrinsic    :: ISO_Fortran_env, only: real64 
implicit none 
real(real64),intent(inout) :: flx_est(3),flx_err(3) 

flux_est e flx_err vettori possono cambiare all'interno della subroutine a seconda delle diverse condizioni e Devo aggiornare i loro valori di conseguenza.

Fortran non sembra gradire questa struttura. Posso risolverlo definendo variabili temporanee:

tmp_flx_est=flx_est(i,:) 
tmp_flx_err=flx_err(i,:) 
call combflx_calc(tmp_flx_est,tmp_flx_err) 
flx_est(i,:)=tmp_flx_est 
flx_err(i,:)=tmp_flx_err 

Ma mi sembra un modo abbastanza sciocco per risolverlo.

Come puoi vedere, non sono un esperto di Fortran, quindi qualsiasi aiuto è più che benvenuto.

Answer 1

Un modo è quello di superare una forma presupposta matrice

real(real64),intent(inout) :: flx_est(:),flx_err(:) 

l'altra è per scambiare le dimensioni della matrice, in modo che è possibile passare una sezione contigua della matrice 2D.

call combflx_calc(flx_est(:,i),flx_err(:,i)) 

Il problema è che la dimensione esplicita argomenti fittizi della vostra procedura (var(n)) richiedono array contigui. Gli array di forme presunte possono avere qualche passo.

Answer 2

L'array temporaneo è in fase di creazione perché si passa una matrice a gradini alla subroutine. Gli array di Fortran sono colonne maggiori, quindi l'indice più a sinistra varia più velocemente in un array, o, meglio detto, l'indice più a sinistra è contiguo in memoria e ogni variabile a destra si trova su quelli a sinistra.

Quando si chiama

call combflx_calc(flx_est(i,:),flx_err(i,:)) 

Queste fette sono le matrici del vostro 3-vettore strided dalla lunghezza di lsign. La subroutine si aspetta che le variabili di una singola dimensione contigue nella memoria, che non è la variabile che si passa in esso. Pertanto, è necessario eseguire un temporaneo per il subroutine su cui operare e quindi copiato nuovamente nella sezione dell'array.

La tua "correzione" non cambia questo, semplicemente non ti avvisa più di un temporaneo perché stai utilizzando una variabile creata in modo esplicito piuttosto che il runtime che lo sta facendo per te.

La risposta di Vladimir offre opzioni per evitare il temporaneo, quindi non le duplicherò qui.