code: def test_gh2002_regression():
      # Add a check that broadcast works in situations where only some
      # x-values are compatible with some of the shape arguments.
      x = np.r_[-2:2:101j]
      a = np.r_[-np.ones(50), np.ones(51)]
      expected = [
          stats.truncnorm.pdf(_x, _a, np.inf) for _x, _a in zip(x, a)
      ]
      ans = stats.truncnorm.pdf(x, a, np.inf)
      npt.assert_array_almost_equal(ans, expected)
  
imports: import sys
  import numpy
  import numpy.testing
  import pytest
  from pytest import raises
  from scipy.integrate import IntegrationWarning
  import itertools
  from scipy import stats
  from .common_tests import check_normalization, check_moment, check_mean_expect, check_var_expect, check_skew_expect, check_kurt_expect, check_entropy, check_private_entropy, check_entropy_vect_scale, check_edge_support, check_named_args, check_random_state_property, check_meth_dtype, check_ppf_dtype, check_cmplx_deriv, check_pickling, check_rvs_broadcast, check_freezing, check_munp_expect
  from scipy.stats._distr_params import distcont
  from scipy.stats._distn_infrastructure import rv_continuous_frozen
label_desc: Use image node_id 13 for calling a global function with example usage: test_gh2002_regression() without return types