const unsigned int degree,
const BooleanOperation & bool_op,
const Mapping<dim0, spacedim> &mapping0 =
- (dealii::ReferenceCells::get_hypercube<dim0>()
+ (ReferenceCells::get_hypercube<dim0>()
.template get_default_linear_mapping<dim0, spacedim>()),
const Mapping<dim1, spacedim> &mapping1 =
- (dealii::ReferenceCells::get_hypercube<dim1>()
+ (ReferenceCells::get_hypercube<dim1>()
.template get_default_linear_mapping<dim1, spacedim>()));
/**
const ReferenceCell reference_cell =
ReferenceCell::n_vertices_to_type(dim, vertices.size());
- if (reference_cell == dealii::ReferenceCells::Line)
+ if (reference_cell == ReferenceCells::Line)
// Return the distance between the two vertices
return (vertices[1] - vertices[0]).norm();
- else if (reference_cell == dealii::ReferenceCells::Triangle)
+ else if (reference_cell == ReferenceCells::Triangle)
// Return the longest of the three edges
return std::max({(vertices[1] - vertices[0]).norm(),
(vertices[2] - vertices[1]).norm(),
(vertices[2] - vertices[0]).norm()});
- else if (reference_cell == dealii::ReferenceCells::Quadrilateral)
+ else if (reference_cell == ReferenceCells::Quadrilateral)
// Return the longer one of the two diagonals of the quadrilateral
return std::max({(vertices[3] - vertices[0]).norm(),
(vertices[2] - vertices[1]).norm()});
- else if (reference_cell == dealii::ReferenceCells::Tetrahedron)
+ else if (reference_cell == ReferenceCells::Tetrahedron)
// Return the longest of the six edges of the tetrahedron
return std::max({(vertices[1] - vertices[0]).norm(),
(vertices[2] - vertices[0]).norm(),
(vertices[3] - vertices[0]).norm(),
(vertices[3] - vertices[1]).norm(),
(vertices[3] - vertices[2]).norm()});
- else if (reference_cell == dealii::ReferenceCells::Pyramid)
+ else if (reference_cell == ReferenceCells::Pyramid)
// Return ...
return std::max({// the longest diagonal of the quadrilateral base
// of the pyramid or ...
(vertices[4] - vertices[1]).norm(),
(vertices[4] - vertices[2]).norm(),
(vertices[4] - vertices[3]).norm()});
- else if (reference_cell == dealii::ReferenceCells::Wedge)
+ else if (reference_cell == ReferenceCells::Wedge)
// Return ...
return std::max({// the longest of the 2*3=6 diagonals of the three
// quadrilateral sides of the wedge or ...
(vertices[4] - vertices[3]).norm(),
(vertices[5] - vertices[4]).norm(),
(vertices[5] - vertices[3]).norm()});
- else if (reference_cell == dealii::ReferenceCells::Hexahedron)
+ else if (reference_cell == ReferenceCells::Hexahedron)
// Return the longest of the four diagonals of the hexahedron
return std::max({(vertices[7] - vertices[0]).norm(),
(vertices[6] - vertices[1]).norm(),
info.face_type =
is_mixed_mesh ?
(dcell->face(f)->reference_cell() !=
- dealii::ReferenceCells::get_hypercube<dim - 1>()) :
+ ReferenceCells::get_hypercube<dim - 1>()) :
0;
info.subface_index =
GeometryInfo<dim>::max_children_per_cell;
info.face_type = is_mixed_mesh ?
(cell->face(face_no)->reference_cell() !=
- dealii::ReferenceCells::get_hypercube<dim - 1>()) :
+ ReferenceCells::get_hypercube<dim - 1>()) :
0;
info.subface_index = GeometryInfo<dim>::max_children_per_cell;
face_data_by_cells[my_q].descriptor[hpq * scale].initialize(
quad_face);
reference_cell_types[my_q][hpq] =
- dealii::ReferenceCells::get_hypercube<dim>();
+ ReferenceCells::get_hypercube<dim>();
}
}
}
reorder_face_derivative_indices(
const unsigned int face_no,
const unsigned int index,
- const ReferenceCell reference_cell = dealii::ReferenceCells::Invalid)
+ const ReferenceCell reference_cell = ReferenceCells::Invalid)
{
Assert(index < dim, ExcInternalError());
- if ((reference_cell == dealii::ReferenceCells::Invalid ||
- reference_cell == dealii::ReferenceCells::get_hypercube<dim>()) ==
- false)
+ if ((reference_cell == ReferenceCells::Invalid ||
+ reference_cell == ReferenceCells::get_hypercube<dim>()) == false)
{
return index;
}
quadrature.push_back(*quadrature_hypercube);
else if (reference_cell.is_simplex())
quadrature.push_back(*quadrature_simplex);
- else if (reference_cell == dealii::ReferenceCells::Wedge)
+ else if (reference_cell == ReferenceCells::Wedge)
quadrature.push_back(*quadrature_wedge);
- else if (reference_cell ==
- dealii::ReferenceCells::Pyramid)
+ else if (reference_cell == ReferenceCells::Pyramid)
quadrature.push_back(*quadrature_pyramid);
else
Assert(false, ExcNotImplemented());
// need hypercube quadratures.
if ((dim < 3) ||
(reference_cell.is_hyper_cube() ||
- (reference_cell == dealii::ReferenceCells::Wedge) ||
- (reference_cell == dealii::ReferenceCells::Pyramid)))
+ (reference_cell == ReferenceCells::Wedge) ||
+ (reference_cell == ReferenceCells::Pyramid)))
quadrature.push_back(*quadrature_hypercube);
// In 3d, if the element is for simplex/wedge/pyramid
// cells, then we also need simplex quadratures
if ((dim == 3) &&
(reference_cell.is_simplex() ||
- (reference_cell == dealii::ReferenceCells::Wedge) ||
- (reference_cell == dealii::ReferenceCells::Pyramid)))
+ (reference_cell == ReferenceCells::Wedge) ||
+ (reference_cell == ReferenceCells::Pyramid)))
quadrature.push_back(*quadrature_simplex);
}
finite_elements.emplace_back(
std::make_shared<dealii::hp::FECollection<dim, spacedim>>(
FE_SimplexDGP<dim, spacedim>(1)));
- else if (reference_cell == dealii::ReferenceCells::Wedge)
+ else if (reference_cell == ReferenceCells::Wedge)
finite_elements.emplace_back(
std::make_shared<dealii::hp::FECollection<dim, spacedim>>(
FE_WedgeDGP<dim, spacedim>(1)));
- else if (reference_cell == dealii::ReferenceCells::Pyramid)
+ else if (reference_cell == ReferenceCells::Pyramid)
finite_elements.emplace_back(
std::make_shared<dealii::hp::FECollection<dim, spacedim>>(
FE_PyramidDGP<dim, spacedim>(1)));
this->reference_cells.clear();
for (const auto &i : reference_cells_ui)
this->reference_cells.emplace_back(
- dealii::internal::ReferenceCell::make_reference_cell_from_int(i));
+ dealii::internal::make_reference_cell_from_int(i));
}
quad_dof_identities
[most_dominating_fe_index][other_fe_index]
[cell->quad(q)->reference_cell() ==
- dealii::ReferenceCells::Quadrilateral],
+ ReferenceCells::Quadrilateral],
most_dominating_fe_index_face_no);
for (const auto &identity : identities)
quad_dof_identities
[most_dominating_fe_index][other_fe_index]
[cell->quad(q)->reference_cell() ==
- dealii::ReferenceCells::Quadrilateral],
+ ReferenceCells::Quadrilateral],
most_dominating_fe_index_face_no);
for (const auto &identity : identities)
Quadrature<dim - 1> quadrature(boundary_points, dummy_weights);
q_projector = QProjector<dim>::project_to_all_faces(
- dealii::ReferenceCells::get_hypercube<dim>(), quadrature);
+ ReferenceCells::get_hypercube<dim>(), quadrature);
}
for (const auto &cell : tria.active_cell_iterators())