aboutsummaryrefslogtreecommitdiffhomepage
path: root/Eigen/src/Core/products
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2010-11-05 14:36:34 +0100
committerGravatar Gael Guennebaud <g.gael@free.fr>2010-11-05 14:36:34 +0100
commit572b5585e3119a96250a0b4e4c51ad61b41137fb (patch)
tree1787ee33b050d961f809827a27042a2c52a9ac67 /Eigen/src/Core/products
parent0e30c4ae3f35523302f1449f67a2be714e30beb8 (diff)
fix Eigen's trsv for complexes
Diffstat (limited to 'Eigen/src/Core/products')
-rw-r--r--Eigen/src/Core/products/TriangularSolverVector.h4
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);
}
}
}