From 572b5585e3119a96250a0b4e4c51ad61b41137fb Mon Sep 17 00:00:00 2001 From: Gael Guennebaud Date: Fri, 5 Nov 2010 14:36:34 +0100 Subject: fix Eigen's trsv for complexes --- Eigen/src/Core/products/TriangularSolverVector.h | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'Eigen/src/Core/products') 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_vector0) rhs[i] -= (cjLhs.row(i).segment(s,k).transpose().cwiseProduct(Map >(rhs+s,k))).sum(); - + if(!(Mode & UnitDiag)) - rhs[i] /= lhs(i,i); + rhs[i] /= cjLhs(i,i); } } } -- cgit v1.2.3