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

debugged zoo of rotation operations and transformations 

all angles now in radians

introduced a rudimentary check in math_init to complain (IO_error) about broken transformations (e.g. quat --> R --> quat)
This commit is contained in:
Philip Eisenlohr 2010-05-06 14:07:21 +00:00
parent 4954228b3a
commit e4067f2d26
3 changed files with 380 additions and 252 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
code/IO.f90
View File

@ -1085,6 +1085,15 @@ endfunction
case (666) case (666)
msg = 'Memory leak detected' msg = 'Memory leak detected'
case (670)
msg = 'math_check: quat -> axisAngle -> quat failed'
case (671)
msg = 'math_check: quat -> R -> quat failed'
case (672)
msg = 'math_check: quat -> euler -> quat failed'
case (673)
msg = 'math_check: R -> euler -> R failed'
case (700) case (700)
msg = 'Singular matrix in stress iteration' msg = 'Singular matrix in stress iteration'
case default case default

2
code/lattice.f90
View File

@ -866,7 +866,7 @@ function lattice_initializeStructure(struct,CoverA)
lattice_tn(:,i,myStructure) = tn(:,i) lattice_tn(:,i,myStructure) = tn(:,i)
lattice_Stwin(:,:,i,myStructure) = math_tensorproduct(td(:,i),tn(:,i)) lattice_Stwin(:,:,i,myStructure) = math_tensorproduct(td(:,i),tn(:,i))
lattice_Stwin_v(:,i,myStructure) = math_Mandel33to6(math_symmetric3x3(lattice_Stwin(:,:,i,myStructure))) lattice_Stwin_v(:,i,myStructure) = math_Mandel33to6(math_symmetric3x3(lattice_Stwin(:,:,i,myStructure)))
lattice_Qtwin(:,:,i,myStructure) = math_RodrigToR(tn(:,i),180.0_pReal*inRad) lattice_Qtwin(:,:,i,myStructure) = math_AxisAngleToR(tn(:,i),180.0_pReal*inRad)
lattice_shearTwin(i,myStructure) = ts(i) lattice_shearTwin(i,myStructure) = ts(i)
enddo enddo
lattice_NslipSystem(1:lattice_maxNslipFamily,myStructure) = myNslipSystem ! number of slip systems in each family lattice_NslipSystem(1:lattice_maxNslipFamily,myStructure) = myNslipSystem ! number of slip systems in each family

567
code/math.f90
View File

File diff suppressed because it is too large Load Diff