From f633c67fe4647d12c26edb0241aec6ddd6692def Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Thu, 2 Jan 2020 20:49:02 +0100 Subject: [PATCH] polishing --- src/math.f90 | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/src/math.f90 b/src/math.f90 index 0b06c9186..d6d87d30a 100644 --- a/src/math.f90 +++ b/src/math.f90 @@ -411,9 +411,9 @@ pure function math_exp33(A,n) endif do i = 1, order - invFac = invFac/real(i,pReal) ! invfac = 1/i! + invFac = invFac/real(i,pReal) ! invfac = 1/(i!) 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 end function math_exp33 @@ -489,8 +489,8 @@ function math_invSym3333(A) real(pReal), dimension(6*(64+2)) :: work logical :: error external :: & - dgetrf, & - dgetri + dgetrf, & + dgetri temp66 = math_sym3333to66(A) 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 integer :: ierr external :: & - dgetrf, & - dgetri + dgetrf, & + dgetri invA = A call dgetrf(size(A,1),size(A,1),invA,size(A,1),ipiv,ierr)