anchor.diff:  
           subroutine sirius_set_lo_radial_function(ia, idxlo, f)&
               &bind(C, name="sirius_set_lo_radial_function")
  +            integer,                 intent(in)  :: ia
  +            integer,                 intent(in)  :: idxlo
  +            real(8),                 intent(in)  :: f
  +        end subroutine
  +
  +        subroutine sirius_set_lo_radial_function_derivative(ia, idxlo, f)&
  +            &bind(C, name="sirius_set_lo_radial_function_derivative")
               integer,                 intent(in)  :: ia
               integer,                 intent(in)  :: idxlo
               real(8),                 intent(in)  :: f
  
positive.code:             integer,                 intent(in) :: idxlo
              integer,                 intent(in) :: order
              real(8),                 intent(in) :: enu
          end subroutine
  
          subroutine sirius_set_aw_surface_derivative(ia, l, io, dm, deriv)&
  
negative.code:             &bind(C, name="sirius_generate_radial_integrals")
          end subroutine
  
          !subroutine sirius_get_aw_deriv_radial_function(ia, l, io, dfdr)&
          !    &bind(C, name="sirius_get_aw_deriv_radial_function")
          !    integer,                 intent(in)  :: ia
          !    integer,                 intent(in)  :: l
          !    integer,                 intent(in)  :: io
          !    real(8),                 intent(out) :: dfdr
          !end subroutine
  
          !subroutine sirius_get_lo_deriv_radial_function(ia, idxlo, dfdr)&
          !    &bind(C, name="sirius_get_lo_deriv_radial_function")
          !    integer,                 intent(in)  :: ia
          !    integer,                 intent(in)  :: idxlo
          !    real(8),                 intent(out) :: dfdr
          !end subroutine
  
          subroutine sirius_get_aw_radial_function(ia, l, io, f)&
              &bind(C, name="sirius_get_aw_radial_function")
  