<h3>Specific improvements</h3>
<ol>
+ <li> Fixed: FE_TraceQ now provides unit support points.
+ <br>
+ (Martin Kronbichler, 2016/09/08)
+ </li>
+
<li> Fixed: Fix MPI_InitFinalize by correctly initializing and destroying
all p4est/libsc related objects by calls to sc_init(), p4est_init(), and
sc_finalize(); compatibility with p4est versions >1.1.
// Initialize face support points
this->unit_face_support_points = fe_q.get_unit_face_support_points();
+
+ // initialize unit support points (this makes it possible to assign initial
+ // values to FE_TraceQ). Note that we simply take the points of fe_q but
+ // skip the last ones which are associated with the interior of FE_Q.
+ this->unit_support_points.resize(this->dofs_per_cell);
+ for (unsigned int i=0; i<this->dofs_per_cell; ++i)
+ this->unit_support_points[i] = fe_q.get_unit_support_points()[i];
+
// Initialize constraint matrices
this->interface_constraints = fe_q.constraints();
}
#include <deal.II/fe/fe_dgq.h>
#include <deal.II/fe/fe_dgp.h>
#include <deal.II/fe/fe_face.h>
+#include <deal.II/fe/fe_trace.h>
#include <deal.II/fe/fe_nedelec.h>
#include <fstream>
{
check1 (FE_FaceQ<dim>(2));
check2 (FE_FaceQ<dim>(2), 0);
+ check1 (FE_TraceQ<dim>(2));
+ check2 (FE_TraceQ<dim>(2), 0);
}
DEAL::1 0.500000
DEAL::2 1.00000
DEAL::dim=2, face=ok
+DEAL::dim=2, cell=ok
+DEAL::dim=2, face=ok
+DEAL::0 0.00000 0.00000
+DEAL::1 1.00000 0.00000
+DEAL::2 0.00000 1.00000
+DEAL::3 1.00000 1.00000
+DEAL::4 0.00000 0.500000
+DEAL::5 1.00000 0.500000
+DEAL::6 0.500000 0.00000
+DEAL::7 0.500000 1.00000
+DEAL::dim=2, cell=ok
+DEAL::0 0.00000
+DEAL::1 1.00000
+DEAL::2 0.500000
+DEAL::dim=2, face=ok
DEAL::dim=3, cell=ok
DEAL::dim=3, face=ok
DEAL::dim=3, cell=ok
DEAL::7 0.500000 1.00000
DEAL::8 1.00000 1.00000
DEAL::dim=3, face=ok
+DEAL::dim=3, cell=ok
+DEAL::dim=3, face=ok
+DEAL::0 0.00000 0.00000 0.00000
+DEAL::1 1.00000 0.00000 0.00000
+DEAL::2 0.00000 1.00000 0.00000
+DEAL::3 1.00000 1.00000 0.00000
+DEAL::4 0.00000 0.00000 1.00000
+DEAL::5 1.00000 0.00000 1.00000
+DEAL::6 0.00000 1.00000 1.00000
+DEAL::7 1.00000 1.00000 1.00000
+DEAL::8 0.00000 0.500000 0.00000
+DEAL::9 1.00000 0.500000 0.00000
+DEAL::10 0.500000 0.00000 0.00000
+DEAL::11 0.500000 1.00000 0.00000
+DEAL::12 0.00000 0.500000 1.00000
+DEAL::13 1.00000 0.500000 1.00000
+DEAL::14 0.500000 0.00000 1.00000
+DEAL::15 0.500000 1.00000 1.00000
+DEAL::16 0.00000 0.00000 0.500000
+DEAL::17 1.00000 0.00000 0.500000
+DEAL::18 0.00000 1.00000 0.500000
+DEAL::19 1.00000 1.00000 0.500000
+DEAL::20 0.00000 0.500000 0.500000
+DEAL::21 1.00000 0.500000 0.500000
+DEAL::22 0.500000 0.00000 0.500000
+DEAL::23 0.500000 1.00000 0.500000
+DEAL::24 0.500000 0.500000 0.00000
+DEAL::25 0.500000 0.500000 1.00000
+DEAL::dim=3, cell=ok
+DEAL::0 0.00000 0.00000
+DEAL::1 1.00000 0.00000
+DEAL::2 0.00000 1.00000
+DEAL::3 1.00000 1.00000
+DEAL::4 0.00000 0.500000
+DEAL::5 1.00000 0.500000
+DEAL::6 0.500000 0.00000
+DEAL::7 0.500000 1.00000
+DEAL::8 0.500000 0.500000
+DEAL::dim=3, face=ok