src/meshTools/triSurface/surfaceLocation/surfaceLocation.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   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 "surfaceLocation.H"
  27 #include "triSurface.H"
  28
  29 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  30
  31 Foam::vector Foam::surfaceLocation::normal(const triSurface& s) const
  32 {
  33     const vectorField& n = s.faceNormals();
  34
  35     if (elementType_ == triPointRef::NONE)
  36     {
  37         return n[index()];
  38     }
  39     else if (elementType_ == triPointRef::EDGE)
  40     {
  41         const labelList& eFaces = s.edgeFaces()[index()];
  42
  43         if (eFaces.size() == 1)
  44         {
  45             return n[eFaces[0]];
  46         }
  47         else
  48         {
  49             vector edgeNormal(vector::zero);
  50
  51             forAll(eFaces, i)
  52             {
  53                 edgeNormal += n[eFaces[i]];
  54             }
  55             return edgeNormal/(mag(edgeNormal) + VSMALL);
  56         }
  57     }
  58     else
  59     {
  60         return s.pointNormals()[index()];
  61     }
  62 }
  63
  64
  65 // * * * * * * * * * * * * * * * Friend Operators  * * * * * * * * * * * * * //
  66
  67 void Foam::surfaceLocation::write(Ostream& os, const triSurface& s) const
  68 {
  69     if (elementType_ == triPointRef::NONE)
  70     {
  71         os  << "trianglecoords:" << s[index()].tri(s.points());
  72     }
  73     else if (elementType() == triPointRef::EDGE)
  74     {
  75         const edge& e = s.edges()[index()];
  76
  77         os  << "edgecoords:" << e.line(s.localPoints());
  78     }
  79     else
  80     {
  81         os  << "pointcoord:" << s.localPoints()[index()];
  82     }
  83 }
  84
  85
  86 Foam::Istream& Foam::operator>>(Istream& is, surfaceLocation& sl)
  87 {
  88     label elType;
  89     is  >> static_cast<pointIndexHit&>(sl)
  90         >> elType >> sl.triangle_;
  91     sl.elementType_ = triPointRef::proxType(elType);
  92     return is;
  93 }
  94
  95
  96 Foam::Ostream& Foam::operator<<(Ostream& os, const surfaceLocation& sl)
  97 {
  98     return os
  99         << static_cast<const pointIndexHit&>(sl)
 100         << token::SPACE << label(sl.elementType_)
 101         << token::SPACE << sl.triangle_;
 102 }
 103
 104
 105 Foam::Ostream& Foam::operator<<
 106 (
 107     Ostream& os,
 108     const InfoProxy<surfaceLocation>& ip
 109 )
 110 {
 111     const surfaceLocation& sl = ip.t_;
 112
 113     if (sl.elementType() == triPointRef::NONE)
 114     {
 115         os  << "coord:" << sl.rawPoint()
 116             << " inside triangle:" << sl.index()
 117             << " excludeTriangle:" << sl.triangle();
 118     }
 119     else if (sl.elementType() == triPointRef::EDGE)
 120     {
 121         os  << "coord:" << sl.rawPoint()
 122             << " on edge:" << sl.index()
 123             << " excludeTriangle:" << sl.triangle();
 124     }
 125     else
 126     {
 127         os  << "coord:" << sl.rawPoint()
 128             << " on point:" << sl.index()
 129             << " excludeTriangle:" << sl.triangle();
 130     }
 131
 132     if (sl.hit())
 133     {
 134         os  << " (hit)";
 135     }
 136     else
 137     {
 138         os  << " (miss)";
 139     }
 140
 141     return os;
 142 }
 143
 144
 145 // ************************************************************************* //