#ifndef  body_INC
#define  body_INC
#include <Eigen/Dense>
#include <iostream>
#include "types.h"
#define R_HEIGHT 2.5e-3            /* measurement noise of height measurement */
#define IMU_NOISE 800e-6            /* IMU process noise */
#define IMU_RANDOMWALK 50e-6      /* Bias process noise */
//#define DOCLAMP
#define MAXHEIGHT -1.5 /* Notion of max and min is reversed due to z pointing down */
#define MINHEIGHT -0.3            /*  */
#define FULLSTATE /* Uncomment to estimate quaternion */


using Eigen::Matrix;
using Eigen::Vector3d;
using Eigen::Quaterniond;

/*
 * =====================================================================================
 *        Class:  Body
 *  Description:  State representation of body.
 * =====================================================================================
 */
class Body
{
    public:
        /* ====================  LIFECYCLE     ======================================= */
        Body (){ gravity_world << 0.,0.,9.80655; };                             /* constructor */

        /* ====================  ACCESSORS     ======================================= */
#ifdef FULLSTATE
        Matrix<double,13,1> asVector();
        void asVector(const Matrix<double,13,1> &m);
        void dx( const Matrix<double,13,1> &del);
        Matrix<double,13,13> F(const Vector3d &ang, const Quaterniond &q, double dt);
        Matrix<double,1,13> H();
        Matrix<double,13,13> Q(double dt);
        Matrix<double,1,1> S(const Matrix<double,13,13> &Pxx);
#else
        Matrix<double,9,9> F(const Vector3d &ang, const Quaterniond &q, double dt);
        Matrix<double,9,1> asVector();
        void asVector(const Matrix<double,9,1> &m);
        void dx( const Matrix<double,9,1> &del);
        Matrix<double,1,9> H();
        Matrix<double,9,9> Q(double dt);
        Matrix<double,1,1> S(const Matrix<double,9,9> &Pxx);
#endif
        Vector3d ned();
        Vector3d vel();

        /* ====================  MUTATORS      ======================================= */
        void accelerometer_bias( const Vector3d &b);
        void clamp();
        void pos( const UTM &utm);
        void vel(const Matrix<double,3,1> &v);

        /* ====================  OPERATORS     ======================================= */
        Matrix<double,1,1> h();
        void motionModel ( const Vector3d &acc, const Vector3d &ang, 
                const Quaterniond &q, const double dt);
        Matrix<double,1,1> R();
        Matrix<double,3,3> skewSymmetric(const Vector3d &x);
        Matrix<double,4,4> omega(const Vector3d &x);
        void unicsv();

    protected:
        /* ====================  METHODS       ======================================= */

        /* ====================  DATA MEMBERS  ======================================= */

    private:
        /* ====================  METHODS       ======================================= */

        /* ====================  DATA MEMBERS  ======================================= */
#ifdef FULLSTATE
        Matrix<double,13,1> X;
#else
        Matrix<double,9,1> X;
#endif
        char utm_c;
        int utm_i;
        Vector3d gravity_world;

}; /* -----  end of class Body  ----- */


#endif   /* ----- #ifndef body_INC  ----- */