FE_Q<dim> fe(fe_degree);
DoFHandler<dim> dof(tria);
dof.distribute_dofs(fe);
- ConstraintMatrix constraints;
+ AffineConstraints<double> constraints;
constraints.close();
do_test<dim, fe_degree, double, fe_degree + 1>(dof, constraints);
FE_Q<dim> fe(fe_degree);
DoFHandler<dim> dof(tria);
dof.distribute_dofs(fe);
- ConstraintMatrix constraints;
+ AffineConstraints<double> constraints;
VectorTools::interpolate_boundary_values(dof,
0,
Functions::ZeroFunction<dim>(),
FE_Q<dim> fe(fe_degree);
DoFHandler<dim> dof(tria);
dof.distribute_dofs(fe);
- ConstraintMatrix constraints;
+ AffineConstraints<double> constraints;
constraints.close();
do_test<dim, fe_degree, double, fe_degree + 1>(dof, constraints);
DoFHandler<dim> dof_handler;
DoFHandler<dim> velocity_dof_handler;
- ConstraintMatrix constraints;
+ AffineConstraints<double> constraints;
BlockSparsityPattern sparsity_pattern;
BlockSparseMatrix<double> system_matrix;
std::vector<SymmetricTensor<2, dim>> symgrad_phi_u(dofs_per_cell);
- std::vector<ConstraintMatrix> boundary_constraints(
+ std::vector<AffineConstraints<double>> boundary_constraints(
triangulation.n_levels());
- std::vector<ConstraintMatrix> boundary_interface_constraints(
+ std::vector<AffineConstraints<double>> boundary_interface_constraints(
triangulation.n_levels());
for (unsigned int level = 0; level < triangulation.n_levels(); ++level)
{
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_out.h>
-#include <deal.II/lac/constraint_matrix.h>
#include <deal.II/lac/precondition.h>
#include <deal.II/lac/solver_cg.h>
#include <deal.II/lac/sparse_matrix.h>
TestPoly<2> pol(deg);
// Project solution onto BR FE field
- ConstraintMatrix hn_constraints;
+ AffineConstraints<double> hn_constraints;
hn_constraints.clear();
DoFTools::make_hanging_node_constraints(*dof_handler, hn_constraints);
hn_constraints.close();
solution_q.reinit(dof_handler->n_dofs());
deformation.reinit(dof_handler_def->n_dofs());
- ConstraintMatrix hn_constraints_def;
+ AffineConstraints<double> hn_constraints_def;
hn_constraints_def.clear();
DoFTools::make_hanging_node_constraints(*dof_handler_def, hn_constraints_def);
hn_constraints_def.close();
#include <deal.II/fe/mapping_q1.h>
#include <deal.II/grid/grid_generator.h>
-#include <deal.II/grid/tria_boundary_lib.h>
#include <deal.II/lac/vector.h>
std::cout << "n_cells: " << dof_handler.get_triangulation().n_active_cells()
<< std::endl;
- ConstraintMatrix constraints;
+ AffineConstraints<double> constraints;
constraints.clear();
dealii::DoFTools::make_hanging_node_constraints(dof_handler, constraints);
constraints.close();
std::cout << "...start plotting shape function" << std::endl;
std::cout << "Patches for output: " << patches << std::endl;
- ConstraintMatrix constraints;
+ AffineConstraints<double> constraints;
constraints.clear();
dealii::DoFTools::make_hanging_node_constraints(dof_handler, constraints);
constraints.close();
Vector<double> localized_solution;
PETScWrappers::MPI::Vector system_rhs;
- ConstraintMatrix constraints;
+ AffineConstraints<double> constraints;
MPI_Comm mpi_communicator;
const unsigned int n_mpi_processes;
deallog << "...start plotting shape function" << std::endl;
deallog << "Patches for output: " << patches << std::endl;
- ConstraintMatrix constraints;
+ AffineConstraints<double> constraints;
constraints.clear();
dealii::DoFTools::make_hanging_node_constraints(dof_handler, constraints);
constraints.close();
dof_handler.distribute_dofs(fe_collection);
- ConstraintMatrix constraints;
+ AffineConstraints<double> constraints;
DoFTools::make_hanging_node_constraints(dof_handler, constraints);
constraints.close();
* Main class
*/
template <int dim>
-class MeshWorkerConstraintMatrixTest
+class MeshWorkerAffineConstraintsTest
{
public:
- MeshWorkerConstraintMatrixTest(const FiniteElement<dim> &fe);
- ~MeshWorkerConstraintMatrixTest();
+ MeshWorkerAffineConstraintsTest(const FiniteElement<dim> &fe);
+ ~MeshWorkerAffineConstraintsTest();
void
run();
template <int dim>
-MeshWorkerConstraintMatrixTest<dim>::MeshWorkerConstraintMatrixTest(
+MeshWorkerAffineConstraintsTest<dim>::MeshWorkerAffineConstraintsTest(
const FiniteElement<dim> &fe)
: mapping()
, dof_handler(triangulation)
template <int dim>
-MeshWorkerConstraintMatrixTest<dim>::~MeshWorkerConstraintMatrixTest()
+MeshWorkerAffineConstraintsTest<dim>::~MeshWorkerAffineConstraintsTest()
{
dof_handler.clear();
}
template <int dim>
void
-MeshWorkerConstraintMatrixTest<dim>::setup_system()
+MeshWorkerAffineConstraintsTest<dim>::setup_system()
{
dof_handler.distribute_dofs(fe);
*/
template <int dim>
void
-MeshWorkerConstraintMatrixTest<dim>::assemble_system_MeshWorker()
+MeshWorkerAffineConstraintsTest<dim>::assemble_system_MeshWorker()
{
MeshWorker::IntegrationInfoBox<dim> info_box;
*/
template <int dim>
void
-MeshWorkerConstraintMatrixTest<dim>::assemble_MeshWorker()
+MeshWorkerAffineConstraintsTest<dim>::assemble_MeshWorker()
{
MeshWorker::IntegrationInfoBox<dim> info_box;
template <int dim>
void
-MeshWorkerConstraintMatrixTest<dim>::createInhomConstraints()
+MeshWorkerAffineConstraintsTest<dim>::createInhomConstraints()
{
this->constraintsInhom.clear();
// boundary constraints
template <int dim>
void
-MeshWorkerConstraintMatrixTest<dim>::run()
+MeshWorkerAffineConstraintsTest<dim>::run()
{
setup_system();
createInhomConstraints();
FE_Q<2> fe(1);
deallog.push(fe.get_name());
- MeshWorkerConstraintMatrixTest<2> test(fe);
+ MeshWorkerAffineConstraintsTest<2> test(fe);
test.run();
deallog.pop();
FE_DGQ<2> fe2(1);
deallog.push(fe2.get_name());
- MeshWorkerConstraintMatrixTest<2> test2(fe2);
+ MeshWorkerAffineConstraintsTest<2> test2(fe2);
test2.run();
deallog.pop();
}
DoFHandler<dim> dof(tria);
dof.distribute_dofs(fe);
dof.distribute_mg_dofs();
- ConstraintMatrix constraints;
+ AffineConstraints<double> constraints;
constraints.close();
MGConstrainedDoFs mg_constrained_dofs;
(++dof_handler.begin_active())->set_active_fe_index(1);
dof_handler.distribute_dofs(fe);
- ConstraintMatrix constraints;
+ AffineConstraints<double> constraints;
DoFTools::make_hanging_node_constraints(dof_handler, constraints);
constraints.close();
#include <deal.II/grid/grid_in.h>
#include <deal.II/grid/grid_tools.h>
#include <deal.II/grid/tria.h>
-#include <deal.II/grid/tria_boundary_lib.h>
#include <deal.II/lac/affine_constraints.h>
#include <deal.II/lac/block_sparse_matrix.h>