def METHOD_NAME(ax, theta, lik):
  dim = theta.shape[1]
  if (dim == 1):
    return
  for i in range(dim):
    for j in range(i + 1, dim):
      if lik:
        ax[i, j].scatter(
            theta[:, j], theta[:, i], marker='o', s=3, alpha=0.5, c=lik)
      else:
        ax[i, j].plot(theta[:, j], theta[:, i], '.', markersize=1)
      ax[i, j].set_xticklabels([])
      ax[i, j].set_yticklabels([])
      ax[i, j].grid(b=True, which='both')