223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

some improvements on HDF5 test, polishing math

This commit is contained in:
Martin Diehl 2013-11-22 10:31:52 +00:00
parent e8de8f69a8
commit c9ce8d2c52
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
code/math.f90
View File

@ -2818,6 +2818,7 @@ end function math_rotate_forward3333
#ifdef Spectral #ifdef Spectral
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief calculates curl field using differentation in Fourier space !> @brief calculates curl field using differentation in Fourier space
!> @todo enable odd resolution
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
function math_curlFFT(geomdim,field) function math_curlFFT(geomdim,field)
use IO, only: & use IO, only: &
@ -2951,6 +2952,7 @@ end function math_curlFFT
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief calculates gradient field using differentation in Fourier space !> @brief calculates gradient field using differentation in Fourier space
!> @todo enable odd resolution
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
function math_gradFFT(geomdim,field) function math_gradFFT(geomdim,field)
use IO, only: & use IO, only: &
@ -3079,6 +3081,7 @@ end function math_gradFFT
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief calculates divergence field using integration in Fourier space !> @brief calculates divergence field using integration in Fourier space
!> @todo enable odd resolution
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
function math_divergenceFFT(geomdim,field) function math_divergenceFFT(geomdim,field)
use IO, only: & use IO, only: &
@ -3268,7 +3271,7 @@ function math_divergenceFDM(geomdim,order,field)
contains contains
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief ! small helper functions for indexing CAREFUL, index and location runs from !> @brief small helper functions for indexing. CAREFUL: index and location runs from
! 0 to N-1 (python style) ! 0 to N-1 (python style)
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
pure function periodic_location(idx,res) pure function periodic_location(idx,res)
@ -3284,7 +3287,7 @@ function math_divergenceFDM(geomdim,order,field)
end function periodic_location end function periodic_location
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief ! small helper functions for indexing CAREFUL, index and location runs from !> @brief small helper functions for indexing CAREFUL: index and location runs from
! 0 to N-1 (python style) ! 0 to N-1 (python style)
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
integer(pInt) pure function periodic_index(location,res) integer(pInt) pure function periodic_index(location,res)