anchor.diff:              &bind(C, name="sirius_write_json_output")
           end subroutine
   
  +        subroutine sirius_ground_state_print_info()&
  +            &bind(C, name="sirius_ground_state_print_info")
  +        end subroutine
  +
           subroutine sirius_density_initialize_aux(rhoit, rhomt, magit, magmt)&
               &bind(C, name="sirius_density_initialize")
               use, intrinsic :: ISO_C_BINDING
  
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)&
  