From 2e441003586ff2a856cc4f07dc05d4500226c73c Mon Sep 17 00:00:00 2001 From: Guido Kanschat Date: Mon, 24 Oct 2005 14:51:36 +0000 Subject: [PATCH] Wolfgang's code removed again git-svn-id: https://svn.dealii.org/trunk@11654 0785d39b-7218-0410-832d-ea1e28bc413d --- .../source/fe/fe_raviart_thomas_nodal.cc | 24 +++++++++---------- 1 file changed, 11 insertions(+), 13 deletions(-) diff --git a/deal.II/deal.II/source/fe/fe_raviart_thomas_nodal.cc b/deal.II/deal.II/source/fe/fe_raviart_thomas_nodal.cc index 93628517a5..3c978fa9d2 100644 --- a/deal.II/deal.II/source/fe/fe_raviart_thomas_nodal.cc +++ b/deal.II/deal.II/source/fe/fe_raviart_thomas_nodal.cc @@ -377,20 +377,18 @@ FE_RaviartThomasNodal::initialize_support_points (const unsigned int deg) for (unsigned int d=0;d - quadrature (dim == 1 ? - QAnisotropic(high) : - (dim == 2 ? - QAnisotropic(((d==0) ? low : high), - ((d==1) ? low : high)) : - QAnisotropic(((d==0) ? low : high), - ((d==1) ? low : high), - ((d==2) ? low : high)))); + QAnisotropic* quadrature; + if (dim == 1) quadrature = new QAnisotropic(high); + if (dim == 2) quadrature = new QAnisotropic(((d==0) ? low : high), + ((d==1) ? low : high)); + if (dim == 3) quadrature = new QAnisotropic(((d==0) ? low : high), + ((d==1) ? low : high), + ((d==2) ? low : high)); + Assert(dim<=3, ExcNotImplemented()); - for (unsigned int k=0;kgeneralized_support_points[current++] = quadrature.point(k); + for (unsigned int k=0;kn_quadrature_points;++k) + this->generalized_support_points[current++] = quadrature->point(k); + delete quadrature; } Assert (current == this->dofs_per_cell, ExcInternalError()); } -- 2.39.5