pytorch - 你如何在pytorch中生成正定矩阵？

Question

我试图用随机生成的协方差矩阵定义多元高斯分布:

psi = torch.zeros(512).normal_(0., 1.).requires_grad_(True)

# Generate random matrix
Sigma_k = torch.rand(512, 512)
# Make it symmetric positive
Sigma_k = Sigma_k * Sigma_k.t()
# Make it definite
Sigma_k.add_(0.001, torch.eye(512)).requires_grad_(True)

multivariate_normal.MultivariateNormal(psi, Sigma_k)

但我最终得到了一个异常(exception):

RuntimeError: Lapack Error in potrf : the leading minor of order 2 is not positive definite at /Users/soumith/mc3build/conda-bld/pytorch_1549597882250/work/aten/src/TH/generic/THTensorLapack.cpp:658

生成正定方阵的正确方法是什么？

Answer 1

在“确定”中
tensor.add()不变tensor ，但只返回更改后的版本。
您要使用 tensor.add_()