Skip to content
Snippets Groups Projects
Commit 11c2edfa authored by Antoine Cyril David Hoffmann's avatar Antoine Cyril David Hoffmann :seedling:
Browse files

Correction of the normalization of the c2r fourier transform

parent 23f28976
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 6 additions and 3 deletions
src/fourier_mod.F90
+ 6
− 3
View file @ 11c2edfa
...@@ -92,7 +92,8 @@ MODULE fourier ...@@ -92,7 +92,8 @@ MODULE fourier
ENDDO ENDDO
call fftw_mpi_execute_dft_c2r(planb, cmpx_data_f, real_data_f) call fftw_mpi_execute_dft_c2r(planb, cmpx_data_f, real_data_f)
call fftw_mpi_execute_dft_c2r(planb, cmpx_data_g, real_data_g) call fftw_mpi_execute_dft_c2r(planb, cmpx_data_g, real_data_g)
bracket_sum_r = bracket_sum_r + real_data_f*inv_Ny*inv_Nx * real_data_g*inv_Ny*inv_Nx ! bracket_sum_r = bracket_sum_r + real_data_f*inv_Ny*inv_Nx * real_data_g*inv_Ny*inv_Nx
bracket_sum_r = bracket_sum_r + real_data_f * real_data_g*inv_Ny*inv_Nx
! Second term -df/dy x dg/dx ! Second term -df/dy x dg/dx
DO ikx = ikxs, ikxe DO ikx = ikxs, ikxe
DO iky = ikys, ikye DO iky = ikys, ikye
...@@ -104,7 +105,8 @@ MODULE fourier ...@@ -104,7 +105,8 @@ MODULE fourier
ENDDO ENDDO
call fftw_mpi_execute_dft_c2r(planb, cmpx_data_f, real_data_f) call fftw_mpi_execute_dft_c2r(planb, cmpx_data_f, real_data_f)
call fftw_mpi_execute_dft_c2r(planb, cmpx_data_g, real_data_g) call fftw_mpi_execute_dft_c2r(planb, cmpx_data_g, real_data_g)
bracket_sum_r = bracket_sum_r - real_data_f*inv_Ny*inv_Nx * real_data_g*inv_Ny*inv_Nx ! bracket_sum_r = bracket_sum_r - real_data_f*inv_Ny*inv_Nx * real_data_g*inv_Ny*inv_Nx
bracket_sum_r = bracket_sum_r - real_data_f * real_data_g*inv_Ny*inv_Nx
END SUBROUTINE poisson_bracket_and_sum END SUBROUTINE poisson_bracket_and_sum
!!! Compute the poisson bracket of [F,G] to real space !!! Compute the poisson bracket of [F,G] to real space
...@@ -122,7 +124,8 @@ SUBROUTINE convolve_and_add( F_, G_) ...@@ -122,7 +124,8 @@ SUBROUTINE convolve_and_add( F_, G_)
ENDDO ENDDO
call fftw_mpi_execute_dft_c2r(planb, cmpx_data_f, real_data_f) call fftw_mpi_execute_dft_c2r(planb, cmpx_data_f, real_data_f)
call fftw_mpi_execute_dft_c2r(planb, cmpx_data_g, real_data_g) call fftw_mpi_execute_dft_c2r(planb, cmpx_data_g, real_data_g)
bracket_sum_r = bracket_sum_r + real_data_f*inv_Ny*inv_Nx * real_data_g*inv_Ny*inv_Nx ! bracket_sum_r = bracket_sum_r + real_data_f*inv_Ny*inv_Nx * real_data_g*inv_Ny*inv_Nx
bracket_sum_r = bracket_sum_r + real_data_f * real_data_g*inv_Ny*inv_Nx
END SUBROUTINE convolve_and_add END SUBROUTINE convolve_and_add
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment