From c5d7c9f0def693e77797e6e43cb1fd2e7f8ae2fd Mon Sep 17 00:00:00 2001 From: Gael Guennebaud Date: Thu, 7 Jan 2010 21:15:32 +0100 Subject: remove the Triangular suffix to Upper, Lower, UnitLower, etc, and remove the respective bit flags --- Eigen/src/Core/SolveTriangular.h | 56 ++++++++++++++++++++-------------------- 1 file changed, 28 insertions(+), 28 deletions(-) (limited to 'Eigen/src/Core/SolveTriangular.h') diff --git a/Eigen/src/Core/SolveTriangular.h b/Eigen/src/Core/SolveTriangular.h index 9dc019d17..73cf77387 100644 --- a/Eigen/src/Core/SolveTriangular.h +++ b/Eigen/src/Core/SolveTriangular.h @@ -57,7 +57,7 @@ struct ei_triangular_solver_selector LhsProductTraits; typedef typename LhsProductTraits::ExtractType ActualLhsType; enum { - IsLowerTriangular = ((Mode&LowerTriangularBit)==LowerTriangularBit) + IsLower = ((Mode&Lower)==Lower) }; static void run(const Lhs& lhs, Rhs& other) { @@ -65,20 +65,20 @@ struct ei_triangular_solver_selector0; - IsLowerTriangular ? pi+=PanelWidth : pi-=PanelWidth) + for(int pi=IsLower ? 0 : size; + IsLower ? pi0; + IsLower ? pi+=PanelWidth : pi-=PanelWidth) { - int actualPanelWidth = std::min(IsLowerTriangular ? size - pi : pi, PanelWidth); + int actualPanelWidth = std::min(IsLower ? size - pi : pi, PanelWidth); - int r = IsLowerTriangular ? pi : size - pi; // remaining size + int r = IsLower ? pi : size - pi; // remaining size if (r > 0) { // let's directly call the low level product function because: // 1 - it is faster to compile // 2 - it is slighlty faster at runtime - int startRow = IsLowerTriangular ? pi : pi-actualPanelWidth; - int startCol = IsLowerTriangular ? 0 : pi; + int startRow = IsLower ? pi : pi-actualPanelWidth; + int startCol = IsLower ? 0 : pi; VectorBlock target(other,startRow,actualPanelWidth); ei_cache_friendly_product_rowmajor_times_vector( @@ -89,12 +89,12 @@ struct ei_triangular_solver_selector0) other.coeffRef(i) -= (lhs.row(i).segment(s,k).transpose().cwiseProduct(other.segment(s,k))).sum(); - if(!(Mode & UnitDiagBit)) + if(!(Mode & UnitDiag)) other.coeffRef(i) /= lhs.coeff(i,i); } } @@ -111,7 +111,7 @@ struct ei_triangular_solver_selector::size, - IsLowerTriangular = ((Mode&LowerTriangularBit)==LowerTriangularBit) + IsLower = ((Mode&Lower)==Lower) }; static void run(const Lhs& lhs, Rhs& other) @@ -120,26 +120,26 @@ struct ei_triangular_solver_selector0; - IsLowerTriangular ? pi+=PanelWidth : pi-=PanelWidth) + for(int pi=IsLower ? 0 : size; + IsLower ? pi0; + IsLower ? pi+=PanelWidth : pi-=PanelWidth) { - int actualPanelWidth = std::min(IsLowerTriangular ? size - pi : pi, PanelWidth); - int startBlock = IsLowerTriangular ? pi : pi-actualPanelWidth; - int endBlock = IsLowerTriangular ? pi + actualPanelWidth : 0; + int actualPanelWidth = std::min(IsLower ? size - pi : pi, PanelWidth); + int startBlock = IsLower ? pi : pi-actualPanelWidth; + int endBlock = IsLower ? pi + actualPanelWidth : 0; for(int k=0; k0) other.segment(s,r) -= other.coeffRef(i) * Block(lhs, s, i, r, 1); } - int r = IsLowerTriangular ? size - endBlock : startBlock; // remaining size + int r = IsLower ? size - endBlock : startBlock; // remaining size if (r > 0) { // let's directly call the low level product function because: @@ -198,16 +198,16 @@ struct ei_triangular_solver_unroller; template struct ei_triangular_solver_unroller { enum { - IsLowerTriangular = ((Mode&LowerTriangularBit)==LowerTriangularBit), - I = IsLowerTriangular ? Index : Size - Index - 1, - S = IsLowerTriangular ? 0 : I+1 + IsLower = ((Mode&Lower)==Lower), + I = IsLower ? Index : Size - Index - 1, + S = IsLower ? 0 : I+1 }; static void run(const Lhs& lhs, Rhs& rhs) { if (Index>0) rhs.coeffRef(I) -= ((lhs.row(I).template segment(S).transpose()).cwiseProduct(rhs.template segment(S))).sum(); - if(!(Mode & UnitDiagBit)) + if(!(Mode & UnitDiag)) rhs.coeffRef(I) /= lhs.coeff(I,I); ei_triangular_solver_unroller::run(lhs,rhs); @@ -245,8 +245,8 @@ void TriangularView::solveInPlace(const MatrixBase::Flags & RowMajorBit && OtherDerived::IsVectorAtCompileTime }; typedef typename ei_meta_if