From e6cf63aaf4b934467e481354409e4b171f0e3d74 Mon Sep 17 00:00:00 2001 From: Ralf Hartmann Date: Mon, 19 Mar 2001 09:54:06 +0000 Subject: [PATCH] Two fixes. git-svn-id: https://svn.dealii.org/trunk@4241 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/mapping_q.cc | 23 ++++++++++++++++++++--- 1 file changed, 20 insertions(+), 3 deletions(-) diff --git a/deal.II/deal.II/source/fe/mapping_q.cc b/deal.II/deal.II/source/fe/mapping_q.cc index 5a887dd409..223effc706 100644 --- a/deal.II/deal.II/source/fe/mapping_q.cc +++ b/deal.II/deal.II/source/fe/mapping_q.cc @@ -442,11 +442,18 @@ MappingQ::set_laplace_on_quad_vector(std::vector > &loq } else { - // not precomputed, then do so now + // not precomputed, then do so now if (dim==2) compute_laplace_vector(loqvs); - else - Assert(false, ExcNotImplemented()); + + // for dim==3 don't throw an + // ExcNotImplemented here to + // allow the creating of that + // MappingQ<3> object. But an + // ExcLaplaceVectorNotSet + // assertion is thrown when the + // apply_laplace_vector + // function is called. } // the sum of weights of the points @@ -986,6 +993,16 @@ MappingQ<3>::add_quad_support_points(const Triangulation<3>::cell_iterator &cell // separately if (lines_at_boundary>0) { + // call of function + // apply_laplace_vector + // increases size of b + // about 1. There + // resize b for the + // case the mentioned + // function was already + // called. + b.resize(4*degree); + // b is of size // 4*degree, make sure // that this is the -- 2.39.5