diff options
author | Gael Guennebaud <g.gael@free.fr> | 2010-11-05 14:36:34 +0100 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2010-11-05 14:36:34 +0100 |
commit | 572b5585e3119a96250a0b4e4c51ad61b41137fb (patch) | |
tree | 1787ee33b050d961f809827a27042a2c52a9ac67 /Eigen | |
parent | 0e30c4ae3f35523302f1449f67a2be714e30beb8 (diff) |
fix Eigen's trsv for complexes
Diffstat (limited to 'Eigen')
-rw-r--r-- | Eigen/src/Core/products/TriangularSolverVector.h | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/Eigen/src/Core/products/TriangularSolverVector.h b/Eigen/src/Core/products/TriangularSolverVector.h index db1e7f029..16184abb5 100644 --- a/Eigen/src/Core/products/TriangularSolverVector.h +++ b/Eigen/src/Core/products/TriangularSolverVector.h @@ -85,9 +85,9 @@ struct triangular_solve_vector<LhsScalar, RhsScalar, Index, OnTheLeft, Mode, Con Index s = IsLower ? pi : i+1; if (k>0) rhs[i] -= (cjLhs.row(i).segment(s,k).transpose().cwiseProduct(Map<Matrix<RhsScalar,Dynamic,1> >(rhs+s,k))).sum(); - + if(!(Mode & UnitDiag)) - rhs[i] /= lhs(i,i); + rhs[i] /= cjLhs(i,i); } } } |