src/finiteVolume/interpolation/interpolation/interpolationCellPointFace/interpolationCellPointFace.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     option) any later version.
  15
  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
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 Class
  26     Foam::interpolationCellPointFace
  27
  28 Description
  29     Foam::interpolationCellPointFace
  30
  31 \*---------------------------------------------------------------------------*/
  32
  33 #ifndef interpolationCellPointFace_H
  34 #define interpolationCellPointFace_H
  35
  36 #include "interpolation.H"
  37
  38 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  39
  40 namespace Foam
  41 {
  42
  43 /*---------------------------------------------------------------------------*\
  44                            Class interpolationCellPointFace Declaration
  45 \*---------------------------------------------------------------------------*/
  46
  47 template<class Type>
  48 class interpolationCellPointFace
  49 :
  50     public interpolation<Type>
  51 {
  52     // Private data
  53
  54         //- Interpolated volfield
  55         const GeometricField<Type, pointPatchField, pointMesh> psip_;
  56
  57         //- Linearly interpolated volfield
  58         const GeometricField<Type, fvsPatchField, surfaceMesh> psis_;
  59
  60         bool findTet
  61         (
  62             const vector& position,
  63             const label nFace,
  64             vector tetPoints[],
  65             label tetLabelCandidate[],
  66             label tetPointLabels[],
  67             scalar phi[],
  68             scalar phiCandidate[],
  69             label& closestFace,
  70             scalar& minDistance
  71         ) const;
  72
  73         bool findTriangle
  74         (
  75             const vector& position,
  76             const label nFace,
  77             label tetPointLabels[],
  78             scalar phi[]
  79         ) const;
  80
  81
  82 public:
  83
  84     //- Runtime type information
  85     TypeName("cellPointFace");
  86
  87
  88     // Constructors
  89
  90         //- Construct from components
  91         interpolationCellPointFace
  92         (
  93             const GeometricField<Type, fvPatchField, volMesh>& psi
  94         );
  95
  96
  97     // Member Functions
  98
  99         //- Interpolate field to the given point in the given cell
 100         Type interpolate
 101         (
 102             const vector& position,
 103             const label nCell,
 104             const label facei = -1
 105         ) const;
 106 };
 107
 108
 109 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 110
 111 } // End namespace Foam
 112
 113 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 114
 115 #ifdef NoRepository
 116 #   include "interpolationCellPointFace.C"
 117 #endif
 118
 119 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 120
 121 #endif
 122
 123 // ************************************************************************* //