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

following standard: Output first

This commit is contained in:
Martin Diehl 2019-09-20 18:28:46 -07:00
parent 53283d5c01
commit 97127f9034
2 changed files with 2 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/homogenization_mech_RGC.f90
View File

@ -816,7 +816,7 @@ module procedure mech_RGC_updateState
integer :: i,j,iBase integer :: i,j,iBase
logical :: error logical :: error
call math_invert33(matmul(transpose(avgF),avgF),invC,detF,error) call math_invert33(invC,detF,error,matmul(transpose(avgF),avgF))
surfaceCorrection = 0.0_pReal surfaceCorrection = 0.0_pReal
do iBase = 1,3 do iBase = 1,3

3
src/math.f90
View File

@ -484,9 +484,8 @@ end function math_inv33
!> @brief Cramer inversion of 33 matrix (subroutine) !> @brief Cramer inversion of 33 matrix (subroutine)
!> @details Direct Cramer inversion of matrix A. Also returns determinant !> @details Direct Cramer inversion of matrix A. Also returns determinant
! Returns an error if not possible, i.e. if determinant is close to zero ! Returns an error if not possible, i.e. if determinant is close to zero
! ToDo: Output arguments should be first
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
pure subroutine math_invert33(A, InvA, DetA, error) pure subroutine math_invert33(InvA, DetA, error, A)
logical, intent(out) :: error logical, intent(out) :: error
real(pReal),dimension(3,3),intent(in) :: A real(pReal),dimension(3,3),intent(in) :: A