16 #ifndef dealii_derivative_form_h
17 #define dealii_derivative_form_h
57 template <
int order,
int dim,
int spacedim,
typename Number =
double>
112 template <
typename OtherNumber>
175 <<
"Invalid DerivativeForm index " << arg1);
197 template <
int order,
int dim,
int spacedim,
typename Number>
202 ExcMessage(
"Only allowed for forms with dim==spacedim."));
204 for (
unsigned int j = 0; j < dim; ++j)
210 template <
int order,
int dim,
int spacedim,
typename Number>
214 for (
unsigned int j = 0; j < spacedim; ++j)
220 template <
int order,
int dim,
int spacedim,
typename Number>
228 for (
unsigned int j = 0; j < dim; ++j)
235 template <
int order,
int dim,
int spacedim,
typename Number>
240 for (
unsigned int j = 0; j < spacedim; ++j)
247 template <
int order,
int dim,
int spacedim,
typename Number>
252 Assert((1 == spacedim) && (order == 1),
253 ExcMessage(
"Only allowed for spacedim==1 and order==1."));
262 template <
int order,
int dim,
int spacedim,
typename Number>
263 template <
typename OtherNumber>
268 for (
unsigned int j = 0; j < spacedim; ++j)
275 template <
int order,
int dim,
int spacedim,
typename Number>
286 template <
int order,
int dim,
int spacedim,
typename Number>
289 const unsigned int i)
const
298 template <
int order,
int dim,
int spacedim,
typename Number>
302 Assert((1 == spacedim) && (order == 1),
310 template <
int order,
int dim,
int spacedim,
typename Number>
319 for (
unsigned int j = 0; j < dim; ++j)
327 template <
int order,
int dim,
int spacedim,
typename Number>
334 for (
unsigned int i = 0; i < spacedim; ++i)
335 for (
unsigned int j = 0; j < dim; ++j)
336 tt[j][i] = (*
this)[i][j];
343 template <
int order,
int dim,
int spacedim,
typename Number>
350 for (
unsigned int i = 0; i < spacedim; ++i)
351 for (
unsigned int j = 0; j < dim; ++j)
352 dest[i][j] = (*
this)[i] *
T[j];
359 template <
int order,
int dim,
int spacedim,
typename Number>
364 for (
unsigned int i = 0; i < spacedim; ++i)
365 sum_of_squares += tensor[i].norm_square();
371 template <
int order,
int dim,
int spacedim,
typename Number>
387 for (
unsigned int i = 0; i < dim; ++i)
388 for (
unsigned int j = 0; j < dim; ++j)
389 G[i][j] = DF_t[i] * DF_t[j];
397 template <
int order,
int dim,
int spacedim,
typename Number>
411 for (
unsigned int i = 0; i < dim; ++i)
412 for (
unsigned int j = 0; j < dim; ++j)
413 G[i][j] = DF_t[i] * DF_t[j];
415 return (this->times_T_t(
invert(G)));
420 template <
int order,
int dim,
int spacedim,
typename Number>
452 template <
int spacedim,
int dim,
typename Number1,
typename Number2>
458 for (
unsigned int i = 0; i < spacedim; ++i)
459 dest[i] = grad_F[i] * d_x;
474 template <
int spacedim,
int dim,
typename Number1,
typename Number2>
484 for (
unsigned int i = 0; i < dim; ++i)
503 template <
int dim,
typename Number1,
typename Number2>
509 for (
unsigned int i = 0; i < dim; ++i)
524 template <
int spacedim,
540 for (
unsigned int i = 0; i < n_components; ++i)
563 template <
int spacedim,
int dim,
typename Number1,
typename Number2>
570 for (
unsigned int i = 0; i < spacedim; ++i)
584 template <
int dim,
int spacedim,
typename Number>
constexpr DEAL_II_HOST Number determinant(const SymmetricTensor< 2, dim, Number > &)
constexpr DEAL_II_HOST SymmetricTensor< 2, dim, Number > invert(const SymmetricTensor< 2, dim, Number > &t)
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define Assert(cond, exc)
static ::ExceptionBase & ExcInvalidTensorIndex(int arg1)
#define AssertIndexRange(index, range)
#define DeclException1(Exception1, type1, outsequence)
static ::ExceptionBase & ExcMessage(std::string arg1)
::VectorizedArray< Number, width > sqrt(const ::VectorizedArray< Number, width > &)
typename internal::ProductTypeImpl< typename std::decay< T >::type, typename std::decay< U >::type >::type type
constexpr DEAL_II_HOST Number determinant(const SymmetricTensor< 2, dim, Number > &)