// Copyright (c) 1995-2009 Nokia Corporation and/or its subsidiary(-ies).

// All rights reserved.

// This component and the accompanying materials are made available

// under the terms of the License "Eclipse Public License v1.0"

// which accompanies this distribution, and is available

// at the URL "http://www.eclipse.org/legal/epl-v10.html".

//

// Initial Contributors:

// Nokia Corporation - initial contribution.

//

// Contributors:

//

// Description:

// e32test\math\t_r32.cpp

// T_R32.CPP - Test routines for TReal32

// Also note that these tests do not generally include testing of special values.  This is done

// in T_R96 tests and conversions are tested thoroughly, so explicit tests are unnecessary here.

// Overview:

// Test functionality of operations on 32bit real numbers.

// API Information:

// TReal32. 

// Details:

// - Test the conversion from TReal to TReal32 is as expected.

// - Check addition, subtraction and multiplication of 32-bit floating point

// numbers are as expected.

// - Check division of 32-bit floating-point numbers and verify that it is

// panicked when divided by zero.

// - Test arithmetic exceptions are raised for

// - overflow error during addition, subtraction.

// - overflow, underflow errors during multiplication.

// - overflow, underflow, divided by zero errors during division.

// - overflow, underflow, invalid operation errors during conversion

// from double to float.

// - Check unary operator, equalities and inequalities operators, pre/post 

// increment, decrement operators with TReal32 are as expected.

// Platforms/Drives/Compatibility:

// All 

// Assumptions/Requirement/Pre-requisites:

// Failures and causes:

// Base Port information:

// 

//

#include "t_math.h"

#include "t_real32.h"

#if defined(__VC32__)

// Some symbols generated by the VC++ compiler for floating point stuff.

extern "C" {

GLDEF_D TInt _adj_fdiv_m32;

}

#endif

// Data for tests from T_R32DTA.cpp 

GLREF_D TReal32 addInput[];

GLREF_D TReal32 subInput[];

GLREF_D TReal32 multInput[];

GLREF_D TReal32 divInput[];

GLREF_D TReal32 unaryInput[];

GLREF_D TReal32 incDecInput[];

GLREF_D TInt sizeAdd;

GLREF_D TInt sizeSub;

GLREF_D TInt sizeMult;

GLREF_D TInt sizeDiv;

GLREF_D TInt sizeUnary;

GLREF_D TInt sizeIncDec;

#if defined (__WINS__) || defined (__X86__)

// Functions from EMGCC32.CPP

GLREF_C TReal32 __addsf3(TReal32 a1,TReal32 a2);

GLREF_C TReal32 __subsf3(TReal32 a1,TReal32 a2);

GLREF_C TReal32 __mulsf3(TReal32 a1,TReal32 a2);

GLREF_C TReal32 __divsf3(TReal32 a1,TReal32 a2);

GLREF_C TReal32 __truncdfsf2(TReal64 a1);

#endif

GLDEF_D TReal32 NaNTReal32;

GLDEF_D TReal32 posInfTReal32;

GLDEF_D TReal32 negInfTReal32;

GLDEF_D TReal NaNTReal;

GLDEF_D TReal posInfTReal;

GLDEF_D TReal negInfTReal;

enum TOrder

	{

	ELessThan,

	EEqual,

	EGreaterThan

	};

LOCAL_D RTest test(_L("T_R32"));

LOCAL_C void initSpecialValues()

//

// Initialise special values

//

	{

	SReal32 *p32=(SReal32*)&NaNTReal32;

	p32->sign=0;

	p32->exp=KTReal32SpecialExponent;

	p32->man=0x7fffff;

	p32=(SReal32*)&posInfTReal32;

	p32->sign=0;

	p32->exp=KTReal32SpecialExponent;

	p32->man=0;

	p32=(SReal32*)&negInfTReal32;

	p32->sign=1;

	p32->exp=KTReal32SpecialExponent;

	p32->man=0;

	SReal64 *p64=(SReal64*)&NaNTReal;

	p64->sign=0;

	p64->exp=KTReal64SpecialExponent;

	p64->lsm=0xffffffffu;

	p64->msm=0xfffff;

	p64=(SReal64*)&posInfTReal;

	p64->sign=0;

	p64->exp=KTReal64SpecialExponent;

	p64->lsm=0;

	p64->msm=0;

	p64=(SReal64*)&negInfTReal;

	p64->sign=1;

	p64->exp=KTReal64SpecialExponent;

	p64->lsm=0;

	p64->msm=0;

	}

