s_spline_matrix_banded__get_element Function

  function s_spline_matrix_banded__get_element( self, i, j ) result( a_ij )
    class(sll_t_spline_matrix_banded), intent(inout) :: self
    integer                          , intent(in   ) :: i
    integer                          , intent(in   ) :: j
    real(wp)                                         :: a_ij

    SLL_ASSERT( max( 1, j-self%ku ) <= i )
    SLL_ASSERT( i <= min( self%n, j+self%kl ) )

    a_ij = self%q(self%kl+self%ku+1+i-j,j)

  end function s_spline_matrix_banded__get_element