code: def sparse_hellinger(ind1, data1, ind2, data2):
      aux_inds, aux_data = sparse_mul(ind1, data1, ind2, data2)
      result = 0.0
      norm1 = np.sum(data1)
      norm2 = np.sum(data2)
      sqrt_norm_prod = np.sqrt(norm1 * norm2)
  
      for val in aux_data:
          result += np.sqrt(val)
  
      if norm1 == 0.0 and norm2 == 0.0:
          return 0.0
      elif norm1 == 0.0 or norm2 == 0.0:
          return 1.0
      elif result > sqrt_norm_prod:
          return 0.0
      else:
          return np.sqrt(1.0 - (result / sqrt_norm_prod))
  
imports: from __future__ import print_function
  import locale
  import numpy
  import numba
  from pynndescent.utils import norm, tau_rand
  from pynndescent.distances import kantorovich, jensen_shannon_divergence, symmetric_kl_divergence
label_desc: Use image node_id 35 for calling a global function with example usage: sparse_hellinger(ind1, data1, ind2, data2) and returns: int, int, int, np