From: hartmann Date: Thu, 23 Sep 2004 10:45:39 +0000 (+0000) Subject: Extension of get_interpolation_matrix to arbitrary source fe's. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=5f14f43bc30ee70043ad5deb3ac8513e42d17b0e;p=dealii-svn.git Extension of get_interpolation_matrix to arbitrary source fe's. git-svn-id: https://svn.dealii.org/trunk@9674 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/fe/fe_dgp_monomial.cc b/deal.II/deal.II/source/fe/fe_dgp_monomial.cc index 32defe27b4..b1f5c9e177 100644 --- a/deal.II/deal.II/source/fe/fe_dgp_monomial.cc +++ b/deal.II/deal.II/source/fe/fe_dgp_monomial.cc @@ -102,7 +102,7 @@ namespace p[i](1)=points3d[start_index3d[k]+i][1]; p[i](2)=points3d[start_index3d[k]+i][2]; } - } + } } @@ -176,34 +176,49 @@ FE_DGPMonomial::clone() const template void FE_DGPMonomial:: -get_interpolation_matrix (const FiniteElementBase &x_source_fe, +get_interpolation_matrix (const FiniteElementBase &source_fe, FullMatrix &interpolation_matrix) const -{ - // this is only implemented, if the - // source FE is also a - // DGPMonomial element - AssertThrow ((x_source_fe.get_name().find ("FE_DGPMonomial<") == 0) - || - (dynamic_cast*>(&x_source_fe) != 0), - typename FiniteElementBase:: - ExcInterpolationNotImplemented()); - - // ok, source is a Q element, so - // we will be able to do the work - const FE_DGPMonomial &source_fe - = dynamic_cast&>(x_source_fe); - - const unsigned int m=interpolation_matrix.m(); - const unsigned int n=interpolation_matrix.n(); - Assert (m == this->dofs_per_cell, ExcDimensionMismatch (m, this->dofs_per_cell)); - Assert (n == source_fe.dofs_per_cell, ExcDimensionMismatch (n, source_fe.dofs_per_cell)); - - const unsigned int min_mn= - interpolation_matrix.m() *source_dgp_monomial + = dynamic_cast *>(&source_fe); + + if (source_dgp_monomial) + { + // ok, source_fe is a DGP_Monomial + // element. Then, the interpolation + // matrix is simple + const unsigned int m=interpolation_matrix.m(); + const unsigned int n=interpolation_matrix.n(); + Assert (m == this->dofs_per_cell, ExcDimensionMismatch (m, this->dofs_per_cell)); + Assert (n == source_dgp_monomial->dofs_per_cell, + ExcDimensionMismatch (n, source_dgp_monomial->dofs_per_cell)); + + const unsigned int min_mn= + interpolation_matrix.m() > unit_points(this->dofs_per_cell); + generate_unit_points(this->degree, unit_points); + + FullMatrix source_fe_matrix(unit_points.size(), source_fe.dofs_per_cell); + for (unsigned int j=0; j this_matrix(this->dofs_per_cell, this->dofs_per_cell); + for (unsigned int j=0; jdofs_per_cell; ++j) + for (unsigned int k=0; kpoly_space.compute_value (j, unit_points[k]); + + this_matrix.gauss_jordan(); + + this_matrix.mmult(interpolation_matrix, source_fe_matrix); + } }