for (const auto &cell : dof_handler.active_cell_iterators())
{
fe_values.reinit(cell);
+
local_matrix = 0;
local_rhs = 0;
for (const auto &cell : dof_handler.active_cell_iterators())
{
fe_values.reinit(cell);
+
local_matrix = 0;
local_preconditioner_matrix = 0;
local_rhs = 0;
for (const auto &cell : dof_handler.active_cell_iterators())
{
fe_values.reinit(cell);
+
cell_matrix = 0;
cell_rhs = 0;
for (const auto &cell : dof_handler.active_cell_iterators())
if (cell->is_locally_owned())
{
+ fe_values.reinit(cell);
+
cell_matrix = 0.;
cell_rhs = 0.;
- fe_values.reinit(cell);
-
for (unsigned int q_point = 0; q_point < n_q_points; ++q_point)
{
const double rhs_value =
// coefficient value at each quadrature point.
for (const auto &cell : dof_handler.active_cell_iterators())
{
+ fe_values.reinit(cell);
+
cell_matrix = 0.;
cell_rhs = 0.;
- fe_values.reinit(cell);
-
for (const unsigned int q_index : fe_values.quadrature_point_indices())
{
const double current_coefficient =
for (const auto &cell : dof_handler.active_cell_iterators())
{
+ fe_values.reinit(cell);
+
cell_matrix = 0;
cell_rhs = 0;
- fe_values.reinit(cell);
-
for (const unsigned int q_index : fe_values.quadrature_point_indices())
{
const double current_coefficient =
// previous examples, so we only comment on the things that have changed.
for (const auto &cell : dof_handler.active_cell_iterators())
{
+ fe_values.reinit(cell);
+
cell_matrix = 0.;
cell_rhs = 0.;
- fe_values.reinit(cell);
-
right_hand_side.value_list(fe_values.get_quadrature_points(),
rhs_values);
// Now we can begin with the loop over all cells:
for (const auto &cell : dof_handler.active_cell_iterators())
{
+ fe_values.reinit(cell);
+
cell_matrix = 0;
cell_rhs = 0;
- fe_values.reinit(cell);
-
// Next we get the values of the coefficients at the quadrature
// points. Likewise for the right hand side:
lambda.value_list(fe_values.get_quadrature_points(), lambda_values);