anchor.diff:  real(C_DOUBLE), intent(in) :: J0
   end subroutine
   end interface
  +
   call sirius_set_atom_type_hubbard_aux(handler,label,l,n,occ,U,J,alpha,beta,J0)
   end subroutine sirius_set_atom_type_hubbard
   
  
positive.code: real(C_DOUBLE), intent(out) :: stress_tensor
  end subroutine
  end interface
  call sirius_get_stress_tensor_aux(handler,label,stress_tensor)
  end subroutine sirius_get_stress_tensor
  
  
negative.code: real(C_DOUBLE), intent(in) :: f_rg
  end subroutine
  end interface
  
  call sirius_set_periodic_function_ptr_aux(handler,label,f_mt,f_rg)
  end subroutine sirius_set_periodic_function_ptr
  
  