code: def testGEVSampleMultidimensionalVar(self):
      loc = np.array([2.0, 4.0, 5.0], dtype=self._dtype)
      scale = np.array([1.0, 0.8, 0.5], dtype=self._dtype)
      conc = np.array([0.2], dtype=self._dtype)
      gev_dist = stats.genextreme(-conc, loc=loc, scale=scale)
      n = int(1e6)
  
      dist = gev.GeneralizedExtremeValue(
          loc=self.make_tensor(loc),
          scale=self.make_tensor(scale),
          concentration=self.make_tensor(conc),
          validate_args=True,
      )
  
      samples = dist.sample(n, seed=test_util.test_seed())
      sample_values = self.evaluate(samples)
      self.assertAllClose(
          gev_dist.var(), sample_values.var(axis=0), rtol=0.03, atol=0
      )
  
imports: import numpy
  from scipy import stats
  import tensorflow.compat.v1
  import tensorflow.compat.v2
  from tensorflow_probability.python.distributions import gev
  from tensorflow_probability.python.internal import test_util
  from tensorflow_probability.python.math import gradient
label_desc: Use image node_id 14 for calling the _GEVTest obj's underlying member method code with example usage: obj.testGEVSampleMultidimensionalVar() without return types