From: guido Date: Sat, 28 May 2005 22:56:37 +0000 (+0000) Subject: fix node values in 2d and accomodate for Wolfgangs project_to_all_faces in 3d X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=6048c49d3e5d121270f1fffb917013b349fb5ac3;p=dealii-svn.git fix node values in 2d and accomodate for Wolfgangs project_to_all_faces in 3d git-svn-id: https://svn.dealii.org/trunk@10772 0785d39b-7218-0410-832d-ea1e28bc413d --- 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 8f57d435a2..4dad76b7ab 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 @@ -28,6 +28,7 @@ # include #endif +#include template FE_RaviartThomasNodal::FE_RaviartThomasNodal (const unsigned int deg) @@ -125,7 +126,7 @@ FE_RaviartThomasNodal::get_dpo_vector (const unsigned int deg) // the element is face-based and we have // (deg+1)^(dim-1) DoFs per face unsigned int dofs_per_face = 1; - for (unsigned int d=0; d::initialize_unit_support_points (const unsigned int d for (unsigned int k=0;kdofs_per_face;++k) this->unit_face_support_points[k] = face_points.point(k); Quadrature faces = QProjector::project_to_all_faces(face_points); - for (unsigned int k=0;kdofs_per_face*GeometryInfo::faces_per_cell;++k) this->unit_support_points[k] = faces.point(k); - current = faces.n_quadrature_points; + current = this->dofs_per_face*GeometryInfo::faces_per_cell; } + + if (deg==0) return; // In the interior, we need // anisotropic Gauss quadratures, // different for each direction. @@ -271,8 +275,8 @@ FE_RaviartThomasNodal::initialize_node_matrix () for (unsigned int i=0;ishape_value_component( i, this->unit_support_points[current], - GeometryInfo< dim >::unit_normal_direction[face]) - * GeometryInfo< dim >::unit_normal_orientation[face]; + GeometryInfo< dim >::unit_normal_direction[face]); +// * GeometryInfo< dim >::unit_normal_orientation[face]; ++current; } // Interior degrees of freedom in each direction