From ac6c92630c57a915fa03409d34d9fff73685a492 Mon Sep 17 00:00:00 2001 From: wolf Date: Fri, 13 Sep 2002 22:27:41 +0000 Subject: [PATCH] Fix embedding matrices for 2d. git-svn-id: https://svn.dealii.org/trunk@6424 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/fe_nedelec_2d.cc | 42 +++++++++++----------- 1 file changed, 20 insertions(+), 22 deletions(-) diff --git a/deal.II/deal.II/source/fe/fe_nedelec_2d.cc b/deal.II/deal.II/source/fe/fe_nedelec_2d.cc index d24a486a27..d55ffd5905 100644 --- a/deal.II/deal.II/source/fe/fe_nedelec_2d.cc +++ b/deal.II/deal.II/source/fe/fe_nedelec_2d.cc @@ -23,48 +23,46 @@ // the child cell are obtained from the degrees of freedom on the // mother cell // -// TODO: [Anna] check whether the following paragraph is correct. if so, then please multiply the values in the four following matrices by two - // note the following: since the shape functions themselves and not // only the gradients are transformed using the mapping object from // the unit cell to the real cell, the actual values of the function // on the real cell is degree of freedom times value of the shape -// function on the unit cell times Jacobian. Thus, what has the DoF -// value 1 on the mother cell must have the DoF value 2 on the child -// cell since the latter is smaller by a (linear scaling) factor of -// two. +// function on the unit cell times inverse Jacobian. Thus, what has +// the DoF value 1 on the mother cell must have the DoF value 1/2 on +// the child cell since the latter is smaller by a (linear scaling) +// factor of two. namespace FE_Nedelec_2d { static const double q1_into_q1_refined_0[] = { - 1., 0, 0, 0, - 0, 0.5,0, 0.5, - 0.5, 0, 0.5,0, - 0, 0, 0, 1 + .5, 0, 0 , 0, + 0, 0.25,0, 0.25, + 0.25, 0, 0.25,0, + 0, 0, 0, .5 }; static const double q1_into_q1_refined_1[] = { - 1., 0., 0., 0., - 0., 1., 0., 0., - 0.5, 0., 0.5, 0., - 0., 0.5, 0., 0.5, + .5, 0., 0., 0., + 0., .5, 0., 0., + 0.25, 0., 0.25, 0., + 0., 0.25, 0., 0.25, }; static const double q1_into_q1_refined_2[] = { - 0.5, 0., 0.5, 0., - 0., 1., 0., 0., - 0., 0., 1., 0., - 0., 0.5, 0., 0.5, + 0.25, 0., 0.25, 0., + 0., .5, 0., 0., + 0., 0., .5, 0., + 0., 0.25, 0., 0.25, }; static const double q1_into_q1_refined_3[] = { - 0.5, 0., 0.5, 0., - 0., 0.5, 0., 0.5, - 0., 0., 1., 0., - 0., 0., 0., 1., + 0.25, 0., 0.25, 0., + 0., 0.25, 0., 0.25, + 0., 0., .5, 0., + 0., 0., 0., .5, }; }; // namespace FE_Nedelec_2d -- 2.39.5