From c8ea8198cdcfbea6d2a0f0bb2fd43adfaf4c697c Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Wed, 27 Mar 2013 15:49:58 +0000
Subject: [PATCH 1/2] STYLE: Updated solver header documentation

---
 .../solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C  | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/applications/solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C b/applications/solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C
index b58df5402e..67df5b61e0 100644
--- a/applications/solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C
+++ b/applications/solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C
@@ -22,7 +22,7 @@ License
     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
 Application
-    buoyantSimpleRadiationFoam
+    buoyantSimpleFoam
 
 Description
     Steady-state solver for buoyant, turbulent flow of compressible fluids,

From f303199ffa31af619dc2dc17399f72d7371c0108 Mon Sep 17 00:00:00 2001
From: laurence <laurence>
Date: Wed, 27 Mar 2013 17:14:57 +0000
Subject: [PATCH 2/2] BUG: PatchTools::sortedEdgeFaces(): vector now normalised
 before dot product

---
 .../PatchTools/PatchToolsSortEdges.C             | 16 ++++++++++------
 1 file changed, 10 insertions(+), 6 deletions(-)

diff --git a/src/OpenFOAM/meshes/primitiveMesh/PatchTools/PatchToolsSortEdges.C b/src/OpenFOAM/meshes/primitiveMesh/PatchTools/PatchToolsSortEdges.C
index 54306d7230..0877a50b5f 100644
--- a/src/OpenFOAM/meshes/primitiveMesh/PatchTools/PatchToolsSortEdges.C
+++ b/src/OpenFOAM/meshes/primitiveMesh/PatchTools/PatchToolsSortEdges.C
@@ -77,10 +77,12 @@ Foam::PatchTools::sortedEdgeFaces
             {
                 if (f0[fpI] != e.start())
                 {
-                    const vector faceEdgeDir = localPoints[f0[fpI]] - edgePt;
-                    const scalar angle = faceEdgeDir & e2;
+                    vector faceEdgeDir = localPoints[f0[fpI]] - edgePt;
+                    faceEdgeDir /= mag(faceEdgeDir) + VSMALL;
 
-                    if (angle < maxAngle)
+                    const scalar angle = e2 & faceEdgeDir;
+
+                    if (mag(angle) < maxAngle)
                     {
                         maxAngle = angle;
                         maxAngleEdgeDir = faceEdgeDir;
@@ -114,10 +116,12 @@ Foam::PatchTools::sortedEdgeFaces
                 {
                     if (f[fpI] != e.start())
                     {
-                        const vector faceEdgeDir = localPoints[f[fpI]] - edgePt;
-                        const scalar angle = faceEdgeDir & e2;
+                        vector faceEdgeDir = localPoints[f[fpI]] - edgePt;
+                        faceEdgeDir /= mag(faceEdgeDir) + VSMALL;
 
-                        if (angle < maxAngle)
+                        const scalar angle = e2 & faceEdgeDir;
+
+                        if (mag(angle) < maxAngle)
                         {
                             maxAngle = angle;
                             maxAngleEdgeDir = faceEdgeDir;