ENH: boundBox: Add function to return the faces.

src/OpenFOAM/meshes/boundBox/boundBox.C

@@ -163,6 +163,49 @@ Foam::tmp<Foam::pointField> Foam::boundBox::points() const
 }
 
 
+Foam::faceList Foam::boundBox::faces()
+{
+    faceList faces(6);
+
+    forAll(faces, fI)
+    {
+        faces[fI].setSize(4);
+    }
+
+    faces[0][0] = 0;
+    faces[0][1] = 1;
+    faces[0][2] = 2;
+    faces[0][3] = 3;
+
+    faces[1][0] = 2;
+    faces[1][1] = 6;
+    faces[1][2] = 7;
+    faces[1][3] = 3;
+
+    faces[2][0] = 0;
+    faces[2][1] = 4;
+    faces[2][2] = 5;
+    faces[2][3] = 1;
+
+    faces[3][0] = 4;
+    faces[3][1] = 7;
+    faces[3][2] = 6;
+    faces[3][3] = 5;
+
+    faces[4][0] = 3;
+    faces[4][1] = 7;
+    faces[4][2] = 4;
+    faces[4][3] = 0;
+
+    faces[5][0] = 1;
+    faces[5][1] = 5;
+    faces[5][2] = 6;
+    faces[5][3] = 2;
+
+    return faces;
+}
+
+
 void Foam::boundBox::inflate(const scalar s)
 {
     vector ext = vector::one*s*mag();

src/OpenFOAM/meshes/boundBox/boundBox.H

@@ -33,6 +33,7 @@ Description
 #define boundBox_H
 
 #include "pointField.H"
+#include "faceList.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -160,6 +161,9 @@ public:
             //- Return corner points in an order corresponding to a 'hex' cell
             tmp<pointField> points() const;
 
+            //- Return faces with correct point order
+            static faceList faces();
+
 
         // Manipulate