expression N() corrected

AMDiS/src/expressions/coords_expr.hpp

@@ -176,9 +176,9 @@ namespace AMDiS
     {
       typedef WorldVector<double> value_type;
       mutable WorldVector<double> normal;
-      int side;
+      int boundary;
 
-      Normals(int side_) : side(side_) {}
+      Normals(int boundary_) : boundary(boundary_) {}
 
       template<typename List>
       void insertFeSpaces(List& feSpaces) const {}
@@ -187,19 +187,25 @@ namespace AMDiS
 
       template<typename OT>
       void initElement(OT* ot, const ElInfo* elInfo,
-	      SubAssembler* subAssembler, Quadrature *quad, 
-		    const BasisFunction *basisFct = NULL)
+		       SubAssembler* subAssembler, Quadrature *quad, 
+		       const BasisFunction *basisFct = NULL)
       {
-	elInfo->getNormal(side, normal);
+	int dim = elInfo->getMesh()->getDim();
+	for (int side = 0; side < dim+1; ++side) {
+	  if (elInfo->getBoundary(side) == boundary) {
+	    elInfo->getNormal(side, normal);
+	    break;
+	  }
+	}
       }
 
 
       template<typename OT>
       void initElement(OT* ot, const ElInfo* smallElInfo, const ElInfo* largeElInfo,
-	      SubAssembler* subAssembler, Quadrature *quad, 
-		    const BasisFunction *basisFct = NULL)
+		       SubAssembler* subAssembler, Quadrature *quad, 
+		       const BasisFunction *basisFct = NULL)
       {
-	smallElInfo->getNormal(side, normal);
+	initElement(ot, smallElInfo, subAssembler, quad, basisFct);
       }
 
       inline value_type operator()(const int& iq) const { return normal; }