Il mio problema è semplice: ho una funzione a 3 variabili (non fattorizzabile), di cui calcolarne la derivata seconda, e dovrei eseguire una fft solo in 2 direzioni, x ed y, per poi usare uno schema di Padè al 2° ordine per z. Le subroutine sono già scritte e testate su funzioni monodimensionali, la mia domanda è: come faccio a fare una fft su una sola variabile alla volta? In generale, come faccio a far girare una subroutine su una sola delle 3 variabili della funzione? Ho provato con dei semplici cicli 'do', ad esempio:
do iz=1, Nz
do iy=1, Ny
call fft(f,Nx,fhat)
end do
end do
ma il programma, nonostante giri senza alcun errore, mi restituisce un vettore vuoto. Aiutatemi, vi prego
