16#ifndef dealii_mesh_worker_local_results_h
17#define dealii_mesh_worker_local_results_h
221 template <
typename number>
355 template <
typename MatrixType>
367 template <
typename MatrixType>
387 template <
typename StreamType>
401 std::vector<number>
J;
407 std::vector<BlockVector<number>>
R;
413 std::vector<MatrixBlock<FullMatrix<number>>>
M1;
421 std::vector<MatrixBlock<FullMatrix<number>>>
M2;
431 template <
typename number>
439 template <
typename number>
447 template <
typename number>
448 template <
typename MatrixType>
454 M1.resize(matrices.
size());
456 M2.resize(matrices.
size());
457 for (
unsigned int i = 0; i < matrices.
size(); ++i)
459 const unsigned int row = matrices.
block(i).
row;
473 template <
typename number>
474 template <
typename MatrixType>
480 M1.resize(matrices.
size());
482 M2.resize(matrices.
size());
483 for (
unsigned int i = 0; i < matrices.
size(); ++i)
486 const unsigned int row = o[o.
min_level()].row;
487 const unsigned int col = o[o.
min_level()].column;
500 template <
typename number>
508 for (
unsigned int i = 0; i < n; ++i)
521 template <
typename number>
524 const unsigned int nv)
526 quadrature_data.reinit(np, nv);
530 template <
typename number>
538 template <
typename number>
546 template <
typename number>
554 template <
typename number>
558 return quadrature_data.n_rows();
562 template <
typename number>
566 return quadrature_data.n_cols();
570 template <
typename number>
579 template <
typename number>
588 template <
typename number>
602 template <
typename number>
605 const unsigned int i)
607 return quadrature_data(k, i);
611 template <
typename number>
615 return quadrature_data;
619 template <
typename number>
628 template <
typename number>
637 template <
typename number>
651 template <
typename number>
654 const unsigned int i)
const
656 return quadrature_data(k, i);
660 template <
typename number>
661 template <
typename StreamType>
665 os <<
"J: " << J.size() << std::endl;
666 os <<
"R: " << R.size() << std::endl;
667 for (
unsigned int i = 0; i < R.size(); ++i)
669 os <<
" " << R[i].n_blocks() <<
" -";
670 for (
unsigned int j = 0; j < R[i].n_blocks(); ++j)
671 os <<
' ' << R[i].block(j).size();
674 os <<
"M: " << M1.size() <<
" face " << M2.size() << std::endl;
675 for (
unsigned int i = 0; i < M1.size(); ++i)
677 os <<
" " << M1[i].row <<
',' << M1[i].column <<
' '
678 << M1[i].matrix.m() <<
'x' << M1[i].matrix.n();
680 os <<
" face " << M2[i].row <<
',' << M2[i].column <<
' '
681 << M2[i].matrix.m() <<
'x' << M2[i].matrix.n();
unsigned int min_level() const
const value_type & block(size_type i) const
unsigned int size() const
unsigned int size() const
Number of stored data objects.
const value_type & block(size_type i) const
MatrixBlock< FullMatrix< number > > & matrix(const unsigned int i, const bool external=false)
unsigned int n_quadrature_values() const
BlockVector< number > & vector(const unsigned int i)
unsigned int n_vectors() const
number & quadrature_value(const unsigned int k, const unsigned int i)
number value(const unsigned int i) const
std::vector< BlockVector< number > > R
std::size_t memory_consumption() const
const BlockVector< number > & vector(const unsigned int i) const
const MatrixBlock< FullMatrix< number > > & matrix(const unsigned int i, const bool external=false) const
number quadrature_value(const unsigned int k, const unsigned int i) const
Table< 2, number > quadrature_data
void initialize_matrices(const MatrixBlockVector< MatrixType > &matrices, bool both)
void initialize_matrices(const MGMatrixBlockVector< MatrixType > &matrices, bool both)
std::vector< MatrixBlock< FullMatrix< number > > > M2
void reinit(const BlockIndices &local_sizes)
void initialize_matrices(const unsigned int n, bool both)
std::vector< MatrixBlock< FullMatrix< number > > > M1
unsigned int n_values() const
void initialize_quadrature(const unsigned int np, const unsigned int nv)
Table< 2, number > & quadrature_values()
void initialize_numbers(const unsigned int n)
void initialize_vectors(const unsigned int n)
unsigned int n_matrices() const
void print_debug(StreamType &os) const
unsigned int n_quadrature_points() const
number & value(const unsigned int i)
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define AssertIndexRange(index, range)