MappedDGDivergence_MappedVector2D Subroutine

  1. SELF_MappedVector_2D.f90
  2. SELF_MappedVector_2D
  3. MappedDGDivergence_MappedVector2D

public subroutine MappedDGDivergence_MappedVector2D(this, df)

Arguments

TypeIntentOptionalAttributesName
class(MappedVector2D), intent(inout) :: this
type(c_ptr), intent(out) :: df

Calls

proc~~mappeddgdivergence_mappedvector2d~~CallsGraph proc~mappeddgdivergence_mappedvector2d MappedDGDivergence_MappedVector2D interface~contravariantprojection_2d_gpu ContravariantProjection_2D_gpu proc~mappeddgdivergence_mappedvector2d->interface~contravariantprojection_2d_gpu interface~divergence_2d_gpu Divergence_2D_gpu proc~mappeddgdivergence_mappedvector2d->interface~divergence_2d_gpu interface~dg_boundarycontribution_2d_gpu DG_BoundaryContribution_2D_gpu proc~mappeddgdivergence_mappedvector2d->interface~dg_boundarycontribution_2d_gpu interface~jacobianweight_2d_gpu JacobianWeight_2D_gpu proc~mappeddgdivergence_mappedvector2d->interface~jacobianweight_2d_gpu

Contents

Source Code

MappedDGDivergence_MappedVector2D

Source Code

  subroutine MappedDGDivergence_MappedVector2D(this,df)
    implicit none
    class(MappedVector2D),intent(inout) :: this
    type(c_ptr),intent(out) :: df

    ! Contravariant projection
    call ContravariantProjection_2D_gpu(this%interior_gpu, &
                                        this%geometry%dsdx%interior_gpu,this%interp%N,this%nvar,this%nelem)

    call Divergence_2D_gpu(this%interior_gpu,df,this%interp%dgMatrix_gpu, &
                           this%interp%N,this%nvar,this%nelem)

    ! Boundary terms --> TO DO : problem here when nvar > 1
    call DG_BoundaryContribution_2D_gpu(this%interp%bmatrix_gpu,this%interp%qweights_gpu, &
                                        this%boundarynormal_gpu,df,this%interp%N,this%nvar,this%nelem)

    call JacobianWeight_2D_gpu(df,this%geometry%J%interior_gpu,this%interp%N,this%nVar,this%nelem)

  endsubroutine MappedDGDivergence_MappedVector2D