anchor.diff:              &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_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")
  
positive.code: 
          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
  
negative.code:             integer,                 intent(in) :: idxlo
              integer,                 intent(in) :: order
              real(8),                 intent(in) :: enu
          end subroutine
  
          subroutine sirius_get_aw_surface_derivative(ia, l, io, dm, deriv)&
              &bind(C, name="sirius_get_aw_surface_derivative")
              integer,                 intent(in)  :: ia
              integer,                 intent(in)  :: l
              integer,                 intent(in)  :: io
              integer,                 intent(in)  :: dm
              real(8),                 intent(out) :: deriv
          end subroutine
  
          subroutine sirius_set_aw_surface_derivative(ia, l, io, dm, deriv)&
  