From 9dd55c02cb32eda9840bc604217a5a145621fec8 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 21 Nov 2011 13:59:42 +0000
Subject: [PATCH] ENH: Added further ouput to volPointInterp test app
---
.../Test-volPointInterpolation.C | 27 ++++++++++++++++---
1 file changed, 23 insertions(+), 4 deletions(-)
diff --git a/applications/test/volPointInterpolation/Test-volPointInterpolation.C b/applications/test/volPointInterpolation/Test-volPointInterpolation.C
index 1bbcd6da5e4..96ffeeac377 100644
--- a/applications/test/volPointInterpolation/Test-volPointInterpolation.C
+++ b/applications/test/volPointInterpolation/Test-volPointInterpolation.C
@@ -33,11 +33,10 @@ Application
int main(int argc, char *argv[])
{
+ #include "setRootCase.H"
-# include "setRootCase.H"
-
-# include "createTime.H"
-# include "createMesh.H"
+ #include "createTime.H"
+ #include "createMesh.H"
Info<< "Reading field p\n" << endl;
volScalarField p
@@ -67,9 +66,29 @@ int main(int argc, char *argv[])
mesh
);
+ const pointMesh& pMesh = pointMesh::New(mesh);
+ const pointBoundaryMesh& pbm = pMesh.boundary();
+
+ Info<< "pointMesh boundary" << nl;
+ forAll(pbm, patchI)
+ {
+ Info<< "patch=" << pbm[patchI].name()
+ << ", type=" << pbm[patchI].type()
+ << ", coupled=" << pbm[patchI].coupled()
+ << endl;
+ }
+
const volPointInterpolation& pInterp = volPointInterpolation::New(mesh);
+
pointScalarField pp(pInterp.interpolate(p));
+ Info<< pp.name() << " boundary" << endl;
+ forAll(pp.boundaryField(), patchI)
+ {
+ Info<< pbm[patchI].name() << " coupled="
+ << pp.boundaryField()[patchI].coupled()<< endl;
+ }
+
pp.write();
pointVectorField pU(pInterp.interpolate(U));
--
GitLab