│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
1066#ifdef DUNE_ISTL_WITH_CHECKING
│ │ │ │ -
1067 if (i<0 || i>=
N()) DUNE_THROW(
ISTLError,
"row index out of range");
│ │ │ │ -
1068 if (j<0 || i>=
M()) DUNE_THROW(
ISTLError,
"column index out of range");
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
1029 template<
typename B>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1035 if constexpr (IsNumber<B>())
│ │ │ │ +
1036 ostr << rowidx <<
" " << colidx <<
" " << entry << std::endl;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1039 for (
auto row=entry.begin(); row != entry.end(); ++row, ++rowidx) {
│ │ │ │ +
│ │ │ │ +
1041 for (
auto col = row->begin();
col != row->end(); ++
col, ++coli)
│ │ │ │ +
1042 ostr<< rowidx<<
" "<<coli<<
" "<<*
col<<std::endl;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1048 template<
typename V>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1050 const std::integral_constant<int,1>&,
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1053 ostr<<Simd::lane(lane,entry)<<std::endl;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1057 template<
typename V>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1059 const std::integral_constant<int,0>&,
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1062 using namespace MatrixMarketImpl;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1065 const int isnumeric = mm_numeric_type<Simd::Scalar<typename V::block_type>>::is_numeric;
│ │ │ │ +
1066 typedef typename V::const_iterator VIter;
│ │ │ │ +
│ │ │ │ +
1068 for(VIter i=vector.begin(); i != vector.end(); ++i)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1071 std::integral_constant<int,isnumeric>(),
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
1075 template<
typename T,
typename A>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1078 return vector.size();
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1081 template<
typename T,
typename A,
int i>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1084 return vector.size()*i;
│ │ │ │ +
│ │ │ │
│ │ │ │
│ │ │ │ -
1087 template<
class T,
class A>
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
1091 using real_type =
typename FieldTraits<field_type>::real_type;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
Helper functions for determining the vector/matrix block level.
│ │ │ │ -
This file implements a vector space as a tensor product of a given vector space. The number of compon...
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
1088 template<
typename V>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1090 const std::integral_constant<int,0>&)
│ │ │ │ +
│ │ │ │ +
1092 using namespace MatrixMarketImpl;
│ │ │ │ +
1093 typedef typename V::field_type field_type;
│ │ │ │ +
│ │ │ │ +
1095 ostr<<
countEntries(vector)<<
" "<<Simd::lanes<field_type>()<<std::endl;
│ │ │ │ +
1096 const int isnumeric = mm_numeric_type<Simd::Scalar<V>>::is_numeric;
│ │ │ │ +
1097 for(
size_t l=0;l<Simd::lanes<field_type>(); ++l){
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1103 template<
typename M>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1106 const std::integral_constant<int,1>&)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1112 typedef typename M::const_iterator riterator;
│ │ │ │ +
1113 typedef typename M::ConstColIterator citerator;
│ │ │ │ +
1114 for(riterator row=matrix.begin(); row != matrix.end(); ++row)
│ │ │ │ +
1115 for(citerator
col = row->begin();
col != row->end(); ++
col)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1125 template<
typename M>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1129 using namespace MatrixMarketImpl;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1132 mm_header_printer<M>::print(ostr);
│ │ │ │ +
1133 mm_block_structure_header<M>::print(ostr,matrix);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1150 template<
typename M>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1152 std::string filename,
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1156 std::string rfilename;
│ │ │ │ +
│ │ │ │ +
1158 if (extension !=
"") {
│ │ │ │ +
1159 rfilename = pureFilename + extension;
│ │ │ │ +
1160 file.open(rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
1162 DUNE_THROW(IOError,
"Could not open file for storage: " << rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1166 rfilename = pureFilename +
".mm";
│ │ │ │ +
1167 file.open(rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
1169 DUNE_THROW(IOError,
"Could not open file for storage: " << rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1172 file.setf(std::ios::scientific,std::ios::floatfield);
│ │ │ │ +
│ │ │ │ +
1174 file.precision(prec);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1194 template<
typename M,
typename G,
typename L>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1196 std::string filename,
│ │ │ │ +
│ │ │ │ +
1198 bool storeIndices=
true,
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1205 std::string rfilename;
│ │ │ │ +
│ │ │ │ +
1207 if (extension !=
"") {
│ │ │ │ +
1208 rfilename = pureFilename +
"_" + std::to_string(rank) + extension;
│ │ │ │ +
1209 file.open(rfilename.c_str());
│ │ │ │ +
1210 dverb<< rfilename <<std::endl;
│ │ │ │ +
│ │ │ │ +
1212 DUNE_THROW(IOError,
"Could not open file for storage: " << rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1216 rfilename = pureFilename +
"_" + std::to_string(rank) +
".mm";
│ │ │ │ +
1217 file.open(rfilename.c_str());
│ │ │ │ +
1218 dverb<< rfilename <<std::endl;
│ │ │ │ +
│ │ │ │ +
1220 DUNE_THROW(IOError,
"Could not open file for storage: " << rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
1222 file.setf(std::ios::scientific,std::ios::floatfield);
│ │ │ │ +
│ │ │ │ +
1224 file.precision(prec);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1232 rfilename = pureFilename +
"_" + std::to_string(rank) +
".idx";
│ │ │ │ +
1233 file.open(rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
1235 DUNE_THROW(IOError,
"Could not open file for storage: " << rfilename.c_str());
│ │ │ │ +
1236 file.setf(std::ios::scientific,std::ios::floatfield);
│ │ │ │ +
│ │ │ │ +
1238 typedef typename IndexSet::const_iterator Iterator;
│ │ │ │ +
1239 for(Iterator iter = comm.
indexSet().begin();
│ │ │ │ +
1240 iter != comm.
indexSet().end(); ++iter) {
│ │ │ │ +
1241 file << iter->global()<<
" "<<(std::size_t)iter->local()<<
" "
│ │ │ │ +
1242 <<(int)iter->local().attribute()<<
" "<<(int)iter->local().isPublic()<<std::endl;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1245 file<<
"neighbours:";
│ │ │ │ +
1246 const std::set<int>& neighbours=comm.
remoteIndices().getNeighbours();
│ │ │ │ +
1247 typedef std::set<int>::const_iterator SIter;
│ │ │ │ +
1248 for(SIter neighbour=neighbours.begin(); neighbour != neighbours.end(); ++neighbour) {
│ │ │ │ +
1249 file<<
" "<< *neighbour;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1268 template<
typename M,
typename G,
typename L>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1270 const std::string& filename,
│ │ │ │ +
│ │ │ │ +
1272 bool readIndices=
true)
│ │ │ │ +
│ │ │ │ +
1274 using namespace MatrixMarketImpl;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1277 typedef typename LocalIndexT::Attribute Attribute;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1282 std::string rfilename;
│ │ │ │ +
│ │ │ │ +
1284 if (extension !=
"") {
│ │ │ │ +
1285 rfilename = pureFilename +
"_" + std::to_string(rank) + extension;
│ │ │ │ +
1286 file.open(rfilename.c_str(), std::ios::in);
│ │ │ │ +
1287 dverb<< rfilename <<std::endl;
│ │ │ │ +
│ │ │ │ +
1289 DUNE_THROW(IOError,
"Could not open file: " << rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1293 rfilename = pureFilename +
"_" + std::to_string(rank) +
".mm";
│ │ │ │ +
1294 file.open(rfilename.c_str(), std::ios::in);
│ │ │ │ +
│ │ │ │ +
1296 rfilename = pureFilename +
"_" + std::to_string(rank) +
".mtx";
│ │ │ │ +
1297 file.open(rfilename.c_str(), std::ios::in);
│ │ │ │ +
1298 dverb<< rfilename <<std::endl;
│ │ │ │ +
│ │ │ │ +
1300 DUNE_THROW(IOError,
"Could not open file: " << rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1311 IndexSet& pis=comm.pis;
│ │ │ │ +
1312 rfilename = pureFilename +
"_" + std::to_string(rank) +
".idx";
│ │ │ │ +
1313 file.open(rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
1315 DUNE_THROW(IOError,
"Could not open file: " << rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
1317 DUNE_THROW(InvalidIndexSetState,
"Index set is not empty!");
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1320 while(!file.eof() && file.peek()!=
'n') {
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1329 pis.add(g,LocalIndexT(l,Attribute(c),b));
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1337 if(s!=
"neighbours:")
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1340 while(!file.eof()) {
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1346 comm.ri.setNeighbours(nb);
│ │ │ │ +
│ │ │ │ +
1348 comm.ri.template rebuild<false>();
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1363 template<
typename M>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1365 const std::string& filename)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1368 std::string rfilename;
│ │ │ │ +
│ │ │ │ +
1370 if (extension !=
"") {
│ │ │ │ +
1371 rfilename = pureFilename + extension;
│ │ │ │ +
1372 file.open(rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
1374 DUNE_THROW(IOError,
"Could not open file: " << rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
1378 rfilename = pureFilename +
".mm";
│ │ │ │ +
1379 file.open(rfilename.c_str(), std::ios::in);
│ │ │ │ +
│ │ │ │ +
1381 rfilename = pureFilename +
".mtx";
│ │ │ │ +
1382 file.open(rfilename.c_str(), std::ios::in);
│ │ │ │ +
│ │ │ │ +
1384 DUNE_THROW(IOError,
"Could not open file: " << rfilename.c_str());
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
This file implements a vector space as a tensor product of a given vector space. The number of compon...
│ │ │ │ +
Implementation of the BCRSMatrix class.
│ │ │ │ +
Some handy generic functions for ISTL matrices.
│ │ │ │ +
Classes providing communication interfaces for overlapping Schwarz methods.
│ │ │ │ +
Col col
Definition matrixmatrix.hh:351
│ │ │ │ +
auto countNonZeros(const M &, typename std::enable_if_t< Dune::IsNumber< M >::value > *sfinae=nullptr)
Get the number of nonzero fields in the matrix.
Definition matrixutils.hh:119
│ │ │ │ +
void readMatrixMarket(Dune::BlockVector< T, A > &vector, std::istream &istr)
Reads a BlockVector from a matrix market file.
Definition matrixmarket.hh:930
│ │ │ │ +
void storeMatrixMarket(const M &matrix, std::string filename, int prec=default_precision)
Stores a parallel matrix/vector in matrix market format in a file.
Definition matrixmarket.hh:1151
│ │ │ │ +
void loadMatrixMarket(M &matrix, const std::string &filename, OwnerOverlapCopyCommunication< G, L > &comm, bool readIndices=true)
Load a parallel matrix/vector stored in matrix market format.
Definition matrixmarket.hh:1269
│ │ │ │ +
std::size_t countEntries(const BlockVector< T, A > &vector)
Definition matrixmarket.hh:1076
│ │ │ │ +
void writeMatrixMarket(const V &vector, std::ostream &ostr, const std::integral_constant< int, 0 > &)
Definition matrixmarket.hh:1089
│ │ │ │ +
void mm_print_vector_entry(const V &entry, std::ostream &ostr, const std::integral_constant< int, 1 > &, size_t lane)
Definition matrixmarket.hh:1049
│ │ │ │ +
static const int default_precision
Definition matrixmarket.hh:1138
│ │ │ │ +
void mm_read_vector_entries(Dune::BlockVector< T, A > &vector, std::size_t size, std::istream &istr, size_t lane)
Definition matrixmarket.hh:900
│ │ │ │ +
void mm_read_header(std::size_t &rows, std::size_t &cols, MatrixMarketImpl::MMHeader &header, std::istream &istr, bool isVector)
Definition matrixmarket.hh:871
│ │ │ │ +
void mm_print_entry(const B &entry, std::size_t rowidx, std::size_t colidx, std::ostream &ostr)
Definition matrixmarket.hh:1030
│ │ │ │ +
│ │ │ │
Definition allocator.hh:11
│ │ │ │ +
std::tuple< std::size_t, std::size_t, std::size_t > calculateNNZ(std::size_t rows, std::size_t cols, std::size_t entries, const MMHeader &header)
Definition matrixmarket.hh:547
│ │ │ │ +
bool operator<(const IndexData< T > &i1, const IndexData< T > &i2)
LessThan operator.
Definition matrixmarket.hh:630
│ │ │ │ +
LineType
Definition matrixmarket.hh:296
│ │ │ │ +
@ DATA
Definition matrixmarket.hh:296
│ │ │ │ +
@ MM_HEADER
Definition matrixmarket.hh:296
│ │ │ │ +
@ MM_ISTLSTRUCT
Definition matrixmarket.hh:296
│ │ │ │ +
bool readMatrixMarketBanner(std::istream &file, MMHeader &mmHeader)
Definition matrixmarket.hh:353
│ │ │ │ +
void readSparseEntries(Dune::BCRSMatrix< T, A > &matrix, std::istream &file, std::size_t entries, const MMHeader &mmHeader, const D &)
Definition matrixmarket.hh:765
│ │ │ │ +
MM_TYPE
Definition matrixmarket.hh:299
│ │ │ │ +
@ array_type
Definition matrixmarket.hh:299
│ │ │ │ +
@ coordinate_type
Definition matrixmarket.hh:299
│ │ │ │ +
@ unknown_type
Definition matrixmarket.hh:299
│ │ │ │ +
std::istream & operator>>(std::istream &is, NumericWrapper< T > &num)
Definition matrixmarket.hh:614
│ │ │ │ +
void skipComments(std::istream &file)
Definition matrixmarket.hh:339
│ │ │ │ +
bool lineFeed(std::istream &file)
Definition matrixmarket.hh:315
│ │ │ │ +
@ MM_MAX_LINE_LENGTH
Definition matrixmarket.hh:297
│ │ │ │ +
MM_STRUCTURE
Definition matrixmarket.hh:303
│ │ │ │ +
@ skew_symmetric
Definition matrixmarket.hh:303
│ │ │ │ +
@ general
Definition matrixmarket.hh:303
│ │ │ │ +
@ hermitian
Definition matrixmarket.hh:303
│ │ │ │ +
@ unknown_structure
Definition matrixmarket.hh:303
│ │ │ │ +
@ symmetric
Definition matrixmarket.hh:303
│ │ │ │ +
MM_CTYPE
Definition matrixmarket.hh:301
│ │ │ │ +
@ unknown_ctype
Definition matrixmarket.hh:301
│ │ │ │ +
@ pattern
Definition matrixmarket.hh:301
│ │ │ │ +
@ complex_type
Definition matrixmarket.hh:301
│ │ │ │ +
@ double_type
Definition matrixmarket.hh:301
│ │ │ │ +
@ integer_type
Definition matrixmarket.hh:301
│ │ │ │ +
std::enable_if_t<!is_complex< T >::value, T > conj(const T &r)
Definition matrixmarket.hh:733
│ │ │ │ +
std::tuple< std::string, std::string > splitFilename(const std::string &filename)
Definition matrixmarket.hh:852
│ │ │ │ +
A sparse block matrix with compressed row storage.
Definition bcrsmatrix.hh:466
│ │ │ │ +
Iterator begin()
Get iterator to first row.
Definition bcrsmatrix.hh:671
│ │ │ │ +
Iterator end()
Get iterator to one beyond last row.
Definition bcrsmatrix.hh:677
│ │ │ │ +
CreateIterator createend()
get create iterator pointing to one after the last block
Definition bcrsmatrix.hh:1100
│ │ │ │ +
size_type M() const
number of columns (counted in blocks)
Definition bcrsmatrix.hh:2007
│ │ │ │ +
CreateIterator createbegin()
get initial create iterator
Definition bcrsmatrix.hh:1094
│ │ │ │ +
size_type N() const
number of rows (counted in blocks)
Definition bcrsmatrix.hh:2001
│ │ │ │ +
void setBuildMode(BuildMode bm)
Sets the build mode of the matrix.
Definition bcrsmatrix.hh:830
│ │ │ │ +
void setSize(size_type rows, size_type columns, size_type nnz=0)
Set the size of the matrix.
Definition bcrsmatrix.hh:858
│ │ │ │
A vector of blocks with memory management.
Definition bvector.hh:392
│ │ │ │ -
derive error class from the base class in common
Definition istlexception.hh:19
│ │ │ │ -
A Vector of blocks with different blocksizes.
Definition matrix.hh:44
│ │ │ │ -
Imp::BlockVectorWindow< B, A > window_type
Definition matrix.hh:70
│ │ │ │ -
BlockVector< B, A > block_type
Same as value_type, here for historical reasons.
Definition matrix.hh:67
│ │ │ │ -
DenseMatrixBase & operator=(const DenseMatrixBase &a)
assignment
Definition matrix.hh:182
│ │ │ │ -
typename Imp::BlockTraits< B >::field_type field_type
export the type representing the field
Definition matrix.hh:50
│ │ │ │ -
Iterator beforeBegin() const
Definition matrix.hh:385
│ │ │ │ -
void resize(size_type rows, size_type columns)
same effect as constructor with same argument
Definition matrix.hh:153
│ │ │ │ -
const window_type const_reference
Definition matrix.hh:74
│ │ │ │ -
DenseMatrixBase(size_type rows, size_type columns)
Definition matrix.hh:95
│ │ │ │ -
Iterator end()
end Iterator
Definition matrix.hh:371
│ │ │ │ -
reference operator[](size_type i)
random access to blocks
Definition matrix.hh:241
│ │ │ │ -
BlockVector< B, A > value_type
Type of the elements of the outer vector, i.e., dynamic vectors of B.
Definition matrix.hh:63
│ │ │ │ -
Iterator find(size_type i)
random access returning iterator (end if not contained)
Definition matrix.hh:391
│ │ │ │ -
size_type N() const
number of blocks in the vector (are of variable size here)
Definition matrix.hh:539
│ │ │ │ -
ConstIterator beforeEnd() const
Definition matrix.hh:525
│ │ │ │ -
ConstIterator end() const
end ConstIterator
Definition matrix.hh:518
│ │ │ │ -
window_type reference
Definition matrix.hh:72
│ │ │ │ -
ConstIterator rend() const
end ConstIterator
Definition matrix.hh:531
│ │ │ │ -
Iterator beforeEnd()
Definition matrix.hh:378
│ │ │ │ -
ConstIterator begin() const
begin ConstIterator
Definition matrix.hh:512
│ │ │ │ -
A allocator_type
export the allocator type
Definition matrix.hh:53
│ │ │ │ -
DenseMatrixBase()
Definition matrix.hh:82
│ │ │ │ -
ConstIterator find(size_type i) const
random access returning iterator (end if not contained)
Definition matrix.hh:397
│ │ │ │ -
A::size_type size_type
The size type for the index access.
Definition matrix.hh:56
│ │ │ │ -
~DenseMatrixBase()
free dynamic memory
Definition matrix.hh:142
│ │ │ │ -
Iterator begin()
begin Iterator
Definition matrix.hh:365
│ │ │ │ -
DenseMatrixBase(const DenseMatrixBase &a)
copy constructor, has copy semantics
Definition matrix.hh:116
│ │ │ │ -
Iterator class for sequential access.
Definition matrix.hh:263
│ │ │ │ -
Iterator & operator--()
prefix decrement
Definition matrix.hh:308
│ │ │ │ -
size_type index() const
Definition matrix.hh:352
│ │ │ │ -
Iterator(Iterator &other)=default
│ │ │ │ -
Iterator(Iterator &&other)=default
│ │ │ │ -
bool operator!=(const Iterator &it) const
inequality
Definition matrix.hh:322
│ │ │ │ -
Iterator & operator=(Iterator &&other)
Move assignment.
Definition matrix.hh:282
│ │ │ │ -
Iterator & operator++()
prefix increment
Definition matrix.hh:300
│ │ │ │ -
Iterator()
constructor, no arguments
Definition matrix.hh:266
│ │ │ │ -
window_type & operator*() const
dereferencing
Definition matrix.hh:340
│ │ │ │ -
bool operator==(const Iterator &it) const
equality
Definition matrix.hh:316
│ │ │ │ -
Iterator & operator=(Iterator &other)
Copy assignment.
Definition matrix.hh:291
│ │ │ │ -
Iterator(B *data, size_type columns, size_type _i)
constructor
Definition matrix.hh:276
│ │ │ │ -
window_type * operator->() const
arrow
Definition matrix.hh:346
│ │ │ │ -
ConstIterator class for sequential access.
Definition matrix.hh:404
│ │ │ │ -
const window_type * operator->() const
arrow
Definition matrix.hh:487
│ │ │ │ -
const window_type & operator*() const
dereferencing
Definition matrix.hh:481
│ │ │ │ -
ConstIterator & operator++()
prefix increment
Definition matrix.hh:441
│ │ │ │ -
ConstIterator(const B *data, size_type columns, size_type _i)
constructor from pointer
Definition matrix.hh:414
│ │ │ │ -
ConstIterator & operator--()
prefix decrement
Definition matrix.hh:449
│ │ │ │ -
ConstIterator(const Iterator &it)
constructor from non_const iterator
Definition matrix.hh:420
│ │ │ │ -
bool operator!=(const ConstIterator &it) const
inequality
Definition matrix.hh:463
│ │ │ │ -
ConstIterator()
constructor
Definition matrix.hh:407
│ │ │ │ -
bool operator==(const ConstIterator &it) const
equality
Definition matrix.hh:457
│ │ │ │ -
size_type index() const
Definition matrix.hh:493
│ │ │ │ -
ConstIterator & operator=(Iterator &&other)
Definition matrix.hh:424
│ │ │ │ -
ConstIterator & operator=(Iterator &other)
Definition matrix.hh:432
│ │ │ │ +
void resize(size_type size)
Resize the vector.
Definition bvector.hh:496
│ │ │ │ +
typename Imp::BlockTraits< B >::field_type field_type
export the type representing the field
Definition bvector.hh:398
│ │ │ │
A generic dynamic dense matrix.
Definition matrix.hh:561
│ │ │ │ -
size_type cols_
Number of columns of the matrix.
Definition matrix.hh:1084
│ │ │ │ -
FieldTraits< ft >::real_type infinity_norm() const
infinity norm (row sum norm, how to generalize for blocks?)
Definition matrix.hh:988
│ │ │ │ -
A allocator_type
Export the allocator.
Definition matrix.hh:571
│ │ │ │ -
FieldTraits< ft >::real_type infinity_norm_real() const
simplified infinity norm (uses Manhattan norm for complex values)
Definition matrix.hh:1007
│ │ │ │ -
void usmhv(const field_type &alpha, const X &x, Y &y) const
y += alpha A^H x
Definition matrix.hh:952
│ │ │ │ -
void usmv(const field_type &alpha, const X &x, Y &y) const
Definition matrix.hh:850
│ │ │ │ -
A::size_type size_type
Type for indices and sizes.
Definition matrix.hh:577
│ │ │ │ -
MatrixImp::DenseMatrixBase< T, A > data_
Abuse DenseMatrixBase as an engine for a 2d array ISTL-style.
Definition matrix.hh:1077
│ │ │ │ -
Matrix transpose() const
Return the transpose of the matrix.
Definition matrix.hh:741
│ │ │ │ -
void mtv(const X &x, Y &y) const
y = A^T x
Definition matrix.hh:803
│ │ │ │ -
void umv(const X &x, Y &y) const
y += A x
Definition matrix.hh:816
│ │ │ │ -
void mv(const X &x, Y &y) const
y = A x
Definition matrix.hh:784
│ │ │ │ -
MatrixImp::DenseMatrixBase< T, A >::ConstIterator ConstRowIterator
Const iterator over the matrix rows.
Definition matrix.hh:586
│ │ │ │ -
void setSize(size_type rows, size_type cols)
Change the matrix size.
Definition matrix.hh:604
│ │ │ │ -
RowIterator beforeBegin()
Definition matrix.hh:630
│ │ │ │ -
RowIterator beforeEnd()
Definition matrix.hh:623
│ │ │ │ -
Matrix()
Create empty matrix.
Definition matrix.hh:592
│ │ │ │ -
Matrix & operator-=(const Matrix &b)
Subtract the entries of another matrix from this one.
Definition matrix.hh:731
│ │ │ │ -
FieldTraits< field_type >::real_type frobenius_norm2() const
square of frobenius norm, need for block recursion
Definition matrix.hh:976
│ │ │ │ -
row_type::iterator ColIterator
Iterator for the entries of each row.
Definition matrix.hh:583
│ │ │ │ -
ConstRowIterator beforeEnd() const
Definition matrix.hh:649
│ │ │ │ -
Matrix & operator=(const field_type &t)
Assignment from scalar.
Definition matrix.hh:662
│ │ │ │ -
RowIterator end()
Get iterator to one beyond last row.
Definition matrix.hh:616
│ │ │ │ -
const row_type operator[](size_type row) const
The const index operator.
Definition matrix.hh:680
│ │ │ │ -
ConstRowIterator end() const
Get const iterator to one beyond last row.
Definition matrix.hh:642
│ │ │ │ -
friend Y operator*(const Matrix< T > &m, const X &vec)
Generic matrix-vector multiplication.
Definition matrix.hh:766
│ │ │ │ -
Matrix< T > & operator*=(const field_type &scalar)
Multiplication with a scalar.
Definition matrix.hh:701
│ │ │ │ -
row_type operator[](size_type row)
The index operator.
Definition matrix.hh:669
│ │ │ │ -
void mmv(const X &x, Y &y) const
y -= A x
Definition matrix.hh:833
│ │ │ │ -
Matrix & operator+=(const Matrix &b)
Add the entries of another matrix to this one.
Definition matrix.hh:717
│ │ │ │ -
ConstRowIterator begin() const
Get const iterator to first row.
Definition matrix.hh:636
│ │ │ │ -
void mmhv(const X &x, Y &y) const
y -= A^H x
Definition matrix.hh:935
│ │ │ │ -
RowIterator begin()
Get iterator to first row.
Definition matrix.hh:610
│ │ │ │ -
typename Imp::BlockTraits< T >::field_type field_type
Export the type representing the underlying field.
Definition matrix.hh:565
│ │ │ │ -
row_type::const_iterator ConstColIterator
Const iterator for the entries of each row.
Definition matrix.hh:589
│ │ │ │ -
size_type M() const
Return the number of columns.
Definition matrix.hh:696
│ │ │ │ -
T block_type
Export the type representing the components.
Definition matrix.hh:568
│ │ │ │ -
bool exists(size_type i, size_type j) const
return true if (i,j) is in pattern
Definition matrix.hh:1064
│ │ │ │ -
Matrix< T > & operator/=(const field_type &scalar)
Division by a scalar.
Definition matrix.hh:707
│ │ │ │ -
friend Matrix< T > operator*(const Matrix< T > &m1, const Matrix< T > &m2)
Generic matrix multiplication.
Definition matrix.hh:751
│ │ │ │ -
void umtv(const X &x, Y &y) const
y += A^T x
Definition matrix.hh:867
│ │ │ │ -
void mmtv(const X &x, Y &y) const
y -= A^T x
Definition matrix.hh:884
│ │ │ │ -
MatrixImp::DenseMatrixBase< T, A >::window_type row_type
The type implementing a matrix row.
Definition matrix.hh:574
│ │ │ │ -
FieldTraits< field_type >::real_type frobenius_norm() const
frobenius norm: sqrt(sum over squared values of entries)
Definition matrix.hh:970
│ │ │ │ -
void usmtv(const field_type &alpha, const X &x, Y &y) const
y += alpha A^T x
Definition matrix.hh:901
│ │ │ │ -
void umhv(const X &x, Y &y) const
y += A^H x
Definition matrix.hh:918
│ │ │ │ -
size_type N() const
Return the number of rows.
Definition matrix.hh:691
│ │ │ │ -
ConstRowIterator beforeBegin() const
Definition matrix.hh:656
│ │ │ │ -
Matrix(size_type rows, size_type cols)
Create uninitialized matrix of size rows x cols.
Definition matrix.hh:597
│ │ │ │ -
MatrixImp::DenseMatrixBase< T, A >::Iterator RowIterator
Iterator over the matrix rows.
Definition matrix.hh:580
│ │ │ │ -
typename Matrix< T, A >::field_type field_type
Definition matrix.hh:1090
│ │ │ │ -
typename FieldTraits< field_type >::real_type real_type
Definition matrix.hh:1091
│ │ │ │ +
Helper metaprogram to get the matrix market string representation of the numeric type.
Definition matrixmarket.hh:76
│ │ │ │ +
@ is_numeric
Whether T is a supported numeric type.
Definition matrixmarket.hh:81
│ │ │ │ +
static std::string str()
Definition matrixmarket.hh:95
│ │ │ │ +
static std::string str()
Definition matrixmarket.hh:111
│ │ │ │ +
static std::string str()
Definition matrixmarket.hh:127
│ │ │ │ +
static std::string str()
Definition matrixmarket.hh:143
│ │ │ │ +
static std::string str()
Definition matrixmarket.hh:159
│ │ │ │ +
Meta program to write the correct Matrix Market header.
Definition matrixmarket.hh:174
│ │ │ │ +
static void print(std::ostream &os)
Definition matrixmarket.hh:179
│ │ │ │ +
static void print(std::ostream &os)
Definition matrixmarket.hh:189
│ │ │ │ +
static void print(std::ostream &os)
Definition matrixmarket.hh:199
│ │ │ │ +
static void print(std::ostream &os)
Definition matrixmarket.hh:209
│ │ │ │ +
Metaprogram for writing the ISTL block structure header.
Definition matrixmarket.hh:225
│ │ │ │ +
static void print(std::ostream &os, const M &)
Definition matrixmarket.hh:233
│ │ │ │ +
BlockVector< T, A > M
Definition matrixmarket.hh:230
│ │ │ │ +
BlockVector< FieldVector< T, i >, A > M
Definition matrixmarket.hh:243
│ │ │ │ +
static void print(std::ostream &os, const M &)
Definition matrixmarket.hh:245
│ │ │ │ +
BCRSMatrix< T, A > M
Definition matrixmarket.hh:255
│ │ │ │ +
static void print(std::ostream &os, const M &)
Definition matrixmarket.hh:258
│ │ │ │ +
BCRSMatrix< FieldMatrix< T, i, j >, A > M
Definition matrixmarket.hh:268
│ │ │ │ +
static void print(std::ostream &os, const M &)
Definition matrixmarket.hh:270
│ │ │ │ +
static void print(std::ostream &os, const M &m)
Definition matrixmarket.hh:283
│ │ │ │ +
FieldMatrix< T, i, j > M
Definition matrixmarket.hh:281
│ │ │ │ +
static void print(std::ostream &os, const M &m)
Definition matrixmarket.hh:292
│ │ │ │ +
FieldVector< T, i > M
Definition matrixmarket.hh:290
│ │ │ │ +
Definition matrixmarket.hh:306
│ │ │ │ +
MM_STRUCTURE structure
Definition matrixmarket.hh:312
│ │ │ │ +
MM_TYPE type
Definition matrixmarket.hh:310
│ │ │ │ +
MMHeader()
Definition matrixmarket.hh:307
│ │ │ │ +
MM_CTYPE ctype
Definition matrixmarket.hh:311
│ │ │ │ +
Definition matrixmarket.hh:578
│ │ │ │ +
std::size_t index
Definition matrixmarket.hh:579
│ │ │ │ +
a wrapper class of numeric values.
Definition matrixmarket.hh:595
│ │ │ │ +
T number
Definition matrixmarket.hh:596
│ │ │ │ +
Utility class for marking the pattern type of the MatrixMarket matrices.
Definition matrixmarket.hh:607
│ │ │ │ +
Functor to the data values of the matrix.
Definition matrixmarket.hh:677
│ │ │ │ +
void operator()(const std::vector< std::set< IndexData< D > > > &rows, BCRSMatrix< T > &matrix)
Sets the matrix values.
Definition matrixmarket.hh:684
│ │ │ │ +
void operator()(const std::vector< std::set< IndexData< D > > > &rows, BCRSMatrix< FieldMatrix< T, brows, bcols > > &matrix)
Sets the matrix values.
Definition matrixmarket.hh:702
│ │ │ │ +
void operator()(const std::vector< std::set< IndexData< PatternDummy > > > &rows, M &matrix)
Definition matrixmarket.hh:723
│ │ │ │ +
Definition matrixmarket.hh:728
│ │ │ │ +
Definition matrixmarket.hh:744
│ │ │ │ +
Definition matrixmarket.hh:868
│ │ │ │ +
Definition matrixutils.hh:27
│ │ │ │ +
Test whether a type is an ISTL Matrix.
Definition matrixutils.hh:504
│ │ │ │ +
A class setting up standard communication for a two-valued attribute set with owner/overlap/copy sema...
Definition owneroverlapcopy.hh:174
│ │ │ │ +
const ParallelIndexSet & indexSet() const
Get the underlying parallel index set.
Definition owneroverlapcopy.hh:462
│ │ │ │ +
const Communication< MPI_Comm > & communicator() const
Definition owneroverlapcopy.hh:299
│ │ │ │ +
const RemoteIndices & remoteIndices() const
Get the underlying remote indices.
Definition owneroverlapcopy.hh:471
│ │ │ │ +
Dune::ParallelIndexSet< GlobalIdType, LI, 512 > ParallelIndexSet
The type of the parallel index set.
Definition owneroverlapcopy.hh:449
│ │ │ │
1.9.8
│ │ │ │
│ │ │ │
│ │ │ │ ├── html2text {}
│ │ │ │ │ @@ -1,19 +1,35 @@
│ │ │ │ │ dune-istl 2.10
│ │ │ │ │ Loading...
│ │ │ │ │ Searching...
│ │ │ │ │ No Matches
│ │ │ │ │ * _d_u_n_e
│ │ │ │ │ * _i_s_t_l
│ │ │ │ │ -_N_a_m_e_s_p_a_c_e_s
│ │ │ │ │ -bccsmatrix.hh File Reference
│ │ │ │ │ -#include