-/*******************************************************************************
+/*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Header: FGRotation.h
Author: Jon Berndt
--------------------------------------------------------------------------------
12/02/98 JSB Created
-********************************************************************************
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
COMMENTS, REFERENCES, and NOTES
-********************************************************************************
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[1] Cooke, Zyda, Pratt, and McGhee, "NPSNET: Flight Simulation Dynamic Modeling
Using Quaternions", Presence, Vol. 1, No. 4, pp. 404-420 Naval Postgraduate
School, January 1994
The order of rotations used in this class corresponds to a 3-2-1 sequence,
or Y-P-R, or Z-Y-X, if you prefer.
-********************************************************************************
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
SENTRY
-*******************************************************************************/
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
#ifndef FGROTATION_H
#define FGROTATION_H
-/*******************************************************************************
+/*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
INCLUDES
-*******************************************************************************/
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
#ifdef FGFS
# include <simgear/compiler.h>
-# ifdef FG_HAVE_STD_INCLUDES
-# include <cmath>
-# else
+# include <math.h>
+#else
+# if defined (sgi) && !defined(__GNUC__)
# include <math.h>
+# else
+# include <cmath>
# endif
-#else
-# include <cmath>
#endif
#include "FGModel.h"
-#include "FGMatrix.h"
+#include "FGMatrix33.h"
+#include "FGColumnVector3.h"
+#include "FGColumnVector4.h"
-/*******************************************************************************
-CLASS DECLARATION
-*******************************************************************************/
+#define ID_ROTATION "$Id$"
-#pragma warn -8026
+/*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+CLASS DECLARATION
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*/
class FGRotation : public FGModel
{
- FGColumnVector vPQR;
- FGColumnVector vPQRdot;
- FGColumnVector vMoments;
- FGColumnVector vEuler;
-
- float Ixx, Iyy, Izz, Ixz;
- float dt;
-
- void GetState(void);
-
public:
FGRotation(FGFDMExec*);
- ~FGRotation(void);
+ ~FGRotation();
bool Run(void);
- inline FGColumnVector GetPQR(void) {return vPQR;}
- inline FGColumnVector GetPQRdot(void) {return vPQRdot;}
- inline FGColumnVector GetEuler(void) {return vEuler;}
- inline void SetPQR(FGColumnVector tt) {vPQR = tt;}
- inline void SetEuler(FGColumnVector tt) {vEuler = tt;}
- inline float Getphi(void) {return vEuler(1);}
- inline float Gettht(void) {return vEuler(2);}
- inline float Getpsi(void) {return vEuler(3);}
+ inline FGColumnVector3& GetPQR(void) {return vPQR;}
+ inline double GetPQR(int axis) {return vPQR(axis);}
+ inline FGColumnVector3& GetAeroPQR(void) {return vAeroPQR;}
+ inline double GetAeroPQR(int axis) {return vAeroPQR(axis);}
+ inline FGColumnVector3& GetPQRdot(void) {return vPQRdot;}
+ inline double GetPQRdot(int idx) {return vPQRdot(idx);}
+ inline FGColumnVector3& GetEuler(void) {return vEuler;}
+ inline double GetEuler(int axis) {return vEuler(axis);}
+ inline FGColumnVector3& GetEulerRates(void) { return vEulerRates; }
+ inline double GetEulerRates(int axis) { return vEulerRates(axis); }
+ inline void SetPQR(FGColumnVector3 tt) {vPQR = tt;}
+ inline void SetPQR(double p, double q, double r) {vPQR(eP)=p;
+ vPQR(eQ)=q;
+ vPQR(eR)=r;}
+ inline void SetEuler(FGColumnVector3 tt) {vEuler = tt;}
+
+ inline double Getphi(void) {return vEuler(1);}
+ inline double Gettht(void) {return vEuler(2);}
+ inline double Getpsi(void) {return vEuler(3);}
+
+ inline double GetCosphi(void) {return cPhi;}
+ inline double GetCostht(void) {return cTht;}
+ inline double GetCospsi(void) {return cPsi;}
+
+ inline double GetSinphi(void) {return sPhi;}
+ inline double GetSintht(void) {return sTht;}
+ inline double GetSinpsi(void) {return sPsi;}
+
+private:
+ FGColumnVector3 vPQR;
+ FGColumnVector3 vAeroPQR;
+ FGColumnVector3 vPQRdot;
+ FGColumnVector3 vMoments;
+ FGColumnVector3 vEuler;
+ FGColumnVector3 vEulerRates;
+ FGColumnVector3 vlastPQRdot;
+
+ double cTht,sTht;
+ double cPhi,sPhi;
+ double cPsi,sPsi;
+
+ double Ixx, Iyy, Izz, Ixz;
+ double dt;
+
+ void GetState(void);
+
+ void Debug(int from);
};
-#pragma warn .8026
-/******************************************************************************/
+//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#endif
+