diff options
Diffstat (limited to 'sci-libs/openfoam-solvers/files/OpenFOAM-1.5-ggi.patch')
| -rw-r--r-- | sci-libs/openfoam-solvers/files/OpenFOAM-1.5-ggi.patch | 233 |
1 files changed, 233 insertions, 0 deletions
diff --git a/sci-libs/openfoam-solvers/files/OpenFOAM-1.5-ggi.patch b/sci-libs/openfoam-solvers/files/OpenFOAM-1.5-ggi.patch new file mode 100644 index 000000000..155905018 --- /dev/null +++ b/sci-libs/openfoam-solvers/files/OpenFOAM-1.5-ggi.patch @@ -0,0 +1,233 @@ +--- OpenFOAM-1.5/src/OpenFOAM/meshes/polyMesh/polyPatches/basic/coupled/coupledPolyPatch.H 2009-03-12 16:57:20.000000000 +0100 ++++ OpenFOAM-1.5-dev/src/OpenFOAM/meshes/polyMesh/polyPatches/basic/coupled/coupledPolyPatch.H 2009-01-21 13:25:35.000000000 +0100 +@@ -54,9 +54,14 @@ + : + public polyPatch + { +- // Private data ++protected: ++ ++ // Protected data ++ ++ // Alowing derived classes to create their own transformation ++ // HJ, 24/Oct/2008 + +- //- offset (distance) vector from one side of the couple to the other ++ //- Offset (distance) vector from one side of the couple to the other + mutable vectorField separation_; + + //- Face transformation tensor +--- OpenFOAM-1.5/src/OpenFOAM/meshes/primitiveMesh/PrimitivePatch/PrimitivePatch.H 2009-03-12 16:57:20.000000000 +0100 ++++ OpenFOAM-1.5-dev/src/OpenFOAM/meshes/primitiveMesh/PrimitivePatch/PrimitivePatch.H 2009-01-21 13:25:37.000000000 +0100 +@@ -173,6 +173,11 @@ + mutable Field<PointType>* pointNormalsPtr_; + + ++ // Static Private Data ++ ++ static const bool nSquaredProjection_; ++ ++ + // Private Member Functions + + //- Calculate edges of the patch +@@ -412,6 +417,22 @@ + labelHashSet* setPtr = NULL + ) const; + ++ //- Write VTK patch, HJ, 14/Jan/2009 ++ static void writeVTK ++ ( ++ const fileName& name, ++ const FaceListType& faces, ++ const Field<PointType>& points ++ ); ++ ++ //- Write VTK patch normals, HJ, 14/Jan/2009 ++ static void writeVTKNormals ++ ( ++ const fileName& name, ++ const FaceListType& faces, ++ const Field<PointType>& points ++ ); ++ + + // Edit + +--- OpenFOAM-1.5/src/OpenFOAM/meshes/primitiveMesh/PrimitivePatch/PrimitivePatchCheck.C 2009-03-12 16:57:20.000000000 +0100 ++++ OpenFOAM-1.5-dev/src/OpenFOAM/meshes/primitiveMesh/PrimitivePatch/PrimitivePatchCheck.C 2009-01-21 13:25:37.000000000 +0100 +@@ -30,6 +30,7 @@ + #include "PrimitivePatch.H" + #include "Map.H" + #include "ListOps.H" ++#include "OFstream.H" + + // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +@@ -331,6 +332,166 @@ + } + + ++template ++< ++ class Face, ++ template<class> class FaceList, ++ class PointField, ++ class PointType ++> ++void PrimitivePatch<Face, FaceList, PointField, PointType>::writeVTK ++( ++ const fileName& name, ++ const FaceListType& faces, ++ const Field<PointType>& points ++) ++{ ++ // Write patch and points into VTK ++ OFstream mps(name + ".vtk"); ++ ++ mps << "# vtk DataFile Version 2.0" << nl ++ << name << ".vtk" << nl ++ << "ASCII" << nl ++ << "DATASET POLYDATA" << nl ++ << "POINTS " << points.size() << " float" << nl; ++ ++ // Write points ++ List<float> mlpBuffer(3*points.size()); ++ ++ label counter = 0; ++ forAll (points, i) ++ { ++ mlpBuffer[counter++] = float(points[i].x()); ++ mlpBuffer[counter++] = float(points[i].y()); ++ mlpBuffer[counter++] = float(points[i].z()); ++ } ++ ++ forAll (mlpBuffer, i) ++ { ++ mps << mlpBuffer[i] << ' '; ++ ++ if (i > 0 && (i % 10) == 0) ++ { ++ mps << nl; ++ } ++ } ++ ++ // Write faces ++ label nFaceVerts = 0; ++ ++ forAll (faces, faceI) ++ { ++ nFaceVerts += faces[faceI].size() + 1; ++ } ++ labelList mlfBuffer(nFaceVerts); ++ ++ counter = 0; ++ forAll (faces, faceI) ++ { ++ const Face& f = faces[faceI]; ++ ++ mlfBuffer[counter++] = f.size(); ++ ++ forAll (f, fpI) ++ { ++ mlfBuffer[counter++] = f[fpI]; ++ } ++ } ++ mps << nl; ++ ++ mps << "POLYGONS " << faces.size() << ' ' << nFaceVerts << endl; ++ ++ forAll (mlfBuffer, i) ++ { ++ mps << mlfBuffer[i] << ' '; ++ ++ if (i > 0 && (i % 10) == 0) ++ { ++ mps << nl; ++ } ++ } ++ mps << nl; ++} ++ ++ ++template ++< ++ class Face, ++ template<class> class FaceList, ++ class PointField, ++ class PointType ++> ++void PrimitivePatch<Face, FaceList, PointField, PointType>::writeVTKNormals ++( ++ const fileName& name, ++ const FaceListType& faces, ++ const Field<PointType>& points ++) ++{ ++ // Write patch and points into VTK ++ OFstream mps(name + ".vtk"); ++ ++ mps << "# vtk DataFile Version 2.0" << nl ++ << name << ".vtk" << nl ++ << "ASCII" << nl ++ << "DATASET POLYDATA" << nl ++ << "POINTS " << faces.size() << " float" << nl; ++ ++ // Write points ++ List<float> mlPointBuffer(3*faces.size()); ++ ++ label counter = 0; ++ forAll (faces, i) ++ { ++ const vector c = faces[i].centre(points); ++ ++ mlPointBuffer[counter++] = float(c.x()); ++ mlPointBuffer[counter++] = float(c.y()); ++ mlPointBuffer[counter++] = float(c.z()); ++ } ++ ++ forAll (mlPointBuffer, i) ++ { ++ mps << mlPointBuffer[i] << ' '; ++ ++ if (i > 0 && (i % 10) == 0) ++ { ++ mps << nl; ++ } ++ } ++ mps << nl; ++ ++ // Write normals ++ mps << "POINT_DATA " << faces.size() << nl ++ << "FIELD attributes " << 1 << nl ++ << "normals" << " 3 " ++ << faces.size() << " float" << nl; ++ ++ List<float> mlNormalBuffer(3*faces.size()); ++ ++ counter = 0; ++ forAll (faces, i) ++ { ++ const vector n = faces[i].normal(points); ++ ++ mlNormalBuffer[counter++] = float(n.x()); ++ mlNormalBuffer[counter++] = float(n.y()); ++ mlNormalBuffer[counter++] = float(n.z()); ++ } ++ ++ forAll (mlNormalBuffer, i) ++ { ++ mps << mlNormalBuffer[i] << ' '; ++ ++ if (i > 0 && (i % 10) == 0) ++ { ++ mps << nl; ++ } ++ } ++ mps << nl; ++} ++ ++ + // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + + } // End namespace Foam |