LOCAL_C void testConvert()

//

//	Conversion tests

//

	{

	TRealX f;

	TReal input[]=

		{

		KMaxTReal32inTReal,KMinTReal32inTReal,-KMaxTReal32inTReal,-KMinTReal32inTReal,

		KMaxTReal32inTReal,KMinTReal32inTReal,-KMaxTReal32inTReal,-KMinTReal32inTReal,

		3.4027E+38,1.1755E-38,-3.4027E+38,-1.1755E-38,

		0.0,64.5,-64.5,1.54E+18,-1.54E+18,4.72E-22,-4.72E-22,

		posInfTReal,negInfTReal,KNegZeroTReal,

		1.4E-45,-1.4E-45,2E-41,-2E-41,1E-38,-1E-38		

		};

	TReal32 expect[]=

		{

		KMaxTReal32,KMinTReal32,-KMaxTReal32,-KMinTReal32,

		KMaxTReal32,KMinTReal32,-KMaxTReal32,-KMinTReal32,

		3.4027E+38f,1.17550E-38f,-3.40270E+38f,-1.17550E-38f,

		0.0f,64.5f,-64.5f,1.54E+18f,-1.54E+18f,4.72E-22f,-4.72E-22f,

		posInfTReal32,negInfTReal32,KNegZeroTReal32,

		1.4E-45f,-1.4E-45f,2E-41f,-2E-41f,1E-38f,-1E-38f				

		};

	TInt size=sizeof(input)/sizeof(TReal);

	for (TInt ii=0; ii<size; ii++)

		{

		f=TRealX(expect[ii]);

		test(f==TRealX(TReal32(input[ii])));

		}

	// NaN

//	TReal a=NaNTReal;

	TReal32 b=NaNTReal32;

	f=TRealX(b);

//	test(f!=TRealX(TReal32(a)));

	test(f.IsNaN());

	// See EON Software Defects Bug Report no. HA-287

	// There is a bug in MSDev compiler which means comparing TReal32's directly

	// does not always work, hence...

	/*

	test(BitTest(TReal32(3.40270E+38),3.40270E+38f));	// this works

	// (BitTest() checks for all 32 bits being identical

	TReal32 a=TReal32(3.40270E+38);

	TReal32 b=3.40270E+38f;

	TReal64 c=3.40270E+38;

	TReal32 d=TReal32(c);

	test(a==b);									// this works

	test(d==b);									// this works

	test(TRealX(TReal32(c))==TRealX(b));		// this works

	test(TReal64(TReal32(c))==TReal64(b));		// this panics

	test(TReal32(c)==b);						// this panics					  

 	test(TReal32(3.40270E+38)==3.40270E+38f);	// this panics

	// As expected, all these work fine under ARM.

	*/

	}

LOCAL_C void testAdd()

//

//	Addition tests

//

	{

	TReal32 f,g,h;

	TRealX ff,gg,hh;

	for (TInt ii=0; ii<sizeAdd-1; ii++)

		{										  

		f=addInput[ii];

		g=addInput[ii+1];

		ff=TRealX(f);

		gg=TRealX(g);

		// Test commute

		test(f+g == g+f);

		// Test PC real addition using fp-hardware same as TRealX addition

		test(TRealX(f+g)==TRealX(TReal32(ff+gg)));

		test(TRealX(g+f)==TRealX(TReal32(ff+gg)));

		// Test hex-encoded constants for TReal32s generated on PC using fp-hardware same as 

		// TRealX addition

		test(TRealX(*(TReal32*)&addArray[ii])==TRealX(f+g));

		test(TRealX(*(TReal32*)&addArray[ii])==TRealX(g+f));

		// similarly to tests above ...

		h=g;

		hh=gg;

		hh+=ff;

		test(TRealX(h+=f)==TRealX(TReal32(hh)));

		test(TRealX(h)==TRealX(TReal32(hh)));

		test(TRealX(*(TReal32*)&addArray[ii])==TRealX(h));

		//

		h=f;

		hh=ff;

		hh+=gg;

		test(TRealX(h+=g)==TRealX(TReal32(hh)));

		test(h==TReal32(hh));

		test(TRealX(*(TReal32*)&addArray[ii])==TRealX(h));

		}

	}

LOCAL_C void testSubt()

