|
|
1
2
对于标题中的公式(以前为logdet(AA^T)):一些随机As的det(AA^T)可以简单地为0。 该函数将失败,因为计算日志(0)无效。 请注意,理论上det(AA^T)不能为负,因为AA^T是 positive semi-definite matrix (这意味着所有特征值都是非负的,并且意味着det>=0)。 对于代码中的公式(logdet(1+AA^T))
你可能应该使用
从其 documentation : 计算数组行列式的符号和(自然)对数。 如果数组的行列式非常小或非常大,那么对det的调用可能会溢出或下溢。 这个例程对此类问题更为稳健,因为它计算的是行列式的对数,而不是行列式本身。” |