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

some comments (talking to Ricardo)

This commit is contained in:
Martin Diehl 2011-02-25 16:24:23 +00:00
parent a6d1e5c911
commit 138e9c43eb
2 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
code/mpie_spectral.f90
View File

@ -283,7 +283,7 @@ program mpie_spectral
xi(1) = real(k_s(1), pReal)/geomdimension(1) xi(1) = real(k_s(1), pReal)/geomdimension(1)
if (any(xi /= 0.0_pReal)) then if (any(xi /= 0.0_pReal)) then
do l = 1,3; do m = 1,3 do l = 1,3; do m = 1,3
xinormdyad(l,m) = xi(l)*xi(m)/sum(xi**2) xinormdyad(l,m) = xi(l)*xi(m)/sum(xi**2) ! Unit sphere, unit vectors in Fourier space
enddo; enddo enddo; enddo
temp33_Real = math_inv3x3(math_mul3333xx33(c0, xinormdyad)) temp33_Real = math_inv3x3(math_mul3333xx33(c0, xinormdyad))
else else
@ -511,7 +511,7 @@ program mpie_spectral
temp33_Complex(m,n) = sum(gamma_hat(1,1,1,m,n,:,:) *(workfft(i*2-1,j,k,:,:)& temp33_Complex(m,n) = sum(gamma_hat(1,1,1,m,n,:,:) *(workfft(i*2-1,j,k,:,:)&
+workfft(i*2 ,j,k,:,:)*img)) +workfft(i*2 ,j,k,:,:)*img))
enddo; enddo enddo; enddo
workfft(i*2-1,j,k,:,:) = real (temp33_Complex) workfft(i*2-1,j,k,:,:) = real (temp33_Complex) ! change of strain
workfft(i*2 ,j,k,:,:) = aimag(temp33_Complex) workfft(i*2 ,j,k,:,:) = aimag(temp33_Complex)
enddo; enddo; enddo enddo; enddo; enddo
else !use precalculated gamma-operator else !use precalculated gamma-operator
@ -520,7 +520,7 @@ program mpie_spectral
temp33_Complex(m,n) = sum(gamma_hat(i,j,k, m,n,:,:) *(workfft(i*2-1,j,k,:,:)& temp33_Complex(m,n) = sum(gamma_hat(i,j,k, m,n,:,:) *(workfft(i*2-1,j,k,:,:)&
+ workfft(i*2 ,j,k,:,:)*img)) + workfft(i*2 ,j,k,:,:)*img))
enddo; enddo enddo; enddo
workfft(i*2-1,j,k,:,:) = real (temp33_Complex) workfft(i*2-1,j,k,:,:) = real (temp33_Complex) ! change of strain
workfft(i*2 ,j,k,:,:) = aimag(temp33_Complex) workfft(i*2 ,j,k,:,:) = aimag(temp33_Complex)
enddo; enddo; enddo enddo; enddo; enddo
endif endif

2
code/prec_single.f90
View File

@ -9,7 +9,7 @@
integer, parameter :: pReal = selected_real_kind(6,37) ! 6 significant digits, up to 1e+-37 integer, parameter :: pReal = selected_real_kind(6,37) ! 6 significant digits, up to 1e+-37
integer, parameter :: pInt = selected_int_kind(9) ! up to +- 1e9 integer, parameter :: pInt = selected_int_kind(9) ! up to +- 1e9
integer, parameter :: pLongInt = 8 ! should be 64bit integer, parameter :: pLongInt = 8 ! should be 64bit
real(pReal), parameter :: tol_math_check = 1.0e-6_pReal real(pReal), parameter :: tol_math_check = 1.0e-5_pReal
real(pReal), parameter :: tol_gravityNodePos = 1.0e-36_pReal real(pReal), parameter :: tol_gravityNodePos = 1.0e-36_pReal
type :: p_vec type :: p_vec