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

added public statements to fftw functions and constants after whole math module is private now

This commit is contained in:
Martin Diehl 2013-06-06 09:25:27 +00:00
parent 622e2dcf15
commit 06b977ff61
1 changed files with 13 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
code/math.f90
View File

@ -31,7 +31,7 @@ module math
pInt pInt
implicit none implicit none
private ! because FFTW is included in math.f90 private
real(pReal), parameter, public :: PI = 3.14159265358979323846264338327950288419716939937510_pReal !< ratio of a circle's circumference to its diameter real(pReal), parameter, public :: PI = 3.14159265358979323846264338327950288419716939937510_pReal !< ratio of a circle's circumference to its diameter
real(pReal), parameter, public :: INDEG = 180.0_pReal/PI !< conversion from radian into degree real(pReal), parameter, public :: INDEG = 180.0_pReal/PI !< conversion from radian into degree
real(pReal), parameter, public :: INRAD = PI/180.0_pReal !< conversion from degree into radian real(pReal), parameter, public :: INRAD = PI/180.0_pReal !< conversion from degree into radian
@ -225,6 +225,17 @@ real(pReal), dimension(4,36), parameter, private :: &
math_rotate_forward3333 math_rotate_forward3333
#ifdef Spectral #ifdef Spectral
public :: & public :: &
fftw_set_timelimit, &
fftw_plan_dft_3d, &
fftw_plan_many_dft_r2c, &
fftw_plan_many_dft_c2r, &
fftw_plan_with_nthreads, &
fftw_init_threads, &
fftw_alloc_complex, &
fftw_execute_dft, &
fftw_execute_dft_r2c, &
fftw_execute_dft_c2r, &
fftw_destroy_plan, &
math_curlFFT, & math_curlFFT, &
math_gradFFT, & math_gradFFT, &
math_divergenceFFT, & math_divergenceFFT, &
@ -346,12 +357,11 @@ subroutine math_init
endif endif
! +++ check rotation sense of q and R +++ ! +++ check rotation sense of q and R +++
q = math_qRand() ! random quaternion q = math_qRand() ! random quaternion
call halton(3_pInt,v) ! random vector call halton(3_pInt,v) ! random vector
R = math_qToR(q) R = math_qToR(q)
if (any(abs(math_mul33x3(R,v) - math_qRot(q,v)) > tol_math_check)) then if (any(abs(math_mul33x3(R,v) - math_qRot(q,v)) > tol_math_check)) then
write(6,'(a,4(f8.3,x))') 'q',q write(6,'(a,4(f8.3,1x))') 'q',q
call IO_error(409_pInt) call IO_error(409_pInt)
endif endif