setboundarycondition_LinearEuler2D Subroutine

public subroutine setboundarycondition_LinearEuler2D(this)

Boundary conditions are set to periodic boundary conditions

Arguments

TypeIntentOptionalAttributesName
class(LinearEuler2D), intent(inout) :: this

Calls

proc~~setboundarycondition_lineareuler2d~~CallsGraph proc~setboundarycondition_lineareuler2d setboundarycondition_LinearEuler2D hipmemcpy hipmemcpy proc~setboundarycondition_lineareuler2d->hipmemcpy gpucheck gpucheck proc~setboundarycondition_lineareuler2d->gpucheck interface~setboundarycondition_lineareuler2d_gpu setboundarycondition_LinearEuler2D_gpu proc~setboundarycondition_lineareuler2d->interface~setboundarycondition_lineareuler2d_gpu

Contents


Source Code

  subroutine setboundarycondition_LinearEuler2D(this)
    !! Boundary conditions are set to periodic boundary conditions
    implicit none
    class(LinearEuler2D),intent(inout) :: this
    ! local
    integer :: i,iEl,j,e2,bcid
    real(prec) :: x(1:2)

    if(this%prescribed_bcs_enabled) then
      call gpuCheck(hipMemcpy(c_loc(this%solution%extboundary), &
                              this%solution%extboundary_gpu,sizeof(this%solution%extboundary), &
                              hipMemcpyDeviceToHost))

      ! Prescribed boundaries are still done on the GPU
      do iEl = 1,this%solution%nElem ! Loop over all elements
        do j = 1,4 ! Loop over all sides

          bcid = this%mesh%sideInfo(5,j,iEl) ! Boundary Condition ID
          e2 = this%mesh%sideInfo(3,j,iEl) ! Neighboring Element ID

          if(e2 == 0) then
            if(bcid == SELF_BC_PRESCRIBED) then

              do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
                x = this%geometry%x%boundary(i,j,iEl,1,1:2)

                this%solution%extBoundary(i,j,iEl,1:this%nvar) = &
                  this%hbc2d_Prescribed(x,this%t)
              enddo

            endif
          endif

        enddo
      enddo

      call gpuCheck(hipMemcpy(this%solution%extBoundary_gpu, &
                              c_loc(this%solution%extBoundary), &
                              sizeof(this%solution%extBoundary), &
                              hipMemcpyHostToDevice))
    endif
    call setboundarycondition_LinearEuler2D_gpu(this%solution%extboundary_gpu, &
                                                this%solution%boundary_gpu, &
                                                this%mesh%sideInfo_gpu, &
                                                this%geometry%nhat%boundary_gpu, &
                                                this%solution%interp%N, &
                                                this%solution%nelem,this%solution%nvar)

  endsubroutine setboundarycondition_LinearEuler2D