1 #ifndef _RHEOLEF_INIT_EXPR_QUADRATURE_H
2 #define _RHEOLEF_INIT_EXPR_QUADRATURE_H
30 #include "rheolef/form_expr_variational.h"
32 namespace rheolef {
namespace details {
34 template<
class T,
class M>
44 if (iopt.
get_order() == std::numeric_limits<quadrature_option::size_type>::max()) {
45 size_type k = X.get_constitution().degree_max();
51 template<
class T,
class M>
62 if (new_iopt.
get_order() == std::numeric_limits<size_type>::max()) {
63 size_type k1 = X.get_constitution().degree_max();
64 size_type k2 = Y.get_constitution().degree_max();
66 if (omega_K.get_background_geo().sizes().ownership_by_variant[
reference_element::q].dis_size() != 0 ||
67 omega_K.get_background_geo().sizes().ownership_by_variant[
reference_element::P].dis_size() != 0 ||
68 omega_K.get_background_geo().sizes().ownership_by_variant[
reference_element::H].dis_size() != 0) {
73 if (quad_order >= n_derivative) quad_order -= n_derivative;
89 #endif // _RHEOLEF_INIT_EXPR_QUADRATURE_H