From: hartmann Date: Wed, 25 Jul 2001 19:48:53 +0000 (+0000) Subject: Simplify computation of laplace vector. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=03e3d49654940c95450e75cdd38828631c298b5f;p=dealii-svn.git Simplify computation of laplace vector. git-svn-id: https://svn.dealii.org/trunk@4860 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/fe/mapping_q.cc b/deal.II/deal.II/source/fe/mapping_q.cc index 2e93e12171..ae159135f5 100644 --- a/deal.II/deal.II/source/fe/mapping_q.cc +++ b/deal.II/deal.II/source/fe/mapping_q.cc @@ -620,51 +620,18 @@ MappingQ::compute_laplace_vector(std::vector > &lvs) co // S:=S_1*T S_1.mmult(S_1_T,T); - // compute the inner - // unit_support_points - std::vector > inner_unit_support_points(n_inner); - const double step = 1./degree; - const unsigned int z_end=(dim==3) ? degree : 2; - unsigned int iall=0; - for (unsigned int iz=1; iz &p=inner_unit_support_points[iall]; - p(0)=ix*step; - p(1)=iy*step; - if (dim==3) - p(2)=iz*step; - } - Assert(iall==n_inner, ExcInternalError()); - - // Compute the shape values at - // the inner - // unit_support_points - InternalData support_data(n_shape_functions); - support_data.shape_values.resize(n_shape_functions * n_inner); - - compute_shapes(inner_unit_support_points, support_data); - // Resize and initialize the // lvs lvs.resize(n_inner); - for (unsigned int unit_point=0; unit_point &lv=lvs[unit_point]; + std::vector &lv=lvs[i]; for (unsigned int k=0; k