src/OpenFOAM/coordinateSystems/coordinateRotation/EulerCoordinateRotation.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright held by original author
   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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "EulerCoordinateRotation.H"
  28
  29 #include "Switch.H"
  30 #include "mathematicalConstants.H"
  31 #include "dictionary.H"
  32 #include "addToRunTimeSelectionTable.H"
  33
  34 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  35
  36 namespace Foam
  37 {
  38     defineTypeNameAndDebug(EulerCoordinateRotation, 0);
  39     addToRunTimeSelectionTable
  40     (
  41         coordinateRotation,
  42         EulerCoordinateRotation,
  43         dictionary
  44     );
  45 }
  46
  47 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  48
  49 void Foam::EulerCoordinateRotation::calcTransform
  50 (
  51     const scalar phiAngle,
  52     const scalar thetaAngle,
  53     const scalar psiAngle,
  54     const bool inDegrees
  55 )
  56 {
  57     scalar phi = phiAngle;
  58     scalar theta = thetaAngle;
  59     scalar psi = psiAngle;
  60
  61     if (inDegrees)
  62     {
  63         phi   *= mathematicalConstant::pi/180.0;
  64         theta *= mathematicalConstant::pi/180.0;
  65         psi   *= mathematicalConstant::pi/180.0;
  66     }
  67
  68     tensor::operator=
  69     (
  70         tensor
  71         (
  72             cos(phi)*cos(psi) - sin(phi)*sin(psi)*cos(theta),
  73             -sin(phi)*cos(psi)*cos(theta) - cos(phi)*sin(psi),
  74             sin(phi)*sin(theta),
  75
  76             cos(phi)*sin(psi)*cos(theta) + sin(phi)*cos(psi),
  77             cos(phi)*cos(psi)*cos(theta) - sin(phi)*sin(psi),
  78             -cos(phi)*sin(theta),
  79
  80             sin(psi)*sin(theta),
  81             cos(psi)*sin(theta),
  82             cos(theta)
  83         )
  84     );
  85 }
  86
  87
  88 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  89
  90 Foam::EulerCoordinateRotation::EulerCoordinateRotation()
  91 :
  92     coordinateRotation()
  93 {}
  94
  95
  96 Foam::EulerCoordinateRotation::EulerCoordinateRotation
  97 (
  98     const vector& phiThetaPsi,
  99     const bool inDegrees
 100 )
 101 :
 102     coordinateRotation()
 103 {
 104     calcTransform
 105     (
 106         phiThetaPsi.component(vector::X),
 107         phiThetaPsi.component(vector::Y),
 108         phiThetaPsi.component(vector::Z),
 109         inDegrees
 110     );
 111 }
 112
 113
 114 Foam::EulerCoordinateRotation::EulerCoordinateRotation
 115 (
 116     const scalar phiAngle,
 117     const scalar thetaAngle,
 118     const scalar psiAngle,
 119     const bool inDegrees
 120 )
 121 :
 122     coordinateRotation()
 123 {
 124     calcTransform(phiAngle, thetaAngle, psiAngle, inDegrees);
 125 }
 126
 127
 128 Foam::EulerCoordinateRotation::EulerCoordinateRotation
 129 (
 130     const dictionary& dict
 131 )
 132 :
 133     coordinateRotation()
 134 {
 135     vector rotation(dict.lookup("rotation"));
 136
 137     calcTransform
 138     (
 139         rotation.component(vector::X),
 140         rotation.component(vector::Y),
 141         rotation.component(vector::Z),
 142         dict.lookupOrDefault<Switch>("degrees", true)
 143     );
 144 }
 145
 146
 147 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //