code: def fisher_rao_dist(point_a, point_b, n_draws):
      return (
          2
          * gs.sqrt(n_draws)
          * gs.arccos(gs.sum(gs.sqrt(point_a * point_b), axis=-1))
      )
  
imports: import pytest
  import geomstats.backend
  from geomstats.test.vectorization import generate_vectorization_data
  from geomstats.test_cases.geometry.base import LevelSetTestCase
  from geomstats.test_cases.geometry.riemannian_metric import RiemannianMetricTestCase
  from geomstats.test_cases.information_geometry.base import InformationManifoldMixinTestCase
label_desc: Use image node_id 1 for calling a global function with example usage: fisher_rao_dist(point_a, point_b, n_draws) and returns: unknown