mg_comm_init Subroutine

public subroutine mg_comm_init(mg, comm)

Uses

  • proc~~mg_comm_init~3~~UsesGraph proc~mg_comm_init~3 mg_comm_init mpi mpi proc~mg_comm_init~3->mpi

Prolong from a parent to a child with index offset dix. This method could sometimes work better than the default prolongation, which does not take the variation in epsilon into account Initialize MPI if needed, and store MPI information

Arguments

Type IntentOptional Attributes Name
type(mg_t), intent(inout) :: mg
integer, intent(in), optional :: comm

MPI communicator (default: MPI_COMM_WORLD)

Calls

proc~~mg_comm_init~3~~CallsGraph proc~mg_comm_init~3 mg_comm_init mpi_comm_rank mpi_comm_rank proc~mg_comm_init~3->mpi_comm_rank mpi_comm_size mpi_comm_size proc~mg_comm_init~3->mpi_comm_size mpi_init mpi_init proc~mg_comm_init~3->mpi_init mpi_initialized mpi_initialized proc~mg_comm_init~3->mpi_initialized

Called by

proc~~mg_comm_init~3~~CalledByGraph proc~mg_comm_init~3 mg_comm_init proc~mg_setup_multigrid mg_setup_multigrid proc~mg_setup_multigrid->proc~mg_comm_init~3 proc~improve_initial_condition improve_initial_condition proc~improve_initial_condition->proc~mg_setup_multigrid proc~main main proc~main->proc~mg_setup_multigrid proc~main->proc~improve_initial_condition program~agile agile program~agile->proc~main

Source Code

  subroutine mg_comm_init(mg, comm)
    use mpi
    type(mg_t), intent(inout)     :: mg
    !> MPI communicator (default: MPI_COMM_WORLD)
    integer, intent(in), optional :: comm
    integer                       :: ierr
    logical                       :: initialized

    call mpi_initialized(initialized, ierr)
    if (.not. initialized) then
       call mpi_init(ierr)
    end if

    if (present(comm)) then
       mg%comm = comm
    else
       mg%comm = MPI_COMM_WORLD
    end if

    call mpi_comm_rank(mg%comm, mg%my_rank, ierr)
    call mpi_comm_size(mg%comm, mg%n_cpu, ierr)
  end subroutine mg_comm_init