s_spline_matrix_banded__set_element Subroutine

  subroutine s_spline_matrix_banded__set_element( self, i, j, a_ij )
    class(sll_t_spline_matrix_banded), intent(inout) :: self
    integer                          , intent(in   ) :: i
    integer                          , intent(in   ) :: j
    real(wp)                         , intent(in   ) :: a_ij

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

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

  end subroutine s_spline_matrix_banded__set_element