DualNumber.h

浏览该文件的文档.

/************************************************************************
* @file DualNumber.h
*
* @brief 对偶数类
*
* @details 对偶数类
*
* @author dixu
*
* @version 1.0
*
* @date 2014-9-11
*
* @license 北京华科软科技有限公司
*
*************************************************************************/
 
#ifndef BIMHOMEBASE_DUALNUMBER_H
#define BIMHOMEBASE_DUALNUMBER_H
 
#include <cmath>
 
namespace Base
{
    class DualNumber
    {
    public:
        DualNumber()
        {
        }
 
        DualNumber(double re, double du = 0.0) : m_re(re), m_du(du)
        {
        }
 
        DualNumber operator - () const
        {
            return DualNumber(-m_re, -m_du);
        }
 
        double getRE() const
        {
            return m_re;
        }
 
        double getDU() const
        {
            return m_du;
        }
 
        void setRE(double dRe)
        {
            this->m_re = dRe;
        }
 
        void setDU(double dDu)
        {
            this->m_du = dDu;
        }
 
    protected:
        double m_re = 0.0;      
        double m_du = 0.0;      
    };
 
    /* 对偶数基本运算 */
    inline DualNumber operator + (DualNumber a, DualNumber b)
    {
        return DualNumber(a.getRE() + b.getRE(), a.getDU() + b.getDU());
    }
 
    inline DualNumber operator + (DualNumber a, double b)
    {
        return DualNumber(a.getRE() + b, a.getDU());
    }
 
    inline DualNumber operator + (double a, DualNumber b)
    {
        return DualNumber(a + b.getRE(), b.getDU());
    }
 
    inline DualNumber operator - (DualNumber a, DualNumber b)
    {
        return DualNumber(a.getRE() - b.getRE(), a.getDU() - b.getDU());
    }
 
    inline DualNumber operator - (DualNumber a, double b)
    {
        return DualNumber(a.getRE() - b, a.getDU());
    }
 
    inline DualNumber operator - (double a, DualNumber b)
    {
        return DualNumber(a - b.getRE(), -b.getDU());
    }
 
    inline DualNumber operator * (DualNumber a, DualNumber b)
    {
        return DualNumber(a.getRE() * b.getRE(), a.getRE() * b.getDU() + a.getDU() * b.getRE());
    }
 
    inline DualNumber operator * (double a, DualNumber b)
    {
        return DualNumber(a * b.getRE(), a * b.getDU());
    }
 
    inline DualNumber operator * (DualNumber a, double b)
    {
        return DualNumber(a.getRE() * b, a.getDU() * b);
    }
 
    inline DualNumber operator / (DualNumber a, DualNumber b)
    {
        return DualNumber(a.getRE() / b.getRE(), (a.getDU() * b.getRE() - a.getRE() * b.getDU()) / (b.getRE() * b.getRE()));
    }
 
    inline DualNumber operator / (DualNumber a, double b)
    {
        return DualNumber(a.getRE() / b, a.getDU() / b);
    }
 
    inline DualNumber pow(DualNumber a, double pw)
    {
        return DualNumber(std::pow(a.getRE(), pw), pw * std::pow(a.getRE(), pw - 1.0) * a.getDU());
    }
} //namespace
 
 
#endif