message: Setting configuration mask for FaceElementTransformations in GetSharedFaceTransformations
diff: @@ -2404,6 +2404,11 @@ GetSharedFaceTransformations(int sf, bool fill2)
  bool is_slave = Nonconforming() && IsSlaveFace(face_info);
  bool is_ghost = Nonconforming() && FaceNo >= GetNumFaces();
  + int mask = 0;
  + FaceElemTr.SetConfigurationMask(0);
  + FaceElemTr.Elem1 = NULL;
  + FaceElemTr.Elem2 = NULL;
  +
  NCFaceInfo* nc_info = NULL;
  if (is_slave) { nc_info = &nc_faces_info[face_info.NCFace]; }
  @@ -2415,6 +2420,7 @@ GetSharedFaceTransformations(int sf, bool fill2)
  FaceElemTr.Elem1No = face_info.Elem1No;
  GetElementTransformation(FaceElemTr.Elem1No, &Transformation);
  FaceElemTr.Elem1 = &Transformation;
  + mask += 1;
  // setup the transformation for the second (neighbor) element
  if (fill2)
  @@ -2422,6 +2428,7 @@ GetSharedFaceTransformations(int sf, bool fill2)
  FaceElemTr.Elem2No = -1 - face_info.Elem2No;
  GetFaceNbrElementTransformation(FaceElemTr.Elem2No, &Transformation2);
  FaceElemTr.Elem2 = &Transformation2;
  + mask += 2;
  }
  else
  {
  @@ -2433,6 +2440,7 @@ GetSharedFaceTransformations(int sf, bool fill2)
  {
  GetFaceTransformation(FaceNo, &FaceElemTr);
  // NOTE: The above call overwrites FaceElemTr.Loc1
  + mask += 16;
  }
  else
  {
  @@ -2443,12 +2451,14 @@ GetSharedFaceTransformations(int sf, bool fill2)
  int elem_type = GetElementType(face_info.Elem1No);
  GetLocalFaceTransformation(face_type, elem_type, FaceElemTr.Loc1.Transf,
  face_info.Elem1Inf);
  + mask += 4;
  if (fill2)
  {
  elem_type = face_nbr_elements[FaceElemTr.Elem2No]->GetType();
  GetLocalFaceTransformation(face_type, elem_type, FaceElemTr.Loc2.Transf,
  face_info.Elem2Inf);
  + mask += 8;
  }
  // adjust Loc1 or Loc2 of the master face if this is a slave face
  @@ -2479,6 +2489,8 @@ GetSharedFaceTransformations(int sf, bool fill2)
  GetGhostFaceTransformation(&FaceElemTr, face_type, face_geom);
  }
  + FaceElemTr.SetConfigurationMask(mask);
  +
  return &FaceElemTr;
  }
  