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

[skip ci] indentation hick-up

This commit is contained in:
Martin Diehl 2019-09-23 22:13:52 -07:00
parent ffcb7a0a24
commit 116f43f401
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/mesh_marc.f90
View File

@ -888,11 +888,11 @@ subroutine buildCells(thisMesh,elem,connectivity_elem)
call results_closeJobFile call results_closeJobFile
#endif #endif
contains contains
!-------------------------------------------------------------------------------------------------- !------------------------------------------------------------------------------------------------
!> @brief count unique rows (same rows need to be stored consequtively) !> @brief count unique rows (same rows need to be stored consecutively)
!-------------------------------------------------------------------------------------------------- !------------------------------------------------------------------------------------------------
pure function uniqueRows(A) result(u) pure function uniqueRows(A) result(u)
integer, dimension(:,:), intent(in) :: A !< array, rows need to be sorted integer, dimension(:,:), intent(in) :: A !< array, rows need to be sorted
@ -1223,7 +1223,7 @@ subroutine mesh_build_ipAreas
forall(n = 1:FE_NcellnodesPerCellface(c)) & forall(n = 1:FE_NcellnodesPerCellface(c)) &
nodePos(1:3,n) = mesh_cellnode(1:3,mesh_cell(theMesh%elem%cellface(n,f),i,e)) nodePos(1:3,n) = mesh_cellnode(1:3,mesh_cell(theMesh%elem%cellface(n,f),i,e))
normal = math_cross(nodePos(1:3,2) - nodePos(1:3,1), & normal = math_cross(nodePos(1:3,2) - nodePos(1:3,1), &
nodePos(1:3,3) - nodePos(1:3,1)) nodePos(1:3,3) - nodePos(1:3,1))
mesh_ipArea(f,i,e) = norm2(normal) mesh_ipArea(f,i,e) = norm2(normal)
mesh_ipAreaNormal(1:3,f,i,e) = normal / norm2(normal) ! ensure unit length of area normal mesh_ipAreaNormal(1:3,f,i,e) = normal / norm2(normal) ! ensure unit length of area normal
enddo enddo
@ -1242,7 +1242,7 @@ subroutine mesh_build_ipAreas
forall(n = 1:FE_NcellnodesPerCellface(c)) & forall(n = 1:FE_NcellnodesPerCellface(c)) &
normals(1:3,n) = 0.5_pReal & normals(1:3,n) = 0.5_pReal &
* math_cross(nodePos(1:3,1+mod(n ,m)) - nodePos(1:3,n), & * math_cross(nodePos(1:3,1+mod(n ,m)) - nodePos(1:3,n), &
nodePos(1:3,1+mod(n+1,m)) - nodePos(1:3,n)) nodePos(1:3,1+mod(n+1,m)) - nodePos(1:3,n))
normal = 0.5_pReal * sum(normals,2) normal = 0.5_pReal * sum(normals,2)
mesh_ipArea(f,i,e) = norm2(normal) mesh_ipArea(f,i,e) = norm2(normal)
mesh_ipAreaNormal(1:3,f,i,e) = normal / norm2(normal) mesh_ipAreaNormal(1:3,f,i,e) = normal / norm2(normal)