faster computation of sigmas in NoiseModel
parent
b0f7d3498b
commit
1cb7b223ae
|
@ -143,8 +143,8 @@ bool Gaussian::equals(const Base& expected, double tol) const {
|
||||||
|
|
||||||
/* ************************************************************************* */
|
/* ************************************************************************* */
|
||||||
Vector Gaussian::sigmas() const {
|
Vector Gaussian::sigmas() const {
|
||||||
// TODO(frank): can this be done faster?
|
Matrix Rinv = thisR().inverse(); // inverse of triangular matrix is fast
|
||||||
return Vector((thisR().transpose() * thisR()).inverse().diagonal()).cwiseSqrt();
|
return Vector((Rinv * Rinv.transpose()).diagonal()).cwiseSqrt();
|
||||||
}
|
}
|
||||||
|
|
||||||
/* ************************************************************************* */
|
/* ************************************************************************* */
|
||||||
|
|
Loading…
Reference in New Issue