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

better use rotation class functionality

This commit is contained in:
Martin Diehl 2019-09-19 23:32:19 -07:00
parent 70b864e745
commit 27732b1c0d
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/material.f90
View File

@ -693,6 +693,7 @@ subroutine material_parseTexture
integer, dimension(:), allocatable :: chunkPos integer, dimension(:), allocatable :: chunkPos
real(pReal), dimension(3,3) :: texture_transformation ! maps texture to microstructure coordinate system real(pReal), dimension(3,3) :: texture_transformation ! maps texture to microstructure coordinate system
real(pReal), dimension(3) :: Eulers ! Euler angles in degrees from file real(pReal), dimension(3) :: Eulers ! Euler angles in degrees from file
type(rotation) :: transformation
do t=1, size(config_texture) do t=1, size(config_texture)
if (config_texture(t)%countKeys('(gauss)') /= 1) call IO_error(147,ext_msg='count((gauss)) != 1') if (config_texture(t)%countKeys('(gauss)') /= 1) call IO_error(147,ext_msg='count((gauss)) != 1')
@ -740,7 +741,8 @@ subroutine material_parseTexture
end select end select
enddo enddo
if(dNeq(math_det33(texture_transformation),1.0_pReal)) call IO_error(157,t) if(dNeq(math_det33(texture_transformation),1.0_pReal)) call IO_error(157,t)
call texture_Eulers(t)%fromRotationMatrix(matmul(texture_Eulers(t)%asRotationMatrix(),texture_transformation)) call transformation%fromRotationMatrix(texture_transformation)
texture_Eulers(t) = texture_Eulers(t) * transformation
endif endif
enddo enddo