ax_helm_template.f90

Go to the documentation of this file.

 
module ax_helm_template
  use ax_product, only : ax_t, ax_helm_allocate, register_ax_helm
  use ax_helm, only : ax_helm_t
  use coefs, only : coef_t
  use mesh, only : mesh_t
  use num_types, only : rp
  use space, only : space_t
  implicit none
  private
 
  type, extends(ax_helm_t) :: ax_helm_template_t
   contains
     procedure, nopass :: compute => ax_helm_template_compute
  end type ax_helm_template_t
 
  public :: ax_helm_template_register_types
 
contains
 
  subroutine ax_helm_template_compute(w, u, coef, msh, Xh)
    type(mesh_t), intent(in) :: msh
    type(space_t), intent(in) :: Xh
    type(coef_t), intent(in) :: coef
    real(kind=rp), intent(inout) :: w(xh%lx, xh%ly, xh%lz, msh%nelv)
    real(kind=rp), intent(in) :: u(xh%lx, xh%ly, xh%lz, msh%nelv)
 
    ! TODO: Replace the identity operation with the desired Ax product.
    w = u
  end subroutine ax_helm_template_compute
 
  subroutine ax_helm_template_register_types()
    procedure(ax_helm_allocate), pointer :: allocator
    allocator => ax_helm_template_allocate
    call register_ax_helm('ax_helm_template', allocator)
  end subroutine ax_helm_template_register_types
 
  subroutine ax_helm_template_allocate(obj)
    class(ax_t), allocatable, intent(inout) :: obj
    allocate(ax_helm_template_t :: obj)
  end subroutine ax_helm_template_allocate
end module ax_helm_template