00001 /*---------------------------------------------------------------------------*\ 00002 ========= | 00003 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox 00004 \\ / O peration | 00005 \\ / A nd | Copyright (C) 1991-2010 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 00013 the Free Software Foundation, either version 3 of the License, or 00014 (at your 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, see <http://www.gnu.org/licenses/>. 00023 00024 \*---------------------------------------------------------------------------*/ 00025 00026 #include <OpenFOAM/error.H> 00027 00028 #include "lineDivide.H" 00029 #include "curvedEdge.H" 00030 00031 00032 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // 00033 00034 Foam::lineDivide::lineDivide 00035 ( 00036 const curvedEdge& cedge, 00037 const label ndiv, 00038 const scalar xratio 00039 ) 00040 : 00041 points_(ndiv + 1), 00042 divisions_(ndiv + 1) 00043 { 00044 divisions_[0] = 0.0; 00045 divisions_[ndiv] = 1.0; 00046 00047 // calculate the spacing 00048 if (xratio == 1.0) 00049 { 00050 for (label i=1; i < ndiv; i++) 00051 { 00052 divisions_[i] = scalar(i)/ndiv; 00053 } 00054 } 00055 else 00056 { 00057 for (label i=1; i < ndiv; i++) 00058 { 00059 divisions_[i] = (1.0 - pow(xratio, i))/(1.0 - pow(xratio, ndiv)); 00060 } 00061 } 00062 00063 // calculate the points 00064 for (label i=0; i <= ndiv; i++) 00065 { 00066 points_[i] = cedge.position(divisions_[i]); 00067 } 00068 } 00069 00070 00071 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // 00072 00073 const Foam::pointField& Foam::lineDivide::points() const 00074 { 00075 return points_; 00076 } 00077 00078 00079 const Foam::scalarList& Foam::lineDivide::lambdaDivisions() const 00080 { 00081 return divisions_; 00082 } 00083 00084 00085 // ************************ vim: set sw=4 sts=4 et: ************************ //
