From 6981516d9d8cf4e67dcbd7f79a9d93c143d29681 Mon Sep 17 00:00:00 2001 From: Markus Buerg Date: Thu, 17 Mar 2011 19:57:18 +0000 Subject: [PATCH] Fix bug in get_face_interpolation_matrix. git-svn-id: https://svn.dealii.org/trunk@23483 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/source/fe/fe_nedelec.cc | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/deal.II/source/fe/fe_nedelec.cc b/deal.II/source/fe/fe_nedelec.cc index 09fed4bf79..d1b96aa9b5 100644 --- a/deal.II/source/fe/fe_nedelec.cc +++ b/deal.II/source/fe/fe_nedelec.cc @@ -2491,13 +2491,14 @@ const for (unsigned int j = 0; j < deg; ++j) { interpolation_matrix - (i + (j + GeometryInfo<2>::lines_per_cell) * source_fe.degree, - i + (j + GeometryInfo<2>::lines_per_cell) * this->degree) + ((i + GeometryInfo::lines_per_face) * (source_fe.degree - 1) + + j + GeometryInfo::lines_per_face, + (i + GeometryInfo::lines_per_face) * deg + j + + GeometryInfo::lines_per_face) = 1.0; interpolation_matrix - ((i * (source_fe.degree - 1) - + GeometryInfo<2>::lines_per_cell) * source_fe.degree + j, - (i * deg + GeometryInfo<2>::lines_per_cell) * this->degree) + (i + (j + GeometryInfo::lines_per_face) * source_fe.degree, + i + (j + GeometryInfo::lines_per_face) * this->degree) = 1.0; } } -- 2.39.5