of view, etc. The present implementation is only a rudimentary hack.
-Review the restruction matrices. I'm not really sure about their
+Review the restriction matrices. I'm not really sure about their
meaning and how they are defined, so they may be wrong for linear
elements and they are not implemented at all for quadratic ones.
+Remove the workaround with the BoundaryHelper class which was
+ introduced to work around a problem in the egcs-19980803 snapshot
+ (previous ones worked flawless)
+
+
+
+
DEAL:
/*---------------------------- boundary-function.h ---------------------------*/
#include <grid/point.h>
+#include <grid/geometry_info.h>
+
+
+
+/**
+ * Workaround for a bug in egcs snapshot 1998/08/03.
+ */
+template <int dim> struct BoundaryHelper;
+template <> struct BoundaryHelper<2> {
+ typedef const Point<2> *PointArray[GeometryInfo<2>::vertices_per_face];
+};
+
+
/**
* This class is used to represent a boundary to a triangulation.
* Typedef an array of the needed number
* of old points.
*/
- typedef const Point<dim>* PointArray[1<<(dim-1)];
-
+ typedef typename BoundaryHelper<dim>::PointArray PointArray;
+
/**
* This function calculates the position
* of the new vertex.
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