]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
use identity for embedding
authorguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 24 May 2005 13:07:00 +0000 (13:07 +0000)
committerguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 24 May 2005 13:07:00 +0000 (13:07 +0000)
git-svn-id: https://svn.dealii.org/trunk@10725 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_dgp_nonparametric.cc

index c0911a2231a57da1fcad885cab66a43378c4ba9a..43af89dcf05944035bf80869386ddbd0d6a280e2 100644 (file)
@@ -37,20 +37,13 @@ FE_DGPNonparametric<dim>::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<GeometryInfo<dim>::children_per_cell; ++c)
-//        this->prolongation[c].fill (Matrices::embedding[degree][c]);
-//    else
-//      for (unsigned int i=0; i<GeometryInfo<dim>::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<GeometryInfo<dim>::children_per_cell;++i)
-    this->prolongation[i].reinit(0, 0);
-
+    {
+      this->prolongation[i].reinit(n_dofs, n_dofs);
+      for (unsigned int j=0;j<n_dofs;++j)
+       this->prolongation[i](j,j) = 1.;
+    }
                                    // restriction can be defined
                                    // through projection for
                                    // discontinuous elements, but is

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.