code: def numerical_jvp(f, primals, tangents, eps=EPS):
      delta = scalar_mul(tangents, eps)
      f_pos = f(*add(primals, delta))
      f_neg = f(*sub(primals, delta))
      return scalar_mul(safe_sub(f_pos, f_neg), 0.5 / eps)
  
imports: from functools import partial
  import operator
  from jax._src import api
  from jax._src import config
  from jax._src import dtypes
  from jax._src.tree_util import tree_map, tree_reduce
  import numpy
label_desc: Use image node_id 14 for calling a global function with example usage: numerical_jvp(f, primals, tangents, eps) and returns: scalar_mul