From a15b55649c61819d81785f2805c4b9515b4d65dd Mon Sep 17 00:00:00 2001 From: kanschat Date: Wed, 7 Dec 2011 17:05:37 +0000 Subject: [PATCH] remove unneeded specializations for codim 1 git-svn-id: https://svn.dealii.org/trunk@24796 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/source/fe/fe.cc | 51 +++++++---------------------------------- 1 file changed, 8 insertions(+), 43 deletions(-) diff --git a/deal.II/source/fe/fe.cc b/deal.II/source/fe/fe.cc index a61a2cb5de..e8db6a6bfe 100644 --- a/deal.II/source/fe/fe.cc +++ b/deal.II/source/fe/fe.cc @@ -927,38 +927,6 @@ FiniteElement::memory_consumption () const } -template<> -void -FiniteElement<1,2>::compute_2nd ( - const Mapping<1,2> &, - const Triangulation<1,2>::cell_iterator &, - const unsigned int, - Mapping<1,2>::InternalDataBase &, - InternalDataBase &, - FEValuesData<1,2> &) const -{ - - Assert(false, ExcNotImplemented()); -} - - - -template<> -void -FiniteElement<2,3>::compute_2nd ( - const Mapping<2,3> &, - const Triangulation<2,3>::cell_iterator &, - const unsigned int, - Mapping<2,3>::InternalDataBase &, - InternalDataBase &, - FEValuesData<2,3> &) const -{ - - Assert(false, ExcNotImplemented()); -} - - - template void FiniteElement::compute_2nd ( @@ -969,6 +937,7 @@ FiniteElement::compute_2nd ( InternalDataBase &fe_internal, FEValuesData &data) const { + AssertDimension(dim,spacedim); Assert ((fe_internal.update_each | fe_internal.update_once) & update_hessians, ExcInternalError()); @@ -999,8 +968,8 @@ FiniteElement::compute_2nd ( // direction (first index) and at // all q-points (second index) std::vector > > - diff_quot (dim, std::vector > (n_q_points)); - std::vector > diff_quot2 (n_q_points); + diff_quot (spacedim, std::vector > (n_q_points)); + std::vector > diff_quot2 (n_q_points); // for all nonzero components of // all shape functions at all @@ -1028,7 +997,7 @@ FiniteElement::compute_2nd ( // so find out about // the actual component // if necessary - Tensor<1,dim> right, left; + Tensor<1,spacedim> right, left; if (is_primitive(shape_index)) { right = fe_internal.differences[d1]->shape_grad(shape_index, q+offset); @@ -1077,7 +1046,7 @@ FiniteElement::compute_2nd ( // derivative from a // symmetric difference // approximation - for (unsigned int d=0; d::compute_2nd ( // unit cell, so transform it // to something on the real // cell - for (unsigned int d=0; d