deallog.depth_console(0);
deallog.threshold_double(1.e-10);
- check ("gauss", 2, QGauss2<1>());
- check ("gauss", 2, QGauss2<2>());
- check ("gauss", 2, QGauss2<3>());
+ check ("gauss", 2, QGauss<1>(2));
+ check ("gauss", 2, QGauss<2>(2));
+ check ("gauss", 2, QGauss<3>(2));
- check ("gauss", 2, QGauss2<3>());
- check ("gauss", 6, QGauss6<3>());
+ check ("gauss", 2, QGauss<3>(2));
+ check ("gauss", 6, QGauss<3>(6));
check ("gauss", 10, QGauss<3>(10));
check ("weddle", 0, QWeddle<2>());
{
Vector<float> estimated_error_per_cell (this->triangulation->n_active_cells());
KellyErrorEstimator<dim>::estimate (this->dof_handler,
- QGauss3<dim-1>(),
+ QGauss<dim-1>(3),
typename FunctionMap<dim>::type(),
this->solution,
estimated_error_per_cell);
{
MappingQ<dim> mapping(4);
KellyErrorEstimator<dim>::estimate (mapping, primal_solver.dof_handler,
- QGauss3<dim-1>(),
+ QGauss<dim-1>(3),
typename FunctionMap<dim>::type(),
primal_solver.solution,
error_indicators);
u);
KellyErrorEstimator<3>::estimate (dof_handler,
- QGauss2<2>(),
+ QGauss<2>(2),
FunctionMap<3>::type(),
u,
e);
const Function<dim> &function,
Vector<double> &vec,
const bool enforce_zero_boundary = false,
- const Quadrature<dim-1> & = QGauss2<dim-1>(),
+ const Quadrature<dim-1> & = QGauss<dim-1>(2),
const bool project_to_boundary_first = false)
{
Assert (dof.get_fe().n_components() == function.n_components,
function_map[0] = &coefficient;
hp::QCollection<dim-1> face_quadrature;
- face_quadrature.push_back (QGauss6<dim-1>());
+ face_quadrature.push_back (QGauss<dim-1>(6));
std::vector<unsigned int> dof_to_boundary_mapping;
DoFTools::map_dof_to_boundary_indices (dof,
mapping.push_back (MappingQ<dim>(3));
hp::QCollection<dim> quadrature;
- quadrature.push_back (QGauss6<dim>());
+ quadrature.push_back (QGauss<dim>(6));
// create sparsity pattern. note
// that different components should