// Distributed under the Boost Software License, Version 1.0. (See
// accompanying file LICENSE_1_0.txt or copy at
// http://www.boost.org/LICENSE_1_0.txt)
// (C) Copyright 2007 Anthony Williams
#include "fixed.hpp"
#include <assert.h>
#include <iostream>
#include <iomanip>

bool valuesAreWithinEpsilon(double t,double u,double epsilon)
{
    double const maxValue=std::max<double>(t,u);
    double const minValue=std::min<double>(t,u);
    epsilon=abs(epsilon);
    bool const res=(maxValue>=0 && (maxValue-epsilon)<=minValue)
        || (minValue<=0 && (minValue+epsilon)>=maxValue);

    if(!res)
    {
        std::cout<<"assertion failed (a==b), a="<<std::setprecision(20)<<std::dec<<t<<", b="<<u<<", e="<<epsilon<<std::endl;
    }
    return res;
}

#define TEST_ASSERT_EQUALS_FP(a,b,e) (assert(valuesAreWithinEpsilon((a),(b),(e))))
#define TEST_ASSERT_EQUALS(a,b) (assert((a)==(b)))

int main()
{
    TEST_ASSERT_EQUALS_FP(sin(fixed(0)),sin(0.0),1e-8);
    TEST_ASSERT_EQUALS_FP(sin(fixed_half_pi),sin(fixed_half_pi.as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(sin(fixed_quarter_pi),sin(fixed_half_pi.as_double()/2),1e-8);
    TEST_ASSERT_EQUALS_FP(sin(fixed_pi),sin(fixed_pi.as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(sin(3*fixed_quarter_pi),sin(3*fixed_half_pi.as_double()/2),1e-8);
    TEST_ASSERT_EQUALS_FP(sin(-fixed_half_pi),sin(-fixed_half_pi.as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(sin(-fixed_quarter_pi),sin(-fixed_half_pi.as_double()/2),1e-8);
    TEST_ASSERT_EQUALS_FP(sin(-fixed_pi),sin(-fixed_pi.as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(sin(-3*fixed_quarter_pi),sin(-3*fixed_half_pi.as_double()/2),1e-8);
    TEST_ASSERT_EQUALS_FP(cos(fixed(0)),cos(0.0),1e-8);
    TEST_ASSERT_EQUALS_FP(cos(fixed_half_pi),cos(fixed_half_pi.as_double()),1e-7);
    TEST_ASSERT_EQUALS_FP(cos(fixed_quarter_pi),cos(fixed_half_pi.as_double()/2),1e-8);
    TEST_ASSERT_EQUALS_FP(cos(-fixed_half_pi),cos(-fixed_half_pi.as_double()),1e-7);
    TEST_ASSERT_EQUALS_FP(cos(-fixed_quarter_pi),cos(-fixed_half_pi.as_double()/2),1e-8);
    TEST_ASSERT_EQUALS_FP(cos(-fixed_pi),cos(-fixed_pi.as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(cos(fixed_pi),cos(fixed_pi.as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(cos(3*fixed_quarter_pi),cos(3*fixed_half_pi.as_double()/2),1e-8);
    TEST_ASSERT_EQUALS_FP(cos(-3*fixed_quarter_pi),cos(-3*fixed_half_pi.as_double()/2),1e-8);
    TEST_ASSERT_EQUALS_FP(tan(fixed(0)),tan(0.0),1e-8);
    TEST_ASSERT_EQUALS_FP(tan(fixed_quarter_pi),tan(fixed_half_pi.as_double()/2),1e-7);
    TEST_ASSERT_EQUALS_FP(tan(-fixed_quarter_pi),tan(-fixed_half_pi.as_double()/2),1e-7);
    TEST_ASSERT_EQUALS_FP(tan(-fixed_pi),tan(-fixed_pi.as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(tan(fixed_pi),tan(fixed_pi.as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(tan(3*fixed_quarter_pi),tan(3*fixed_half_pi.as_double()/2),1e-7);
    TEST_ASSERT_EQUALS_FP(tan(-3*fixed_quarter_pi),tan(-3*fixed_half_pi.as_double()/2),1e-7);
    TEST_ASSERT_EQUALS_FP(fixed(1.5*2000000)/2000000.0,1.5,1e-8);
    TEST_ASSERT_EQUALS_FP(90*fixed_half_pi,90*fixed_half_pi.as_double(),1e-8);
    TEST_ASSERT_EQUALS_FP(45*fixed_half_pi,45*fixed_half_pi.as_double(),1e-8);
    TEST_ASSERT_EQUALS_FP(180*fixed_half_pi,180*fixed_half_pi.as_double(),1e-8);
    TEST_ASSERT_EQUALS_FP(90*fixed_half_pi/90.0,90*fixed_half_pi.as_double()/90.0,1e-8);
    TEST_ASSERT_EQUALS_FP(45*fixed_half_pi/90.0,45*fixed_half_pi.as_double()/90.0,1e-8);
    TEST_ASSERT_EQUALS_FP(180*fixed_half_pi/90.0,180*fixed_half_pi.as_double()/90.0,1e-8);
    TEST_ASSERT_EQUALS_FP(fixed(-180.3).floor(),-181.0,1e-8);
    TEST_ASSERT_EQUALS_FP(fixed(-180.0).floor(),-180.0,1e-8);
    TEST_ASSERT_EQUALS_FP(fixed(-180.3).ceil(),-180.0,1e-8);
    TEST_ASSERT_EQUALS_FP(fixed(10000)/fixed(0.9),10000.0/fixed(0.9).as_double(),1e-7);
    TEST_ASSERT_EQUALS_FP(fixed(0.9)/fixed(10000),fixed(0.9).as_double()/10000.0,1e-7);
    TEST_ASSERT_EQUALS_FP((fixed(10000)/fixed(0.9))*fixed(0.9),10000.0,1e-7);
    TEST_ASSERT_EQUALS_FP(fixed(10000)/fixed(9),10000.0/9.0,1e-8);
    TEST_ASSERT_EQUALS_FP(fixed(0.09),0.09,1e-8);
    TEST_ASSERT_EQUALS_FP(std::arg(std::polar<fixed>(1,0)),0.0,1e-8);
    TEST_ASSERT_EQUALS_FP(std::arg(std::polar<fixed>(1,fixed_quarter_pi)),fixed_quarter_pi,1e-8);
    TEST_ASSERT_EQUALS_FP(std::arg(std::polar<fixed>(1,fixed_half_pi)),fixed_half_pi,1e-7);
    TEST_ASSERT_EQUALS_FP(std::arg(std::polar<fixed>(1,3*fixed_quarter_pi)),3*fixed_quarter_pi,1e-8);
    TEST_ASSERT_EQUALS_FP(std::arg(std::polar<fixed>(3,0)),0.0,1e-8);
    TEST_ASSERT_EQUALS_FP(std::arg(std::polar<fixed>(3,fixed_quarter_pi)),fixed_quarter_pi,1e-8);
    TEST_ASSERT_EQUALS_FP(std::arg(std::polar<fixed>(3,fixed_half_pi)),fixed_half_pi,1e-7);
    TEST_ASSERT_EQUALS_FP(std::arg(std::polar<fixed>(3,3*fixed_quarter_pi)),3*fixed_quarter_pi,1e-8);
    TEST_ASSERT_EQUALS_FP(exp(fixed(0)),1.0,1e-7);
    TEST_ASSERT_EQUALS_FP(exp(fixed(1)),exp(1.0),1e-7);
    TEST_ASSERT_EQUALS_FP(exp(fixed(10)),exp(10.0),1e-4);
    TEST_ASSERT_EQUALS_FP(exp(fixed(15.7)),exp(fixed(15.7).as_double()),2e-2);
    TEST_ASSERT_EQUALS_FP(exp(fixed(0.00013)),exp(0.00013),1e-7);
    TEST_ASSERT_EQUALS_FP(exp(fixed(fixed::internal(),1)),exp(fixed(fixed::internal(),1).as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(exp(fixed(fixed::internal(),0x7fffffff)),exp(fixed(fixed::internal(),0x7fffffff).as_double()),1e-5);
    TEST_ASSERT_EQUALS_FP(exp(fixed(-1)),exp(-1.0),1e-7);
    TEST_ASSERT_EQUALS_FP(exp(fixed(-10)),exp(-10.0),1e-4);
    TEST_ASSERT_EQUALS_FP(exp(fixed(-15.7)),exp(fixed(-15.7).as_double()),2e-2);
    TEST_ASSERT_EQUALS_FP(exp(fixed(-0.00013)),exp(-0.00013),1e-7);
    TEST_ASSERT_EQUALS_FP(exp(fixed(fixed::internal(),-1)),exp(fixed(fixed::internal(),-1).as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(exp(fixed(fixed::internal(),-0x7fffffff)),exp(fixed(fixed::internal(),-0x7fffffff).as_double()),1e-5);
    TEST_ASSERT_EQUALS_FP(log(fixed_max),log(fixed_max.as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(log(fixed(0.25)),log(0.25),1e-8);
    TEST_ASSERT_EQUALS_FP(log(fixed(0.5)),log(0.5),1e-8);
    TEST_ASSERT_EQUALS_FP(log(fixed(1)),log(1.0),1e-8);
    TEST_ASSERT_EQUALS_FP(log(fixed(2)),log(2.0),1e-8);
    TEST_ASSERT_EQUALS_FP(log(fixed(4)),log(4.0),1e-8);
    TEST_ASSERT_EQUALS_FP(log(fixed(5)),log(5.0),1e-8);
    TEST_ASSERT_EQUALS_FP(log(fixed(10)),log(10.0),1e-8);
    TEST_ASSERT_EQUALS_FP(log(fixed(15.7)),log(fixed(15.7).as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(log(fixed(0.00013)),log(0.00013),1e-4);
    TEST_ASSERT_EQUALS_FP(log(fixed(fixed::internal(),1)),log(fixed(fixed::internal(),1).as_double()),1e-8);
    TEST_ASSERT_EQUALS_FP(log(fixed(fixed::internal(),0x7fffffff)),log(fixed(fixed::internal(),0x7fffffff).as_double()),1e-8);
    TEST_ASSERT_EQUALS(fixed(5)%fixed(4),1);
    TEST_ASSERT_EQUALS(fixed(6)%fixed(4),2);
    TEST_ASSERT_EQUALS(fixed(6.5)%fixed(4),2.5);
}