src/meshLibrary/cartesianMesh/cartesianMeshExtractor/cartesianMeshExtractorPointsAndAddressing.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | cfMesh: A library for mesh generation
   4    \\    /   O peration     |
   5     \\  /    A nd           | Author: Franjo Juretic (franjo.juretic@c-fields.com)
   6      \\/     M anipulation  | Copyright (C) Creative Fields, Ltd.
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of cfMesh.
  10
  11     cfMesh 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 3 of the License, or (at your
  14     option) any later version.
  15
  16     cfMesh 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 cfMesh.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Description
  25
  26 \*---------------------------------------------------------------------------*/
  27
  28 #include "cartesianMeshExtractor.H"
  29 #include "meshOctree.H"
  30
  31 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  32
  33 namespace Foam
  34 {
  35
  36 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  37 // Private member functions
  38
  39 void cartesianMeshExtractor::createPointsAndAddressing()
  40 {
  41     Info << "Creating octree vertices" << endl;
  42
  43     Info << "Octree nodes " << octreeCheck_.numberOfNodes() << endl;
  44
  45     //- set the size of the point field
  46     pointFieldPMG& points = mesh_.points();
  47     points.setSize(octreeCheck_.numberOfNodes());
  48
  49     //- store vertices into the pointField
  50     const pointField& octreePoints = octreeCheck_.octreePoints();
  51
  52     forAll(points, pointI)
  53         points[pointI] = octreePoints[pointI];
  54
  55     Info << "Finished creating octree vertices" << endl;
  56 }
  57
  58 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  59
  60 } // End namespace Foam
  61
  62 // ************************************************************************* //