src/OpenFOAM/coordinateSystems/coordinateRotation/axisCoordinateRotation.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 "axisCoordinateRotation.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(axisCoordinateRotation, 0);
  39     addToRunTimeSelectionTable
  40     (
  41         coordinateRotation,
  42         axisCoordinateRotation,
  43         dictionary
  44     );
  45 }
  46
  47 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  48
  49 void Foam::axisCoordinateRotation::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(theta)*cos(psi),
  73             sin(phi)*sin(theta)*cos(psi) - cos(phi)*sin(psi),
  74             cos(phi)*sin(theta)*cos(psi) + sin(phi)*sin(psi),
  75
  76             cos(theta)*sin(psi),
  77             sin(phi)*sin(theta)*sin(psi) + cos(phi)*cos(psi),
  78             cos(phi)*sin(theta)*sin(psi) - sin(phi)*cos(psi),
  79
  80             -sin(theta),
  81             sin(phi)*cos(theta),
  82             cos(phi)*cos(theta)
  83         )
  84     );
  85 }
  86
  87
  88 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  89
  90 Foam::axisCoordinateRotation::axisCoordinateRotation
  91 (
  92     const scalar phiAngle,
  93     const scalar thetaAngle,
  94     const scalar psiAngle,
  95     const bool inDegrees
  96 )
  97 :
  98     coordinateRotation()
  99 {
 100     calcTransform(phiAngle, thetaAngle, psiAngle, inDegrees);
 101 }
 102
 103
 104 Foam::axisCoordinateRotation::axisCoordinateRotation
 105 (
 106     const dictionary& dict
 107 )
 108 :
 109     coordinateRotation()
 110 {
 111     scalar phi = readScalar(dict.lookup("phi"));
 112     scalar theta = readScalar(dict.lookup("theta"));
 113     scalar psi = readScalar(dict.lookup("psi"));
 114
 115     bool inDegrees = true;
 116     if (dict.found("degrees"))
 117     {
 118         inDegrees = Switch(dict.lookup("degrees"));
 119     }
 120
 121     calcTransform
 122     (
 123         phi,
 124         theta,
 125         psi,
 126         dict.lookupOrDefault<Switch>("degrees", true)
 127     );
 128 }
 129
 130
 131 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //