Init_RefinementIndicator2D Subroutine

  1. SELF_RefinementIndicator_2D.f90
  2. SELF_RefinementIndicator_2D
  3. Init_RefinementIndicator2D

public subroutine Init_RefinementIndicator2D(this, interp, nElem, refineThreshold, coarsenThreshold)

Arguments

TypeIntentOptionalAttributesName
class(RefinementIndicator2D), intent(out) :: this
type(Lagrange), intent(in), target:: interp
integer, intent(in) :: nElem
real(kind=prec), intent(in) :: refineThreshold
real(kind=prec), intent(in) :: coarsenThreshold

Calls

proc~~init_refinementindicator2d~~CallsGraph proc~init_refinementindicator2d Init_RefinementIndicator2D proc~init_refinementindicator2d_t Init_RefinementIndicator2D_t proc~init_refinementindicator2d->proc~init_refinementindicator2d_t interface~hipmalloc hipMalloc proc~init_refinementindicator2d->interface~hipmalloc proc~gpucheck gpuCheck proc~init_refinementindicator2d->proc~gpucheck proc~buildmodaltransform BuildModalTransform proc~init_refinementindicator2d_t->proc~buildmodaltransform proc~normalizedlegendre NormalizedLegendre proc~buildmodaltransform->proc~normalizedlegendre proc~invertmatrix InvertMatrix proc~buildmodaltransform->proc~invertmatrix

Contents

Source Code

Init_RefinementIndicator2D

Source Code

  subroutine Init_RefinementIndicator2D(this,interp,nElem,refineThreshold,coarsenThreshold)
    implicit none
    class(RefinementIndicator2D),intent(out) :: this
    type(Lagrange),intent(in),target :: interp
    integer,intent(in) :: nElem
    real(prec),intent(in) :: refineThreshold
    real(prec),intent(in) :: coarsenThreshold

    ! Reuse the portable host initialization (allocation + modal transform build).
    call Init_RefinementIndicator2D_t(this,interp,nElem,refineThreshold,coarsenThreshold)

    ! The device kernel uses per-thread scratch sized to (N+1) <= 16 (AMR2D_MAXNP in
    ! SELF_Refinement.cpp). Fail early rather than overrun that bound on the device.
    if(interp%N+1 > 16) then
      print*,__FILE__,':',__LINE__, &
        ' : Error : GPU refinement indicator supports degree N <= 15 '// &
        '(raise AMR2D_MAXNP in SELF_Refinement.cpp to extend).'
      stop 1
    endif

    call gpuCheck(hipMalloc(this%Pmodal_gpu,sizeof(this%Pmodal)))
    call gpuCheck(hipMalloc(this%indicator_gpu,sizeof(this%indicator)))
    call gpuCheck(hipMalloc(this%flag_gpu,sizeof(this%flag)))

    call this%UpdateDevice()

  endsubroutine Init_RefinementIndicator2D