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

cleaned

This commit is contained in:
Martin Diehl 2015-09-28 17:47:00 +00:00
parent b15eb5d787
commit 0f0b4eaa6b
1 changed files with 6 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
code/DAMASK_spectral_utilities.f90
View File

@ -386,7 +386,6 @@ end subroutine utilities_init
!> @details Sets the current reference stiffness to the stiffness given as an argument. !> @details Sets the current reference stiffness to the stiffness given as an argument.
!> If the gamma operator is precalculated, it is calculated with this stiffness. !> If the gamma operator is precalculated, it is calculated with this stiffness.
!> In case of a on-the-fly calculation, only the reference stiffness is updated. !> In case of a on-the-fly calculation, only the reference stiffness is updated.
!> The gamma operator is filtered depening on the filter selected in numerics.
!> Also writes out the current reference stiffness for restart. !> Also writes out the current reference stiffness for restart.
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine utilities_updateGamma(C,saveReference) subroutine utilities_updateGamma(C,saveReference)
@ -455,10 +454,9 @@ subroutine utilities_FFTtensorForward()
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! applying filter ! applying filter
do k = 1_pInt, grid3; do j = 1_pInt, grid(2); do i = 1_pInt,grid1Red forall(k = 1_pInt:grid3, j = 1_pInt:grid(2), i = 1_pInt:grid1Red) &
tensorField_fourier(1:3,1:3,i,j,k) = utilities_getFilter(xi2nd(1:3,i,j,k))* & tensorField_fourier(1:3,1:3,i,j,k) = utilities_getFilter(xi2nd(1:3,i,j,k))* &
tensorField_fourier(1:3,1:3,i,j,k) tensorField_fourier(1:3,1:3,i,j,k)
enddo; enddo; enddo
end subroutine utilities_FFTtensorForward end subroutine utilities_FFTtensorForward
@ -493,10 +491,9 @@ subroutine utilities_FFTscalarForward()
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! applying filter ! applying filter
do k = 1_pInt, grid3; do j = 1_pInt, grid(2); do i = 1_pInt,grid1Red forall(k = 1_pInt:grid3, j = 1_pInt:grid(2), i = 1_pInt:grid1Red) &
scalarField_fourier(i,j,k) = utilities_getFilter(xi2nd(1:3,i,j,k))* & scalarField_fourier(i,j,k) = utilities_getFilter(xi2nd(1:3,i,j,k))* &
scalarField_fourier(i,j,k) scalarField_fourier(i,j,k)
enddo; enddo; enddo
end subroutine utilities_FFTscalarForward end subroutine utilities_FFTscalarForward
@ -532,10 +529,9 @@ subroutine utilities_FFTvectorForward()
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! applying filter ! applying filter
do k = 1_pInt, grid3; do j = 1_pInt, grid(2); do i = 1_pInt,grid1Red forall(k = 1_pInt:grid3, j = 1_pInt:grid(2), i = 1_pInt:grid1Red) &
vectorField_fourier(1:3,i,j,k) = utilities_getFilter(xi2nd(1:3,i,j,k))* & vectorField_fourier(1:3,i,j,k) = utilities_getFilter(xi2nd(1:3,i,j,k))* &
vectorField_fourier(1:3,i,j,k) vectorField_fourier(1:3,i,j,k)
enddo; enddo; enddo
end subroutine utilities_FFTvectorForward end subroutine utilities_FFTvectorForward
@ -1123,7 +1119,7 @@ end function utilities_forwardField
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief calculates filter for fourier convolution depending on type given in numerics.config !> @brief calculates filter for fourier convolution depending on type given in numerics.config
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
complex(pReal) pure function utilities_getFilter(k) pure function utilities_getFilter(k)
use math, only: & use math, only: &
PI PI
use mesh, only: & use mesh, only: &
@ -1131,6 +1127,7 @@ complex(pReal) pure function utilities_getFilter(k)
implicit none implicit none
real(pReal), intent(in), dimension(3) :: k !< indices of frequency real(pReal), intent(in), dimension(3) :: k !< indices of frequency
complex(pReal) :: utilities_getFilter
select case (spectral_filter_ID) select case (spectral_filter_ID)
case (FILTER_NONE_ID) ! default, no weighting case (FILTER_NONE_ID) ! default, no weighting