src/foam/coordinateSystems/coordinateRotation/axisCoordinateRotation.C

   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 #include "axisCoordinateRotation.H"
  27
  28 #include "Switch.H"
  29 #include "mathematicalConstants.H"
  30 #include "dictionary.H"
  31 #include "addToRunTimeSelectionTable.H"
  32
  33 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  34
  35 namespace Foam
  36 {
  37     defineTypeNameAndDebug(axisCoordinateRotation, 0);
  38     addToRunTimeSelectionTable
  39     (
  40         coordinateRotation,
  41         axisCoordinateRotation,
  42         dictionary
  43     );
  44 }
  45
  46 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  47
  48 void Foam::axisCoordinateRotation::calcTransform
  49 (
  50     const scalar phiAngle,
  51     const scalar thetaAngle,
  52     const scalar psiAngle,
  53     const bool inDegrees
  54 )
  55 {
  56     scalar phi   = phiAngle;
  57     scalar theta = thetaAngle;
  58     scalar psi   = psiAngle;
  59
  60     if (inDegrees)
  61     {
  62         phi *= mathematicalConstant::pi/180.0;
  63         theta *= mathematicalConstant::pi/180.0;
  64         psi *= mathematicalConstant::pi/180.0;
  65     }
  66
  67     tensor::operator=
  68     (
  69         tensor
  70         (
  71             cos(theta)*cos(psi),
  72             sin(phi)*sin(theta)*cos(psi) - cos(phi)*sin(psi),
  73             cos(phi)*sin(theta)*cos(psi) + sin(phi)*sin(psi),
  74
  75             cos(theta)*sin(psi),
  76             sin(phi)*sin(theta)*sin(psi) + cos(phi)*cos(psi),
  77             cos(phi)*sin(theta)*sin(psi) - sin(phi)*cos(psi),
  78
  79             -sin(theta),
  80             sin(phi)*cos(theta),
  81             cos(phi)*cos(theta)
  82         )
  83     );
  84 }
  85
  86
  87 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  88
  89 Foam::axisCoordinateRotation::axisCoordinateRotation
  90 (
  91     const scalar phiAngle,
  92     const scalar thetaAngle,
  93     const scalar psiAngle,
  94     const bool inDegrees
  95 )
  96 :
  97     coordinateRotation()
  98 {
  99     calcTransform(phiAngle, thetaAngle, psiAngle, inDegrees);
 100 }
 101
 102
 103 Foam::axisCoordinateRotation::axisCoordinateRotation
 104 (
 105     const dictionary& dict
 106 )
 107 :
 108     coordinateRotation()
 109 {
 110     scalar phi = readScalar(dict.lookup("phi"));
 111     scalar theta = readScalar(dict.lookup("theta"));
 112     scalar psi = readScalar(dict.lookup("psi"));
 113
 114     calcTransform
 115     (
 116         phi,
 117         theta,
 118         psi,
 119         dict.lookupOrDefault<Switch>("degrees", true)
 120     );
 121 }
 122
 123
 124 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //