code: def test_skewnorm_pdf_gh16038():
      rng = np.random.default_rng(0)
      x, a = -np.inf, 0
      npt.assert_equal(stats.skewnorm.pdf(x, a), stats.norm.pdf(x))
      x, a = rng.random(size=(3, 3)), rng.random(size=(3, 3))
      mask = rng.random(size=(3, 3)) < 0.5
      a[mask] = 0
      x_norm = x[mask]
      res = stats.skewnorm.pdf(x, a)
      npt.assert_equal(res[mask], stats.norm.pdf(x_norm))
      npt.assert_equal(
          res[~mask], stats.skewnorm.pdf(x[~mask], a[~mask])
      )
  
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 43 for calling a global function with example usage: test_skewnorm_pdf_gh16038() without return types