1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
7 -------------------------------------------------------------------------------
9 This file is part of OpenFOAM.
11 OpenFOAM is free software: you can redistribute it and/or modify it
12 under the terms of the GNU General Public License as published by
13 the Free Software Foundation, either version 3 of the License, or
14 (at your option) any later version.
16 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 You should have received a copy of the GNU General Public License
22 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
24 \*---------------------------------------------------------------------------*/
26 #include "volPointInterpolation.H"
27 #include "sampledCuttingPlane.H"
28 #include "volFieldsFwd.H"
29 #include "pointFields.H"
30 #include "volPointInterpolation.H"
32 // * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
35 Foam::tmp<Foam::Field<Type> >
36 Foam::sampledCuttingPlane::sampleField
38 const GeometricField<Type, fvPatchField, volMesh>& vField
41 return tmp<Field<Type> >(new Field<Type>(vField, surface().meshCells()));
46 Foam::tmp<Foam::Field<Type> >
47 Foam::sampledCuttingPlane::interpolateField
49 const interpolation<Type>& interpolator
52 // Get fields to sample. Assume volPointInterpolation!
53 const GeometricField<Type, fvPatchField, volMesh>& volFld =
56 if (subMeshPtr_.valid())
58 tmp<GeometricField<Type, fvPatchField, volMesh> > tvolSubFld =
59 subMeshPtr_().interpolate(volFld);
61 const GeometricField<Type, fvPatchField, volMesh>& volSubFld =
64 tmp<GeometricField<Type, pointPatchField, pointMesh> > tpointSubFld =
65 volPointInterpolation::New(volSubFld.mesh()).interpolate(volSubFld);
68 return surface().interpolate
72 ? pointAverage(tpointSubFld())()
80 tmp<GeometricField<Type, pointPatchField, pointMesh> > tpointFld
82 volPointInterpolation::New(volFld.mesh()).interpolate(volFld)
86 return surface().interpolate
90 ? pointAverage(tpointFld())()
99 // ************************************************************************* //