sll_m_bsplines_base.F90

Go to the documentation of this file.

 
 
module sll_m_bsplines_base
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#include "sll_assert.h"
 
   use sll_m_working_precision, only: f64
 
   implicit none
 
   public :: &
      sll_c_bsplines
 
   private
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
   integer, parameter :: wp = f64
 
   type, abstract :: sll_c_bsplines
 
      integer :: degree
      logical :: periodic
      logical :: uniform
      logical :: radial
      integer :: ncells
      integer :: nbasis
 
      real(wp) :: xmin
      real(wp) :: xmax
 
      real(wp), allocatable :: knots(:) ! Only used by non-uniform B-splines
 
   contains
      procedure(i_fun_find_cell), deferred :: find_cell
      procedure(i_sub_eval_basis), deferred :: eval_basis
      procedure(i_sub_eval_deriv), deferred :: eval_deriv
      procedure(i_sub_eval_basis_and_n_derivs), deferred :: eval_basis_and_n_derivs
      procedure(i_sub_free), deferred :: free
 
   end type sll_c_bsplines
 
   abstract interface
 
      !---------------------------------------------------------------------------
      !---------------------------------------------------------------------------
      sll_pure function i_fun_find_cell(self, x) result(icell)
         import sll_c_bsplines, wp
         class(sll_c_bsplines), intent(in) :: self
         real(wp), intent(in) :: x
         integer :: icell
      end function i_fun_find_cell
 
      !---------------------------------------------------------------------------
      !---------------------------------------------------------------------------
      sll_pure subroutine i_sub_eval_basis(self, x, values, jmin)
         import sll_c_bsplines, wp
         class(sll_c_bsplines), intent(in) :: self
         real(wp), intent(in) :: x
         real(wp), intent(out) :: values(:)
         integer, intent(out) :: jmin
      end subroutine i_sub_eval_basis
 
      !---------------------------------------------------------------------------
      !---------------------------------------------------------------------------
      sll_pure subroutine i_sub_eval_deriv(self, x, derivs, jmin)
         import sll_c_bsplines, wp
         class(sll_c_bsplines), intent(in) :: self
         real(wp), intent(in) :: x
         real(wp), intent(out) :: derivs(:)
         integer, intent(out) :: jmin
      end subroutine i_sub_eval_deriv
 
      !---------------------------------------------------------------------------
      !---------------------------------------------------------------------------
      sll_pure subroutine i_sub_eval_basis_and_n_derivs(self, x, n, derivs, jmin)
         import sll_c_bsplines, wp
         class(sll_c_bsplines), intent(in) :: self
         real(wp), intent(in) :: x
         integer, intent(in) :: n
         real(wp), intent(out) :: derivs(:, :)
         integer, intent(out) :: jmin
      end subroutine i_sub_eval_basis_and_n_derivs
 
      !---------------------------------------------------------------------------
      !---------------------------------------------------------------------------
      subroutine i_sub_free(self)
         import sll_c_bsplines
         class(sll_c_bsplines), intent(inout) :: self
      end subroutine i_sub_free
 
   end interface
 
end module sll_m_bsplines_base