From c94265160977dd4278a8ebc12c18a0c78af05998 Mon Sep 17 00:00:00 2001 From: guido Date: Tue, 24 May 2005 13:07:00 +0000 Subject: [PATCH] use identity for embedding git-svn-id: https://svn.dealii.org/trunk@10725 0785d39b-7218-0410-832d-ea1e28bc413d --- .../deal.II/source/fe/fe_dgp_nonparametric.cc | 19 ++++++------------- 1 file changed, 6 insertions(+), 13 deletions(-) diff --git a/deal.II/deal.II/source/fe/fe_dgp_nonparametric.cc b/deal.II/deal.II/source/fe/fe_dgp_nonparametric.cc index c0911a2231..43af89dcf0 100644 --- a/deal.II/deal.II/source/fe/fe_dgp_nonparametric.cc +++ b/deal.II/deal.II/source/fe/fe_dgp_nonparametric.cc @@ -37,20 +37,13 @@ FE_DGPNonparametric::FE_DGPNonparametric (const unsigned int degree) degree(degree), polynomial_space (Polynomials::Legendre::generate_complete_basis(degree)) { - // if defined, copy over matrices - // from precomputed arrays -// if ((degree < Matrices::n_embedding_matrices) && -// (Matrices::embedding[degree][0] != 0)) -// for (unsigned int c=0; c::children_per_cell; ++c) -// this->prolongation[c].fill (Matrices::embedding[degree][c]); -// else -// for (unsigned int i=0; i::children_per_cell;++i) -// this->prolongation[i].reinit(0,0); - // since not implemented, set to - // "empty" + const unsigned int n_dofs = this->dofs_per_cell; for (unsigned int i=0;i::children_per_cell;++i) - this->prolongation[i].reinit(0, 0); - + { + this->prolongation[i].reinit(n_dofs, n_dofs); + for (unsigned int j=0;jprolongation[i](j,j) = 1.; + } // restriction can be defined // through projection for // discontinuous elements, but is -- 2.39.5