hp::FEFaceValues<dim,spacedim> x_fe_face_values (mapping_collection,
fe_collection,
face_quadrature_collection,
- update_q_points |
+ update_quadrature_points |
update_normal_vectors);
// have a map that stores normal vectors for each vector-dof tuple we want
hp::FEFaceValues<dim,spacedim> x_fe_face_values (mapping_collection,
fe_collection,
face_quadrature_collection,
- update_q_points |
+ update_quadrature_points |
update_normal_vectors);
// Extract a list that collects all vector components that belong to the
support_quadrature(degree),
fe_values(mapping_q_eulerian.euler_dof_handler->get_fe(),
support_quadrature,
- update_values | update_q_points)
+ update_values | update_quadrature_points)
{}
template <int dim, class VectorType, int spacedim>
QGauss<2> quadrature(3);
MappingQ<3> mapping(2);
FEFaceValues<3> fe_face_values1 (mapping, fe, quadrature,
- update_q_points | update_JxW_values |
+ update_quadrature_points | update_JxW_values |
update_normal_vectors);
FESubfaceValues<3> fe_face_values2 (mapping, fe, quadrature,
- update_q_points | update_JxW_values |
+ update_quadrature_points | update_JxW_values |
update_normal_vectors);
void check_this (Triangulation<3> &tria)
{
Scratch(const FiniteElement<dim> &fe, const Quadrature<dim> &quadrature) :
fe_collection(fe), quadrature_collection(quadrature), x_fe_values(
- fe_collection, quadrature_collection, update_q_points), rhs_values(
+ fe_collection, quadrature_collection, update_quadrature_points), rhs_values(
quadrature_collection.size())
{
}
Scratch(const Scratch &data) :
fe_collection(data.fe_collection), quadrature_collection(
data.quadrature_collection), x_fe_values(fe_collection,
- quadrature_collection, update_q_points), rhs_values(
+ quadrature_collection, update_quadrature_points), rhs_values(
data.rhs_values)
{
}
{
Scratch(const FiniteElement<dim> &fe, const Quadrature<dim> &quadrature) :
fe_collection(fe), quadrature_collection(quadrature), x_fe_values(
- fe_collection, quadrature_collection, update_q_points), rhs_values(
+ fe_collection, quadrature_collection, update_quadrature_points), rhs_values(
quadrature_collection.size())
{
}
Scratch(const Scratch &data) :
fe_collection(data.fe_collection), quadrature_collection(
data.quadrature_collection), x_fe_values(fe_collection,
- quadrature_collection, update_q_points), rhs_values(
+ quadrature_collection, update_quadrature_points), rhs_values(
data.rhs_values)
{
}
FEValues<2> fe_values (dof_handler.get_fe(),
qr, UpdateFlags(update_values |
update_gradients |
- update_q_points |
+ update_quadrature_points |
update_JxW_values));
MatrixTools::create_laplace_matrix (dof_handler, qr, B.block(0,0));
FEValues<dim> fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
Quadrature<2> quadrature_formula(fe.get_unit_support_points());
FEValues<2> fe_values (fe, quadrature_formula,
- update_values | update_gradients | update_JxW_values|update_q_points);
+ update_values | update_gradients | update_JxW_values|update_quadrature_points);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
std::vector<types::global_dof_index> local_dof_indices (dofs_per_cell);
QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (fe, q_face,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
FESubfaceValues<dim> fe_subface_values (fe, q_face,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
Point<dim> n1, n2;
for (typename DoFHandler<dim>::active_cell_iterator
QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
FESubfaceValues<dim> fe_subface_values (mapping, fe, q_face,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
Point<dim> n1, n2;
for (typename DoFHandler<dim>::active_cell_iterator
QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
FESubfaceValues<dim> fe_subface_values (mapping, fe, q_face,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
Point<dim> n1, n2;
for (typename DoFHandler<dim>::active_cell_iterator
QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
FESubfaceValues<dim> fe_subface_values (mapping, fe, q_face,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
Point<dim> n1, n2;
for (typename DoFHandler<dim>::active_cell_iterator
QGauss<2> quadrature(3);
FE_Q<3> fe(1);
FEFaceValues<3> fe_face_values1 (fe, quadrature,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
FEFaceValues<3> fe_face_values2 (fe, quadrature,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
DoFHandler<3> dof_handler (tria);
dof_handler.distribute_dofs (fe);
update_values |
update_gradients |
update_second_derivatives |
- update_q_points |
+ update_quadrature_points |
update_jacobians);
fe_values.reinit (dof_handler.begin_active(), 0);
QGauss<1> quad (6);
FEFaceValues<2> fe_v_face (dof_handler.get_fe (), quad,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_gradients |
update_normal_vectors |
update_JxW_values));
FEFaceValues<2> fe_v_face_n (dof_handler.get_fe (), quad,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_gradients |
update_normal_vectors |
update_JxW_values));
FEValues<dim> x_fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
FEFaceValues<dim> x_fe_face_values (fe, face_quadrature_formula,
- update_values | update_q_points |
+ update_values | update_quadrature_points |
update_normal_vectors | update_JxW_values);
const RightHandSide<dim> right_hand_side;
const UpdateFlags update_flags = update_values
| update_gradients
- | update_q_points
+ | update_quadrature_points
| update_JxW_values;
const UpdateFlags face_update_flags = update_values
- | update_q_points
+ | update_quadrature_points
| update_JxW_values
| update_normal_vectors;
const UpdateFlags update_flags = update_values
| update_gradients
- | update_q_points
+ | update_quadrature_points
| update_JxW_values;
const UpdateFlags face_update_flags = update_values
- | update_q_points
+ | update_quadrature_points
| update_JxW_values
| update_normal_vectors;
assemble_rhs (Vector<double> &rhs) const
{
FEValues<dim> fe_values (*this->fe, *this->quadrature,
- update_values | update_q_points |
+ update_values | update_quadrature_points |
update_JxW_values);
const unsigned int dofs_per_cell = this->fe->dofs_per_cell;
QTrapez<dim> vertex_quadrature;
FEValues<dim> fe_values (dof_handler.get_fe(),
vertex_quadrature,
- update_gradients | update_q_points);
+ update_gradients | update_quadrature_points);
std::vector<Tensor<1,dim> >
solution_gradients (vertex_quadrature.size());
assemble_rhs (Vector<double> &rhs) const
{
FEValues<dim> fe_values (*this->fe, *this->quadrature,
- update_values | update_q_points |
+ update_values | update_quadrature_points |
update_JxW_values);
const unsigned int dofs_per_cell = this->fe->dofs_per_cell;
QGauss<dim> quadrature(4);
FEValues<dim> fe_values (dof_handler.get_fe(), quadrature,
update_gradients |
- update_q_points |
+ update_quadrature_points |
update_JxW_values);
const unsigned int n_q_points = fe_values.n_quadrature_points;
const unsigned int dofs_per_cell = dof_handler.get_fe().dofs_per_cell;
fe_values (fe, quadrature,
update_values |
update_second_derivatives |
- update_q_points |
+ update_quadrature_points |
update_JxW_values),
right_hand_side (&right_hand_side)
{
QGauss<dim> quadrature_formula(4);
FEValues<dim> fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
FEValues<dim> fe_values (fe, quadrature,
update_values | update_gradients |
update_second_derivatives |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
FEValues<dim> neighbor_fe_values (fe, quadrature,
update_gradients);
QGauss<dim> quadrature_formula(4);
FEValues<dim> fe_values (dof_handler.get_fe(), quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int n_q_points = quadrature_formula.size();
FEValues<dim> fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
FEValues<dim> fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
FEValues<dim> x_fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
FEValues<dim> x_fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
FEValues<dim> x_fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
FEValues<dim> x_fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
FEValues<dim> x_fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
FEValues<dim> x_fe_values (*fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
FEFaceValues<dim> x_fe_face_values (*fe, face_quadrature_formula,
- update_values | update_q_points |
+ update_values | update_quadrature_points |
update_normal_vectors | update_JxW_values);
const RightHandSide<dim> right_hand_side;
FEValues<dim> x_fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
FEFaceValues<dim> fe_face_values (fe, q_face,
update_normal_vectors |
- update_q_points |
+ update_quadrature_points |
update_JxW_values);
FESubfaceValues<dim> fe_subface_values (fe, q_face,
update_normal_vectors |
- update_q_points |
+ update_quadrature_points |
update_JxW_values);
double v1=0, v2=0;
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_normal_vectors |
- update_q_points |
+ update_quadrature_points |
update_JxW_values);
FESubfaceValues<dim> fe_subface_values (mapping, fe, q_face,
update_normal_vectors |
- update_q_points |
+ update_quadrature_points |
update_JxW_values);
double v1=0, v2=0;
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_normal_vectors |
- update_q_points |
+ update_quadrature_points |
update_JxW_values);
FESubfaceValues<dim> fe_subface_values (mapping, fe, q_face,
update_normal_vectors |
- update_q_points |
+ update_quadrature_points |
update_JxW_values);
double v1=0, v2=0;
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_normal_vectors |
- update_q_points |
+ update_quadrature_points |
update_JxW_values);
FESubfaceValues<dim> fe_subface_values (mapping, fe, q_face,
update_normal_vectors |
- update_q_points |
+ update_quadrature_points |
update_JxW_values);
double v1=0, v2=0;
if (fe.has_support_points())
{
Quadrature<dim> quad(fe.get_unit_support_points());
- fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_q_points));
+ fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_quadrature_points));
}
for (typename DoFHandler<dim>::active_cell_iterator cell=dof.begin_active();
if (fe.has_support_points())
{
Quadrature<dim> quad(fe.get_unit_support_points());
- fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_q_points));
+ fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_quadrature_points));
}
for (typename DoFHandler<dim>::cell_iterator cell=dof.begin(level);
if (fe.has_support_points())
{
Quadrature<dim> quad(fe.get_unit_support_points());
- fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_q_points));
+ fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_quadrature_points));
}
for (typename DoFHandler<dim>::active_cell_iterator cell=dof.begin_active();
if (fe.has_support_points())
{
Quadrature<dim> quad(fe.get_unit_support_points());
- fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_q_points));
+ fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_quadrature_points));
}
for (typename DoFHandler<dim>::active_cell_iterator cell=dof.begin_active();
if (fe.has_support_points())
{
Quadrature<dim> quad(fe.get_unit_support_points());
- fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_q_points));
+ fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_quadrature_points));
}
for (typename DoFHandler<dim>::active_cell_iterator cell=dof.begin_active();
if (fe.has_support_points())
{
Quadrature<dim> quad(fe.get_unit_support_points());
- fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_q_points));
+ fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_quadrature_points));
}
for (typename DoFHandler<dim>::active_cell_iterator cell=dof.begin_active();
if (fe.has_support_points())
{
Quadrature<dim> quad(fe.get_unit_support_points());
- fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_q_points));
+ fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_quadrature_points));
}
for (typename DoFHandler<dim>::active_cell_iterator cell=dof.begin_active();
if (fe.has_support_points())
{
Quadrature<dim> quad(fe.get_unit_support_points());
- fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_q_points));
+ fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_quadrature_points));
}
for (typename DoFHandler<dim>::active_cell_iterator cell=dof.begin_active();
if (fe.has_support_points())
{
Quadrature<dim> quad(fe.get_unit_support_points());
- fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_q_points));
+ fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_quadrature_points));
}
for (typename DoFHandler<dim>::active_cell_iterator cell=dof.begin_active();
if (fe.has_support_points())
{
Quadrature<dim> quad(fe.get_unit_support_points());
- fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_q_points));
+ fevalues = std_cxx11::shared_ptr<FEValues<dim> >(new FEValues<dim>(fe, quad, update_quadrature_points));
}
for (typename DoFHandler<dim>::active_cell_iterator cell=dof.begin_active();
QGauss<2> quad (6);
FEValues<2> fe_values (mapping, dof_handler->get_fe (), quad,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_gradients |
update_JxW_values|
update_contravariant_transformation));
hp::FEValues<dim> hp_fe_values (fe_collection,
quadrature_collection,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const RightHandSide<dim> rhs_function;
QTrapez<dim> vertex_quadrature;
hp::FEValues<dim> fe_values (dof_handler.get_fe(),
vertex_quadrature,
- update_gradients | update_q_points);
+ update_gradients | update_quadrature_points);
std::vector<Tensor<1,dim> >
solution_gradients (vertex_quadrature.size());
assemble_rhs (Vector<double> &rhs) const
{
hp::FEValues<dim> fe_values (*this->fe, *this->quadrature,
- update_values | update_q_points |
+ update_values | update_quadrature_points |
update_JxW_values);
const unsigned int dofs_per_cell = (*this->fe)[0].dofs_per_cell;
QGauss<dim> quadrature(4);
hp::FEValues<dim> fe_values (dof_handler.get_fe(), quadrature,
update_gradients |
- update_q_points |
+ update_quadrature_points |
update_JxW_values);
const unsigned int n_q_points = fe_values.get_present_fe_values().n_quadrature_points;
const unsigned int dofs_per_cell = dof_handler.get_fe().dofs_per_cell;
fe_values (fe, quadrature,
update_values |
update_second_derivatives |
- update_q_points |
+ update_quadrature_points |
update_JxW_values),
right_hand_side (&right_hand_side)
{
hp::QCollection<dim> quadrature_formula(QGauss<dim>(4));
hp::FEValues<dim> fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe[0].dofs_per_cell;
const unsigned int n_q_points = quadrature_formula[0].size();
hp::FEValues<dim> fe_values (fe, quadrature,
update_values | update_gradients |
update_second_derivatives |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
hp::FEValues<dim> neighbor_fe_values (fe, quadrature,
update_gradients);
hp::QCollection<dim> quadrature_formula(QGauss<dim>(4));
hp::FEValues<dim> fe_values (dof_handler.get_fe(), quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int n_q_points = quadrature_formula[0].size();
QGauss<2> quad (6);
FEValues<2> fe_values (mapping, dof_handler->get_fe (), quad,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_JxW_values));
const unsigned int n_q_points = quad.size();
QGauss<2> quad (6);
FEValues<2> fe_values (mapping, dof_handler.get_fe (), quad,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_gradients |
update_JxW_values|
update_contravariant_transformation));
QGauss<2> quad (3);
FEValues<2> fe_values (dof_handler->get_fe (), quad,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_gradients |
update_JxW_values));
Vector<double> &rhs_vector,
const Function<dim> *const coefficient = 0)
{
- UpdateFlags update_flags = UpdateFlags(update_values | update_JxW_values | update_q_points);
+ UpdateFlags update_flags = UpdateFlags(update_values | update_JxW_values | update_quadrature_points);
if (coefficient != 0)
- update_flags = UpdateFlags (update_flags | update_q_points);
+ update_flags = UpdateFlags (update_flags | update_quadrature_points);
FEValues<dim> fe_values (mapping, dof.get_fe(), q, update_flags);
rhs_vector = 0;
UpdateFlags update_flags = UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_JxW_values);
FEValues<dim> fe_values (mapping, fe, quadrature, update_flags);
QTrapez<1> q_trapez;
const unsigned int div=2;
QIterated<dim> q(q_trapez, div);
- FEValues<dim> fe(element, q, update_values|update_gradients|update_q_points);
+ FEValues<dim> fe(element, q, update_values|update_gradients|update_quadrature_points);
fe.reinit(c);
for (unsigned int q_point=0; q_point< q.size(); ++q_point)
QGauss<2> quad (3);
FEValues<2> fe_values (dof_handler->get_fe (), quad,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_gradients |
update_JxW_values));
Vector<double> &rhs_vector,
const Function<dim> *const coefficient = 0)
{
- UpdateFlags update_flags = UpdateFlags(update_values | update_JxW_values | update_q_points);
+ UpdateFlags update_flags = UpdateFlags(update_values | update_JxW_values | update_quadrature_points);
if (coefficient != 0)
- update_flags = UpdateFlags (update_flags | update_q_points);
+ update_flags = UpdateFlags (update_flags | update_quadrature_points);
FEValues<dim> fe_values (mapping, dof.get_fe(), q, update_flags);
rhs_vector = 0;
UpdateFlags update_flags = UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_JxW_values);
FEValues<dim> fe_values (mapping, fe, quadrature, update_flags);
Assert (q.size() == 1, ExcInternalError());
FEFaceValues<dim> fe_values(mapping, fe, q,
- UpdateFlags(update_q_points
+ UpdateFlags(update_quadrature_points
| update_JxW_values
| update_values
| update_gradients
Assert (q.size() == 1, ExcInternalError());
FESubfaceValues<dim> fe_values(mapping, fe, q,
- UpdateFlags(update_q_points
+ UpdateFlags(update_quadrature_points
| update_JxW_values
| update_values
| update_gradients
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
- update_values | update_gradients | update_q_points);
+ update_values | update_gradients | update_quadrature_points);
fe_values.reinit (dof.begin_active());
// let the FEValues object compute the
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
- update_values | update_gradients | update_q_points);
+ update_values | update_gradients | update_quadrature_points);
fe_values.reinit (dof.begin_active());
// let the FEValues object compute the
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
- update_values | update_gradients | update_q_points);
+ update_values | update_gradients | update_quadrature_points);
fe_values.reinit (dof.begin_active());
// let the FEValues object compute the
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
- update_values | update_gradients | update_q_points);
+ update_values | update_gradients | update_quadrature_points);
fe_values.reinit (dof.begin_active());
// let the FEValues object compute the
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
- update_values | update_gradients | update_q_points);
+ update_values | update_gradients | update_quadrature_points);
fe_values.reinit (dof.begin_active());
// let the FEValues object compute the
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
- update_values | update_gradients | update_q_points);
+ update_values | update_gradients | update_quadrature_points);
fe_values.reinit (dof.begin_active());
// let the FEValues object compute the
QTrapez<1> q_trapez;
const unsigned int div=2;
QIterated<dim> q(q_trapez, div);
- FEValues<dim> fe(element, q, update_values|update_gradients|update_q_points);
+ FEValues<dim> fe(element, q, update_values|update_gradients|update_quadrature_points);
fe.reinit(c);
for (unsigned int q_point=0; q_point< q.size(); ++q_point)
MappingCartesian<dim> mapping;
FEFaceValues<dim> feface(mapping, fe_ned, quadrature,
- update_values|update_q_points|update_normal_vectors);
+ update_values|update_quadrature_points|update_normal_vectors);
FEFaceValues<dim> feneighbor(mapping, fe_ned, quadrature,
- update_values|update_q_points|update_normal_vectors);
+ update_values|update_quadrature_points|update_normal_vectors);
FESubfaceValues<dim> fesubface(mapping, fe_ned, quadrature,
- update_values|update_q_points|update_normal_vectors);
+ update_values|update_quadrature_points|update_normal_vectors);
for (typename DoFHandler<dim>::active_cell_iterator
std::vector<Vector<double> > shape_values (quadrature.size(),
Vector<double>(dim));
FEValues<dim> fe(fe_ned, quadrature,
- update_values|update_q_points);
+ update_values|update_quadrature_points);
for (typename DoFHandler<dim>::active_cell_iterator
c = dof.begin_active();
QTrapez<1> q_trapez;
const unsigned int div=2;
QIterated<dim> q(q_trapez, div);
- FEValues<dim> fe(element, q, update_values|update_gradients|update_q_points);
+ FEValues<dim> fe(element, q, update_values|update_gradients|update_quadrature_points);
fe.reinit(c);
for (unsigned int q_point=0; q_point< q.size(); ++q_point)
QTrapez<1> q_trapez;
const unsigned int div=10;
QIterated<dim> q(q_trapez, div);
- FEValues<dim> fe(fe_rt, q, update_values|update_gradients|update_q_points);
+ FEValues<dim> fe(fe_rt, q, update_values|update_gradients|update_quadrature_points);
fe.reinit(c);
Assert (fe.get_fe().n_components() == dim, ExcInternalError());
QGauss<2> quad (6);
FEValues<2> fe_values (mapping, dof_handler->get_fe (), quad,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_gradients |
update_JxW_values|
update_contravariant_transformation));
FEFaceValues<2> fe_v_face (dof_handler->get_fe (), quad,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_gradients |
update_normal_vectors |
update_JxW_values));
FEValues<2> fe_v (dof_handler->get_fe (), qproject,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_gradients |
update_JxW_values));
FEValues<2> fe_v_n (dof_handler->get_fe (), qproject,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_gradients |
update_JxW_values));
QGauss<1> quad (6);
FEFaceValues<2> fe_v_face (dof_handler->get_fe (), quad,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_gradients |
update_normal_vectors |
update_JxW_values));
FEFaceValues<2> fe_v_face_n (dof_handler->get_fe (), quad,
UpdateFlags(update_values |
- update_q_points |
+ update_quadrature_points |
update_gradients |
update_normal_vectors |
update_JxW_values));
QTrapez<1> q_trapez;
const unsigned int div=10;
QIterated<dim> q(q_trapez, div);
- FEValues<dim> fe(fe_rt, q, update_values|update_gradients|update_q_points);
+ FEValues<dim> fe(fe_rt, q, update_values|update_gradients|update_quadrature_points);
fe.reinit(c);
for (unsigned int q_point=0; q_point<q.size(); ++q_point)
QTrapez<1> q_trapez;
QIterated<dim-1> q(q_trapez, div);
- FEFaceValues<dim> fe(mapping, finel, q, UpdateFlags(uflags | update_q_points));
- FESubfaceValues<dim> sub(mapping, finel, q, UpdateFlags(uflags | update_q_points));
+ FEFaceValues<dim> fe(mapping, finel, q, UpdateFlags(uflags | update_quadrature_points));
+ FESubfaceValues<dim> sub(mapping, finel, q, UpdateFlags(uflags | update_quadrature_points));
sprintf(fname, "Face%dd-%s", dim, name);
deallog.push(fname);
QTrapez<2> quadrature;
FE_Q<3> fe(1);
FEFaceValues<3> fe_face_values1 (fe, quadrature,
- update_q_points | update_JxW_values |
+ update_quadrature_points | update_JxW_values |
update_normal_vectors);
FESubfaceValues<3> fe_face_values2 (fe, quadrature,
- update_q_points | update_JxW_values |
+ update_quadrature_points | update_JxW_values |
update_normal_vectors);
DoFHandler<3> dof_handler (tria);
QTrapez<2> quadrature;
FE_Q<3> fe(1);
FEFaceValues<3> fe_face_values1 (fe, quadrature,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
FEFaceValues<3> fe_face_values2 (fe, quadrature,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
DoFHandler<3> dof_handler (tria);
dof_handler.distribute_dofs (fe);
QTrapez<2> quadrature;
FE_Q<3> fe(1);
FEFaceValues<3> fe_face_values1 (fe, quadrature,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
FEFaceValues<3> fe_face_values2 (fe, quadrature,
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
DoFHandler<3> dof_handler (tria);
dof_handler.distribute_dofs (fe);
hp::FEValues<dim> hp_fe_values (fe_collection,
quadrature_collection,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const RightHandSide<dim> rhs_function;
hp::FEValues<dim> hp_fe_values (fe_collection,
quadrature_collection,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const RightHandSide<dim> rhs_function;
const UpdateFlags update_flags = update_values
| update_gradients
- | update_q_points
+ | update_quadrature_points
| update_JxW_values;
const UpdateFlags face_update_flags = update_values
- | update_q_points
+ | update_quadrature_points
| update_JxW_values
| update_normal_vectors;
const UpdateFlags update_flags = update_values
| update_gradients
- | update_q_points
+ | update_quadrature_points
| update_JxW_values;
const UpdateFlags face_update_flags = update_values
- | update_q_points
+ | update_quadrature_points
| update_JxW_values
| update_normal_vectors;
assemble_rhs (Vector<double> &rhs) const
{
hp::FEValues<dim> fe_values (*this->fe, *this->quadrature,
- update_values | update_q_points |
+ update_values | update_quadrature_points |
update_JxW_values);
const unsigned int dofs_per_cell = (*this->fe)[0].dofs_per_cell;
hp::FEValues<dim> x_fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe[0].dofs_per_cell;
const unsigned int n_q_points = quadrature_formula[0].size();
hp::FEValues<dim> x_fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe[0].dofs_per_cell;
const unsigned int n_q_points = quadrature_formula[0].size();
hp::FEValues<dim> x_fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe[0].dofs_per_cell;
const unsigned int n_q_points = quadrature_formula[0].size();
hp::FEValues<dim> x_fe_values (*fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
hp::FEFaceValues<dim> x_fe_face_values (*fe, face_quadrature_formula,
- update_values | update_q_points |
+ update_values | update_quadrature_points |
update_normal_vectors | update_JxW_values);
const RightHandSide<dim> right_hand_side;
hp::FEValues<dim> x_fe_values (fe, quadrature_formula,
update_values | update_gradients |
- update_q_points | update_JxW_values);
+ update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe[0].dofs_per_cell;
const unsigned int n_q_points = quadrature_formula[0].size();
FEValues<dim> fe (mapping, element, quadrature,
update_values | update_gradients
- | update_q_points | update_JxW_values);
+ | update_quadrature_points | update_JxW_values);
std::vector <types::global_dof_index> global_dofs (element.dofs_per_cell);
std::vector <double> function (quadrature.size());
Assert (q.size() == 1, ExcInternalError());
FEFaceValues<dim> fe_values(mapping, fe, q,
- UpdateFlags(update_q_points
+ UpdateFlags(update_quadrature_points
| update_JxW_values
| update_values
| update_gradients
Assert (q.size() == 1, ExcInternalError());
FESubfaceValues<dim> fe_values(mapping, fe, q,
- UpdateFlags(update_q_points
+ UpdateFlags(update_quadrature_points
| update_JxW_values
| update_values
| update_gradients
Assert (q.size() == 1, ExcInternalError());
FEFaceValues<dim> fe_values(mapping, fe, q,
- UpdateFlags(update_q_points
+ UpdateFlags(update_quadrature_points
| update_JxW_values
| update_values
| update_gradients
Assert (q.size() == 1, ExcInternalError());
FESubfaceValues<dim> fe_values(mapping, fe, q,
- UpdateFlags(update_q_points
+ UpdateFlags(update_quadrature_points
| update_JxW_values
| update_values
| update_gradients
QTrapez<1> q_trapez;
QIterated<dim> q(q_trapez, div);
FEValues<dim> fe_values(mapping, fe, q,
- UpdateFlags(update_q_points
+ UpdateFlags(update_quadrature_points
| update_JxW_values));
fe_values.reinit(cell);
const unsigned int nq = (unsigned int) (.01 + std::pow(q.size(), 1./(dim-1)));
FEFaceValues<dim> fe_values(mapping, fe, q,
- UpdateFlags(update_q_points
+ UpdateFlags(update_quadrature_points
| update_JxW_values
| update_normal_vectors));
const unsigned int nq = (unsigned int) (.01 + std::pow(q.size(), 1./(dim-1)));
FESubfaceValues<dim> fe_values(mapping, fe, q,
- UpdateFlags(update_q_points
+ UpdateFlags(update_quadrature_points
| update_normal_vectors));
for (unsigned int face_nr=0;
face_nr < GeometryInfo<dim>::faces_per_cell;
const QTrapez<1> quadrature;
const MappingC1<2> c1_mapping;
FEFaceValues<2> c1_values (c1_mapping, fe, quadrature,
- update_q_points | update_normal_vectors);
+ update_quadrature_points | update_normal_vectors);
// to compare with, also print the
// normal vectors as generated by a
// cubic mapping
const MappingQ<2> q3_mapping(3);
FEFaceValues<2> q3_values (q3_mapping, fe, quadrature,
- update_q_points | update_normal_vectors);
+ update_quadrature_points | update_normal_vectors);
for (DoFHandler<2>::active_cell_iterator cell=dof_handler.begin_active();
cell!=dof_handler.end(); ++cell)
Vector<double> solution;
MappingQGeneric<dim> mapping(1);
QMidpoint<dim> q_midpoint;
- FEValues<dim> fe_values(mapping, fe, q_midpoint, update_q_points);
+ FEValues<dim> fe_values(mapping, fe, q_midpoint, update_quadrature_points);
dof_handler.distribute_dofs (fe);
solution.reinit(dof_handler.n_dofs());
FEValues<dim> fe (mapping, element, quadrature,
update_values
- | update_q_points | update_JxW_values);
+ | update_quadrature_points | update_JxW_values);
std::vector <types::global_dof_index> global_dofs (element.dofs_per_cell);
std::vector <double> function (quadrature.size());
FEValues<dim> fe2 (mapping, element, quadrature,
update_values | update_gradients
| update_second_derivatives
- | update_q_points | update_JxW_values);
+ | update_quadrature_points | update_JxW_values);
double l2 = 0.;
double h1 = 0.;
const std::vector<Point<dim-1> > &
unit_support_points = fe.get_unit_face_support_points();
Quadrature<dim-1> quadrature(unit_support_points);
- FEFaceValues<dim, dim> fe_face_values(fe, quadrature, update_q_points);
+ FEFaceValues<dim, dim> fe_face_values(fe, quadrature, update_quadrature_points);
typename DoFHandler<dim>::active_cell_iterator
cell = dof.begin_active(),
endc = dof.end();
const std::vector<Point<dim-1> > &
unit_support_points = fe.get_unit_face_support_points();
Quadrature<dim-1> quadrature(unit_support_points);
- FEFaceValues<dim, dim> fe_face_values(fe, quadrature, update_q_points);
+ FEFaceValues<dim, dim> fe_face_values(fe, quadrature, update_quadrature_points);
typename DoFHandler<dim>::active_cell_iterator
cell = dof.begin_active(),
endc = dof.end();