cout << " Assembling matrices..." << endl;
UpdateFlags update_flags = UpdateFlags(update_q_points | update_gradients |
- update_jacobians | update_JxW_values);
+ update_JxW_values);
ProblemBase<dim>::FunctionMap dirichlet_bc;
dirichlet_bc[0] = boundary_values;
cout << " Assembling matrices..." << endl;
UpdateFlags update_flags = UpdateFlags(update_q_points | update_gradients |
- update_jacobians | update_JxW_values);
+ update_JxW_values);
ProblemBase<dim>::FunctionMap dirichlet_bc;
dirichlet_bc[0] = solution_function;
dirichlet_bc[0] = boundary_values;
assemble (equation, quadrature, fe,
UpdateFlags(update_gradients | update_JxW_values |
- update_jacobians | update_q_points),
+ update_q_points),
dirichlet_bc, *boundary);
cout << " Solving..." << endl;
FEFaceValues<dim> fe_face_values_cell (fe, quadrature,
UpdateFlags(update_gradients |
update_JxW_values |
- update_jacobians |
update_normal_vectors));
- FEFaceValues<dim> fe_face_values_neighbor (fe, quadrature,
- UpdateFlags(update_gradients |
- update_jacobians));
- FESubfaceValues<dim> fe_subface_values (fe, quadrature,
- UpdateFlags(update_gradients |
- update_jacobians));
+ FEFaceValues<dim> fe_face_values_neighbor (fe, quadrature, update_gradients);
+ FESubfaceValues<dim> fe_subface_values (fe, quadrature, update_gradients);
// loop over all cells
const active_cell_iterator endc = dof.end();
dSMatrix &matrix,
const Function<dim> * const a) {
dVector dummy; // no entries, should give an error if accessed
- UpdateFlags update_flags = UpdateFlags(update_jacobians |
- update_JxW_values);
+ UpdateFlags update_flags = update_JxW_values;
if (a != 0)
update_flags = UpdateFlags (update_flags | update_q_points);
const AssemblerData<dim> data (dof,
dVector &rhs_vector,
const Function<dim> * const a) {
UpdateFlags update_flags = UpdateFlags(update_q_points |
- update_jacobians |
update_JxW_values);
const AssemblerData<dim> data (dof,
true, true,
const Function<dim> * const a) {
dVector dummy; // no entries, should give an error if accessed
UpdateFlags update_flags = UpdateFlags(update_gradients |
- update_jacobians |
update_JxW_values);
if (a != 0)
update_flags = UpdateFlags(update_flags | update_q_points);
const Function<dim> * const a) {
UpdateFlags update_flags = UpdateFlags(update_q_points |
update_gradients |
- update_jacobians |
update_JxW_values);
const AssemblerData<dim> data (dof,
true, true,
difference.reinit (dof.get_tria().n_active_cells());
UpdateFlags update_flags = UpdateFlags (update_q_points |
- update_jacobians |
update_JxW_values);
if ((norm==H1_seminorm) || (norm==H1_norm))
update_flags = UpdateFlags (update_flags | update_gradients);
cout << " Assembling matrices..." << endl;
UpdateFlags update_flags = UpdateFlags(update_q_points | update_gradients |
- update_jacobians | update_JxW_values);
+ update_JxW_values);
ProblemBase<dim>::FunctionMap dirichlet_bc;
dirichlet_bc[0] = boundary_values;
cout << " Assembling matrices..." << endl;
UpdateFlags update_flags = UpdateFlags(update_q_points | update_gradients |
- update_jacobians | update_JxW_values);
+ update_JxW_values);
ProblemBase<dim>::FunctionMap dirichlet_bc;
dirichlet_bc[0] = solution_function;
dirichlet_bc[0] = boundary_values;
assemble (equation, quadrature, fe,
UpdateFlags(update_gradients | update_JxW_values |
- update_jacobians | update_q_points),
+ update_q_points),
dirichlet_bc, *boundary);
cout << " Solving..." << endl;