//

// Subtraction tests

//

	{

	TReal32 f,g,h;

	TRealX ff,gg,hh;

	for (TInt ii=0; ii<sizeSub-1; ii++)

		{

		f=subInput[ii];

		g=subInput[ii+1];

		ff=TRealX(f);

		gg=TRealX(g);

		//

// This test fails on GCC (with -O1 switch). The reason is that

// comparing two intermediate floats is unpredictable.

// See http://www.parashift.com/c++-faq-lite/newbie.html#faq-29.18

#ifndef __GCC32__	

		test(f-g == -(g-f));

#endif

		//

		test(TRealX(f-g)==TRealX(TReal32(ff-gg)));

		test(TRealX(g-f)==TRealX(TReal32(gg-ff)));

		test(TRealX(*(TReal32*)&subArray[ii])==TRealX(f-g));

		test(TRealX(*(TReal32*)&subArray[ii])==TRealX(-(g-f)));

		//

		h=g;

		hh=gg;

		hh-=ff;

		test(TRealX(h-=f)==TRealX(TReal32(hh)));

		test(TRealX(h)==TRealX(TReal32(hh)));

		test(TRealX(*(TReal32*)&subArray[ii])==TRealX(-h));

		//

		h=f;

		hh=ff;

		hh-=gg;

		test(TRealX(h-=g)==TRealX(TReal32(hh)));

		test(TRealX(h)==TRealX(TReal32(hh)));

		test(TRealX(*(TReal32*)&subArray[ii])==TRealX(h));

		}

	}

LOCAL_C void testMult()

//

//	Multiplication test

//

	{

	TReal32 f,g,h;

	TRealX ff,gg,hh;

	for (TInt ii=0; ii<sizeMult-1; ii++)

		{

		f=multInput[ii];

		g=multInput[ii+1];

		ff=TRealX(f);

		gg=TRealX(g);

		//

		test(f*g == g*f);

		//

		test(TRealX(f*g)==TRealX(TReal32(ff*gg)));

		test(TRealX(g*f)==TRealX(TReal32(gg*ff)));

		test(TRealX(*(TReal32*)&multArray[ii])==TRealX(f*g));

		test(TRealX(*(TReal32*)&multArray[ii])==TRealX(g*f));

		//

		h=f;		

		hh=ff;

		hh*=gg;

		test(TRealX(h*=g)==TRealX(TReal32(hh)));

		test(TRealX(h)==TRealX(TReal32(hh)));

		test(TRealX(*(TReal32*)&multArray[ii])==TRealX(h));

		//

		h=g;

		hh=gg;

		hh*=ff;

		test(TRealX(h*=f)==TRealX(TReal32(hh)));

		test(TRealX(h)==TRealX(TReal32(hh)));

		test(TRealX(*(TReal32*)&multArray[ii])==TRealX(h));

		}

	}

LOCAL_C void testDiv()

//

//	Division test

//

	{

	TReal32 f,g,h;

	TRealX ff,gg,hh;

	TInt count=0;

	// Panic: Divide by Zero

	// f=1.0;

	// g=0.0;

	// f/=g;

	for (TInt ii=0; ii<sizeDiv-1; ii++)

		{

		f=divInput[ii];

		g=divInput[ii+1];

		ff=TRealX(f);

		gg=TRealX(g);

		if (g!=0.0)

			{

			test(TRealX(f/g)==TRealX(TReal32(ff/gg)));

			test(TRealX(*(TReal32*)&divArray[count])==TRealX(f/g));

			//

			h=f;

			hh=ff;

			hh/=gg;

			test(TRealX(h/=g)==TRealX(TReal32(hh)));

			test(TRealX(h)==TRealX(TReal32(hh)));

			test(TRealX(*(TReal32*)&divArray[count])==TRealX(h));

			++count;

			}

		if (f!=0.0)

			{

			test(TRealX(g/f)==TRealX(TReal32(gg/ff)));

			h=g;

			hh=gg;

			hh/=ff;

			test(TRealX(h/=f)==TRealX(TReal32(hh)));

			test(h==TReal32(hh));

			}

		};

	//Additional test

	f=3.9999f;

	g=KMinTReal32;

	ff=TRealX(f);

	gg=TRealX(g);

	test(TRealX(f/g)==TRealX(TReal32(ff/gg)));

	h=f;

	hh=ff;

	hh/=gg;

	test(TRealX(h/=g)==TRealX(TReal32(hh)));

	test(TRealX(h)==TRealX(TReal32(hh)));

	}

