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

connectivity is not needed

This commit is contained in:
Martin Diehl 2019-06-07 10:34:02 +02:00
parent 33ed1fa06c
commit ca97c3db18
1 changed files with 4 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
src/mesh_grid.f90
View File

@ -166,7 +166,7 @@ subroutine mesh_init(ip,el)
allocate(FEsolving_execIP(2_pInt,theMesh%nElems), source=1_pInt) ! parallel loop bounds set to comprise from first IP... allocate(FEsolving_execIP(2_pInt,theMesh%nElems), source=1_pInt) ! parallel loop bounds set to comprise from first IP...
forall (j = 1_pInt:theMesh%nElems) FEsolving_execIP(2,j) = theMesh%elem%nIPs ! ...up to own IP count for each element forall (j = 1_pInt:theMesh%nElems) FEsolving_execIP(2,j) = theMesh%elem%nIPs ! ...up to own IP count for each element
call discretization_init(mesh_element(3,:),mesh_element(4,:),& call discretization_init(mesh_element(1,:),mesh_element(2,:),&
reshape(mesh_ipCoordinates,[3,grid(1)*grid(2)*grid3]),& reshape(mesh_ipCoordinates,[3,grid(1)*grid(2)*grid3]),&
mesh_node0) mesh_node0)
@ -370,23 +370,12 @@ subroutine mesh_spectral_build_elements()
e, & e, &
elemOffset elemOffset
allocate(mesh_element (4_pInt+8_pInt,theMesh%nElems), source = 0_pInt) allocate(mesh_element (2,theMesh%nElems), source = 0_pInt)
elemOffset = product(grid(1:2))*grid3Offset elemOffset = product(grid(1:2))*grid3Offset
do e=1, theMesh%nElems do e=1, theMesh%nElems
mesh_element( 1,e) = -1_pInt ! DEPRECATED mesh_element( 1,e) = mesh_homogenizationAt(e)
mesh_element( 2,e) = -1_pInt ! DEPRECATED mesh_element( 2,e) = microGlobal(e+elemOffset)
mesh_element( 3,e) = mesh_homogenizationAt(e)
mesh_element( 4,e) = microGlobal(e+elemOffset) ! microstructure
mesh_element( 5,e) = e + (e-1_pInt)/grid(1) + &
((e-1_pInt)/(grid(1)*grid(2)))*(grid(1)+1_pInt) ! base node
mesh_element( 6,e) = mesh_element(5,e) + 1_pInt
mesh_element( 7,e) = mesh_element(5,e) + grid(1) + 2_pInt
mesh_element( 8,e) = mesh_element(5,e) + grid(1) + 1_pInt
mesh_element( 9,e) = mesh_element(5,e) +(grid(1) + 1_pInt) * (grid(2) + 1_pInt) ! second floor base node
mesh_element(10,e) = mesh_element(9,e) + 1_pInt
mesh_element(11,e) = mesh_element(9,e) + grid(1) + 2_pInt
mesh_element(12,e) = mesh_element(9,e) + grid(1) + 1_pInt
enddo enddo
end subroutine mesh_spectral_build_elements end subroutine mesh_spectral_build_elements