matrix_for_bc);
}
- rhs.compress (Add);
+ rhs.compress (VectorOperation::add);
// Now that we have the right linear system, we solve it using the CG
// method with a simple Jacobi preconditioner:
StokesSystem<dim> (stokes_fe));
stokes_matrix.compress();
- stokes_rhs.compress(Add);
+ stokes_rhs.compress(VectorOperation::add);
rebuild_stokes_matrix = false;
Assembly::CopyData::
TemperatureRHS<dim> (temperature_fe));
- temperature_rhs.compress(Add);
+ temperature_rhs.compress(VectorOperation::add);
}
cell->get_dof_indices(local_dof_indices);
const std::vector<Point<dim> > &q_points = fe_v.get_quadrature_points();
- const std::vector<Point<dim> > &normals = fe_v.get_cell_normal_vectors();
+ const std::vector<Point<dim> > &normals = fe_v.get_normal_vectors();
wind.vector_value_list(q_points, cell_wind);
// We then form the integral over the current cell for all degrees of
std::vector<Vector<double> > singular_cell_wind( singular_quadrature.size(),
Vector<double>(dim) );
- const std::vector<Point<dim> > &singular_normals = fe_v_singular.get_cell_normal_vectors();
+ const std::vector<Point<dim> > &singular_normals = fe_v_singular.get_normal_vectors();
const std::vector<Point<dim> > &singular_q_points = fe_v_singular.get_quadrature_points();
wind.vector_value_list(singular_q_points, singular_cell_wind);
fe_v.reinit(cell);
const std::vector<Point<dim> > &q_points = fe_v.get_quadrature_points();
- const std::vector<Point<dim> > &normals = fe_v.get_cell_normal_vectors();
+ const std::vector<Point<dim> > &normals = fe_v.get_normal_vectors();
cell->get_dof_indices(dofs);
fe_v.get_function_values(phi, local_phi);
#include <deal.II/lac/solver_gmres.h>
#include <deal.II/lac/sparse_ilu.h>
#include <deal.II/lac/sparse_direct.h>
+#include <deal.II/lac/constraint_matrix.h>
#include <deal.II/grid/tria.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/dofs/dof_accessor.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/dofs/dof_renumbering.h>
-#include <deal.II/dofs/dof_constraints.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_values.h>
#include <deal.II/base/timer.h>
#include <deal.II/base/work_stream.h>
-#include <deal.II/dofs/dof_constraints.h>
#include <deal.II/dofs/dof_renumbering.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/lac/precondition_selector.h>
#include <deal.II/lac/solver_cg.h>
#include <deal.II/lac/sparse_direct.h>
+#include <deal.II/lac/constraint_matrix.h>
#include <deal.II/numerics/data_out.h>
#include <deal.II/numerics/vector_tools.h>