std::vector<std_cxx11::shared_ptr<dealii::FEValues<dim> > >
fe_values (current_data.quadrature.size());
UpdateFlags update_flags_feval =
- (update_flags & update_inverse_jacobians ? update_jacobians : update_default) |
- (update_flags & update_jacobian_grads ? update_jacobian_grads : update_default) |
- (update_flags & update_quadrature_points ? update_quadrature_points : update_default);
+ ((update_flags & update_inverse_jacobians) ? update_jacobians : update_default) |
+ ((update_flags & update_jacobian_grads) ? update_jacobian_grads : update_default) |
+ ((update_flags & update_quadrature_points) ? update_quadrature_points : update_default);
// resize the fields that have fixed size or for which we know
// something from an earlier loop
const unsigned int n_q_points = mapping_data.quadrature_points.size();
- std::vector<double> values(fe_internal.update_each & update_values ? this->dofs_per_cell : 0);
- std::vector<Tensor<1,dim> > grads(fe_internal.update_each & update_gradients ? this->dofs_per_cell : 0);
- std::vector<Tensor<2,dim> > grad_grads(fe_internal.update_each & update_hessians ? this->dofs_per_cell : 0);
+ std::vector<double> values((fe_internal.update_each & update_values) ? this->dofs_per_cell : 0);
+ std::vector<Tensor<1,dim> > grads((fe_internal.update_each & update_gradients) ? this->dofs_per_cell : 0);
+ std::vector<Tensor<2,dim> > grad_grads((fe_internal.update_each & update_hessians) ? this->dofs_per_cell : 0);
std::vector<Tensor<3,dim> > empty_vector_of_3rd_order_tensors;
std::vector<Tensor<4,dim> > empty_vector_of_4th_order_tensors;
const unsigned int n_q_points = mapping_data.quadrature_points.size();
- std::vector<double> values(fe_internal.update_each & update_values ? this->dofs_per_cell : 0);
- std::vector<Tensor<1,dim> > grads(fe_internal.update_each & update_gradients ? this->dofs_per_cell : 0);
- std::vector<Tensor<2,dim> > grad_grads(fe_internal.update_each & update_hessians ? this->dofs_per_cell : 0);
+ std::vector<double> values((fe_internal.update_each & update_values) ? this->dofs_per_cell : 0);
+ std::vector<Tensor<1,dim> > grads((fe_internal.update_each & update_gradients) ? this->dofs_per_cell : 0);
+ std::vector<Tensor<2,dim> > grad_grads((fe_internal.update_each & update_hessians) ? this->dofs_per_cell : 0);
std::vector<Tensor<3,dim> > empty_vector_of_3rd_order_tensors;
std::vector<Tensor<4,dim> > empty_vector_of_4th_order_tensors;
const unsigned int n_q_points = mapping_data.quadrature_points.size();
- std::vector<double> values(fe_internal.update_each & update_values ? this->dofs_per_cell : 0);
- std::vector<Tensor<1,dim> > grads(fe_internal.update_each & update_gradients ? this->dofs_per_cell : 0);
- std::vector<Tensor<2,dim> > grad_grads(fe_internal.update_each & update_hessians ? this->dofs_per_cell : 0);
+ std::vector<double> values((fe_internal.update_each & update_values) ? this->dofs_per_cell : 0);
+ std::vector<Tensor<1,dim> > grads((fe_internal.update_each & update_gradients) ? this->dofs_per_cell : 0);
+ std::vector<Tensor<2,dim> > grad_grads((fe_internal.update_each & update_hessians) ? this->dofs_per_cell : 0);
std::vector<Tensor<3,dim> > empty_vector_of_3rd_order_tensors;
std::vector<Tensor<4,dim> > empty_vector_of_4th_order_tensors;