3 #ifndef DUNE_MATRIX_UTILS_HH
4 #define DUNE_MATRIX_UTILS_HH
9 #include <dune/common/typetraits.hh>
10 #include <dune/common/static_assert.hh>
11 #include <dune/common/fmatrix.hh>
12 #include <dune/common/diagonalmatrix.hh>
13 #include <dune/common/unused.hh>
21 template<
typename B,
typename A>
24 template<
typename K,
int n,
int m>
27 template<
class T,
class A>
47 static typename M::size_type
count(
const M& matrix)
49 typedef typename M::ConstRowIterator RowIterator;
51 RowIterator endRow = matrix.end();
52 typename M::size_type nonZeros = 0;
54 for(RowIterator
row = matrix.begin();
row != endRow; ++
row) {
55 typedef typename M::ConstColIterator Entry;
56 Entry endEntry =
row->end();
57 for(Entry entry =
row->begin(); entry != endEntry; ++entry) {
66 struct NonZeroCounter<1>
69 static typename M::size_type
count(
const M& matrix)
71 return matrix.N()*matrix.M();
81 template<
class Matrix, std::
size_t blocklevel, std::
size_t l=blocklevel>
90 DUNE_UNUSED_PARAMETER(mat);
91 #ifdef DUNE_ISTL_WITH_CHECKING
95 Entry diagonal =
row->find(
row.index());
96 if(diagonal==
row->end())
97 DUNE_THROW(
ISTLError,
"Missing diagonal value in row "<<
row.index()
98 <<
" at block recursion level "<<l-blocklevel);
106 template<
class Matrix, std::
size_t l>
114 DUNE_THROW(
ISTLError,
"Missing diagonal value in row "<<
row.index()
115 <<
" at block recursion level "<<l);
120 template<
typename T1,
typename T2,
typename T3,
typename T4,
typename T5,
121 typename T6,
typename T7,
typename T8,
typename T9>
124 template<
typename T1,
typename T2,
typename T3,
typename T4,
typename T5,
125 typename T6,
typename T7,
typename T8,
typename T9, std::size_t blocklevel, std::size_t l>
137 #ifdef DUNE_ISTL_WITH_CHECKING
168 template<
class G,
class M>
169 bool operator()(
const std::pair<G,M>& p1,
const std::pair<G,M>& p2)
171 return p1.first<p2.first;
176 template<
class M,
class C>
179 typedef typename C::ParallelIndexSet::const_iterator IIter;
180 typedef typename C::OwnerSet OwnerSet;
181 typedef typename C::ParallelIndexSet::GlobalIndex GlobalIndex;
185 for(IIter idx=ooc.indexSet().begin(), eidx=ooc.indexSet().end();
187 gmax=std::max(gmax,idx->global());
189 gmax=ooc.communicator().max(gmax);
190 ooc.buildGlobalLookup();
192 for(IIter idx=ooc.indexSet().begin(), eidx=ooc.indexSet().end();
194 if(OwnerSet::contains(idx->local().attribute()))
196 typedef typename M::block_type Block;
198 std::set<std::pair<GlobalIndex,Block>,CompPair> entries;
201 typedef typename M::ConstColIterator CIter;
202 for(CIter c=mat[idx->local()].begin(), cend=mat[idx->local()].end();
204 const typename C::ParallelIndexSet::IndexPair* pair
205 =ooc.globalLookup().pair(c.index());
207 entries.insert(std::make_pair(pair->global(), *c));
211 GlobalIndex rowidx = idx->global();
212 GlobalIndex cur=std::numeric_limits<GlobalIndex>::max();
214 cur=ooc.communicator().min(rowidx);
217 typedef typename std::set<std::pair<GlobalIndex,Block>,CompPair>::iterator SIter;
218 for(SIter s=entries.begin(), send=entries.end(); s!=send; ++s)
219 os<<idx->global()<<
" "<<s->first<<
" "<<s->second<<std::endl;
225 ooc.freeGlobalLookup();
227 GlobalIndex cur=std::numeric_limits<GlobalIndex>::max();
228 while(cur!=ooc.communicator().min(cur)) ;
232 struct MatrixDimension
236 template<
typename B,
typename TA>
258 if (A.j.get()[k]==c) {
294 std::vector<size_type> coldims(A.
M(),
295 std::numeric_limits<size_type>::max());
300 if (coldims[
col.index()]==std::numeric_limits<size_type>::max())
304 for (
typename std::vector<size_type>::iterator it=coldims.begin();
305 it!=coldims.end(); ++it)
315 template<
typename B,
int n,
int m,
typename TA>
340 template<
typename K,
int n,
int m>
367 template<
typename K,
int n,
int m,
typename TA>
394 template<
typename K,
int n>
421 template<
typename K,
int n>
474 template<
typename T,
typename A>
static size_type coldim(const Matrix &A, size_type c)
Definition: matrixutils.hh:252
static size_type rowdim(const Matrix &A)
Definition: matrixutils.hh:280
A::size_type size_type
Type for indices and sizes.
Definition: matrix.hh:41
static size_type coldim(const Matrix &, size_type)
Definition: matrixutils.hh:405
row_type::const_iterator ConstColIterator
Const iterator for the entries of each row.
Definition: matrix.hh:53
int countNonZeros(const M &matrix)
Get the number of nonzero fields in the matrix.
Definition: matrixutils.hh:155
FieldMatrix< K, n, m > Matrix
Definition: matrixutils.hh:343
void printGlobalSparseMatrix(const M &mat, C &ooc, std::ostream &os)
Definition: matrixutils.hh:177
size_type M() const
number of columns (counted in blocks)
Definition: bcrsmatrix.hh:1856
static size_type coldim(const ThisMatrix &A)
Definition: matrixutils.hh:388
static size_type rowdim(const Matrix &, size_type)
Definition: matrixutils.hh:321
A::size_type size_type
The type for the index access and the size.
Definition: bcrsmatrix.hh:448
static size_type rowdim(const Matrix &, size_type)
Definition: matrixutils.hh:346
Matrix::size_type size_type
Definition: matrixutils.hh:241
ScaledIdentityMatrix< K, n > Matrix
Definition: matrixutils.hh:424
Definition: matrixutils.hh:486
static size_type coldim(const Matrix &A)
Definition: matrixutils.hh:335
static size_type rowdim(const Matrix &, size_type)
Definition: matrixutils.hh:427
Matrix::size_type size_type
Definition: matrixutils.hh:344
static void check(const Matrix &mat)
Definition: matrixutils.hh:109
Check whether the a matrix has diagonal values on blocklevel recursion levels.
Definition: matrixutils.hh:82
B * getptr()
get pointer
Definition: bvector.hh:1051
A sparse block matrix with compressed row storage.
Definition: bcrsmatrix.hh:413
Matrix< FieldMatrix< K, n, m >, TA > ThisMatrix
Definition: matrixutils.hh:370
Iterator end()
Get iterator to one beyond last row.
Definition: bcrsmatrix.hh:632
std::size_t size_type
The type used for the index access and size operations.
Definition: scaledidmatrix.hh:41
A generic dynamic dense matrix.
Definition: matrix.hh:24
static size_type rowdim(const ThisMatrix &, size_type)
Definition: matrixutils.hh:373
static size_type rowdim(const Matrix &, size_type)
Definition: matrixutils.hh:400
static size_type rowdim(const ThisMatrix &A)
Definition: matrixutils.hh:383
size_type M() const
Return the number of columns.
Definition: matrix.hh:165
Definition: bcrsmatrix.hh:72
RowIterator end()
Get iterator to one beyond last row.
Definition: matrix.hh:85
size_type N() const
number of rows (counted in blocks)
Definition: bcrsmatrix.hh:1850
static size_type rowdim(const Matrix &)
Definition: matrixutils.hh:437
Col col
Definition: matrixmatrix.hh:347
Matrix::size_type size_type
Definition: matrixutils.hh:398
BCRSMatrix< FieldMatrix< B, n, m >,TA > Matrix
Definition: matrixutils.hh:318
Dune::CheckIfDiagonalPresent< MultiTypeBlockMatrix< T1, T2, T3, T4, T5, T6, T7, T8, T9 >, blocklevel, l >::Matrix MultiTypeBlockMatrix< T1, T2, T3, T4, T5, T6, T7, T8, T9 > Matrix
Definition: matrixutils.hh:129
static size_type rowdim(const Matrix &A, size_type i)
Definition: matrixutils.hh:243
static size_type coldim(const Matrix &)
Definition: matrixutils.hh:415
Matrix & A
Definition: matrixmatrix.hh:216
static size_type coldim(const Matrix &A)
Definition: matrixutils.hh:287
iterator class for sequential access
Definition: basearray.hh:594
bool operator()(const T *l, const T *r)
Definition: matrixutils.hh:488
VariableBlockVector< T, A >::ConstIterator ConstRowIterator
Const iterator over the matrix rows.
Definition: matrix.hh:50
size_type * getindexptr()
get pointer
Definition: bvector.hh:1057
static size_type coldim(const Matrix &)
Definition: matrixutils.hh:361
This file implements a quadratic matrix of fixed size which is a multiple of the identity.
Matrix::block_type block_type
Definition: matrixutils.hh:240
A multiple of the identity matrix of static size.
Definition: scaledidmatrix.hh:27
Test whether a type is an ISTL Matrix.
Definition: matrixutils.hh:452
Row row
Definition: matrixmatrix.hh:345
Matrix::size_type size_type
Definition: matrixutils.hh:319
RowIterator begin()
Get iterator to first row.
Definition: matrix.hh:79
BCRSMatrix< B, TA > Matrix
Definition: matrixutils.hh:239
Iterator begin()
Get iterator to first row.
Definition: bcrsmatrix.hh:626
size_type N() const
Return the number of rows.
Definition: matrix.hh:160
std::size_t count
Definition: matrixmatrix.hh:215
B block_type
export the type representing the components
Definition: bcrsmatrix.hh:439
derive error class from the base class in common
Definition: istlexception.hh:16
static size_type coldim(const Matrix &)
Definition: matrixutils.hh:442
True if T is an ISTL matrix.
Definition: matrixutils.hh:458
static size_type coldim(const Matrix &, size_type)
Definition: matrixutils.hh:351
Matrix::size_type size_type
Definition: matrixutils.hh:425
static size_type coldim(const Matrix &, size_type)
Definition: matrixutils.hh:326
static size_type rowdim(const Matrix &A)
Definition: matrixutils.hh:331
static size_type rowdim(const Matrix &)
Definition: matrixutils.hh:410
DiagonalMatrix< K, n > Matrix
Definition: matrixutils.hh:397
static void check(const Matrix &mat)
Check whether the a matrix has diagonal values on blocklevel recursion levels.
Definition: matrixutils.hh:88
Matrix & mat
Definition: matrixmatrix.hh:343
Dune::CheckIfDiagonalPresent< MultiTypeBlockMatrix< T1, T2, T3, T4, T5, T6, T7, T8, T9 >, blocklevel, l >::check static void check(const Matrix &mat)
Check whether the a matrix has diagonal values on blocklevel recursion levels.
Definition: matrixutils.hh:135
Definition: matrixutils.hh:122
static size_type coldim(const ThisMatrix &, size_type)
Definition: matrixutils.hh:378
size_type getsize() const
get size
Definition: bvector.hh:1074
static size_type rowdim(const Matrix &)
Definition: matrixutils.hh:356
ThisMatrix::size_type size_type
Definition: matrixutils.hh:371
Definition: matrixutils.hh:25
static size_type coldim(const Matrix &, size_type)
Definition: matrixutils.hh:432
Iterator access to matrix rows
Definition: bcrsmatrix.hh:528