code: def _erf_complex(ctx, z):
      z2 = ctx.square_exp_arg(z, -1)
      # z2 = -z**2
      v = (2 / ctx.sqrt(ctx.pi)) * z * ctx.hyp1f1((1, 2), (3, 2), z2)
      if not ctx._re(z):
          v = ctx._im(v) * ctx.j
      return v
  
imports: from .functions import defun, defun_wrapped
label_desc: Use image node_id 1 for calling a global function with example usage: _erf_complex(ctx, z) and returns: v