OpenFOAM logo
Open Source CFD Toolkit

interpolationCellPointI.H

Go to the documentation of this file.
00001 /*---------------------------------------------------------------------------*\
00002   =========                 |
00003   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
00004    \\    /   O peration     |
00005     \\  /    A nd           | Copyright (C) 1991-2005 OpenCFD Ltd.
00006      \\/     M anipulation  |
00007 -------------------------------------------------------------------------------
00008 License
00009     This file is part of OpenFOAM.
00010 
00011     OpenFOAM is free software; you can redistribute it and/or modify it
00012     under the terms of the GNU General Public License as published by the
00013     Free Software Foundation; either version 2 of the License, or (at your
00014     option) any later version.
00015 
00016     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
00017     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00018     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
00019     for more details.
00020 
00021     You should have received a copy of the GNU General Public License
00022     along with OpenFOAM; if not, write to the Free Software Foundation,
00023     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
00024 
00025 \*---------------------------------------------------------------------------*/
00026 
00027 namespace Foam
00028 {
00029 
00030 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
00031 
00032 template<class Type>
00033 inline Type interpolationCellPoint<Type>::interpolate
00034 (
00035     const cellPointWeight& cpw
00036 ) const
00037 {
00038     const FixedList<scalar, 4>& weights = cpw.weights();
00039     const FixedList<label, 3>& faceVertices = cpw.faceVertices();
00040 
00041     Type t = psip_[faceVertices[0]]*weights[0];
00042     t += psip_[faceVertices[1]]*weights[1];
00043     t += psip_[faceVertices[2]]*weights[2];
00044     t += this->psi_[cpw.cell()]*weights[3];
00045 
00046     return t;
00047 }
00048 
00049 
00050 template<class Type>
00051 inline Type interpolationCellPoint<Type>::interpolate
00052 (
00053     const vector& position,
00054     const label celli,
00055     const label facei
00056 ) const
00057 {
00058     return interpolate(cellPointWeight(this->pMesh_, position, celli, facei));
00059 }
00060 
00061 
00062 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00063 
00064 } // End namespace Foam
00065 
00066 // ************************************************************************* //
For further information go to www.openfoam.org