src/foam/primitives/Vector/VectorTemplateI.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend 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     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  27
  28 namespace Foam
  29 {
  30
  31 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  32
  33 // Construct null
  34 template <class Cmpt>
  35 inline Vector<Cmpt>::Vector()
  36 {}
  37
  38
  39 // Construct given VectorSpace
  40 template <class Cmpt>
  41 inline Vector<Cmpt>::Vector(const VectorSpace<Vector<Cmpt>, Cmpt, 3>& vs)
  42 :
  43     VectorSpace<Vector<Cmpt>, Cmpt, 3>(vs)
  44 {}
  45
  46
  47 // Construct given three Cmpts
  48 template <class Cmpt>
  49 inline Vector<Cmpt>::Vector(const Cmpt& vx, const Cmpt& vy, const Cmpt& vz)
  50 {
  51     this->v_[X] = vx;
  52     this->v_[Y] = vy;
  53     this->v_[Z] = vz;
  54 }
  55
  56
  57 // Construct from Istream
  58 template <class Cmpt>
  59 inline Vector<Cmpt>::Vector(Istream& is)
  60 :
  61     VectorSpace<Vector<Cmpt>, Cmpt, 3>(is)
  62 {}
  63
  64
  65 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  66
  67 template <class Cmpt>
  68 inline const Cmpt& Vector<Cmpt>::x() const
  69 {
  70     return this->v_[X];
  71 }
  72
  73 template <class Cmpt>
  74 inline const Cmpt& Vector<Cmpt>::y() const
  75 {
  76     return this->v_[Y];
  77 }
  78
  79 template <class Cmpt>
  80 inline const Cmpt& Vector<Cmpt>::z() const
  81 {
  82     return this->v_[Z];
  83 }
  84
  85
  86 template <class Cmpt>
  87 inline Cmpt& Vector<Cmpt>::x()
  88 {
  89     return this->v_[X];
  90 }
  91
  92 template <class Cmpt>
  93 inline Cmpt& Vector<Cmpt>::y()
  94 {
  95     return this->v_[Y];
  96 }
  97
  98 template <class Cmpt>
  99 inline Cmpt& Vector<Cmpt>::z()
 100 {
 101     return this->v_[Z];
 102 }
 103
 104
 105 template <class Cmpt>
 106 inline const Cmpt& Vector<Cmpt>::operator()
 107 (
 108     const direction i
 109 ) const
 110 {
 111     return this->operator[](i);
 112 }
 113
 114
 115 template <class Cmpt>
 116 inline Cmpt& Vector<Cmpt>::operator()
 117 (
 118     const direction i
 119 )
 120 {
 121     return this->operator[](i);
 122 }
 123
 124
 125 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 126
 127 template <class Cmpt>
 128 inline const Vector<Cmpt>& Vector<Cmpt>::centre
 129 (
 130     const Foam::List<Vector<Cmpt> >&
 131 )const
 132 {
 133     return *this;
 134 }
 135
 136
 137 // * * * * * * * * * * * * * * * Global Operators  * * * * * * * * * * * * * //
 138
 139 template <class Cmpt>
 140 inline typename innerProduct<Vector<Cmpt>, Vector<Cmpt> >::type
 141 operator&(const Vector<Cmpt>& v1, const Vector<Cmpt>& v2)
 142 {
 143     return Cmpt(v1.x()*v2.x() + v1.y()*v2.y() + v1.z()*v2.z());
 144 }
 145
 146
 147 template <class Cmpt>
 148 inline Vector<Cmpt> operator^(const Vector<Cmpt>& v1, const Vector<Cmpt>& v2)
 149 {
 150     return Vector<Cmpt>
 151     (
 152         (v1.y()*v2.z() - v1.z()*v2.y()),
 153         (v1.z()*v2.x() - v1.x()*v2.z()),
 154         (v1.x()*v2.y() - v1.y()*v2.x())
 155     );
 156 }
 157
 158
 159 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 160
 161 } // End namespace Foam
 162
 163 // ************************************************************************* //