FreeFOAM The Cross-Platform CFD Toolkit
lineDivide.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 1991-2010 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
13  the Free Software Foundation, either version 3 of the License, or
14  (at your 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, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 #include <OpenFOAM/error.H>
27 
28 #include "lineDivide.H"
29 #include "curvedEdge.H"
30 
31 
32 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
33 
35 (
36  const curvedEdge& cedge,
37  const label ndiv,
38  const scalar xratio
39 )
40 :
41  points_(ndiv + 1),
42  divisions_(ndiv + 1)
43 {
44  divisions_[0] = 0.0;
45  divisions_[ndiv] = 1.0;
46 
47  // calculate the spacing
48  if (xratio == 1.0)
49  {
50  for (label i=1; i < ndiv; i++)
51  {
52  divisions_[i] = scalar(i)/ndiv;
53  }
54  }
55  else
56  {
57  for (label i=1; i < ndiv; i++)
58  {
59  divisions_[i] = (1.0 - pow(xratio, i))/(1.0 - pow(xratio, ndiv));
60  }
61  }
62 
63  // calculate the points
64  for (label i=0; i <= ndiv; i++)
65  {
66  points_[i] = cedge.position(divisions_[i]);
67  }
68 }
69 
70 
71 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
72 
74 {
75  return points_;
76 }
77 
78 
80 {
81  return divisions_;
82 }
83 
84 
85 // ************************ vim: set sw=4 sts=4 et: ************************ //