diff options
Diffstat (limited to 'eigen/Eigen/src/Eigen2Support/TriangularSolver.h')
-rw-r--r-- | eigen/Eigen/src/Eigen2Support/TriangularSolver.h | 42 |
1 files changed, 42 insertions, 0 deletions
diff --git a/eigen/Eigen/src/Eigen2Support/TriangularSolver.h b/eigen/Eigen/src/Eigen2Support/TriangularSolver.h new file mode 100644 index 0000000..ebbeb3b --- /dev/null +++ b/eigen/Eigen/src/Eigen2Support/TriangularSolver.h @@ -0,0 +1,42 @@ +// This file is part of Eigen, a lightweight C++ template library +// for linear algebra. +// +// Copyright (C) 2010 Gael Guennebaud <gael.guennebaud@inria.fr> +// +// This Source Code Form is subject to the terms of the Mozilla +// Public License v. 2.0. If a copy of the MPL was not distributed +// with this file, You can obtain one at http://mozilla.org/MPL/2.0/. + +#ifndef EIGEN_TRIANGULAR_SOLVER2_H +#define EIGEN_TRIANGULAR_SOLVER2_H + +namespace Eigen { + +const unsigned int UnitDiagBit = UnitDiag; +const unsigned int SelfAdjointBit = SelfAdjoint; +const unsigned int UpperTriangularBit = Upper; +const unsigned int LowerTriangularBit = Lower; + +const unsigned int UpperTriangular = Upper; +const unsigned int LowerTriangular = Lower; +const unsigned int UnitUpperTriangular = UnitUpper; +const unsigned int UnitLowerTriangular = UnitLower; + +template<typename ExpressionType, unsigned int Added, unsigned int Removed> +template<typename OtherDerived> +typename ExpressionType::PlainObject +Flagged<ExpressionType,Added,Removed>::solveTriangular(const MatrixBase<OtherDerived>& other) const +{ + return m_matrix.template triangularView<Added>().solve(other.derived()); +} + +template<typename ExpressionType, unsigned int Added, unsigned int Removed> +template<typename OtherDerived> +void Flagged<ExpressionType,Added,Removed>::solveTriangularInPlace(const MatrixBase<OtherDerived>& other) const +{ + m_matrix.template triangularView<Added>().solveInPlace(other.derived()); +} + +} // end namespace Eigen + +#endif // EIGEN_TRIANGULAR_SOLVER2_H |