From: hartmann Date: Fri, 12 Mar 2004 10:16:28 +0000 (+0000) Subject: Implementation of generate_unit_points in 1d. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=6858ae41023c738f85767b17f23109ec7b267605;p=dealii-svn.git Implementation of generate_unit_points in 1d. git-svn-id: https://svn.dealii.org/trunk@8722 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/fe/fe_dgp_monomial.cc b/deal.II/deal.II/source/fe/fe_dgp_monomial.cc index 9178be1a3f..e9773a7817 100644 --- a/deal.II/deal.II/source/fe/fe_dgp_monomial.cc +++ b/deal.II/deal.II/source/fe/fe_dgp_monomial.cc @@ -79,10 +79,13 @@ namespace std::vector > &); template <> - void generate_unit_points (const unsigned int, - std::vector > &) + void generate_unit_points (const unsigned int k, + std::vector > &p) { - Assert(false, ExcNotImplemented()); + Assert(p.size()==k+1, ExcDimensionMismatch(p.size(), k+1)); + const double h = 1./k; + for (unsigned int i=0; i