How to use the harmonica.forward.point_mass.jit_point_mass_spherical function in harmonica

distance_size_ratio,
                stack,
                small_tesseroids,
                radial_discretization,
            )
            # Get total number of point masses, their coordinates and weights
            n_point_masses = n_nodes * n_splits
            tesseroids_to_point_masses(
                small_tesseroids[:n_splits],
                glq_nodes,
                glq_weights,
                point_masses,
                weights,
            )
            # Compute gravitational fields
            jit_point_mass_spherical(
                coordinates[0][m : m + 1],  # slice lon to pass a single element array
                coordinates[1][m : m + 1],  # slice lat to pass a single element array
                coordinates[2][m : m + 1],  # slice rad to pass a single element array
                point_masses[0, :n_point_masses],
                point_masses[1, :n_point_masses],
                point_masses[2, :n_point_masses],
                density[l] * weights[:n_point_masses],
                result[m : m + 1],
                kernel,
            )