15#ifndef dealii_fe_point_evaluation_h
16#define dealii_fe_point_evaluation_h
47 <<
"You are requesting information from an FEPointEvaluationBase "
48 <<
"object for which this kind of information has not been computed. "
49 <<
"What information these objects compute is determined by the update_* "
50 <<
"flags you pass to MappingInfo() in the Constructor. Here, "
51 <<
"the operation you are attempting requires the <" << arg1
52 <<
"> flag to be set, but it was apparently not specified "
53 <<
"upon initialization.");
63 typename Enable =
void>
69 typename ::internal::VectorizedArrayTrait<
78 n_components == spacedim,
90 const unsigned int component,
94 result[component] = vector_entry;
97 static scalar_value_type
103 static scalar_value_type
108 for (
unsigned int c = 0; c < n_components; ++c)
109 result_scalar[c] = result[c].sum();
111 return result_scalar;
119 return result[component].
sum();
124 const unsigned int vector_lane,
127 for (
unsigned int i = 0; i < n_components; ++i)
128 for (
unsigned int d = 0; d < dim; ++d)
131 value[d][i], vector_lane);
136 const unsigned int vector_lane,
139 for (
unsigned int i = 0; i < n_components; ++i)
140 for (
unsigned int d = 0; d < dim; ++d)
142 value[d][i], vector_lane) = result[i][d];
147 const unsigned int vector_lane,
150 for (
unsigned int i = 0; i < n_components; ++i)
151 for (
unsigned int d = 0; d < dim; ++d)
153 value[d][i], vector_lane) = result[i][d];
158 const unsigned int vector_lane)
160 for (
unsigned int i = 0; i < n_components; ++i)
161 for (
unsigned int d = 0; d < spacedim; ++d)
168 const unsigned int vector_lane,
171 for (
unsigned int i = 0; i < n_components; ++i)
172 result[i] =
value[i][vector_lane];
185 const unsigned int vector_lane,
188 for (
unsigned int i = 0; i < n_components; ++i)
189 value[i][vector_lane] = result[i];
203 for (
unsigned int i = 0; i < n_components; ++i)
210 const unsigned int vector_lane,
211 const unsigned int component,
215 vector_lane) += shape_value;
220 const unsigned int vector_lane,
221 const unsigned int component)
229 const unsigned int vector_lane,
230 const unsigned int component,
233 for (
unsigned int d = 0; d < spacedim; ++d)
241 const unsigned int vector_lane,
242 const unsigned int component)
245 for (
unsigned int d = 0; d < spacedim; ++d)
253 template <
int dim,
int spacedim,
typename Number>
259 typename ::internal::VectorizedArrayTrait<
276 result = vector_entry;
279 static scalar_value_type
285 static scalar_value_type
299 const unsigned int vector_lane,
302 for (
unsigned int d = 0; d < dim; ++d)
303 result[d] =
value[d][vector_lane];
316 const unsigned int vector_lane,
319 for (
unsigned int d = 0; d < dim; ++d)
320 value[d][vector_lane] = result[d];
333 const unsigned int vector_lane)
335 for (
unsigned int d = 0; d < spacedim; ++d)
342 const unsigned int vector_lane,
345 result =
value[vector_lane];
358 const unsigned int vector_lane,
361 value[vector_lane] = result;
380 const unsigned int vector_lane,
390 const unsigned int vector_lane,
398 const unsigned int vector_lane,
402 for (
unsigned int d = 0; d < spacedim; ++d)
409 const unsigned int vector_lane,
413 for (
unsigned int d = 0; d < spacedim; ++d)
420 template <
int dim,
typename Number>
425 std::enable_if_t<dim != 1>>
430 typename ::internal::VectorizedArrayTrait<
444 const unsigned int component,
448 result[component] = vector_entry;
451 static scalar_value_type
457 static scalar_value_type
462 for (
unsigned int c = 0; c < dim; ++c)
463 result_scalar[c] = result[c].sum();
465 return result_scalar;
473 return result[component].
sum();
478 const unsigned int vector_lane,
481 for (
unsigned int i = 0; i < dim; ++i)
482 for (
unsigned int d = 0; d < dim; ++d)
485 value[d][i], vector_lane);
490 const unsigned int vector_lane,
493 for (
unsigned int i = 0; i < dim; ++i)
494 for (
unsigned int d = 0; d < dim; ++d)
496 value[d][i], vector_lane) = result[i][d];
501 const unsigned int vector_lane)
503 for (
unsigned int i = 0; i < dim; ++i)
504 for (
unsigned int d = 0; d < dim; ++d)
511 const unsigned int vector_lane,
514 for (
unsigned int i = 0; i < dim; ++i)
515 result[i] =
value[i][vector_lane];
528 const unsigned int vector_lane,
531 for (
unsigned int i = 0; i < dim; ++i)
532 value[i][vector_lane] = result[i];
546 for (
unsigned int i = 0; i < dim; ++i)
553 const unsigned int vector_lane,
554 const unsigned int component,
558 vector_lane) += shape_value;
563 const unsigned int vector_lane,
564 const unsigned int component)
572 const unsigned int vector_lane,
573 const unsigned int component,
576 for (
unsigned int d = 0; d < dim; ++d)
584 const unsigned int vector_lane,
585 const unsigned int component)
588 for (
unsigned int d = 0; d < dim; ++d)
596 template <
int dim,
int spacedim>
599 const unsigned int base_element_number);
601 template <
int dim,
int spacedim>
605 template <
int dim,
int spacedim>
606 std::vector<Polynomials::Polynomial<double>>
619template <
int n_components_,
622 typename Number =
double>
635 using ETT =
typename internal::FEPointEvaluation::
636 EvaluatorTypeTraits<dim, spacedim, n_components, Number>;
642 typename ETT::interface_vectorized_unit_gradient_type;
666 const unsigned int first_selected_component = 0);
690 const unsigned int first_selected_component = 0,
799 JxW(
const unsigned int point_index)
const;
865 setup(
const unsigned int first_selected_component);
872 template <
bool is_face,
bool is_linear>
905 std::vector<Polynomials::Polynomial<double>>
poly;
1041 std::unique_ptr<NonMatching::MappingInfo<dim, spacedim, Number>>
1120template <
int n_components_,
1123 typename Number =
double>
1137 using ETT =
typename internal::FEPointEvaluation::
1138 EvaluatorTypeTraits<dim, spacedim, n_components, Number>;
1145 typename ETT::interface_vectorized_unit_gradient_type;
1168 const unsigned int first_selected_component = 0);
1189 const unsigned int first_selected_component = 0);
1218 reinit(
const unsigned int cell_index);
1232 template <std::
size_t str
ide_view>
1275 template <std::
size_t str
ide_view>
1279 const bool sum_into_values =
false);
1306 const bool sum_into_values =
false);
1334 template <std::
size_t str
ide_view>
1339 const bool sum_into_values =
false);
1370 const bool sum_into_values =
false);
1408 template <
bool is_linear, std::
size_t str
ide_view>
1417 template <
bool is_linear, std::
size_t str
ide_view>
1422 const unsigned int n_shapes,
1423 const unsigned int qb,
1430 template <
bool is_linear, std::
size_t str
ide_view>
1439 template <std::
size_t str
ide_view>
1450 template <
bool is_linear>
1454 const unsigned int n_shapes,
1455 const unsigned int qb,
1465 template <
bool is_linear, std::
size_t str
ide_view>
1470 const bool sum_into_values);
1475 template <
bool do_JxW,
bool is_linear, std::
size_t str
ide_view>
1480 const bool sum_into_values);
1485 template <
bool do_JxW, std::
size_t str
ide_view>
1490 const bool sum_into_values);
1495 template <
bool do_JxW, std::
size_t str
ide_view>
1500 const bool sum_into_values);
1530template <
int n_components_,
1533 typename Number =
double>
1547 using ETT =
typename internal::FEPointEvaluation::
1548 EvaluatorTypeTraits<dim, spacedim, n_components, Number>;
1555 typename ETT::interface_vectorized_unit_gradient_type;
1564 const unsigned int first_selected_component = 0);
1571 reinit(
const unsigned int cell_index,
const unsigned int face_number);
1578 reinit(
const unsigned int face_index);
1591 template <std::
size_t str
ide_view>
1634 template <std::
size_t str
ide_view>
1638 const bool sum_into_values =
false);
1665 const bool sum_into_values =
false);
1689 template <std::
size_t str
ide_view>
1694 const bool sum_into_values =
false);
1721 const bool sum_into_values =
false);
1729 template <
int str
ide_face_dof = VectorizedArrayType::size()>
1740 template <
int str
ide_face_dof = VectorizedArrayType::size()>
1744 const bool sum_into_values =
false);
1778 template <
bool is_linear, std::
size_t str
ide_view>
1784 template <
bool do_JxW,
bool is_linear, std::
size_t str
ide_view>
1789 const bool sum_into_values);
1795 template <
bool is_linear,
int str
ide_face_dof>
1804 template <
bool do_JxW,
bool is_linear,
int str
ide_face_dof>
1809 const bool sum_into_values);
1817template <
int n_components_,
int dim,
int spacedim,
typename Number>
1822 const unsigned int first_selected_component)
1831 std::make_unique<
NonMatching::MappingInfo<dim, spacedim, Number>>(
1840 setup(first_selected_component);
1845template <
int n_components_,
int dim,
int spacedim,
typename Number>
1850 const unsigned int first_selected_component,
1865 setup(first_selected_component);
1870template <
int n_components_,
int dim,
int spacedim,
typename Number>
1893 other.mapping_info_on_the_fly ?
1908template <
int n_components_,
int dim,
int spacedim,
typename Number>
1940template <
int n_components_,
int dim,
int spacedim,
typename Number>
1943 const unsigned int first_selected_component)
1946 fe->n_components() + 1);
1950 bool same_base_element =
true;
1951 unsigned int base_element_number = 0;
1953 unsigned int component = 0;
1954 for (; base_element_number <
fe->n_base_elements(); ++base_element_number)
1955 if (component +
fe->element_multiplicity(base_element_number) >
1956 first_selected_component)
1959 component +
fe->element_multiplicity(base_element_number))
1960 same_base_element =
false;
1965 component +=
fe->element_multiplicity(base_element_number);
1969 *
fe, base_element_number) &&
1977 fe->base_element(base_element_number));
1979 bool is_lexicographic =
true;
1980 for (
unsigned int i = 0; i <
renumber.size(); ++i)
1982 is_lexicographic =
false;
1984 if (is_lexicographic)
2006 const unsigned int component = first_selected_component + d;
2007 for (
unsigned int i = 0; i <
fe->n_dofs_per_cell(); ++i)
2009 const bool is_primitive =
2010 fe->is_primitive() ||
fe->is_primitive(i);
2013 (component ==
fe->system_to_component_index(i).first);
2016 (
fe->get_nonzero_components(i)[component] ==
true);
2026template <
int n_components_,
int dim,
int spacedim,
typename Number>
2027template <
bool is_face,
bool is_linear>
2031 const unsigned int geometry_index =
2032 mapping_info->template compute_geometry_index_offset<is_face>(
2038 mapping_info->get_n_q_points_unvectorized(geometry_index);
2044 const unsigned int n_q_points_before =
n_q_points;
2061 const unsigned int unit_point_offset =
2062 mapping_info->compute_unit_point_index_offset(geometry_index);
2071 const unsigned int data_offset =
2072 mapping_info->compute_data_index_offset(geometry_index);
2073 const unsigned int compressed_data_offset =
2074 mapping_info->compute_compressed_data_index_offset(geometry_index);
2076 const UpdateFlags update_flags_mapping =
2102 const std::size_t n_shapes =
poly.size();
2134 shapes.data() + qb * n_shapes,
2153template <
int n_components_,
int dim,
int spacedim,
typename Number>
2157 Number>::value_type &
2159 const unsigned int point_index)
const
2162 return values[point_index];
2167template <
int n_components_,
int dim,
int spacedim,
typename Number>
2171 Number>::gradient_type &
2173 const unsigned int point_index)
const
2181template <
int n_components_,
int dim,
int spacedim,
typename Number>
2184 const unsigned int point_index)
const
2187 "Only makes sense for a vector field with dim components");
2195template <
int n_components_,
int dim,
int spacedim,
typename Number>
2199 const unsigned int point_index)
2207template <
int n_components_,
int dim,
int spacedim,
typename Number>
2211 const unsigned int point_index)
2219template <
int n_components_,
int dim,
int spacedim,
typename Number>
2222 const Number &
value,
2223 const unsigned int point_index)
2226 "Only makes sense for a vector field with dim components");
2230 for (
unsigned int d = 0; d < dim; ++d)
2236template <
int n_components_,
int dim,
int spacedim,
typename Number>
2239 const unsigned int point_index)
const
2244 ExcFEPointEvaluationAccessToUninitializedMappingField(
2245 "update_jacobians"));
2247 GeometryType::affine ?
2254template <
int n_components_,
int dim,
int spacedim,
typename Number>
2257 const unsigned int point_index)
const
2262 ExcFEPointEvaluationAccessToUninitializedMappingField(
2263 "update_inverse_jacobians"));
2273template <
int n_components_,
int dim,
int spacedim,
typename Number>
2276 const unsigned int point_index)
const
2281 ExcFEPointEvaluationAccessToUninitializedMappingField(
2282 "update_JxW_values"));
2288template <
int n_components_,
int dim,
int spacedim,
typename Number>
2291 const unsigned int point_index)
const
2299template <
int n_components_,
int dim,
int spacedim,
typename Number>
2302 const unsigned int point_index)
const
2307 ExcFEPointEvaluationAccessToUninitializedMappingField(
2308 "update_quadrature_points"));
2314template <
int n_components_,
int dim,
int spacedim,
typename Number>
2317 const unsigned int point_index)
const
2322 for (
unsigned int d = 0; d < dim; ++d)
2330template <
int n_components_,
int dim,
int spacedim,
typename Number>
2340template <
int n_components_,
int dim,
int spacedim,
typename Number>
2344 const unsigned int first_selected_component)
2348 first_selected_component)
2353template <
int n_components_,
int dim,
int spacedim,
typename Number>
2358 const unsigned int first_selected_component)
2363 first_selected_component)
2368template <
int n_components_,
int dim,
int spacedim,
typename Number>
2384template <
int n_components_,
int dim,
int spacedim,
typename Number>
2394template <
int n_components_,
int dim,
int spacedim,
typename Number>
2409 this->
fe_values = std::make_shared<FEValues<dim, spacedim>>(
2413 std::vector<
Point<dim>>(unit_points.begin(), unit_points.end())),
2426template <
int n_components_,
int dim,
int spacedim,
typename Number>
2429 const unsigned int cell_index)
2445 for (
unsigned int d = 0; d < dim; ++d)
2449 this->
fe_values = std::make_shared<FEValues<dim, spacedim>>(
2453 std::vector<
Point<dim>>(unit_points.begin(), unit_points.end())),
2457 this->
mapping_info->get_cell_iterator(this->current_cell_index));
2463template <
int n_components_,
int dim,
int spacedim,
typename Number>
2464template <std::
size_t str
ide_view>
2496template <
int n_components_,
int dim,
int spacedim,
typename Number>
2503 solution_values.
size()),
2509template <
int n_components_,
int dim,
int spacedim,
typename Number>
2510template <std::
size_t str
ide_view>
2515 const bool sum_into_values)
2522template <
int n_components_,
int dim,
int spacedim,
typename Number>
2527 const bool sum_into_values)
2530 solution_values.
size()),
2537template <
int n_components_,
int dim,
int spacedim,
typename Number>
2546 return_value +=
values[point_index] * this->
JxW(point_index);
2548 return ETT::sum_value(return_value);
2553template <
int n_components_,
int dim,
int spacedim,
typename Number>
2560 "Calling this function only makes sense in fully vectorized mode."));
2563 const unsigned int n_filled_lanes =
2565 if (n_filled_lanes == 0)
2568 return n_filled_lanes;
2576template <
int n_components_,
int dim,
int spacedim,
typename Number>
2577template <std::
size_t str
ide_view>
2582 const bool sum_into_values)
2589template <
int n_components_,
int dim,
int spacedim,
typename Number>
2594 const bool sum_into_values)
2597 solution_values.
size()),
2604template <
int n_components_,
int dim,
int spacedim,
typename Number>
2605template <
bool is_linear, std::
size_t str
ide_view>
2610 const unsigned int dofs_per_comp =
2613 for (
unsigned int comp = 0; comp <
n_components; ++comp)
2615 const std::size_t offset =
2620 for (
unsigned int i = 0; i < dofs_per_comp; ++i)
2621 ETT::read_value(solution_values[i + offset],
2627 const unsigned int *renumber_ptr = this->
renumber.data() + offset;
2628 for (
unsigned int i = 0; i < dofs_per_comp; ++i)
2629 ETT::read_value(solution_values[renumber_ptr[i]],
2638template <
int n_components_,
int dim,
int spacedim,
typename Number>
2639template <
bool is_linear, std::
size_t str
ide_view>
2644 const unsigned int n_shapes,
2645 const unsigned int qb,
2651 std::array<vectorized_value_type, dim + 1> result;
2652 if constexpr (is_linear)
2661 stride_view>(solution_values.
data(), this->unit_point_ptr[qb]);
2673 false>(this->
shapes.data() + qb * n_shapes,
2676 gradient[0] = result[0];
2678 gradient[1] = result[1];
2680 gradient[2] = result[2];
2681 value = result[dim];
2685 if constexpr (is_linear)
2692 stride_view>(solution_values.
data(), this->unit_point_ptr[qb]);
2703 this->
shapes.data() + qb * n_shapes,
2711template <
int n_components_,
int dim,
int spacedim,
typename Number>
2712template <
bool is_linear, std::
size_t str
ide_view>
2722 const unsigned int n_shapes = is_linear ? 2 : this->
poly.size();
2724 for (
unsigned int qb = 0; qb < this->
n_q_batches; ++qb)
2730 solution_values, evaluation_flags, n_shapes, qb,
value, gradient);
2734 for (
unsigned int v = 0, offset = qb *
stride;
2745 for (
unsigned int v = 0, offset = qb *
stride;
2750 ETT::set_gradient(gradient, v, unit_gradient);
2760 GeometryType::affine ?
2772template <
int n_components_,
int dim,
int spacedim,
typename Number>
2773template <std::
size_t str
ide_view>
2782 "Not initialized. Please call FEPointEvaluation::reinit()!"));
2784 const std::size_t n_points = this->
fe_values->get_quadrature().size();
2790 for (
unsigned int i = 0; i < this->
fe->n_dofs_per_cell(); ++i)
2795 (this->
fe->is_primitive(i) || this->fe->is_primitive()))
2796 for (
unsigned int qb = 0, q = 0; q < n_points;
2798 for (
unsigned int v = 0;
2801 ETT::access(this->
values[qb],
2806 for (
unsigned int qb = 0, q = 0; q < n_points;
2808 for (
unsigned int v = 0;
2811 ETT::access(this->
values[qb],
2814 this->
fe_values->shape_value_component(i,
2825 this->gradients.end(),
2827 for (
unsigned int i = 0; i < this->
fe->n_dofs_per_cell(); ++i)
2832 (this->
fe->is_primitive(i) || this->fe->is_primitive()))
2833 for (
unsigned int qb = 0, q = 0; q < n_points;
2835 for (
unsigned int v = 0;
2843 for (
unsigned int qb = 0, q = 0; q < n_points;
2845 for (
unsigned int v = 0;
2859template <
int n_components_,
int dim,
int spacedim,
typename Number>
2860template <
bool is_linear>
2864 const unsigned int n_shapes,
2865 const unsigned int qb,
2880 solution_values_vectorized_linear :
2882 this->unit_point_ptr[qb],
2889 this->
shapes.data() + qb * n_shapes,
2892 is_linear ? solution_values_vectorized_linear :
2894 this->unit_point_ptr[qb],
2900template <
int n_components_,
int dim,
int spacedim,
typename Number>
2901template <
bool is_linear, std::
size_t str
ide_view>
2906 const bool sum_into_values)
2909 for (
unsigned int i = 0; i < solution_values.
size(); ++i)
2910 solution_values[i] = 0;
2912 const unsigned int dofs_per_comp =
2915 for (
unsigned int comp = 0; comp <
n_components; ++comp)
2917 const std::size_t offset =
2920 if (is_linear || this->
renumber.empty())
2922 for (
unsigned int i = 0; i < dofs_per_comp; ++i)
2923 if (sum_into_values)
2924 solution_values[i + offset] +=
2925 ETT::sum_value(comp,
2927 *(solution_values_vectorized_linear + i) :
2930 solution_values[i + offset] =
2931 ETT::sum_value(comp,
2933 *(solution_values_vectorized_linear + i) :
2938 const unsigned int *renumber_ptr = this->
renumber.data() + offset;
2939 for (
unsigned int i = 0; i < dofs_per_comp; ++i)
2940 if (sum_into_values)
2941 solution_values[renumber_ptr[i]] +=
2944 solution_values[renumber_ptr[i]] =
2952template <
int n_components_,
int dim,
int spacedim,
typename Number>
2953template <
bool do_JxW,
bool is_linear, std::
size_t str
ide_view>
2958 const bool sum_into_values)
2961 if constexpr (
stride == 1)
2962 if (
const unsigned int n_filled_lanes =
2968 ETT::set_zero_value(this->
values.back(), v);
2971 ETT::set_zero_gradient(this->
gradients.back(), v);
2975 solution_values_vectorized_linear = {};
2978 const unsigned int n_shapes = is_linear ? 2 : this->
poly.size();
2980 const bool cartesian_cell =
2982 const bool affine_cell =
2984 for (
unsigned int qb = 0; qb < this->
n_q_batches; ++qb)
2990 for (
unsigned int v = 0, offset = qb *
stride;
2993 ETT::get_value(
value,
2999 for (
unsigned int v = 0, offset = qb *
stride;
3023 solution_values_vectorized_linear.data());
3028 solution_values_vectorized_linear.
data(),
3034template <
int n_components_,
int dim,
int spacedim,
typename Number>
3035template <
bool do_JxW, std::
size_t str
ide_view>
3040 const bool sum_into_values)
3045 "Not initialized. Please call FEPointEvaluation::reinit()!"));
3046 if (!sum_into_values)
3047 for (
unsigned int i = 0; i < solution_values.
size(); ++i)
3048 solution_values[i] = 0;
3050 const std::size_t n_points = this->
fe_values->get_quadrature().size();
3055 for (
unsigned int i = 0; i < this->
fe->n_dofs_per_cell(); ++i)
3059 (this->
fe->is_primitive(i) || this->fe->is_primitive()))
3060 for (
unsigned int qb = 0, q = 0; q < n_points;
3062 for (
unsigned int v = 0;
3065 solution_values[i] +=
3066 this->
fe_values->shape_value(i, q + v) *
3067 ETT::access(this->
values[qb], v, d) *
3068 (do_JxW ? this->
fe_values->JxW(q + v) : 1.);
3070 for (
unsigned int qb = 0, q = 0; q < n_points;
3072 for (
unsigned int v = 0;
3075 solution_values[i] +=
3076 this->
fe_values->shape_value_component(i, q + v, d) *
3077 ETT::access(this->
values[qb], v, d) *
3078 (do_JxW ? this->
fe_values->JxW(q + v) : 1.);
3085 for (
unsigned int i = 0; i < this->
fe->n_dofs_per_cell(); ++i)
3089 (this->
fe->is_primitive(i) || this->fe->is_primitive()))
3090 for (
unsigned int qb = 0, q = 0; q < n_points;
3092 for (
unsigned int v = 0;
3095 solution_values[i] +=
3097 ETT::access(this->
gradients[qb], v, d) *
3098 (do_JxW ? this->
fe_values->JxW(q + v) : 1.);
3100 for (
unsigned int qb = 0, q = 0; q < n_points;
3102 for (
unsigned int v = 0;
3105 solution_values[i] +=
3106 this->
fe_values->shape_grad_component(i, q + v, d) *
3107 ETT::access(this->
gradients[qb], v, d) *
3108 (do_JxW ? this->
fe_values->JxW(q + v) : 1.);
3115template <
int n_components_,
int dim,
int spacedim,
typename Number>
3116template <
bool do_JxW, std::
size_t str
ide_view>
3121 const bool sum_into_values)
3130 (integration_flags &
3133 if (!sum_into_values)
3134 for (
unsigned int i = 0; i < solution_values.
size(); ++i)
3135 solution_values[i] = 0;
3140 !do_JxW || this->
JxW_ptr !=
nullptr,
3142 "JxW pointer is not set! If you do not want to integrate() use test_and_sum()"));
3162template <
int n_components_,
int dim,
int spacedim,
typename Number>
3165 const unsigned int point_index)
const
3170 ExcFEPointEvaluationAccessToUninitializedMappingField(
3171 "update_normal_vectors"));
3177template <
int n_components_,
int dim,
int spacedim,
typename Number>
3183 const unsigned int point_index)
const
3191template <
int n_components_,
int dim,
int spacedim,
typename Number>
3195 const unsigned int point_index)
3207template <
int n_components_,
int dim,
int spacedim,
typename Number>
3213 const unsigned int first_selected_component)
3217 first_selected_component,
3223template <
int n_components_,
int dim,
int spacedim,
typename Number>
3226 const unsigned int cell_index,
3227 const unsigned int face_number)
3241template <
int n_components_,
int dim,
int spacedim,
typename Number>
3244 const unsigned int face_index)
3259template <
int n_components_,
int dim,
int spacedim,
typename Number>
3260template <std::
size_t str
ide_view>
3267 ExcMessage(
"Object has not been reinitialized!"));
3288template <
int n_components_,
int dim,
int spacedim,
typename Number>
3295 solution_values.
size()),
3301template <
int n_components_,
int dim,
int spacedim,
typename Number>
3302template <
bool is_linear, std::
size_t str
ide_view>
3308 const unsigned int dofs_per_comp =
3313 (is_linear || this->
renumber.empty()))
3314 input = solution_values.
data();
3317 for (
unsigned int comp = 0; comp <
n_components; ++comp)
3319 const std::size_t offset =
3322 if (is_linear || this->
renumber.empty())
3324 for (
unsigned int i = 0; i < dofs_per_comp; ++i)
3326 solution_values[i + offset];
3330 const unsigned int *renumber_ptr = this->
renumber.data() + offset;
3331 for (
unsigned int i = 0; i < dofs_per_comp; ++i)
3333 solution_values[renumber_ptr[i]];
3355template <
int n_components_,
int dim,
int spacedim,
typename Number>
3356template <std::
size_t str
ide_view>
3361 const bool sum_into_values)
3364 ExcMessage(
"Object has not been reinitialized!"));
3370 (integration_flags &
3373 if (!sum_into_values)
3374 for (
unsigned int i = 0; i < solution_values.
size(); ++i)
3375 solution_values[i] = 0;
3393template <
int n_components_,
int dim,
int spacedim,
typename Number>
3398 const bool sum_into_values)
3401 solution_values.
size()),
3408template <
int n_components_,
int dim,
int spacedim,
typename Number>
3409template <std::
size_t str
ide_view>
3414 const bool sum_into_values)
3417 ExcMessage(
"Object has not been reinitialized!"));
3423 (integration_flags &
3426 if (!sum_into_values)
3427 for (
unsigned int i = 0; i < solution_values.
size(); ++i)
3428 solution_values[i] = 0;
3446template <
int n_components_,
int dim,
int spacedim,
typename Number>
3451 const bool sum_into_values)
3454 solution_values.
size()),
3461template <
int n_components_,
int dim,
int spacedim,
typename Number>
3462template <
bool do_JxW,
bool is_linear, std::
size_t str
ide_view>
3467 const bool sum_into_values)
3470 for (
unsigned int i = 0; i < solution_values.
size(); ++i)
3471 solution_values[i] = 0;
3480 (is_linear || this->
renumber.empty()))
3482 if (sum_into_values)
3484 FEFaceNormalEvaluationImpl<dim, is_linear ? 1 : -1,
ScalarNumber>::
3489 solution_values.
data(),
3490 this->current_face_number);
3493 FEFaceNormalEvaluationImpl<dim, is_linear ? 1 : -1,
ScalarNumber>::
3498 solution_values.
data(),
3499 this->current_face_number);
3503 const unsigned int dofs_per_comp_face =
3506 const unsigned int size_input = 3 * dofs_per_comp_face *
n_components;
3510 FEFaceNormalEvaluationImpl<dim, is_linear ? 1 : -1,
ScalarNumber>::
3518 const unsigned int dofs_per_comp =
3521 for (
unsigned int comp = 0; comp <
n_components; ++comp)
3523 const std::size_t offset =
3526 if (is_linear || this->
renumber.empty())
3528 for (
unsigned int i = 0; i < dofs_per_comp; ++i)
3529 if (sum_into_values)
3530 solution_values[i + offset] +=
3531 output[i + comp * dofs_per_comp];
3533 solution_values[i + offset] =
3534 output[i + comp * dofs_per_comp];
3538 const unsigned int *renumber_ptr = this->
renumber.data() + offset;
3539 for (
unsigned int i = 0; i < dofs_per_comp; ++i)
3540 if (sum_into_values)
3541 solution_values[renumber_ptr[i]] +=
3542 output[i + comp * dofs_per_comp];
3544 solution_values[renumber_ptr[i]] =
3545 output[i + comp * dofs_per_comp];
3553template <
int n_components_,
int dim,
int spacedim,
typename Number>
3554template <
int str
ide_face_dof>
3570template <
int n_components_,
int dim,
int spacedim,
typename Number>
3571template <
bool is_linear,
int str
ide_face_dof>
3579 face_dof_values_ptr = face_dof_values;
3582 const unsigned int dofs_per_comp_face =
3584 for (
unsigned int comp = 0; comp <
n_components; ++comp)
3585 for (
unsigned int i = 0; i < 2 * dofs_per_comp_face; ++i)
3586 ETT::read_value(face_dof_values[(i + comp * 3 * dofs_per_comp_face) *
3594 constexpr int stride_face_dof_actual =
3598 const unsigned int n_shapes = is_linear ? 2 : this->
poly.size();
3600 for (
unsigned int qb = 0; qb < this->
n_q_batches; ++qb)
3607 const std::array<vectorized_value_type, dim + 1> interpolated_value =
3614 stride_face_dof_actual>(face_dof_values_ptr,
3625 face_dof_values_ptr);
3627 value = interpolated_value[dim - 1];
3632 gradient[0] = interpolated_value[dim];
3634 gradient[1] = interpolated_value[0];
3636 gradient[2] = interpolated_value[1];
3641 gradient[1] = interpolated_value[dim];
3644 gradient[0] = interpolated_value[1];
3645 gradient[2] = interpolated_value[0];
3648 gradient[0] = interpolated_value[0];
3656 gradient[0] = interpolated_value[0];
3657 gradient[1] = interpolated_value[1];
3658 gradient[2] = interpolated_value[dim];
3673 stride_face_dof_actual>(face_dof_values_ptr,
3683 face_dof_values_ptr);
3688 for (
unsigned int v = 0, offset = qb *
stride;
3699 for (
unsigned int v = 0, offset = qb *
stride;
3704 ETT::set_gradient(gradient, v, unit_gradient);
3714 GeometryType::affine ?
3726template <
int n_components_,
int dim,
int spacedim,
typename Number>
3727template <
int str
ide_face_dof>
3732 const bool sum_into_values)
3746template <
int n_components_,
int dim,
int spacedim,
typename Number>
3747template <
bool do_JxW,
bool is_linear,
int str
ide_face_dof>
3753 const bool sum_into_values)
3756 if constexpr (
stride == 1)
3757 if (
const unsigned int n_filled_lanes =
3763 ETT::set_zero_value(this->
values.back(), v);
3766 ETT::set_zero_gradient(this->
gradients.back(), v);
3771 solution_values_vectorized_linear = {};
3774 const unsigned int n_shapes = is_linear ? 2 : this->
poly.size();
3776 const bool cartesian_cell =
3778 const bool affine_cell =
3780 for (
unsigned int qb = 0; qb < this->
n_q_batches; ++qb)
3786 for (
unsigned int v = 0, offset = qb *
stride;
3789 ETT::get_value(
value,
3795 for (
unsigned int v = 0, offset = qb *
stride;
3815 std::array<vectorized_value_type, 2> value_face = {};
3818 value_face[0] =
value;
3823 value_face[1] = gradient[0];
3825 gradient_in_face[0] = gradient[1];
3827 gradient_in_face[1] = gradient[2];
3832 value_face[1] = gradient[1];
3835 gradient_in_face[0] = gradient[2];
3836 gradient_in_face[1] = gradient[0];
3839 gradient_in_face[0] = gradient[0];
3847 value_face[1] = gradient[2];
3848 gradient_in_face[0] = gradient[0];
3849 gradient_in_face[1] = gradient[1];
3866 is_linear ? solution_values_vectorized_linear.data() :
3867 this->solution_renumbered_vectorized.data(),
3868 this->unit_point_faces_ptr[qb],
3879 is_linear ? solution_values_vectorized_linear.data() :
3880 this->solution_renumbered_vectorized.data(),
3881 this->unit_point_faces_ptr[qb],
3885 const unsigned int dofs_per_comp_face =
3888 for (
unsigned int comp = 0; comp <
n_components; ++comp)
3889 for (
unsigned int i = 0; i < 2 * dofs_per_comp_face; ++i)
3890 if (sum_into_values)
3891 face_dof_values[(i + comp * 3 * dofs_per_comp_face) *
3893 ETT::sum_value(comp,
3895 *(solution_values_vectorized_linear.data() + i) :
3896 this->solution_renumbered_vectorized[i]);
3898 face_dof_values[(i + comp * 3 * dofs_per_comp_face) * stride_face_dof] =
3899 ETT::sum_value(comp,
3901 *(solution_values_vectorized_linear.data() + i) :
3902 this->solution_renumbered_vectorized[i]);
3907template <
int n_components_,
int dim,
int spacedim,
typename Number>
3910 const unsigned int point_index)
const
3915 ExcFEPointEvaluationAccessToUninitializedMappingField(
3916 "update_normal_vectors"));
3920 for (
unsigned int d = 0; d < dim; ++d)
3935template <
int n_components_,
int dim,
int spacedim,
typename Number>
3949template <
int n_components_,
int dim,
int spacedim,
typename Number>
3953 const unsigned int point_index)
value_type * data() const noexcept
typename ETT::vectorized_value_type vectorized_value_type
static constexpr std::size_t n_lanes_internal
void do_integrate_in_face(ScalarNumber *face_dof_values, const EvaluationFlags::EvaluationFlags &integration_flags, const bool sum_into_values)
void integrate_in_face(ScalarNumber *face_dof_values, const EvaluationFlags::EvaluationFlags &integration_flags, const bool sum_into_values=false)
typename ETT::scalar_value_type scalar_value_type
void evaluate(const StridedArrayView< const ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &evaluation_flags)
void do_evaluate_in_face(const ScalarNumber *face_dof_values, const EvaluationFlags::EvaluationFlags &evaluation_flags)
void submit_normal_derivative(const value_type &, const unsigned int point_index)
static constexpr std::size_t stride
void do_evaluate(const StridedArrayView< const ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &evaluation_flags)
static constexpr unsigned int n_components
static constexpr std::size_t n_lanes_user_interface
typename ETT::real_gradient_type gradient_type
void evaluate_in_face(const ScalarNumber *face_dof_values, const EvaluationFlags::EvaluationFlags &evaluation_flags)
void reinit(const unsigned int cell_index, const unsigned int face_number)
void test_and_sum(const StridedArrayView< ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &integration_flags, const bool sum_into_values=false)
typename internal::VectorizedArrayTrait< Number >::value_type ScalarNumber
typename ETT::interface_vectorized_unit_gradient_type interface_vectorized_unit_gradient_type
FEFacePointEvaluation(NonMatching::MappingInfo< dim, spacedim, Number > &mapping_info, const FiniteElement< dim, spacedim > &fe, const bool is_interior=true, const unsigned int first_selected_component=0)
typename ETT::unit_gradient_type unit_gradient_type
typename internal::FEPointEvaluation:: EvaluatorTypeTraits< dim, spacedim, n_components, Number > ETT
const value_type get_normal_derivative(const unsigned int point_index) const
void do_integrate(const StridedArrayView< ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &integration_flags, const bool sum_into_values)
void integrate(const StridedArrayView< ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &integration_flags, const bool sum_into_values=false)
static constexpr unsigned int dimension
Tensor< 1, spacedim, Number > normal_vector(const unsigned int point_index) const
typename ETT::value_type value_type
typename ::internal::VectorizedArrayTrait< Number >::vectorized_value_type VectorizedArrayType
const DerivativeForm< 1, dim, spacedim, Number > * jacobian_ptr
std::unique_ptr< NonMatching::MappingInfo< dim, spacedim, Number > > mapping_info_on_the_fly
std::vector< gradient_type > gradients
unsigned int dofs_per_component
Number get_divergence(const unsigned int point_index) const
unsigned int dofs_per_component_face
Number JxW(const unsigned int point_index) const
const UpdateFlags update_flags
static constexpr std::size_t n_lanes_user_interface
static constexpr std::size_t stride
internal::MatrixFreeFunctions::GeometryType cell_type
std::vector< Polynomials::Polynomial< double > > poly
Point< spacedim, Number > real_point(const unsigned int point_index) const
typename internal::VectorizedArrayTrait< Number >::value_type ScalarNumber
const unsigned int n_q_points_scalar
const Point< dim, VectorizedArrayType > * unit_point_ptr
typename internal::FEPointEvaluation:: EvaluatorTypeTraits< dim, spacedim, n_components, Number > ETT
FEPointEvaluationBase(FEPointEvaluationBase &&other) noexcept
AlignedVector< ScalarNumber > scratch_data_scalar
std_cxx20::ranges::iota_view< unsigned int, unsigned int > quadrature_point_indices() const
const value_type & get_value(const unsigned int point_index) const
const Point< dim - 1, VectorizedArrayType > * unit_point_faces_ptr
std::vector< scalar_value_type > solution_renumbered
SmartPointer< NonMatching::MappingInfo< dim, spacedim, double > > mapping_info
const unsigned int n_q_batches
std::vector< std::array< bool, n_components > > nonzero_shape_function_component
Point< spacedim, Number > quadrature_point(const unsigned int point_index) const
unsigned int n_active_entries_per_quadrature_batch(unsigned int q)
const gradient_type & get_gradient(const unsigned int point_index) const
std::shared_ptr< FEValues< dim, spacedim > > fe_values
FEPointEvaluationBase(const Mapping< dim, spacedim > &mapping, const FiniteElement< dim, spacedim > &fe, const UpdateFlags update_flags, const unsigned int first_selected_component=0)
SmartPointer< const FiniteElement< dim, spacedim > > fe
std::vector< unsigned int > renumber
std::vector< value_type > values
AlignedVector<::ndarray< VectorizedArrayType, 2, dim - 1 > > shapes_faces
Point< dim, Number > unit_point(const unsigned int point_index) const
const unsigned int n_q_points
void submit_divergence(const Number &value, const unsigned int point_index)
bool must_reinitialize_pointers
const Tensor< 1, spacedim, Number > * normal_ptr
DerivativeForm< 1, spacedim, dim, Number > inverse_jacobian(const unsigned int point_index) const
AlignedVector< vectorized_value_type > solution_renumbered_vectorized
static constexpr std::size_t n_lanes_internal
const DerivativeForm< 1, spacedim, dim, Number > * inverse_jacobian_ptr
AlignedVector<::ndarray< VectorizedArrayType, 2, dim > > shapes
typename ETT::value_type value_type
unsigned int current_cell_index
scalar_value_type integrate_value() const
void submit_gradient(const gradient_type &, const unsigned int point_index)
void setup(const unsigned int first_selected_component)
FEPointEvaluationBase(NonMatching::MappingInfo< dim, spacedim, Number > &mapping_info, const FiniteElement< dim, spacedim > &fe, const unsigned int first_selected_component=0, const bool is_interior=true)
typename ETT::scalar_value_type scalar_value_type
static constexpr unsigned int dimension
typename ETT::vectorized_value_type vectorized_value_type
static constexpr unsigned int n_components
unsigned int current_face_number
typename ::internal::VectorizedArrayTrait< Number >::vectorized_value_type VectorizedArrayType
FEPointEvaluationBase(FEPointEvaluationBase &other) noexcept
DerivativeForm< 1, dim, spacedim, Number > jacobian(const unsigned int point_index) const
typename ETT::real_gradient_type gradient_type
const Point< spacedim, Number > * real_point_ptr
unsigned int component_in_base_element
SmartPointer< const Mapping< dim, spacedim > > mapping
internal::MatrixFreeFunctions::ShapeInfo< ScalarNumber > shape_info
void submit_value(const value_type &value, const unsigned int point_index)
typename ETT::interface_vectorized_unit_gradient_type interface_vectorized_unit_gradient_type
static constexpr std::size_t stride
typename ETT::interface_vectorized_unit_gradient_type interface_vectorized_unit_gradient_type
void integrate_slow(const StridedArrayView< ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &integration_flags, const bool sum_into_values)
static constexpr unsigned int dimension
typename ETT::value_type value_type
void do_integrate(const StridedArrayView< ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &integration_flags, const bool sum_into_values)
void evaluate_fast(const StridedArrayView< const ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &evaluation_flags)
typename ETT::scalar_value_type scalar_value_type
typename ETT::unit_gradient_type unit_gradient_type
void evaluate(const StridedArrayView< const ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &evaluation_flags)
const value_type get_normal_derivative(const unsigned int point_index) const
void prepare_evaluate_fast(const StridedArrayView< const ScalarNumber, stride_view > &solution_values)
void test_and_sum(const StridedArrayView< ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &integration_flags, const bool sum_into_values=false)
void integrate(const StridedArrayView< ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &integration_flags, const bool sum_into_values=false)
typename ETT::vectorized_value_type vectorized_value_type
void compute_integrate_fast(const EvaluationFlags::EvaluationFlags &integration_flags, const unsigned int n_shapes, const unsigned int qb, const vectorized_value_type value, const interface_vectorized_unit_gradient_type gradient, vectorized_value_type *solution_values_vectorized_linear)
typename ETT::real_gradient_type gradient_type
FEPointEvaluation(const Mapping< dim, spacedim > &mapping, const FiniteElement< dim, spacedim > &fe, const UpdateFlags update_flags, const unsigned int first_selected_component=0)
typename internal::FEPointEvaluation:: EvaluatorTypeTraits< dim, spacedim, n_components, Number > ETT
void submit_normal_derivative(const value_type &, const unsigned int point_index)
void evaluate_slow(const StridedArrayView< const ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &evaluation_flags)
static constexpr std::size_t n_lanes_internal
static constexpr unsigned int n_components
void finish_integrate_fast(const StridedArrayView< ScalarNumber, stride_view > &solution_values, vectorized_value_type *solution_values_vectorized_linear, const bool sum_into_values)
void internal_reinit_single_cell_state_mapping_info()
void compute_evaluate_fast(const StridedArrayView< const ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &evaluation_flags, const unsigned int n_shapes, const unsigned int qb, vectorized_value_type &value, interface_vectorized_unit_gradient_type &gradient)
void integrate_fast(const StridedArrayView< ScalarNumber, stride_view > &solution_values, const EvaluationFlags::EvaluationFlags &integration_flags, const bool sum_into_values)
static constexpr std::size_t n_lanes_user_interface
typename ::internal::VectorizedArrayTrait< Number >::vectorized_value_type VectorizedArrayType
Tensor< 1, spacedim, Number > normal_vector(const unsigned int point_index) const
typename internal::VectorizedArrayTrait< Number >::value_type ScalarNumber
Abstract base class for mapping classes.
value_type * data() const noexcept
Tensor< rank, dim, Number > sum(const Tensor< rank, dim, Number > &local, const MPI_Comm mpi_communicator)
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
static ::ExceptionBase & ExcNotImplemented()
#define Assert(cond, exc)
#define AssertDimension(dim1, dim2)
static ::ExceptionBase & ExcFEPointEvaluationAccessToUninitializedMappingField(std::string arg1)
#define AssertIndexRange(index, range)
static ::ExceptionBase & ExcNotInitialized()
#define DeclException1(Exception1, type1, outsequence)
static ::ExceptionBase & ExcMessage(std::string arg1)
#define AssertThrow(cond, exc)
TriaIterator< CellAccessor< dim, spacedim > > cell_iterator
const bool IsBlockVector< VectorType >::value
@ update_values
Shape function values.
@ update_normal_vectors
Normal vectors.
@ update_JxW_values
Transformed quadrature weights.
@ update_jacobians
Volume element.
@ update_inverse_jacobians
Volume element.
@ update_gradients
Shape function gradients.
@ update_quadrature_points
Transformed quadrature points.
#define DEAL_II_ASSERT_UNREACHABLE()
MappingQ< dim, spacedim > StaticMappingQ1< dim, spacedim >::mapping
EvaluationFlags
The EvaluationFlags enum.
constexpr T pow(const T base, const int iexp)
std::vector< Polynomials::Polynomial< double > > get_polynomial_space(const FiniteElement< dim, spacedim > &fe)
bool is_fast_path_supported(const FiniteElement< dim, spacedim > &fe, const unsigned int base_element_number)
std::array< typename ProductTypeNoPoint< Number, Number2 >::type, dim+n_values > evaluate_tensor_product_value_and_gradient_shapes(const ::ndarray< Number2, 2, dim > *shapes, const int n_shapes, const Number *values, const std::vector< unsigned int > &renumber={})
void integrate_tensor_product_value_and_gradient(const ::ndarray< Number, 2, dim > *shapes, const unsigned int n_shapes, const Number2 *value, const Tensor< 1, dim, Number2 > &gradient, Number2 *values, const Point< dim, Number > &p, const bool do_add)
void compute_values_of_array_in_pairs(::ndarray< Number, 2, dim > *shapes, const std::vector< Polynomials::Polynomial< double > > &poly, const Point< dim, Number > &p0, const Point< dim, Number > &p1)
ProductTypeNoPoint< Number, Number2 >::type evaluate_tensor_product_value_shapes(const ::ndarray< Number2, 2, dim > *shapes, const int n_shapes, const Number *values, const std::vector< unsigned int > &renumber={})
void compute_values_of_array(::ndarray< Number, 2, dim > *shapes, const std::vector< Polynomials::Polynomial< double > > &poly, const Point< dim, Number > &p, const unsigned int derivative=1)
ProductTypeNoPoint< Number, Number2 >::type evaluate_tensor_product_value_linear(const Number *values, const Point< dim, Number2 > &p)
std::array< typename ProductTypeNoPoint< Number, Number2 >::type, dim+n_values > evaluate_tensor_product_value_and_gradient_linear(const Number *values, const Point< dim, Number2 > &p)
void integrate_tensor_product_value(const ::ndarray< Number, 2, dim > *shapes, const unsigned int n_shapes, const Number2 &value, Number2 *values, const Point< dim, Number > &p, const bool do_add)
static const unsigned int invalid_unsigned_int
boost::integer_range< IncrementableType > iota_view
typename internal::ndarray::HelperArray< T, Ns... >::type ndarray
static void set_zero_gradient(unit_gradient_type &value, const unsigned int vector_lane)
static Tensor< 1, spacedim, ScalarNumber > access(const real_gradient_type &value, const unsigned int vector_lane, const unsigned int component)
static void get_gradient(interface_vectorized_unit_gradient_type &value, const unsigned int vector_lane, const DerivativeForm< 1, dim, n_components, Number > &result)
static void get_value(vectorized_value_type &value, const unsigned int, const vectorized_value_type &result)
static void get_value(vectorized_value_type &value, const unsigned int vector_lane, const scalar_value_type &result)
Tensor< 1, dim, Tensor< 1, dim, VectorizedArrayType > > interface_vectorized_unit_gradient_type
static void set_value(const vectorized_value_type &value, const unsigned int, vectorized_value_type &result)
static ScalarNumber access(const value_type &value, const unsigned int vector_lane, const unsigned int)
Tensor< 1, dim, VectorizedArrayType > vectorized_value_type
typename internal::VectorizedArrayTrait< Number >::value_type ScalarNumber
static void get_gradient(vectorized_unit_gradient_type &value, const unsigned int, const vectorized_unit_gradient_type &result)
static scalar_value_type sum_value(const scalar_value_type &result)
static void access(value_type &value, const unsigned int vector_lane, const unsigned int component, const ScalarNumber &shape_value)
static scalar_value_type sum_value(const scalar_value_type &result)
Tensor< 1, n_components, VectorizedArrayType > vectorized_value_type
static void get_value(vectorized_value_type &value, const unsigned int, const vectorized_value_type &result)
Tensor< 1, dim, Tensor< 1, n_components, VectorizedArrayType > > interface_vectorized_unit_gradient_type
static void set_gradient(const vectorized_unit_gradient_type &value, const unsigned int vector_lane, scalar_unit_gradient_type &result)
static void set_value(const vectorized_value_type &value, const unsigned int vector_lane, scalar_value_type &result)
static void set_zero_value(value_type &value, const unsigned int vector_lane)
Tensor< 2, dim, Number > unit_gradient_type
static Tensor< 1, spacedim, ScalarNumber > access(const real_gradient_type &value, const unsigned int vector_lane, const unsigned int)
static scalar_value_type sum_value(const vectorized_value_type &result)
static Tensor< 1, dim, ScalarNumber > access(const real_gradient_type &value, const unsigned int vector_lane, const unsigned int component)
static void set_value(const vectorized_value_type &value, const unsigned int, vectorized_value_type &result)
static void read_value(const ScalarNumber vector_entry, const unsigned int component, scalar_value_type &result)
static void get_value(vectorized_value_type &value, const unsigned int, const vectorized_value_type &result)
VectorizedArrayType vectorized_value_type
Tensor< 1, n_components, Number > value_type
static void set_zero_value(value_type &value, const unsigned int vector_lane)
Tensor< 1, n_components, ScalarNumber > scalar_value_type
Tensor< 1, dim, VectorizedArrayType > vectorized_unit_gradient_type
static void set_gradient(const interface_vectorized_unit_gradient_type &value, const unsigned int vector_lane, unit_gradient_type &result)
std::conditional_t< n_components==spacedim, Tensor< 2, spacedim, Number >, Tensor< 1, n_components, Tensor< 1, spacedim, Number > > > real_gradient_type
Tensor< 2, dim, VectorizedArrayType > vectorized_unit_gradient_type
Tensor< 1, n_components, Tensor< 1, dim, VectorizedArrayType > > vectorized_unit_gradient_type
static void get_value(vectorized_value_type &value, const unsigned int vector_lane, const scalar_value_type &result)
static void read_value(const ScalarNumber vector_entry, const unsigned int component, scalar_value_type &result)
Tensor< 2, dim, ScalarNumber > scalar_unit_gradient_type
static void set_value(const vectorized_value_type &value, const unsigned int vector_lane, scalar_value_type &result)
Tensor< 1, dim, ScalarNumber > scalar_unit_gradient_type
typename ::internal::VectorizedArrayTrait< Number >::vectorized_value_type VectorizedArrayType
typename internal::VectorizedArrayTrait< Number >::value_type ScalarNumber
static scalar_value_type sum_value(const vectorized_value_type &result)
static void access(real_gradient_type &value, const unsigned int vector_lane, const unsigned int, const Tensor< 1, spacedim, ScalarNumber > &shape_gradient)
vectorized_unit_gradient_type interface_vectorized_unit_gradient_type
static ScalarNumber access(const value_type &value, const unsigned int vector_lane, const unsigned int component)
static ScalarNumber sum_value(const unsigned int component, const vectorized_value_type &result)
static void access(real_gradient_type &value, const unsigned int vector_lane, const unsigned int component, const Tensor< 1, dim, ScalarNumber > &shape_gradient)
static void set_zero_value(value_type &value, const unsigned int vector_lane)
static scalar_value_type sum_value(const vectorized_value_type &result)
Tensor< 1, dim, Number > unit_gradient_type
static void set_gradient(const vectorized_unit_gradient_type &value, const unsigned int, vectorized_unit_gradient_type &result)
static scalar_value_type sum_value(const scalar_value_type &result)
static ScalarNumber sum_value(const unsigned int component, const vectorized_value_type &result)
static ScalarNumber access(const value_type &value, const unsigned int vector_lane, const unsigned int component)
typename ::internal::VectorizedArrayTrait< Number >::vectorized_value_type VectorizedArrayType
Tensor< 1, n_components, Tensor< 1, dim, ScalarNumber > > scalar_unit_gradient_type
static void get_value(vectorized_value_type &value, const unsigned int vector_lane, const scalar_value_type &result)
ScalarNumber scalar_value_type
static void set_gradient(const interface_vectorized_unit_gradient_type &value, const unsigned int vector_lane, unit_gradient_type &result)
Tensor< 1, dim, ScalarNumber > scalar_value_type
static void set_zero_gradient(real_gradient_type &value, const unsigned int vector_lane)
Tensor< 1, spacedim, Number > real_gradient_type
static void set_value(const vectorized_value_type &value, const unsigned int vector_lane, scalar_value_type &result)
static void get_gradient(interface_vectorized_unit_gradient_type &value, const unsigned int vector_lane, const unit_gradient_type &result)
static void get_gradient(interface_vectorized_unit_gradient_type &value, const unsigned int vector_lane, const unit_gradient_type &result)
static void set_value(const vectorized_value_type &value, const unsigned int, vectorized_value_type &result)
static void set_zero_gradient(real_gradient_type &value, const unsigned int vector_lane)
static void access(value_type &value, const unsigned int vector_lane, const unsigned int, const ScalarNumber &shape_value)
static void get_gradient(vectorized_unit_gradient_type &value, const unsigned int vector_lane, const scalar_unit_gradient_type &result)
unit_gradient_type real_gradient_type
static void access(real_gradient_type &value, const unsigned int vector_lane, const unsigned int component, const Tensor< 1, spacedim, ScalarNumber > &shape_gradient)
Tensor< 1, dim, Number > value_type
static ScalarNumber sum_value(const unsigned int, const vectorized_value_type &result)
typename internal::VectorizedArrayTrait< Number >::value_type ScalarNumber
Tensor< 1, n_components, Tensor< 1, dim, Number > > unit_gradient_type
static void access(value_type &value, const unsigned int vector_lane, const unsigned int component, const ScalarNumber &shape_value)
typename ::internal::VectorizedArrayTrait< Number >::vectorized_value_type VectorizedArrayType
static void read_value(const ScalarNumber vector_entry, const unsigned int, scalar_value_type &result)
static constexpr std::size_t width()
static constexpr std::size_t stride()
static value_type & get_from_vectorized(vectorized_value_type &values, unsigned int c)
static value_type & get(value_type &value, unsigned int c)
DEAL_II_HOST constexpr SymmetricTensor< 4, dim, Number > outer_product(const SymmetricTensor< 2, dim, Number > &t1, const SymmetricTensor< 2, dim, Number > &t2)
DEAL_II_HOST constexpr Number trace(const SymmetricTensor< 2, dim2, Number > &)