Eigen  3.3.0
 
Loading...
Searching...
No Matches
Ordering.h
1
2// This file is part of Eigen, a lightweight C++ template library
3// for linear algebra.
4//
5// Copyright (C) 2012 Désiré Nuentsa-Wakam <desire.nuentsa_wakam@inria.fr>
6//
7// This Source Code Form is subject to the terms of the Mozilla
8// Public License v. 2.0. If a copy of the MPL was not distributed
9// with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
10
11#ifndef EIGEN_ORDERING_H
12#define EIGEN_ORDERING_H
13
14namespace Eigen {
15
16#include "Eigen_Colamd.h"
17
18namespace internal {
19
26template<typename MatrixType>
27void ordering_helper_at_plus_a(const MatrixType& A, MatrixType& symmat)
28{
29 MatrixType C;
30 C = A.transpose(); // NOTE: Could be costly
31 for (int i = 0; i < C.rows(); i++)
32 {
33 for (typename MatrixType::InnerIterator it(C, i); it; ++it)
34 it.valueRef() = 0.0;
35 }
36 symmat = C + A;
37}
38
39}
40
41#ifndef EIGEN_MPL2_ONLY
42
51template <typename StorageIndex>
53{
54 public:
56
60 template <typename MatrixType>
61 void operator()(const MatrixType& mat, PermutationType& perm)
62 {
63 // Compute the symmetric pattern
65 internal::ordering_helper_at_plus_a(mat,symm);
66
67 // Call the AMD routine
68 //m_mat.prune(keep_diag());
69 internal::minimum_degree_ordering(symm, perm);
70 }
71
73 template <typename SrcType, unsigned int SrcUpLo>
75 {
77
78 // Call the AMD routine
79 // m_mat.prune(keep_diag()); //Remove the diagonal elements
80 internal::minimum_degree_ordering(C, perm);
81 }
82};
83
84#endif // EIGEN_MPL2_ONLY
85
94template <typename StorageIndex>
96{
97 public:
99
101 template <typename MatrixType>
102 void operator()(const MatrixType& /*mat*/, PermutationType& perm)
103 {
104 perm.resize(0);
105 }
106
107};
108
117template<typename StorageIndex>
119{
120 public:
123
127 template <typename MatrixType>
128 void operator() (const MatrixType& mat, PermutationType& perm)
129 {
130 eigen_assert(mat.isCompressed() && "COLAMDOrdering requires a sparse matrix in compressed mode. Call .makeCompressed() before passing it to COLAMDOrdering");
131
132 StorageIndex m = StorageIndex(mat.rows());
133 StorageIndex n = StorageIndex(mat.cols());
134 StorageIndex nnz = StorageIndex(mat.nonZeros());
135 // Get the recommended value of Alen to be used by colamd
136 StorageIndex Alen = internal::colamd_recommended(nnz, m, n);
137 // Set the default parameters
138 double knobs [COLAMD_KNOBS];
139 StorageIndex stats [COLAMD_STATS];
140 internal::colamd_set_defaults(knobs);
141
142 IndexVector p(n+1), A(Alen);
143 for(StorageIndex i=0; i <= n; i++) p(i) = mat.outerIndexPtr()[i];
144 for(StorageIndex i=0; i < nnz; i++) A(i) = mat.innerIndexPtr()[i];
145 // Call Colamd routine to compute the ordering
146 StorageIndex info = internal::colamd(m, n, Alen, A.data(), p.data(), knobs, stats);
147 EIGEN_UNUSED_VARIABLE(info);
148 eigen_assert( info && "COLAMD failed " );
149
150 perm.resize(n);
151 for (StorageIndex i = 0; i < n; i++) perm.indices()(p(i)) = i;
152 }
153};
154
155} // end namespace Eigen
156
157#endif
Definition: Ordering.h:53
void operator()(const SparseSelfAdjointView< SrcType, SrcUpLo > &mat, PermutationType &perm)
Definition: Ordering.h:74
void operator()(const MatrixType &mat, PermutationType &perm)
Definition: Ordering.h:61
Definition: Ordering.h:119
void operator()(const MatrixType &mat, PermutationType &perm)
Definition: Ordering.h:128
The matrix class, also used for vectors and row-vectors.
Definition: Matrix.h:180
Definition: Ordering.h:96
void operator()(const MatrixType &, PermutationType &perm)
Definition: Ordering.h:102
void resize(Index newSize)
Definition: PermutationMatrix.h:136
Permutation matrix.
Definition: PermutationMatrix.h:309
const IndicesType & indices() const
Definition: PermutationMatrix.h:388
const Scalar * data() const
Definition: PlainObjectBase.h:249
A versatible sparse matrix representation.
Definition: SparseMatrix.h:94
Pseudo expression to manipulate a triangular sparse matrix as a selfadjoint matrix.
Definition: SparseSelfAdjointView.h:45
Namespace containing all symbols from the Eigen library.
Definition: Core:287