#if defined (__WINS__) || defined (__X86__)

LOCAL_C	 void testArithmeticExceptionRaising()

//

// Test that UP_GCC.CPP raise exceptions correctly by calling functions from EMGCC32.CPP which

// are copies of those in UP_GCC.CPP.  To be used in debugger only.

// Added by AnnW, December 1996

//

	{

	TReal32 f,g,h;

	// Addition - possible errors are overflow, argument or none

	// NB no underflow

	f=NaNTReal32;

	h=__addsf3(f,f);	// argument

	f=KMaxTReal32;

	h=__addsf3(f,f);	// overflow

	f=1.0f;

	g=2.0f;

	h=__addsf3(f,g);	// none

	test(h==3.0f);

	// Subtraction - possible errors are overflow, argument or none

	// NB no underflow

	f=NaNTReal32;

	h=__subsf3(f,f);	// argument

	f=KMaxTReal32;

	g=-KMaxTReal32;

	h=__subsf3(f,g);	// overflow

	f=1.0f;

	g=2.0f;

	h=__subsf3(f,g);	// none

	test(h==-1.0f);

	// Multiplication - possible errors are argument, overflow, underflow or none

	f=NaNTReal32;

	h=__mulsf3(f,f);	// argument

	f=KMaxTReal32;

	g=2.0f;

	h=__mulsf3(f,g);	// overflow

	f=minDenormalTReal32;

	g=0.1f;

	h=__mulsf3(f,g);	// underflow

	f=1.0f;

	g=2.0f;

	h=__mulsf3(f,g);	// none

	test(h==2.0f);

	// Division - possible errors are overflow, underflow, divide by zero, argument or none

	f=KMaxTReal32;

	g=0.5f;

	h=__divsf3(f,g);	// overflow

	f=minDenormalTReal32;

	g=10.0f;

	h=__divsf3(f,g);	// underflow

	f=4.0f;

	g=0.0f;

	h=__divsf3(f,g);	// divide by zero

	f=0.0f;

	g=0.0f;

	h=__divsf3(f,g);	// argument

	f=1.0f;

	g=2.0f;

	h=__divsf3(f,g);	// none

	test(h==0.5f);

	// Converting double to float - possible errors are overflow, underflow, invalid operation or none

	TReal64 d;

	d=1.0E+50;

	f=__truncdfsf2(d);	// overflow

	d=1.0E-50;

	f=__truncdfsf2(d);	// underflow

	d=KNaNTReal64;

	f=__truncdfsf2(d);	// invalid operation

	d=4.0;	

	f=__truncdfsf2(d);	// none

	}

#endif

LOCAL_C void testUnary()

//

//	Unary operator tests

//

	{

	TReal32 f;

	TRealX g;

	for (TInt ii=0; ii<sizeUnary-1; ii++)

		{

		f=unaryInput[ii];

		g=TRealX(f);

		test(TRealX(-f)==TRealX(TReal32(-g)));

		test(TRealX(-f)==TRealX(0.0f-f));

		test(TRealX(+f)==TRealX(TReal32(g)));

		test(TRealX(+f)==TRealX(0.0f+f));

		test(TRealX(*(TReal32*)&unaryArray[ii])==TRealX(-f));

		}

	}

LOCAL_C void testEqualities(const TReal& aA, TOrder aOrder, const TReal& aB)

//

//	Test equality/inequality functions on aA and aB

//	aOrder specifies the operand's relative sizes

//

	{

	//	Tautologies

	test((aA>aA) ==FALSE);

	test((aA<aA) ==FALSE);

	test((aA>=aA)==TRUE);

	test((aA<=aA)==TRUE);

	test((aA==aA)==TRUE);

	test((aA!=aA)==FALSE);

	if (aOrder!=EEqual)

		{

		test((aA==aB)==FALSE);

  		test((aA!=aB)==TRUE);

		}

	if (aOrder==ELessThan)

		{

		test((aA<aB) ==TRUE);

		test((aA<=aB)==TRUE);

		test((aA>aB) ==FALSE);

		test((aA>=aB)==FALSE);

		}

	if (aOrder==EEqual)

		{

		test((aA==aB)==TRUE);

		test((aA!=aB)==FALSE);

		test((aA>=aB)==TRUE);