MESH_PARTITIONING.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	mesh_partitioning (mpi_file, nelem, nnode, nparts, conn, w_yn
	Makes mesh partitioning using METIS.

Function/Subroutine Documentation

mesh_partitioning()

subroutine mesh_partitioning	(	character*70	mpi_file,
		integer*4, intent(inout)	nelem,
		integer*4, intent(inout)	nnode,
		integer*4, intent(inout)	nparts,
		integer*4, dimension(nelem,9), intent(in)	conn,
		integer*4	w_yn
	)

Makes mesh partitioning using METIS.

Date

September, 2013

Version

1.0

Warning

Mesh partitioning is made using metis-4.0.3 library. At the moment we have problems in using metis-5.0.2 library.

Todo:

Switch to other mesh partitioner or to ParMetis library.

Parameters

[in]	mpi_file	folder name given in SPEED.input (MPIFILE keyword)
[in]	nelem	number of elements
[in]	nnode	number of mesh nodes
[in]	nparts	number of processors (number of parts of the partition)
[in]	conn	connectivity matrix: i = element number conn(i,1) = Matrial /Block label conn(i,2),...,conn(i,9) = nodes of the hex
[in]	w	control parameter 1 - different weights for the elements 0 - same weight for all elements (for metis-5.0.2 version)
[out]	elemdomain.mpi	file containing elements partitioning

Definition at line 37 of file MESH_PARTITIONING.f90.

 
 
      implicit none
      
 
      character*70 :: mpi_file,u_name
      integer*4 :: ie,i,ic,u_mpi, w_yn, trash
      integer*4 :: ncommon, objval, edgecut
 
      integer*4, intent(inout) :: nelem, nnode, nparts
      integer*4, dimension(nelem,9), intent(in) :: conn
 
      integer*4, dimension(nnode) :: npart
      integer*4, dimension(nelem) :: epart
      integer*4, dimension(nelem) :: vwgt
 
      integer*4, pointer :: vsize(:) => null() 
 
      integer*4, dimension(nelem+1) :: eptr
      integer*4, dimension(nelem*8) :: eind
            
      integer*4 :: options(40) ! see METIS_NOPTIONS in metis.h
      integer*4, parameter :: METIS_OPTION_PTYPE = 1, metis_option_objtype = 2, metis_option_ctype = 3, &
                              metis_option_iptype = 4, metis_option_rtype = 5, metis_option_dbglvl = 6, &
                              metis_option_niter = 7, metis_option_ncuts = 8, metis_option_seed = 9, &
                              metis_option_no2hop = 10, &
                              metis_option_minconn = 11, metis_option_contig = 12, metis_option_compress = 13, &
                              metis_option_ccorder = 14, metis_option_pfactor = 15, metis_option_nseps = 16, &
                              metis_option_ufactor = 17, metis_option_numbering = 18
 
      real*8, pointer :: tpwgts(:) => null() 
 
      eptr(1) = 1
      ic = 1
      do ie = 1, nelem
         eptr(ie+1) = eptr(ie) + 8 
         do i = 2, 9
            eind(ic) = conn(ie,i)
            ic = ic +1
         enddo
      enddo
 
     eptr = eptr -1
     eind = eind - 1
 
     !!!----------   METIS-4.0.3  ------------------------------------!!!
      eind = eind +1
      call metis_partmeshdual(nelem,nnode,eind,3,1,nparts,edgecut,&
                              epart,npart)
      epart = epart - 1
 
     !!!----------   METIS-5.1.0  ------------------------------------!!!
     !ncommon = 4 
     !eptr = eptr + 1
     !eind = eind + 1
     !if(w_yn .eq. 0) then
     !   vwgt = 1
     !else
     !   u_name = 'element.wgt'
     !   u_mpi = 40                                 
     !   open(u_mpi,file=u_name)        
     !   read(u_mpi,*) nelem
     !   
     !   do i = 1, nelem
     !          read(u_mpi,*) trash, vwgt(i)
     !   enddo
     !   close(u_mpi)
     !endif     
 
     !! Setting the Fortran numbering scheme
     !call metis_setdefaultoptions(options)  !INPUT
     !options(metis_option_numbering) = 1
                              
     !call metis_partmeshdual(nelem, nnode, eptr, eind, vwgt, vsize, &
     !                         ncommon, nparts, tpwgts, options, objval, &
     !                         epart, npart)                        
     ! epart = epart - 1 
 
     !!!--------------------------------------------------------------------!!!
  
 
      ! if(w_yn .eq. 0) then
         if(len_trim(mpi_file) .ne. 70) then                                                                                  
            u_name = mpi_file(1:len_trim(mpi_file)) // '/elemdomain.mpi'
         else 
            u_name = 'elemdomain.mpi'     
         endif      
   
         u_mpi = 400                                 
         open(u_mpi,file=u_name)        
         write(u_mpi,*) nelem
         do i = 1, nelem
               write(u_mpi,*) i, epart(i)
         enddo
         close(u_mpi)
 
!       else
!
!         if(len_trim(mpi_file) .ne. 70) then                                                                                  
!            u_name = mpi_file(1:len_trim(mpi_file)) // '/final.part'
!         else 
!            u_name = 'final.part'     
!         endif
!         
!         u_mpi = 400                                 
!         open(u_mpi,file=u_name)        
!         write(u_mpi,*) nelem
!         do i = 1, nelem
!               write(u_mpi,*) i, epart(i)
!         enddo
!         close(u_mpi)
 
      !  endif       
      
      return