From: hartmann Date: Thu, 1 Jul 1999 19:52:31 +0000 (+0000) Subject: implementation of #get(_unit)_support_points# X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=7a035f7f4b108ffdc252f20b484be1345e3b58f3;p=dealii-svn.git implementation of #get(_unit)_support_points# git-svn-id: https://svn.dealii.org/trunk@1524 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/fe/fe_system.cc b/deal.II/deal.II/source/fe/fe_system.cc index fa4a82536a..2ecc79a584 100644 --- a/deal.II/deal.II/source/fe/fe_system.cc +++ b/deal.II/deal.II/source/fe/fe_system.cc @@ -561,40 +561,83 @@ FESystem::shape_grad_grad (const unsigned int i, template -void FESystem::get_unit_support_points (vector > &/*support_points*/) const +void FESystem::get_unit_support_points ( + vector > &unit_support_points) const { - Assert(false, ExcNotImplemented()); -/* - Assert (support_points.size() == total_dofs, - ExcWrongFieldDimension (support_points.size(), total_dofs)); + Assert(unit_support_points.size() == total_dofs, + ExcWrongFieldDimension (unit_support_points.size(), total_dofs)); - vector > base_support_points (base_element->total_dofs); - base_element->get_unit_support_points (base_support_points); - - for (unsigned int i=0; itotal_dofs; ++i) - for (unsigned int n=0; n > base_unit_support_points (base_element(0).total_dofs); + unsigned int component = 0; + for (unsigned int base_el=0 ; base_el > > base_us_points(n_base_elements()); +// for (unsigned int base_el=0 ; base_el -void FESystem::get_support_points (const DoFHandler::cell_iterator &/*cell*/, - vector > &/*support_points*/) const +void FESystem::get_support_points (const DoFHandler::cell_iterator &cell, + vector > &support_points) const { - Assert(false, ExcNotImplemented()); -/* - Assert (support_points.size() == total_dofs, - ExcWrongFieldDimension (support_points.size(), total_dofs)); + Assert(support_points.size() == total_dofs, + ExcWrongFieldDimension (support_points.size(), total_dofs)); - vector > base_support_points (base_element->total_dofs); - base_element->get_support_points (cell, base_support_points); - - for (unsigned int i=0; itotal_dofs; ++i) - for (unsigned int n=0; n > base_support_points (base_element(0).total_dofs); + unsigned int component = 0; + for (unsigned int base_el=0 ; base_el