![]() |
Reference documentation for deal.II version 8.1.0
|
#include <table.h>
Public Member Functions | |
TableBase () | |
TableBase (const TableIndices< N > &sizes) | |
TableBase (const TableBase< N, T > &src) | |
template<typename T2 > | |
TableBase (const TableBase< N, T2 > &src) | |
~TableBase () | |
TableBase< N, T > & | operator= (const TableBase< N, T > &src) |
template<typename T2 > | |
TableBase< N, T > & | operator= (const TableBase< N, T2 > &src) |
bool | operator== (const TableBase< N, T > &T2) const |
void | reset_values () |
void | reinit (const TableIndices< N > &new_size, const bool fast=false) |
unsigned int | size (const unsigned int i) const |
const TableIndices< N > & | size () const |
unsigned int | n_elements () const |
bool | empty () const |
template<typename T2 > | |
void | fill (const T2 *entries) |
void | fill (const T &value) |
std::vector< T >::reference | operator() (const TableIndices< N > &indices) |
std::vector< T >::const_reference | operator() (const TableIndices< N > &indices) const |
void | swap (TableBase< N, T > &v) |
std::size_t | memory_consumption () const |
template<class Archive > | |
void | serialize (Archive &ar, const unsigned int version) |
![]() | |
Subscriptor () | |
Subscriptor (const Subscriptor &) | |
virtual | ~Subscriptor () |
Subscriptor & | operator= (const Subscriptor &) |
void | subscribe (const char *identifier=0) const |
void | unsubscribe (const char *identifier=0) const |
unsigned int | n_subscriptions () const |
void | list_subscribers () const |
DeclException3 (ExcInUse, int, char *, std::string &,<< "Object of class "<< arg2<< " is still used by "<< arg1<< " other objects.\n"<< "(Additional information: "<< arg3<< ")\n"<< "Note the entry in the Frequently Asked Questions of "<< "deal.II (linked to from http://www.dealii.org/) for "<< "more information on what this error means.") | |
DeclException2 (ExcNoSubscriber, char *, char *,<< "No subscriber with identifier \""<< arg2<< "\" did subscribe to this object of class "<< arg1) | |
template<class Archive > | |
void | serialize (Archive &ar, const unsigned int version) |
Protected Member Functions | |
unsigned int | position (const TableIndices< N > &indices) const |
std::vector< T >::reference | el (const TableIndices< N > &indices) |
std::vector< T >::const_reference | el (const TableIndices< N > &indices) const |
std::vector< T >::const_pointer | data () const DEAL_II_DEPRECATED |
Protected Attributes | |
std::vector< T > | values |
TableIndices< N > | table_size |
Friends | |
template<int , typename > | |
class | TableBase |
General class holding an array of objects of templated type in multiple dimensions. If the template parameter indicating the number of dimensions is one, then this is more or less a vector, if it is two then it is a matrix, and so on.
Previously, this data type was emulated in this library by constructs like std::vector<std::vector<T>>
, or even higher nested constructs. However, this has the disadvantage that it is hard to initialize, and most importantly that it is very inefficient if all rows have the same size (which is the usual case), since then the memory for each row is allocated independently, both wasting time and memory. This can be made more efficient by allocating only one chunk of memory for the entire object.
Therefore, this data type was invented. Its implementation is rather straightforward, with two exceptions. The first thing to think about is how to pass the size in each of the coordinate directions to the object; this is done using the TableIndices class. Second, how to access the individual elements. The basic problem here is that we would like to make the number of arguments to be passed to the constructor as well as the access functions dependent on the template parameter N
indicating the number of dimensions. Of course, this is not possible.
The way out of the first problem (and partly the second one as well) is to have a common base class TableBase and a derived class for each value of N
. This derived class has a constructor with the correct number of arguments, namely N
. These then transform their arguments into the data type the base class (this class in fact) uses in the constructor as well as in element access through operator() functions.
The second problem is that we would like to allow access through a sequence of operator[]
calls. This mostly because, as said, this class is a replacement for previous use of nested std::vector
objects, where we had to use the operator[]
access function recursively until we were at the innermost object. Emulating this behavior without losing the ability to do index checks, and in particular without losing performance is possible but nontrivial, and done in the TableBaseAccessors namespace.
In some way, this class is similar to the Tensor class, in that it templatizes on the number of dimensions. However, there are two major differences. The first is that the Tensor class stores only numeric values (as double
s), while the Table class stores arbitrary objects. The second is that the Tensor class has fixed dimensions, also given as a template argument, while this class can handle arbitrary dimensions, which may also be different between different indices.
This has two consequences. First, since the size is not known at compile time, it has to do explicit memory allocating. Second, the layout of individual elements is not known at compile time, so access is slower than for the Tensor class where the number of elements are their location is known at compile time and the compiler can optimize with this knowledge (for example when unrolling loops). On the other hand, this class is of course more flexible, for example when you want a two-dimensional table with the number of rows equal to the number of degrees of freedom on a cell, and the number of columns equal to the number of quadrature points. Both numbers may only be known at run-time, so a flexible table is needed here. Furthermore, you may want to store, say, the gradients of shape functions, so the data type is not a single scalar value, but a tensor itself.
Default constructor. Set all dimensions to zero.
TableBase< N, T >::TableBase | ( | const TableIndices< N > & | sizes | ) |
Constructor. Initialize the array with the given dimensions in each index component.
Copy constructor. Performs a deep copy.
TableBase< N, T >::TableBase | ( | const TableBase< N, T2 > & | src | ) |
Copy constructor. Performs a deep copy from a table object storing some other data type.
TableBase<N,T>& TableBase< N, T >::operator= | ( | const TableBase< N, T > & | src | ) |
Assignment operator. Copy all elements of src
into the matrix. The size is adjusted if needed.
We can't use the other, templatized version since if we don't declare this one, the compiler will happily generate a predefined copy operator which is not what we want.
TableBase<N,T>& TableBase< N, T >::operator= | ( | const TableBase< N, T2 > & | src | ) |
Copy operator. Copy all elements of src
into the array. The size is adjusted if needed.
This function requires that the type T2
is convertible to T
.
bool TableBase< N, T >::operator== | ( | const TableBase< N, T > & | T2 | ) | const |
Test for equality of two tables.
void TableBase< N, T >::reset_values | ( | ) |
Set all entries to their default value (i.e. copy them over with default constructed objects). Do not change the size of the table, though.
void TableBase< N, T >::reinit | ( | const TableIndices< N > & | new_size, |
const bool | fast = false |
||
) |
Set the dimensions of this object to the sizes given in the argument, and newly allocate the required memory. If fast
is set to false
, previous content is deleted, otherwise the memory is not touched.
Size of the table in direction i
.
const TableIndices<N>& TableBase< N, T >::size | ( | ) | const |
Return the sizes of this object in each direction.
Return the number of elements stored in this object, which is the product of the extensions in each dimension.
Return whether the object is empty, i.e. one of the directions is zero. This is equivalent to n_elements()==0
.
void TableBase< N, T >::fill | ( | const T2 * | entries | ) |
Fill array with an array of elements. The input array must be arranged in usual C style, i.e. with the last index running fastest. For two-dimensional tables, this means line by line. No range checking is performed, i.e., it is assumed that the input array entries
contains n_rows()*n_cols() elements, and that the layout refers to the desired shape of this table. The only check we do is that the present array is non-empty.
Note also that the type of the objects of the input array, T2
, must be convertible to the type of the objects of this array.
void TableBase< N, T >::fill | ( | const T & | value | ) |
Fill all table entries with the same value.
std::vector<T>::reference TableBase< N, T >::operator() | ( | const TableIndices< N > & | indices | ) |
Return a read-write reference to the indicated element.
std::vector<T>::const_reference TableBase< N, T >::operator() | ( | const TableIndices< N > & | indices | ) | const |
Return the value of the indicated element as a read-only reference.
We return the requested value as a constant reference rather than by value since this object may hold data types that may be large, and we don't know here whether copying is expensive or not.
Swap the contents of this table and the other table v
. One could do this operation with a temporary variable and copying over the data elements, but this function is significantly more efficient since it only swaps the pointers to the data of the two vectors and therefore does not need to allocate temporary storage and move data around.
This function is analog to the the swap
function of all C++ standard containers. Also, there is a global function swap(u,v)
that simply calls u.swap(v)
, again in analogy to standard functions.
std::size_t TableBase< N, T >::memory_consumption | ( | ) | const |
Determine an estimate for the memory consumption (in bytes) of this object.
void TableBase< N, T >::serialize | ( | Archive & | ar, |
const unsigned int | version | ||
) |
Write or read the data of this object to or from a stream for the purpose of serialization.
|
protected |
Return the position of the indicated element within the array of elements stored one after the other. This function does no index checking.
|
protected |
Return a read-write reference to the indicated element.
This function does no bounds checking and is only to be used internally and in functions already checked.
|
protected |
Return the value of the indicated element as a read-only reference.
This function does no bounds checking and is only to be used internally and in functions already checked.
We return the requested value as a constant reference rather than by value since this object may hold data types that may be large, and we don't know here whether copying is expensive or not.
|
protected |
Direct read-only access to data field. Used by FullMatrix
of the LAC sublibrary (there even with a cast from const), otherwise, keep away!
|
protected |
|
protected |