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

volume calculation was wrong 

could result in negative volumes
This commit is contained in:
Martin Diehl 2019-10-16 20:56:31 +02:00
parent dbe15f88f2
commit 0ae0e23325
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/math.f90
View File

@ -1307,10 +1307,10 @@ real(pReal) pure function math_volTetrahedron(v1,v2,v3,v4)
real(pReal), dimension (3,3) :: m real(pReal), dimension (3,3) :: m
m(1:3,1) = v1-v2 m(1:3,1) = v1-v2
m(1:3,2) = v2-v3 m(1:3,2) = v1-v3
m(1:3,3) = v3-v4 m(1:3,3) = v1-v4
math_volTetrahedron = math_det33(m)/6.0_pReal math_volTetrahedron = abs(math_det33(m))/6.0_pReal
end function math_volTetrahedron end function math_volTetrahedron