fe_values_2.n_quadrature_points;
// Initialize counters and renumbering vectors to zero
- unsigned int first_n_dofs = 0;
- unsigned int second_n_dofs = 0;
-
- unsigned int first_n_quadrature_points = 0;
- unsigned int second_n_quadrature_points = 0;
-
std::vector<unsigned int> first_dofs_map;
std::vector<unsigned int> second_dofs_map;
n_coupling_dofs_ = numbers::invalid_unsigned_int;
first_dofs_map.clear();
second_dofs_map.clear();
- first_n_dofs = first_fe_values->dofs_per_cell;
- second_n_dofs = second_fe_values->dofs_per_cell;
break;
}
case DoFCouplingType::contiguous:
n_coupling_dofs_ =
fe_values_1.dofs_per_cell + fe_values_2.dofs_per_cell;
- first_n_dofs = n_coupling_dofs_;
- second_n_dofs = n_coupling_dofs_;
-
first_dofs_map.resize(n_coupling_dofs_);
second_dofs_map.resize(n_coupling_dofs_);
n_quadrature_points_ =
quadrature_points_1.size() * quadrature_points_2.size();
- first_n_quadrature_points = n_quadrature_points_;
- second_n_quadrature_points = n_quadrature_points_;
-
first_quad_map.resize(n_quadrature_points_);
second_quad_map.resize(n_quadrature_points_);
n_quadrature_points_ = fe_values_1.get_quadrature_points().size();
- first_n_quadrature_points = n_quadrature_points_;
- second_n_quadrature_points = n_quadrature_points_;
-
first_quad_map.clear();
second_quad_map.clear();
n_quadrature_points_ = fe_values_1.get_quadrature_points().size();
- first_n_quadrature_points = n_quadrature_points_;
- second_n_quadrature_points = n_quadrature_points_;
-
first_quad_map.clear();
second_quad_map.clear();
n_quadrature_points_ = fe_values_1.get_quadrature_points().size();
- first_n_quadrature_points = n_quadrature_points_;
- second_n_quadrature_points = n_quadrature_points_;
-
// The first is the id. The second is renumbered.
first_quad_map.clear();
second_quad_map.resize(fe_values_1.get_quadrature_points().size());
}
n_quadrature_points_ = first_quad_map.size();
- first_n_quadrature_points = n_quadrature_points_;
- second_n_quadrature_points = n_quadrature_points_;
break;
}
default: