-//#define WOLFGANG
-
-
template <class DH, class SparsityPattern>
void
DoFTools::make_sparsity_pattern (const DH &dof,
if ((face_constraints(row,i) != 0)
&&
(std::fabs(face_constraints(row,i)) >= 1e-14*abs_sum))
- {
-#ifdef WOLFGANG
- std::cout << " " << slave_dofs[row]
- << " -> " << face_constraints (row,i) << " * "
- << master_dofs[i]
- << std::endl;
-#endif
- constraints.add_entry (slave_dofs[row],
- master_dofs[i],
- face_constraints (row,i));
- }
+ constraints.add_entry (slave_dofs[row],
+ master_dofs[i],
+ face_constraints (row,i));
constraints.set_inhomogeneity (slave_dofs[row], 0.);
}
}
// Add constraints to global constraint
// matrix.
-#ifdef WOLFGANG
- std::cout << "Constraints for cell=" << cell
- << ", face=" << face
- << ", subface=" << c
- << std::endl;
-#endif
-
filter_constraints (master_dofs,
slave_dofs,
*(subface_interpolation_matrices
AssertDimension (slave_dofs.size(),
cell->get_fe().dofs_per_face - dominating_fe.dofs_per_face);
-#ifdef WOLFGANG
- std::cout << "Constraints for cell=" << cell
- << ", face=" << face
- << " (complicated case, mother)"
- << std::endl;
-#endif
filter_constraints (master_dofs,
slave_dofs,
constraint_matrix,
cell->face(face)->child(sf)->get_dof_indices (slave_dofs,
subface_fe_index);
-#ifdef WOLFGANG
- std::cout << "Constraints for cell=" << cell
- << ", face=" << face
- << ", subface=" << sf
- << " (complicated case, children)"
- << std::endl;
-#endif
filter_constraints (master_dofs,
slave_dofs,
constraint_matrix,
// Add constraints to global constraint
// matrix.
-#ifdef WOLFGANG
- std::cout << "p-constraints for cell=" << cell
- << ", face=" << face << std::endl;
-#endif
filter_constraints (master_dofs,
slave_dofs,
*(face_interpolation_matrices