Public Member Functions | |
| SymmetricTensor () | |
| SymmetricTensor (const Tensor< 2, dim, Number > &t) | |
| SymmetricTensor (const Number(&array)[internal::SymmetricTensorAccessors::StorageType< rank, dim, Number >::n_independent_components]) | |
| SymmetricTensor & | operator= (const SymmetricTensor &) |
| SymmetricTensor & | operator= (const Number d) |
| operator Tensor< rank, dim, Number > () const | |
| bool | operator== (const SymmetricTensor &) const |
| bool | operator!= (const SymmetricTensor &) const |
| SymmetricTensor & | operator+= (const SymmetricTensor &) |
| SymmetricTensor & | operator-= (const SymmetricTensor &) |
| SymmetricTensor & | operator*= (const Number factor) |
| SymmetricTensor & | operator/= (const Number factor) |
| SymmetricTensor | operator+ (const SymmetricTensor &s) const |
| SymmetricTensor | operator- (const SymmetricTensor &s) const |
| SymmetricTensor | operator- () const |
| internal::SymmetricTensorAccessors::double_contraction_result < rank, 2, dim, Number >::type | operator* (const SymmetricTensor< 2, dim, Number > &s) const |
| internal::SymmetricTensorAccessors::double_contraction_result < rank, 4, dim, Number >::type | operator* (const SymmetricTensor< 4, dim, Number > &s) const |
| Number & | operator() (const TableIndices< rank > &indices) |
| Number | operator() (const TableIndices< rank > &indices) const |
| internal::SymmetricTensorAccessors::Accessor < rank, dim, true, rank-1, Number > | operator[] (const unsigned int row) const |
| internal::SymmetricTensorAccessors::Accessor < rank, dim, false, rank-1, Number > | operator[] (const unsigned int row) |
| Number | operator[] (const TableIndices< rank > &indices) const |
| Number & | operator[] (const TableIndices< rank > &indices) |
| Number | norm () const |
| void | clear () |
| template<class Archive > | |
| void | serialize (Archive &ar, const unsigned int version) |
Static Public Member Functions | |
| static unsigned int | component_to_unrolled_index (const TableIndices< rank > &indices) |
| static TableIndices< rank > | unrolled_to_component_indices (const unsigned int i) |
| static std::size_t | memory_consumption () |
Static Public Attributes | |
| static const unsigned int | dimension = dim |
| static const unsigned int | n_independent_components |
Private Types | |
| typedef internal::SymmetricTensorAccessors::StorageType < rank, dim, Number > | base_tensor_descriptor |
| typedef base_tensor_descriptor::base_tensor_type | base_tensor_type |
Private Attributes | |
| base_tensor_type | data |
Friends | |
| class | SymmetricTensor |
| template<int dim2, typename Number2 > | |
| Number2 | trace (const SymmetricTensor< 2, dim2, Number2 > &d) |
| template<int dim2, typename Number2 > | |
| Number2 | determinant (const SymmetricTensor< 2, dim2, Number2 > &t) |
| template<int dim2, typename Number2 > | |
| SymmetricTensor< 2, dim2, Number2 > | deviator (const SymmetricTensor< 2, dim2, Number2 > &t) |
| template<int dim2, typename Number2 > | |
| SymmetricTensor< 2, dim2, Number2 > | unit_symmetric_tensor () |
| template<int dim2, typename Number2 > | |
| SymmetricTensor< 4, dim2, Number2 > | deviator_tensor () |
| template<int dim2, typename Number2 > | |
| SymmetricTensor< 4, dim2, Number2 > | identity_tensor () |
| template<int dim2, typename Number2 > | |
| SymmetricTensor< 4, dim2, Number2 > | invert (const SymmetricTensor< 4, dim2, Number2 > &) |
Related Functions | |
(Note that these are not member functions.) | |
| template<int dim, typename Number > | |
| Number | determinant (const SymmetricTensor< 2, dim, Number > &t) |
| template<int dim, typename Number > | |
| double | third_invariant (const SymmetricTensor< 2, dim, Number > &t) |
| template<int dim, typename Number > | |
| Number | trace (const SymmetricTensor< 2, dim, Number > &d) |
| template<int dim, typename Number > | |
| Number | first_invariant (const SymmetricTensor< 2, dim, Number > &t) |
| template<typename Number > | |
| Number | second_invariant (const SymmetricTensor< 2, 1, Number > &) |
| template<typename Number > | |
| Number | second_invariant (const SymmetricTensor< 2, 2, Number > &t) |
| template<typename Number > | |
| Number | second_invariant (const SymmetricTensor< 2, 3, Number > &t) |
| template<int rank, int dim, typename Number > | |
| SymmetricTensor< rank, dim, Number > | transpose (const SymmetricTensor< rank, dim, Number > &t) |
| template<int dim, typename Number > | |
| SymmetricTensor< 2, dim, Number > | deviator (const SymmetricTensor< 2, dim, Number > &t) |
| template<int dim, typename Number > | |
| SymmetricTensor< 4, dim, Number > | invert (const SymmetricTensor< 4, dim, Number > &t) |
| template<> | |
| SymmetricTensor< 4, 3, double > | invert (const SymmetricTensor< 4, 3, double > &t) |
| template<int dim, typename Number > | |
| SymmetricTensor< 4, dim, Number > | outer_product (const SymmetricTensor< 2, dim, Number > &t1, const SymmetricTensor< 2, dim, Number > &t2) |
| template<typename Number > | |
| SymmetricTensor< 2, 1, Number > | symmetrize (const Tensor< 2, 1, Number > &t) |
| template<typename Number > | |
| SymmetricTensor< 2, 2, Number > | symmetrize (const Tensor< 2, 2, Number > &t) |
| template<typename Number > | |
| SymmetricTensor< 2, 3, Number > | symmetrize (const Tensor< 2, 3, Number > &t) |
| template<int rank, int dim, typename Number > | |
| SymmetricTensor< rank, dim, Number > | operator* (const SymmetricTensor< rank, dim, Number > &t, const Number factor) |
| template<int rank, int dim, typename Number > | |
| SymmetricTensor< rank, dim, Number > | operator* (const Number factor, const SymmetricTensor< rank, dim, Number > &t) |
| template<int rank, int dim, typename Number > | |
| SymmetricTensor< rank, dim, Number > | operator/ (const SymmetricTensor< rank, dim, Number > &t, const Number factor) |
| template<int rank, int dim> | |
| SymmetricTensor< rank, dim > | operator* (const SymmetricTensor< rank, dim > &t, const double factor) |
| template<int rank, int dim> | |
| SymmetricTensor< rank, dim > | operator* (const double factor, const SymmetricTensor< rank, dim > &t) |
| template<int rank, int dim> | |
| SymmetricTensor< rank, dim > | operator/ (const SymmetricTensor< rank, dim > &t, const double factor) |
| template<typename Number > | |
| void | double_contract (SymmetricTensor< 2, 1, Number > &tmp, const SymmetricTensor< 4, 1, Number > &t, const SymmetricTensor< 2, 1, Number > &s) |
| template<typename Number > | |
| void | double_contract (SymmetricTensor< 2, 1, Number > &tmp, const SymmetricTensor< 2, 1, Number > &s, const SymmetricTensor< 4, 1, Number > &t) |
| template<typename Number > | |
| void | double_contract (SymmetricTensor< 2, 2, Number > &tmp, const SymmetricTensor< 4, 2, Number > &t, const SymmetricTensor< 2, 2, Number > &s) |
| template<typename Number > | |
| void | double_contract (SymmetricTensor< 2, 2, Number > &tmp, const SymmetricTensor< 2, 2, Number > &s, const SymmetricTensor< 4, 2, Number > &t) |
| template<typename Number > | |
| void | double_contract (SymmetricTensor< 2, 3, Number > &tmp, const SymmetricTensor< 4, 3, Number > &t, const SymmetricTensor< 2, 3, Number > &s) |
| template<typename Number > | |
| void | double_contract (SymmetricTensor< 2, 3, Number > &tmp, const SymmetricTensor< 2, 3, Number > &s, const SymmetricTensor< 4, 3, Number > &t) |
| template<int dim, typename Number > | |
| Tensor< 1, dim, Number > | operator* (const SymmetricTensor< 2, dim, Number > &src1, const Tensor< 1, dim, Number > &src2) |
| template<int dim, typename Number > | |
| std::ostream & | operator<< (std::ostream &out, const SymmetricTensor< 2, dim, Number > &t) |
| template<int dim, typename Number > | |
| std::ostream & | operator<< (std::ostream &out, const SymmetricTensor< 4, dim, Number > &t) |
Provide a class that stores symmetric tensors of rank 2,4,... efficiently, i.e. only store those off-diagonal elements of the full tensor that are not redundant. For example, for symmetric 2x2 tensors, this would be the elements 11, 22, and 12, while the element 21 is equal to the 12 element.
Using this class for symmetric tensors of rank 2 has advantages over matrices in many cases since the dimension is known to the compiler as well as the location of the data. It is therefore possible to produce far more efficient code than for matrices with runtime-dependent dimension. It is also more efficient than using the more general Tensor class, since less elements are stored, and the class automatically makes sure that the tensor represents a symmetric object.
For tensors of higher rank, the savings in storage are even higher. For example for the 3x3x3x3 tensors of rank 4, only 36 instead of the full 81 entries have to be stored.
While the definition of a symmetric rank-2 tensor is obvious, tensors of rank 4 are considered symmetric if they are operators mapping symmetric rank-2 tensors onto symmetric rank-2 tensors. This entails certain symmetry properties on the elements in their 4-dimensional index space, in particular that Cijkl=Cjikl=Cijlk. However, it does not imply the relation Cijkl=Cklij. Consequently, symmetric tensors of rank 4 as understood here are only tensors that map symmetric tensors onto symmetric tensors, but they do not necessarily induce a symmetric scalar product a:C:b=b:C:a or even a positive (semi-)definite form a:C:a, where a,b are symmetric rank-2 tensors and the colon indicates the common double-index contraction that acts as a product for symmetric tensors.
Symmetric tensors are most often used in structural and fluid mechanics, where strains and stresses are usually symmetric tensors, and the stress-strain relationship is given by a symmetric rank-4 tensor.
Note that symmetric tensors only exist with even numbers of indices. In other words, the only objects that you can use are SymmetricTensor<2,dim>, SymmetricTensor<4,dim>, etc, but SymmetricTensor<1,dim> and SymmetricTensor<3,dim> do not exist and their use will most likely lead to compiler errors.
The elements of a tensor t can be accessed using the bracket operator, i.e. for a tensor of rank 4, t[0][1][0][1] accesses the element t0101. This access can be used for both reading and writing (if the tensor is non-constant at least). You may also perform other operations on it, although that may lead to confusing situations because several elements of the tensor are stored at the same location. For example, for a rank-2 tensor that is assumed to be zero at the beginning, writing t[0][1]+=1; t[1][0]+=1; will lead to the same element being increased by one twice, because even though the accesses use different indices, the elements that are accessed are symmetric and therefore stored at the same location. It may therefore be useful in application programs to restrict operations on individual elements to simple reads or writes.
Definition at line 614 of file symmetric_tensor.h.
typedef internal::SymmetricTensorAccessors::StorageType<rank,dim,Number> SymmetricTensor< rank, dim, Number >::base_tensor_descriptor [private] |
A structure that describes properties of the base tensor.
Definition at line 990 of file symmetric_tensor.h.
typedef base_tensor_descriptor::base_tensor_type SymmetricTensor< rank, dim, Number >::base_tensor_type [private] |
Data storage type for a symmetric tensor.
Definition at line 996 of file symmetric_tensor.h.
| SymmetricTensor< rank, dim, Number >::SymmetricTensor | ( | ) |
Default constructor. Creates a zero tensor.
| SymmetricTensor< rank, dim, Number >::SymmetricTensor | ( | const Tensor< 2, dim, Number > & | t ) |
Constructor. Generate a symmetric tensor from a general one. Assumes that t is already symmetric, and in debug mode this is in fact checked. Note that no provision is made to assure that the tensor is symmetric only up to round-off error: if the incoming tensor is not exactly symmetric, then an exception is thrown. If you know that incoming tensor is symmetric only up to round-off, then you may want to call the symmetrize function first. If you aren't sure, it is good practice to check before calling symmetrize.
| SymmetricTensor< rank, dim, Number >::SymmetricTensor | ( | const Number(&) | array[internal::SymmetricTensorAccessors::StorageType< rank, dim, Number >::n_independent_components] ) |
A constructor that creates a symmetric tensor from an array holding its independent elements. Using this constructor assumes that the caller knows the order in which elements are stored in symmetric tensors; its use is therefore discouraged, but if you think you want to use it anyway you can query the order of elements using the unrolled_index() function.
This constructor is currently only implemented for symmetric tensors of rank 2.
The size of the array passed is equal to SymmetricTensor<rank,dim>::n_independent_component; the reason for using the object from the internal namespace is to work around bugs in some older compilers.
| SymmetricTensor& SymmetricTensor< rank, dim, Number >::operator= | ( | const SymmetricTensor< rank, dim, Number > & | ) |
Assignment operator.
| SymmetricTensor& SymmetricTensor< rank, dim, Number >::operator= | ( | const Number | d ) |
This operator assigns a scalar to a tensor. To avoid confusion with what exactly it means to assign a scalar value to a tensor, zero is the only value allowed for d, allowing the intuitive notation t=0 to reset all elements of the tensor to zero.
| SymmetricTensor< rank, dim, Number >::operator Tensor< rank, dim, Number > | ( | ) | const |
Convert the present symmetric tensor into a full tensor with the same elements, but using the different storage scheme of full tensors.
| bool SymmetricTensor< rank, dim, Number >::operator== | ( | const SymmetricTensor< rank, dim, Number > & | ) | const |
Test for equality of two tensors.
| bool SymmetricTensor< rank, dim, Number >::operator!= | ( | const SymmetricTensor< rank, dim, Number > & | ) | const |
Test for inequality of two tensors.
| SymmetricTensor& SymmetricTensor< rank, dim, Number >::operator+= | ( | const SymmetricTensor< rank, dim, Number > & | ) |
Add another tensor.
| SymmetricTensor& SymmetricTensor< rank, dim, Number >::operator-= | ( | const SymmetricTensor< rank, dim, Number > & | ) |
Subtract another tensor.
| SymmetricTensor& SymmetricTensor< rank, dim, Number >::operator*= | ( | const Number | factor ) |
Scale the tensor by factor, i.e. multiply all components by factor.
| SymmetricTensor& SymmetricTensor< rank, dim, Number >::operator/= | ( | const Number | factor ) |
Scale the vector by 1/factor.
| SymmetricTensor SymmetricTensor< rank, dim, Number >::operator+ | ( | const SymmetricTensor< rank, dim, Number > & | s ) | const |
Add two tensors. If possible, you should use operator += instead since this does not need the creation of a temporary.
| SymmetricTensor SymmetricTensor< rank, dim, Number >::operator- | ( | const SymmetricTensor< rank, dim, Number > & | s ) | const |
Subtract two tensors. If possible, you should use operator -= instead since this does not need the creation of a temporary.
| SymmetricTensor SymmetricTensor< rank, dim, Number >::operator- | ( | ) | const |
Unary minus operator. Negate all entries of a tensor.
| internal::SymmetricTensorAccessors::double_contraction_result<rank,2,dim,Number>::type SymmetricTensor< rank, dim, Number >::operator* | ( | const SymmetricTensor< 2, dim, Number > & | s ) | const |
Product between the present symmetric tensor and a tensor of rank 2. For example, if the present object is also a rank-2 tensor, then this is the scalar-product double contraction aijbij over all indices i,j. In this case, the return value evaluates to a single scalar. While it is possible to define other scalar product (and associated induced norms), this one seems to be the most appropriate one.
If the present object is a rank-4 tensor, the the result is a rank-2 tensor, the operation contracts over the last two indices of the present object and the indices of the argument, and the result is a tensor of rank 2.
Note that the multiplication operator for symmetrict tensors is defined to be a double contraction over two indices, while it is defined as a single contraction over only one index for regular Tensor objects. For symmetric tensors it therefore acts in a way that is commonly denoted by a "colon
multiplication" in the mathematica literature.
There are global functions double_contract that do the same work as this operator, but rather than returning the result as a return value, they write it into the first argument to the function.
| internal::SymmetricTensorAccessors::double_contraction_result<rank,4,dim,Number>::type SymmetricTensor< rank, dim, Number >::operator* | ( | const SymmetricTensor< 4, dim, Number > & | s ) | const |
Contraction over two indices of the present object with the rank-4 symmetric tensor given as argument.
| Number& SymmetricTensor< rank, dim, Number >::operator() | ( | const TableIndices< rank > & | indices ) |
Return a read-write reference to the indicated element.
| Number SymmetricTensor< rank, dim, Number >::operator() | ( | const TableIndices< rank > & | 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.
| internal::SymmetricTensorAccessors::Accessor<rank,dim,true,rank-1,Number> SymmetricTensor< rank, dim, Number >::operator[] | ( | const unsigned int | row ) | const |
Access the elements of a row of this symmetric tensor. This function is called for constant tensors.
| internal::SymmetricTensorAccessors::Accessor<rank,dim,false,rank-1,Number> SymmetricTensor< rank, dim, Number >::operator[] | ( | const unsigned int | row ) |
Access the elements of a row of this symmetric tensor. This function is called for non-constant tensors.
| Number SymmetricTensor< rank, dim, Number >::operator[] | ( | const TableIndices< rank > & | indices ) | const |
Access to an element where you specify the entire set of indices.
| Number& SymmetricTensor< rank, dim, Number >::operator[] | ( | const TableIndices< rank > & | indices ) |
Access to an element where you specify the entire set of indices.
| Number SymmetricTensor< rank, dim, Number >::norm | ( | ) | const |
Return the Frobenius-norm of a tensor, i.e. the square root of the sum of squares of all entries. This norm is induced by the scalar product defined above for two symmetric tensors. Note that it includes all entries of the tensor, counting symmetry, not only the unique ones (for example, for rank-2 tensors, this norm includes adding up the squares of upper right as well as lower left entries, not just one of them, although they are equal for symmetric tensors).
| static unsigned int SymmetricTensor< rank, dim, Number >::component_to_unrolled_index | ( | const TableIndices< rank > & | indices ) | [static] |
Tensors can be unrolled by simply pasting all elements into one long vector, but for this an order of elements has to be defined. For symmetric tensors, this function returns which index within the range [0,n_independent_components) the given entry in a symmetric tensor has.
| static TableIndices<rank> SymmetricTensor< rank, dim, Number >::unrolled_to_component_indices | ( | const unsigned int | i ) | [static] |
The opposite of the previous function: given an index
in the unrolled form of the tensor, return what set of indices
(for rank-2 tensors) or
(for rank-4 tensors) corresponds to it.
| void SymmetricTensor< rank, dim, Number >::clear | ( | ) |
Reset all values to zero.
Note that this is partly inconsistent with the semantics of the clear() member functions of the STL and of several other classes within deal.II which not only reset the values of stored elements to zero, but release all memory and return the object into a virginial state. However, since the size of objects of the present type is determined by its template parameters, resizing is not an option, and indeed the state where all elements have a zero value is the state right after construction of such an object.
| static std::size_t SymmetricTensor< rank, dim, Number >::memory_consumption | ( | ) | [static] |
Determine an estimate for the memory consumption (in bytes) of this object.
| void SymmetricTensor< rank, dim, Number >::serialize | ( | Archive & | ar, |
| const unsigned int | version | ||
| ) |
Read or write the data of this object to or from a stream for the purpose of serialization
friend class SymmetricTensor [friend] |
Make all other symmetric tensors friends.
Definition at line 1007 of file symmetric_tensor.h.
| Number2 trace | ( | const SymmetricTensor< 2, dim2, Number2 > & | d ) | [friend] |
Make a few more functions friends.
| Number2 determinant | ( | const SymmetricTensor< 2, dim2, Number2 > & | t ) | [friend] |
| SymmetricTensor<2,dim2,Number2> deviator | ( | const SymmetricTensor< 2, dim2, Number2 > & | t ) | [friend] |
| SymmetricTensor< 2, dim > unit_symmetric_tensor | ( | ) | [friend] |
Return a unit symmetric tensor of rank 2.
Return a unit symmetric tensor of rank 2 for double tensor.
Definition at line 2288 of file symmetric_tensor.h.
| SymmetricTensor< 4, dim > deviator_tensor | ( | ) | [friend] |
Return the tensor of rank 4 that, when multiplied by a symmetric rank 2 tensor t returns the deviator dev t. It is the operator representation of the linear deviator operator.
For every tensor t, there holds the identity deviator(t)==deviator_tensor<dim>()*t, up to numerical round-off. The reason this operator representation is provided is that one sometimes needs to invert operators like identity_tensor<dim>() + delta_t*deviator_tensor<dim>() or similar.
Definition at line 2344 of file symmetric_tensor.h.
| SymmetricTensor< 4, dim > identity_tensor | ( | ) | [friend] |
Returns the fourth-order symmetric identity tensor which maps symmetric second-order tensors to themselves.
Note that this tensor, even though it is the identity, has a somewhat funny form, and in particular does not only consist of zeros and ones. For example, for dim=2, the identity tensor has all zero entries except for id[0][0][0][0]=id[1][1][1][1]=1 and id[0][1][0][1]=id[0][1][1][0]=id[1][0][0][1]=id[1][0][1][0]=1/2. To see why this factor of 1/2 is necessary, consider computing A=Id . B. For the element a_01 we have a_01=id_0100 b_00 + id_0111 b_11 + id_0101 b_01 + id_0110 b_10. On the other hand, we need to have a_01=b_01, and symmetry implies b_01=b_10, leading to a_01=(id_0101+id_0110) b_01, or, again by symmetry, id_0101=id_0110=1/2. Similar considerations hold for the three-dimensional case.
Return the tensor of rank 4 that, when multiplied by a symmetric rank 2 tensor t returns the deviator dev t. It is the operator representation of the linear deviator operator.
Note that this tensor, even though it is the identity, has a somewhat funny form, and in particular does not only consist of zeros and ones. For example, for dim=2, the identity tensor has all zero entries except for id[0][0][0][0]=id[1][1][1][1]=1 and id[0][1][0][1]=id[0][1][1][0]=id[1][0][0][1]=id[1][0][1][0]=1/2. To see why this factor of 1/2 is necessary, consider computing A=Id . B. For the element a_01 we have a_01=id_0100 b_00 + id_0111 b_11 + id_0101 b_01 + id_0110 b_10. On the other hand, we need to have a_01=b_01, and symmetry implies b_01=b_10, leading to a_01=(id_0101+id_0110) b_01, or, again by symmetry, id_0101=id_0110=1/2. Similar considerations hold for the three-dimensional case.
Definition at line 2415 of file symmetric_tensor.h.
| SymmetricTensor<4,dim2,Number2> invert | ( | const SymmetricTensor< 4, dim2, Number2 > & | ) | [friend] |
| Number determinant | ( | const SymmetricTensor< 2, dim, Number > & | t ) | [related] |
Compute the determinant of a tensor or rank 2. The determinant is also commonly referred to as the third invariant of rank-2 tensors.
For a one-dimensional tensor, the determinant equals the only element and is therefore equivalent to the trace.
For greater notational simplicity, there is also a third_invariant function that returns the determinant of a tensor.
Definition at line 2106 of file symmetric_tensor.h.
References Assert, SymmetricTensor< rank, dim, Number >::data, and StandardExceptions::ExcNotImplemented().
Referenced by Tensor< rank_, dim, Number >::determinant(), and SymmetricTensor< rank, dim, Number >::third_invariant().
| double third_invariant | ( | const SymmetricTensor< 2, dim, Number > & | t ) | [related] |
Compute the determinant of a tensor or rank 2. This function therefore computes the same value as the determinant() functions and is only provided for greater notational simplicity (since there are also functions first_invariant and second_invariant).
Definition at line 2142 of file symmetric_tensor.h.
References SymmetricTensor< rank, dim, Number >::determinant().
| Number trace | ( | const SymmetricTensor< 2, dim, Number > & | d ) | [related] |
Compute and return the trace of a tensor of rank 2, i.e. the sum of its diagonal entries. The trace is the first invariant of a rank-2 tensor.
Definition at line 2158 of file symmetric_tensor.h.
References SymmetricTensor< rank, dim, Number >::data.
Referenced by SymmetricTensor< rank, dim, Number >::deviator(), and SymmetricTensor< rank, dim, Number >::first_invariant().
| Number first_invariant | ( | const SymmetricTensor< 2, dim, Number > & | t ) | [related] |
Compute the trace of a tensor or rank 2. This function therefore computes the same value as the trace() functions and is only provided for greater notational simplicity (since there are also functions second_invariant and third_invariant).
Definition at line 2178 of file symmetric_tensor.h.
References SymmetricTensor< rank, dim, Number >::trace().
| Number second_invariant | ( | const SymmetricTensor< 2, 1, Number > & | ) | [related] |
Compute the second invariant of a tensor of rank 2. The second invariant is defined as I2 = 1/2[ (trace sigma)^2 - trace (sigma^2) ].
Definition at line 2193 of file symmetric_tensor.h.
| Number second_invariant | ( | const SymmetricTensor< 2, 2, Number > & | t ) | [related] |
Compute the second invariant of a tensor of rank 2. The second invariant is defined as I2 = 1/2[ (trace sigma)^2 - trace (sigma^2) ].
Definition at line 2209 of file symmetric_tensor.h.
| Number second_invariant | ( | const SymmetricTensor< 2, 3, Number > & | t ) | [related] |
Compute the second invariant of a tensor of rank 2. The second invariant is defined as I2 = 1/2[ (trace sigma)^2 - trace (sigma^2) ].
Definition at line 2225 of file symmetric_tensor.h.
| SymmetricTensor< rank, dim, Number > transpose | ( | const SymmetricTensor< rank, dim, Number > & | t ) | [related] |
Return the transpose of the given symmetric tensor. Since we are working with symmetric objects, the transpose is of course the same as the original tensor. This function mainly exists for compatibility with the Tensor class.
Definition at line 2246 of file symmetric_tensor.h.
| SymmetricTensor< 2, dim, Number > deviator | ( | const SymmetricTensor< 2, dim, Number > & | t ) | [related] |
Compute the deviator of a symmetric tensor, which is defined as dev[s] = s - 1/dim*tr[s]*I, where I is the identity operator. This quantity equals the original tensor minus its contractive or dilative component and refers to the shear in, for example, elasticity.
Definition at line 2265 of file symmetric_tensor.h.
References SymmetricTensor< rank, dim, Number >::data, and SymmetricTensor< rank, dim, Number >::trace().
| SymmetricTensor< 4, dim, Number > invert | ( | const SymmetricTensor< 4, dim, Number > & | t ) | [related] |
Invert a symmetric rank-4 tensor. Since symmetric rank-4 tensors are mappings from and to symmetric rank-2 tensors, they can have an inverse. This function computes it, if it exists, for the case that the dimension equals either 1 or 2.
If a tensor is not invertible, then the result is unspecified, but will likely contain the results of a division by zero or a very small number at the very least.
Definition at line 2485 of file symmetric_tensor.h.
References Assert, SymmetricTensor< rank, dim, Number >::data, and StandardExceptions::ExcNotImplemented().
| SymmetricTensor< 4, 3, double > invert | ( | const SymmetricTensor< 4, 3, double > & | t ) | [related] |
Invert a symmetric rank-4 tensor. Since symmetric rank-4 tensors are mappings from and to symmetric rank-2 tensors, they can have an inverse. This function computes it, if it exists, for the case that the dimension equals 3.
If a tensor is not invertible, then the result is unspecified, but will likely contain the results of a division by zero or a very small number at the very least.
| SymmetricTensor< 4, dim, Number > outer_product | ( | const SymmetricTensor< 2, dim, Number > & | t1, |
| const SymmetricTensor< 2, dim, Number > & | t2 | ||
| ) | [related] |
Return the tensor of rank 4 that is the outer product of the two tensors given as arguments, i.e. the result
satisfies T phi = t1 (t2 : phi) for all symmetric tensors phi.
For example, the deviator tensor can be computed as identity_tensor<dim,Number>() - 1/d*outer_product(unit_symmetric_tensor<dim,Number>(), unit_symmetric_tensor<dim,Number>()), since the (double) contraction with the unit tensor yields the trace of a symmetric tensor.
Definition at line 2601 of file symmetric_tensor.h.
| SymmetricTensor< 2, 1, Number > symmetrize | ( | const Tensor< 2, 1, Number > & | t ) | [related] |
Return the symmetrized version of a full rank-2 tensor, i.e. (t+transpose(t))/2, as a symmetric rank-2 tensor. This is the version for dim==1.
Definition at line 2629 of file symmetric_tensor.h.
| SymmetricTensor< 2, 2, Number > symmetrize | ( | const Tensor< 2, 2, Number > & | t ) | [related] |
Return the symmetrized version of a full rank-2 tensor, i.e. (t+transpose(t))/2, as a symmetric rank-2 tensor. This is the version for dim==2.
Definition at line 2649 of file symmetric_tensor.h.
| SymmetricTensor< 2, 3, Number > symmetrize | ( | const Tensor< 2, 3, Number > & | t ) | [related] |
Return the symmetrized version of a full rank-2 tensor, i.e. (t+transpose(t))/2, as a symmetric rank-2 tensor. This is the version for dim==3.
Definition at line 2669 of file symmetric_tensor.h.
| SymmetricTensor< rank, dim, Number > operator* | ( | const SymmetricTensor< rank, dim, Number > & | t, |
| const Number | factor | ||
| ) | [related] |
Multiplication of a symmetric tensor of general rank with a scalar from the right.
Definition at line 2690 of file symmetric_tensor.h.
| SymmetricTensor< rank, dim, Number > operator* | ( | const Number | factor, |
| const SymmetricTensor< rank, dim, Number > & | t | ||
| ) | [related] |
Multiplication of a symmetric tensor of general rank with a scalar from the left.
Definition at line 2709 of file symmetric_tensor.h.
| SymmetricTensor< rank, dim, Number > operator/ | ( | const SymmetricTensor< rank, dim, Number > & | t, |
| const Number | factor | ||
| ) | [related] |
Division of a symmetric tensor of general rank by a scalar.
Definition at line 2727 of file symmetric_tensor.h.
| SymmetricTensor< rank, dim > operator* | ( | const SymmetricTensor< rank, dim > & | t, |
| const double | factor | ||
| ) | [related] |
Multiplication of a symmetric tensor of general rank with a scalar from the right.
Definition at line 2746 of file symmetric_tensor.h.
| SymmetricTensor< rank, dim > operator* | ( | const double | factor, |
| const SymmetricTensor< rank, dim > & | t | ||
| ) | [related] |
Multiplication of a symmetric tensor of general rank with a scalar from the left.
Definition at line 2765 of file symmetric_tensor.h.
| SymmetricTensor< rank, dim > operator/ | ( | const SymmetricTensor< rank, dim > & | t, |
| const double | factor | ||
| ) | [related] |
Division of a symmetric tensor of general rank by a scalar.
Definition at line 2783 of file symmetric_tensor.h.
| void double_contract | ( | SymmetricTensor< 2, 1, Number > & | tmp, |
| const SymmetricTensor< 4, 1, Number > & | t, | ||
| const SymmetricTensor< 2, 1, Number > & | s | ||
| ) | [related] |
Double contraction between a rank-4 and a rank-2 symmetric tensor, resulting in the symmetric tensor of rank 2 that is given as first argument to this function. This operation is the symmetric tensor analogon of a matrix-vector multiplication.
This function does the same as the member operator* of the SymmetricTensor class. It should not be used, however, since the member operator has knowledge of the actual data storage format and is at least 2 orders of magnitude faster. This function mostly exists for compatibility purposes with the general tensor class.
Definition at line 2812 of file symmetric_tensor.h.
| void double_contract | ( | SymmetricTensor< 2, 1, Number > & | tmp, |
| const SymmetricTensor< 2, 1, Number > & | s, | ||
| const SymmetricTensor< 4, 1, Number > & | t | ||
| ) | [related] |
Double contraction between a rank-4 and a rank-2 symmetric tensor, resulting in the symmetric tensor of rank 2 that is given as first argument to this function. This operation is the symmetric tensor analogon of a matrix-vector multiplication.
This function does the same as the member operator* of the SymmetricTensor class. It should not be used, however, since the member operator has knowledge of the actual data storage format and is at least 2 orders of magnitude faster. This function mostly exists for compatibility purposes with the general tensor class.
Definition at line 2839 of file symmetric_tensor.h.
| void double_contract | ( | SymmetricTensor< 2, 2, Number > & | tmp, |
| const SymmetricTensor< 4, 2, Number > & | t, | ||
| const SymmetricTensor< 2, 2, Number > & | s | ||
| ) | [related] |
Double contraction between a rank-4 and a rank-2 symmetric tensor, resulting in the symmetric tensor of rank 2 that is given as first argument to this function. This operation is the symmetric tensor analogon of a matrix-vector multiplication.
This function does the same as the member operator* of the SymmetricTensor class. It should not be used, however, since the member operator has knowledge of the actual data storage format and is at least 2 orders of magnitude faster. This function mostly exists for compatibility purposes with the general tensor class.
Definition at line 2865 of file symmetric_tensor.h.
| void double_contract | ( | SymmetricTensor< 2, 2, Number > & | tmp, |
| const SymmetricTensor< 2, 2, Number > & | s, | ||
| const SymmetricTensor< 4, 2, Number > & | t | ||
| ) | [related] |
Double contraction between a rank-4 and a rank-2 symmetric tensor, resulting in the symmetric tensor of rank 2 that is given as first argument to this function. This operation is the symmetric tensor analogon of a matrix-vector multiplication.
This function does the same as the member operator* of the SymmetricTensor class. It should not be used, however, since the member operator has knowledge of the actual data storage format and is at least 2 orders of magnitude faster. This function mostly exists for compatibility purposes with the general tensor class.
Definition at line 2898 of file symmetric_tensor.h.
| void double_contract | ( | SymmetricTensor< 2, 3, Number > & | tmp, |
| const SymmetricTensor< 4, 3, Number > & | t, | ||
| const SymmetricTensor< 2, 3, Number > & | s | ||
| ) | [related] |
Double contraction between a rank-4 and a rank-2 symmetric tensor, resulting in the symmetric tensor of rank 2 that is given as first argument to this function. This operation is the symmetric tensor analogon of a matrix-vector multiplication.
This function does the same as the member operator* of the SymmetricTensor class. It should not be used, however, since the member operator has knowledge of the actual data storage format and is at least 2 orders of magnitude faster. This function mostly exists for compatibility purposes with the general tensor class.
Definition at line 2931 of file symmetric_tensor.h.
| void double_contract | ( | SymmetricTensor< 2, 3, Number > & | tmp, |
| const SymmetricTensor< 2, 3, Number > & | s, | ||
| const SymmetricTensor< 4, 3, Number > & | t | ||
| ) | [related] |
Double contraction between a rank-4 and a rank-2 symmetric tensor, resulting in the symmetric tensor of rank 2 that is given as first argument to this function. This operation is the symmetric tensor analogon of a matrix-vector multiplication.
This function does the same as the member operator* of the SymmetricTensor class. It should not be used, however, since the member operator has knowledge of the actual data storage format and is at least 2 orders of magnitude faster. This function mostly exists for compatibility purposes with the general tensor class.
Definition at line 2967 of file symmetric_tensor.h.
| Tensor< 1, dim, Number > operator* | ( | const SymmetricTensor< 2, dim, Number > & | src1, |
| const Tensor< 1, dim, Number > & | src2 | ||
| ) | [related] |
Multiplication operator performing a contraction of the last index of the first argument and the first index of the second argument. This function therefore does the same as the corresponding contract function, but returns the result as a return value, rather than writing it into the reference given as the first argument to the contract function.
Note that for the Tensor class, the multiplication operator only performs a contraction over a single pair of indices. This is in contrast to the multiplication operator for symmetric tensors, which does the double contraction.
Definition at line 3003 of file symmetric_tensor.h.
| std::ostream & operator<< | ( | std::ostream & | out, |
| const SymmetricTensor< 2, dim, Number > & | t | ||
| ) | [related] |
Output operator for symmetric tensors of rank 2. Print the elements consecutively, with a space in between, two spaces between rank 1 subtensors, three between rank 2 and so on. No special amends are made to represents the symmetry in the output, for example by outputting only the unique entries.
Definition at line 3025 of file symmetric_tensor.h.
| std::ostream & operator<< | ( | std::ostream & | out, |
| const SymmetricTensor< 4, dim, Number > & | t | ||
| ) | [related] |
Output operator for symmetric tensors of rank 4. Print the elements consecutively, with a space in between, two spaces between rank 1 subtensors, three between rank 2 and so on. No special amends are made to represents the symmetry in the output, for example by outputting only the unique entries.
Definition at line 3053 of file symmetric_tensor.h.
const unsigned int SymmetricTensor< rank, dim, Number >::dimension = dim [static] |
Provide a way to get the dimension of an object without explicit knowledge of it's data type. Implementation is this way instead of providing a function dimension() because now it is possible to get the dimension at compile time without the expansion and preevaluation of an inlined function; the compiler may therefore produce more efficient code and you may use this value to declare other data types.
Definition at line 634 of file symmetric_tensor.h.
const unsigned int SymmetricTensor< rank, dim, Number >::n_independent_components [static] |
internal::SymmetricTensorAccessors::StorageType<rank,dim,Number>:: n_independent_components
An integer denoting the number of independent components that fully describe a symmetric tensor. In
space dimensions, this number equals
for symmetric tensors of rank 2.
Definition at line 645 of file symmetric_tensor.h.
base_tensor_type SymmetricTensor< rank, dim, Number >::data [private] |
The place where we store the data of the tensor.
Definition at line 1002 of file symmetric_tensor.h.
Referenced by SymmetricTensor< rank, dim, Number >::determinant(), SymmetricTensor< rank, dim, Number >::deviator(), SymmetricTensor< rank, dim, Number >::invert(), and SymmetricTensor< rank, dim, Number >::trace().
documentation generated on Fri Feb 3 2012 06:04:12 by
doxygen
1.7.2