Is there a way to partition a mesh consisting of 20 noded hexahedral elements for parallel processing? I used METIS for partitioning mesh with 8 noded hexahedron elmements, which works fine but i don't know how to extend this to 20 noded element mesh.
My Fortran code is as follows, which works fine for 4 node quad and 8 node hex meshes but when i used 20 noded element it gives error
Segmentation fault
Code
subroutine split_domain(cells,nodes,topol,nprocs, epart)
implicit none
!======================================================== integer,intent(in)::cells,nodes,topol(cells,20),nprocs
integer,intent(out)::epart(cells)
integer::ne, nn
integer::eptr(cells+1),eind(20*cells),npart(nodes)
integer::ncommon,nparts,objval
integer, pointer ::vwgt,vsize,tpwgts,options
integer::j
!========================================================
ne=cells; nn=nodes; eptr(1)=1
do j=1,ne
eind(20*j-19:20*j)=topol(j,:)
eptr(j+1)=eptr(j)+20
enddo
ncommon=1; nparts=nprocs
call METIS_PartMeshDual (ne,nn,eptr,eind,vwgt,vsize,ncommon,nparts,tpwgts,options,objval,epart,npart)
end subroutine split_domain
