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

polishing

This commit is contained in:
Martin Diehl 2020-01-02 20:49:02 +01:00
parent 1ad8464821
commit f633c67fe4
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
src/math.f90
View File

@ -411,9 +411,9 @@ pure function math_exp33(A,n)
endif endif
do i = 1, order do i = 1, order
invFac = invFac/real(i,pReal) ! invfac = 1/i! invFac = invFac/real(i,pReal) ! invfac = 1/(i!)
B = matmul(B,A) B = matmul(B,A)
math_exp33 = math_exp33 + invFac*B ! exp = SUM (A^i)/i! math_exp33 = math_exp33 + invFac*B ! exp = SUM (A^i)/(i!)
enddo enddo
end function math_exp33 end function math_exp33
@ -489,8 +489,8 @@ function math_invSym3333(A)
real(pReal), dimension(6*(64+2)) :: work real(pReal), dimension(6*(64+2)) :: work
logical :: error logical :: error
external :: & external :: &
dgetrf, & dgetrf, &
dgetri dgetri
temp66 = math_sym3333to66(A) temp66 = math_sym3333to66(A)
call dgetrf(6,6,temp66,6,ipiv6,ierr) call dgetrf(6,6,temp66,6,ipiv6,ierr)
@ -519,8 +519,8 @@ subroutine math_invert(InvA, error, A)
real(pReal), dimension(size(A,1)*(64+2)) :: work real(pReal), dimension(size(A,1)*(64+2)) :: work
integer :: ierr integer :: ierr
external :: & external :: &
dgetrf, & dgetrf, &
dgetri dgetri
invA = A invA = A
call dgetrf(size(A,1),size(A,1),invA,size(A,1),ipiv,ierr) call dgetrf(size(A,1),size(A,1),invA,size(A,1),ipiv,ierr)