const Quadrature<dim-1> &q_boundary = QGauss<dim-1>(2),
const bool project_to_boundary_first = false);
+ /**
+ * Declaration of specialization
+ * of the previous function for
+ * 1d. At present, it is not
+ * implemented.
+ *
+ * The default value of the boundary
+ * quadrature formula is an invalid
+ * object since it makes no sense in 1d.
+ */
+ template <class VECTOR>
+ static void project (const DoFHandler<1> &dof,
+ const ConstraintMatrix &constraints,
+ const Quadrature<1> &quadrature,
+ const Function<1> &function,
+ VECTOR &vec,
+ const bool enforce_zero_boundary = false,
+ const Quadrature<0> &q_boundary = *invalid_face_quadrature,
+ const bool project_to_boundary_first = false);
+
/**
* Create a right hand side
* vector. Prior content of the
#include <vector>
#include <cmath>
+
+//TODO[RH]: Use StaticQ1Mapping where appropriate
+
template <int dim>
inline double sqr_point (const Tensor<1,dim> &p)
{
}
+template <class VECTOR>
+void VectorTools::project (const DoFHandler<1> &dof_handler,
+ const ConstraintMatrix &constraints,
+ const Quadrature<1> &quadrature,
+ const Function<1> &function,
+ VECTOR &vec_result,
+ const bool enforce_zero_boundary,
+ const Quadrature<0> &q_boundary,
+ const bool project_to_boundary_first)
+{
+ Assert (DEAL_II_COMPAT_MAPPING, ExcCompatibility("mapping"));
+ static const MappingQ1<1> mapping;
+ project (mapping, dof_handler, constraints, quadrature, function, vec_result,
+ enforce_zero_boundary, q_boundary, project_to_boundary_first);
+}
+
+
#endif