subroutine boundaryflux_advection_diffusion_1d(this)
! this method uses an linear upwind solver for the
! advective flux and the bassi-rebay method for the
! diffusive fluxes
implicit none
class(advection_diffusion_1d),intent(inout) :: this
! Local
integer :: ndof
ndof = this%solution%nelem*this%solution%nvar*2
call boundaryflux_advection_diffusion_1d_gpu(this%solution%boundary_gpu, &
this%solution%extBoundary_gpu,this%solutionGradient%avgBoundary_gpu, &
this%flux%boundarynormal_gpu,this%u,this%nu,ndof)
endsubroutine boundaryflux_advection_diffusion_1d