[-0.7486,  1.3544,  0.1294],
              [ 1.4551,  0.1294,  1.6724]])
      >>> l
      tensor([[ 1.5528,  0.0000,  0.0000],
              [-0.4821,  1.0592,  0.0000],
              [ 0.9371,  0.5487,  0.7023]])
      >>> l @ l.mT
      tensor([[ 2.4112, -0.7486,  1.4551],
              [-0.7486,  1.3544,  0.1294],
              [ 1.4551,  0.1294,  1.6724]])
      >>> a = torch.randn(3, 2, 2) # Example for batched input
      >>> a = a @ a.mT + 1e-03 # make symmetric positive-definite
      >>> l = torch.cholesky(a)
      >>> z = l @ l.mT
      >>> torch.dist(z, a)
      tensor(2.3842e-07)