lagrangian_set_CreateMPIType Subroutine

private impure subroutine lagrangian_set_CreateMPIType(this)

Creates an MPI data type for parallel communication of the Lagrangian objects.

Type Bound

lagrangian_set

Arguments

Type IntentOptional Attributes Name
class(lagrangian_set), intent(inout) :: this

A set of Lagrangian objects

Calls

proc~~lagrangian_set_creatempitype~~CallsGraph proc~lagrangian_set_creatempitype lagrangian_set%lagrangian_set_CreateMPIType SetMPIDataTypeParams SetMPIDataTypeParams proc~lagrangian_set_creatempitype->SetMPIDataTypeParams mpi_type_commit mpi_type_commit proc~lagrangian_set_creatempitype->mpi_type_commit mpi_type_create_resized mpi_type_create_resized proc~lagrangian_set_creatempitype->mpi_type_create_resized mpi_type_create_struct mpi_type_create_struct proc~lagrangian_set_creatempitype->mpi_type_create_struct mpi_type_get_extent mpi_type_get_extent proc~lagrangian_set_creatempitype->mpi_type_get_extent mpi_type_size mpi_type_size proc~lagrangian_set_creatempitype->mpi_type_size

Called by

proc~~lagrangian_set_creatempitype~~CalledByGraph proc~lagrangian_set_creatempitype lagrangian_set%lagrangian_set_CreateMPIType proc~lagrangian_set_init lagrangian_set%lagrangian_set_Init proc~lagrangian_set_init->proc~lagrangian_set_creatempitype proc~marker_set_init marker_set%marker_set_Init proc~marker_set_init->proc~lagrangian_set_creatempitype proc~particle_set_changeparttype particle_set%particle_set_ChangePartType proc~particle_set_changeparttype->proc~lagrangian_set_creatempitype proc~respart_set_init ResPart_set%ResPart_set_Init proc~respart_set_init->proc~lagrangian_set_creatempitype none~initialize~26 marker_set%Initialize proc~respart_set_init->none~initialize~26 none~initialize~25 ResPart_set%Initialize none~initialize~25->proc~respart_set_init none~initialize~26->proc~respart_set_init none~initialize~27 solid_obj%Initialize none~initialize~27->proc~respart_set_init proc~particle_set_prepare particle_set%particle_set_Prepare proc~particle_set_prepare->proc~particle_set_changeparttype proc~cdifs_obj_preparesolver cdifs_obj_PrepareSolver proc~cdifs_obj_preparesolver->none~initialize~25 proc~grans_obj_preparesolver grans_obj_PrepareSolver proc~grans_obj_preparesolver->none~initialize~26 proc~solid_set_init solid_set%solid_set_Init proc~solid_set_init->none~initialize~27 interface~cdifs_obj_preparesolver cdifs_obj%cdifs_obj_PrepareSolver interface~cdifs_obj_preparesolver->proc~cdifs_obj_preparesolver interface~grans_obj_preparesolver grans_obj%grans_obj_PrepareSolver interface~grans_obj_preparesolver->proc~grans_obj_preparesolver

Source Code

    impure subroutine lagrangian_set_CreateMPIType(this)
      !> Creates an MPI data type for parallel communication
      ! of the Lagrangian objects.
      implicit none
      class(lagrangian_set), intent(inout) :: this                             !! A set of Lagrangian objects
      ! Work variables
      integer,                        allocatable :: lengths(:)
      type(MPI_DATATYPE),             allocatable :: types(:)
      integer(kind=MPI_ADDRESS_KIND), allocatable :: displacement(:)
      integer(kind=MPI_ADDRESS_KIND) :: extent,lb
      type(MPI_DATATYPE) :: type_tmp
      integer :: ierr

      associate(mpi=>this%parallel)
        ! Compute types, lengths and displacements
        call this%SetMPIDataTypeParams(types,lengths,displacement)

        ! Create MPI type
        call MPI_TYPE_CREATE_STRUCT(size(lengths),lengths,displacement,types,type_tmp,ierr)
        if (ierr.ne.0) error stop "Problem with MPI_TYPE"

        ! Get actual lower bound and extent of this type
        call MPI_TYPE_GET_EXTENT(type_tmp,lb,extent,ierr)
        if (ierr.ne.0) error stop "Problem with MPI_TYPE"

        ! Adjust upper and lower bounds
        call MPI_TYPE_CREATE_RESIZED(type_tmp,lb,extent,this%MPI_TYPE,ierr)
        call MPI_TYPE_COMMIT(this%MPI_TYPE,ierr)
        if (ierr.ne.0) error stop "Problem with MPI_TYPE"

        ! Get the size of this type
        call MPI_TYPE_SIZE(this%MPI_TYPE,this%MPI_SIZE,ierr)
        if (ierr.ne.0) error stop "Problem with MPI_TYPE"
      end associate

      return
    end subroutine lagrangian_set_CreateMPIType