From b588c259ae1ac7d4a1c2ccddcff32fb749fed5b0 Mon Sep 17 00:00:00 2001 From: bangerth Date: Fri, 12 Aug 2011 23:09:22 +0000 Subject: [PATCH] Make FEFaceValues in 1d also work with MappingCartesian. git-svn-id: https://svn.dealii.org/trunk@24056 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/include/deal.II/base/quadrature.h | 5 +---- deal.II/source/fe/mapping_cartesian.cc | 12 ++++++++++++ 2 files changed, 13 insertions(+), 4 deletions(-) diff --git a/deal.II/include/deal.II/base/quadrature.h b/deal.II/include/deal.II/base/quadrature.h index 1b94cf5c32..a2aba45f49 100644 --- a/deal.II/include/deal.II/base/quadrature.h +++ b/deal.II/include/deal.II/base/quadrature.h @@ -1,7 +1,7 @@ //--------------------------------------------------------------------------- // $Id$ // -// Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2009, 2010 by the deal.II authors +// Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2009, 2010, 2011 by the deal.II authors // // This file is subject to QPL and may not be distributed // without copyright and license information. Please refer @@ -386,7 +386,6 @@ inline const Point & Quadrature::point (const unsigned int i) const { - Assert (dim != 0, ExcNotImplemented()); AssertIndexRange(i, size()); return quadrature_points[i]; } @@ -397,7 +396,6 @@ template double Quadrature::weight (const unsigned int i) const { - Assert (dim != 0, ExcNotImplemented()); AssertIndexRange(i, size()); return weights[i]; } @@ -409,7 +407,6 @@ inline const std::vector > & Quadrature::get_points () const { - Assert (dim > 0, ExcInternalError()); return quadrature_points; } diff --git a/deal.II/source/fe/mapping_cartesian.cc b/deal.II/source/fe/mapping_cartesian.cc index 5fdb56212f..2e5d400b9e 100644 --- a/deal.II/source/fe/mapping_cartesian.cc +++ b/deal.II/source/fe/mapping_cartesian.cc @@ -258,6 +258,18 @@ MappingCartesian::compute_fill (const typename Triangulation + normals[GeometryInfo<1>::faces_per_cell] + = { Point(-1.), + Point( 1.) }; + std::fill (normal_vectors.begin(), + normal_vectors.end(), + normals[face_no]); + break; + } + case 2: { static const Point -- 2.39.5