code: def rmn(rowcov, colcov):
      """
      Generate random draws from a zero-mean matrix-normal distribution.
  
      Parameters
      -----------
      rowcov : np.ndarray
          Row covariance (assumed to be positive definite)
      colcov : np.ndarray
          Column covariance (assumed to be positive definite)
      """
  
      Z = norm.rvs(size=(rowcov.shape[0], colcov.shape[0]))
      return cholesky(rowcov).dot(Z).dot(cholesky(colcov))
  
file_path: brainiak/brainiak/matnormal/utils.py
imports: import tensorflow
  import tensorflow_probability
  from scipy.stats import norm
  from numpy.linalg import cholesky
  import numpy
label_desc: Use image node_id 1 for calling a global function with example usage: rmn(rowcov, colcov) and returns: cholesky