fluxmethod_advection_diffusion_1d Subroutine

  subroutine fluxmethod_advection_diffusion_1d(this)
    implicit none
    class(advection_diffusion_1d),intent(inout) :: this
    ! Local
    integer :: ndof

    ndof = this%solution%nelem*this%solution%nvar*(this%solution%interp%N+1)

    call fluxmethod_advection_diffusion_1d_gpu(this%solution%interior_gpu, &
                                               this%solutiongradient%interior_gpu,this%flux%interior_gpu, &
                                               this%u,this%nu,ndof)

  endsubroutine fluxmethod_advection_diffusion_1d