FE_DGQ<dim> fe(1);
DoFHandler<dim> dof_handler(tria);
Vector<double> solution;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
DataOut<dim> data_out;
dof_handler.distribute_dofs (fe);
void
check ()
{
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
Triangulation<dim> tria;
make_mesh (tria);
FE_Q<dim> element(3);
DoFHandler<dim> dof(tria);
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
dof.distribute_dofs(element);
// test with two different functions: one
FE_Q<dim> element(3);
DoFHandler<dim> dof(tria);
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
dof.distribute_dofs(element);
// test with two different functions: one
DoFHandler<dim> dgq_dof_handler(tria);
Vector<double> q_solution;
Vector<double> dgq_solution;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
DataOut<dim> q_data_out, dgq_data_out;
ConstraintMatrix cm;
cm.close();
void output_results (const unsigned int cycle) const;
Triangulation<dim> triangulation;
- const MappingQ1<dim> mapping;
+ const MappingQGeneric<dim> mapping(1);
FE_DGQ<dim> fe;
DoFHandler<dim> dof_handler;
void test(std::ostream & /*out*/)
{
MyFunction<dim> func;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
parallel::distributed::Triangulation<dim> tr(MPI_COMM_WORLD);
GridGenerator::hyper_cube(tr);
void test(std::ostream & /*out*/)
{
MyFunction<dim> func;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
parallel::distributed::Triangulation<dim> tr(MPI_COMM_WORLD);
GridGenerator::hyper_cube(tr);
FE_Q<dim> fe(1);
deallog << "FE=" << fe.get_name() << std::endl;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
deallog << "Mapping=Q1" << std::endl;
FE_Q<dim> fe(2);
deallog << "FE=" << fe.get_name() << std::endl;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
deallog << "Mapping=Q1" << std::endl;
FE_DGPMonomial<dim> fe(1);
deallog << "FE=" << fe.get_name() << std::endl;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
deallog << "Mapping=Q1" << std::endl;
FE_DGPMonomial<dim> fe(2);
deallog << "FE=" << fe.get_name() << std::endl;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
deallog << "Mapping=Q1" << std::endl;
FE_DGPNonparametric<dim> fe(1);
deallog << "FE=" << fe.get_name() << std::endl;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
deallog << "Mapping=Q1" << std::endl;
FE_DGPNonparametric<dim> fe(2);
deallog << "FE=" << fe.get_name() << std::endl;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
deallog << "Mapping=Q1" << std::endl;
std::vector<Mapping<dim> *> mapping_ptr;
std::vector<std::string> mapping_strings;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
std::string mapping_name = "MappingQ1";
Triangulation<dim> tria;
const QGauss<dim> quadrature (fe.degree+1);\r
const unsigned int n_q_points = quadrature.size ();\r
MappingQ<dim> mapping(1);\r
- //MappingQ1<dim> mapping;\r
+ //MappingQGeneric<dim> mapping(1);\r
std::vector<double> div_v(n_q_points);\r
std::vector<typename FEValuesViews::Vector<dim>::curl_type> curl_v(n_q_points);\r
\r
const unsigned int n_q_points = quadrature.size ();\r
const unsigned int n_face_q_points = face_quadrature.size ();\r
//MappingQ<dim> mapping(2);\r
- MappingQ1<dim> mapping;\r
+ MappingQGeneric<dim> mapping(1);\r
std::vector<double> div_v(n_q_points);\r
std::vector<typename FEValuesViews::Vector<dim>::curl_type> curl_v(n_q_points);\r
std::vector<Tensor<3,dim> > hessians(n_q_points);\r
VectorFunction<dim> fe_function;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (mapping, fe, quadrature,
hp::DoFHandler<dim> dgq_dof_handler(tria);
Vector<double> q_solution;
Vector<double> dgq_solution;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
// refine a few cells
typename Triangulation<dim>::active_cell_iterator cell=tria.begin_active(),
}
hp::DoFHandler<dim> q_dof_handler(tria);
Vector<double> q_solution;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
// refine a few cells
typename Triangulation<dim>::active_cell_iterator cell=tria.begin_active(),
mg.set_zero();
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
MeshWorker::IntegrationInfoBox<dim> info_box;
UpdateFlags update_flags = update_values | update_gradients | update_hessians;
void
test_cochain(const Triangulation<dim> &tr, const FiniteElement<dim> &fe)
{
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
// Initialize DofHandler for a
// block system with local blocks
DoFHandler<dim> dof(tr);
local.cells = true;
local.faces = false;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
MeshWorker::IntegrationInfoBox<dim> info_box;
info_box.initialize_gauss_quadrature(1, 1, 1);
assemble(const DoFHandler<dim> &dof_handler, SparseMatrix<double> &matrix)
{
const FiniteElement<dim> &fe = dof_handler.get_fe();
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
MeshWorker::IntegrationInfoBox<dim> info_box;
const unsigned int n_gauss_points = dof_handler.get_fe().tensor_degree()+1;
MGLevelObject<SparseMatrix<double> > dg_down)
{
const FiniteElement<dim> &fe = dof_handler.get_fe();
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
MeshWorker::IntegrationInfoBox<dim> info_box;
const unsigned int n_gauss_points = dof_handler.get_fe().tensor_degree()+1;
assemble(const DoFHandler<dim> &dof_handler, SparseMatrix<double> &matrix)
{
const FiniteElement<dim> &fe = dof_handler.get_fe();
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
MeshWorker::IntegrationInfoBox<dim> info_box;
const unsigned int n_gauss_points = dof_handler.get_fe().tensor_degree()+1;
MGLevelObject<SparseMatrix<double> > dg_down)
{
const FiniteElement<dim> &fe = dof_handler.get_fe();
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
MeshWorker::IntegrationInfoBox<dim> info_box;
const unsigned int n_gauss_points = dof_handler.get_fe().tensor_degree()+1;
private:
- const MappingQ1<dim> mapping;
+ const MappingQGeneric<dim> mapping(1);
void setup_system ();
local.cells = true;
local.faces = faces;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
MeshWorker::IntegrationInfoBox<dim> info_box;
info_box.initialize_gauss_quadrature(1, 1, 1);
TrilinosWrappers::MPI::Vector vec_rel ( locally_relevant_set);
vec_rel = vec;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
Vector<float> indicators(tr.n_active_cells());
DerivativeApproximation::approximate_gradient (mapping,
dofh,
local.cells = true;
local.faces = faces;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
MeshWorker::IntegrationInfoBox<dim> info_box;
info_box.initialize_gauss_quadrature(1, 1, 1);
void refine_local ();
Triangulation<dim> triangulation;
- const MappingQ1<dim> mapping;
+ const MappingQGeneric<dim> mapping(1);
FESystem<dim> fe;
DoFHandler<dim> mg_dof_handler_renumbered;
typename FunctionMap<dim>::type dirichlet_boundary;
ZeroFunction<dim> homogeneous_dirichlet_bc (1);
dirichlet_boundary[0] = &homogeneous_dirichlet_bc;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
VectorTools::interpolate_boundary_values (mapping,
mg_dof_handler,
dirichlet_boundary,
{
mg_matrices = 0.;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
MeshWorker::IntegrationInfoBox<dim> info_box;
UpdateFlags update_flags = update_values | update_gradients | update_hessians;
info_box.add_update_flags_all(update_flags);
typename FunctionMap<dim>::type dirichlet_boundary;
ZeroFunction<dim> homogeneous_dirichlet_bc (1);
dirichlet_boundary[0] = &homogeneous_dirichlet_bc;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
VectorTools::interpolate_boundary_values (mapping,
mg_dof_handler,
dirichlet_boundary,
{
mg_matrices = 0.;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
MeshWorker::IntegrationInfoBox<dim> info_box;
UpdateFlags update_flags = update_values | update_gradients | update_hessians;
info_box.add_update_flags_all(update_flags);
typename FunctionMap<dim>::type dirichlet_boundary;
ZeroFunction<dim> homogeneous_dirichlet_bc (1);
dirichlet_boundary[0] = &homogeneous_dirichlet_bc;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
VectorTools::interpolate_boundary_values (mapping,
mg_dof_handler,
dirichlet_boundary,
FE_DGQ<dim> fe(2);
DoFHandler<dim> dof_handler(tria);
Vector<double> solution;
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
QMidpoint<dim> q_midpoint;
FEValues<dim> fe_values(mapping, fe, q_midpoint, update_q_points);
dof.distribute_dofs(fe);
QGauss<dim-1> quadrature(degree+2);
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
TestFunction<dim> f(degree-1);
std::map<types::global_dof_index, double> boundary_constraints;
solution.reinit (dof_handler.n_dofs());
system_rhs.reinit (dof_handler.n_dofs());
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
MeshWorker::IntegrationInfoBox<dim> info_box;
UpdateFlags update_flags = update_values | update_gradients;
info_box.add_update_flags_all(update_flags);
solution.reinit (dof_handler.n_dofs());
system_rhs.reinit (dof_handler.n_dofs());
- MappingQ1<dim> mapping;
+ MappingQGeneric<dim> mapping(1);
MeshWorker::IntegrationInfoBox<dim> info_box;
UpdateFlags update_flags = update_values | update_gradients;
info_box.add_update_flags_all(update_flags);