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

safe use of optional argument

This commit is contained in:
Martin Diehl 2019-09-19 14:20:02 -07:00
parent f6c616c84b
commit b41d6182a5
1 changed files with 8 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
src/math.f90
View File

@ -425,12 +425,19 @@ pure function math_exp33(A,n)
real(pReal), dimension(3,3), intent(in) :: A
real(pReal), dimension(3,3) :: B, math_exp33
real(pReal) :: invFac
integer :: order
B = math_I3 ! init
invFac = 1.0_pReal ! 0!
math_exp33 = B ! A^0 = eye2
do i = 1, merge(n,5,present(n))
if (present(n)) then
order = n
else
order = 5
endif
do i = 1, order
invFac = invFac/real(i,pReal) ! invfac = 1/i!
B = matmul(B,A)
math_exp33 = math_exp33 + invFac*B ! exp = SUM (A^i)/i!