/*

* Copyright (c) 2006-2007 Nokia Corporation and/or its subsidiary(-ies). 

* All rights reserved.

* This component and the accompanying materials are made available

* under the terms of "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:   ?Description

*

*/

#include <e32cmn.h>

#include "uiacceltk/huifixmath.h"

static const TInt KFixPi_2 = 0x1921F;

static const TInt KSinMagic = 0x78000;

#define FIX_PT_BITS 16

#define SWAP(x, y) {(x)^=(y);(y)^=(x);(x)^=(y);}

EXPORT_C TInt HuiFixMath::FixSin(TInt aRadians)

	{

    if (aRadians < (-KFixPi*2) || aRadians > (KFixPi*2))

        {

        aRadians%=KFixPi*2;

        }

    if (aRadians < -KFixPi_2)

        {

        aRadians+=2*KFixPi;

        }

    // PI/2 - PI

    if (aRadians > KFixPi_2 && aRadians < KFixPi)

        {

        aRadians = KFixPi-aRadians;

        }

    // PI-PI*3/4

    else if (aRadians >=KFixPi && aRadians < (KFixPi+KFixPi_2))    

        {

        aRadians = -(aRadians-KFixPi);

        }

    // PI3*/4-2PI

    else if (aRadians >=(KFixPi+KFixPi_2) && aRadians <=2*KFixPi)

        {

        aRadians-=KFixPi*2;

        }

    else

        {

        // for PC lint

        }

    TInt64 value = aRadians;

    value = (value*value*value);

    value = (value / KSinMagic)>>FIX_PT_BITS;

    return FixMul(aRadians-value,0xF2E0);

	}

EXPORT_C TInt HuiFixMath::FixCos(TInt aRadians)

	{

	return FixSin(KFixPi_2-aRadians);

	}

EXPORT_C TUint HuiFixMath::FixHypot(TInt aVal1,TInt aVal2)

    {

    TUint64 r;    

    TUint p = Abs(aVal1);

    TUint q = Abs(aVal2);

    if (q == 0)

        {

        return p;

        }

    if (p<q)

        {

        SWAP(p,q)

        }

    for (;;)

        {

        r = ((TUint64(q))<<FIX_PT_BITS)/p;

        r = (r*r) >> FIX_PT_BITS;

        if (r == 0)

            {

            return p;

            }

        r = (r<<FIX_PT_BITS)/(0x40000+r);

        p = p + 2*TUint(((r*p)>>FIX_PT_BITS));

        q = (TUint((q*r))) >> FIX_PT_BITS;

        }

    }