/* ------------------------------------------------------------------
 * Copyright (C) 1998-2009 PacketVideo
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *      http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either
 * express or implied.
 * See the License for the specific language governing permissions
 * and limitations under the License.
 * -------------------------------------------------------------------
 */
/*

 Filename: pv_sine.c

------------------------------------------------------------------------------
 REVISION HISTORY


 Who:                                   Date: MM/DD/YYYY
 Description:

------------------------------------------------------------------------------
 INPUT AND OUTPUT DEFINITIONS

    Int32 x             32-bit integer angle (in Q30) between 0 and pi/2


------------------------------------------------------------------------------
 FUNCTION DESCRIPTION

    Find the sine of a number between 0 and pi/2
------------------------------------------------------------------------------
 REQUIREMENTS


------------------------------------------------------------------------------
 REFERENCES

------------------------------------------------------------------------------
 PSEUDO-CODE

------------------------------------------------------------------------------
*/


/*----------------------------------------------------------------------------
; INCLUDES
----------------------------------------------------------------------------*/

#ifdef AAC_PLUS

#ifdef PARAMETRICSTEREO

#include "pv_audio_type_defs.h"
#include "fxp_mul32.h"
#include "pv_sine.h"

/*----------------------------------------------------------------------------
; MACROS
; Define module specific macros here
----------------------------------------------------------------------------*/


/*----------------------------------------------------------------------------
; DEFINES
; Include all pre-processor statements here. Include conditional
; compile variables also.
----------------------------------------------------------------------------*/

/*----------------------------------------------------------------------------
; LOCAL FUNCTION DEFINITIONS
; Function Prototype declaration
----------------------------------------------------------------------------*/

/*----------------------------------------------------------------------------
; LOCAL STORE/BUFFER/POINTER DEFINITIONS
; Variable declaration - defined here and used outside this module
----------------------------------------------------------------------------*/

#define R_SHIFT     30

#define Q_fmt(x)   (Int32)(x*((Int32)1<<R_SHIFT) + (x>=0?0.5F:-0.5F))

const Int32 sin_table[9] =
{
    Q_fmt(0.00001724684028), Q_fmt(-0.00024606242846),
    Q_fmt(0.00007297328923), Q_fmt(0.00826706596417),
    Q_fmt(0.00003585160465), Q_fmt(-0.16667772526248),
    Q_fmt(0.00000174197440), Q_fmt(0.99999989138797),
    Q_fmt(0.00000000110513)
};


/*----------------------------------------------------------------------------
; EXTERNAL FUNCTION REFERENCES
; Declare functions defined elsewhere and referenced in this module
----------------------------------------------------------------------------*/

/*----------------------------------------------------------------------------
; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES
; Declare variables used in this module but defined elsewhere
----------------------------------------------------------------------------*/

/*----------------------------------------------------------------------------
; FUNCTION CODE
----------------------------------------------------------------------------*/


Int32 pv_sine(Int32 z)
{
    Int32 sine;
    Int32 i;
    const Int32 *pt_table = sin_table;
    Int32 sign = 0;

    if (z < 0)
    {
        z = -z;
        sign = 1;
    }

    if (z > Q_fmt(0.0015))
    {
        sine  = fxp_mul32_Q30(*(pt_table++), z);

        for (i = 7; i != 0; i--)
        {
            sine += *(pt_table++);
            sine  = fxp_mul32_Q30(sine, z);
        }

    }
    else
    {
        sine = z;  /*  better approximation in this range */
    }

    if (sign)
    {
        sine = -sine;
    }

    return sine;
}



Int32 pv_cosine(Int32 z)
{
    Int32 cosine;

    if (z < 0)
    {
        z = -z;     /* sign does not play a role in cosine */
    }

    if (z > Q_fmt(0.0015))
    {
        z = Q_fmt(1.57079632679490) - z;   /* pi/2 - z */

        cosine  = pv_sine(z);
    }
    else
    {   /*  better approximation in this range  */
        cosine = Q_fmt(0.99999999906868) - (fxp_mul32_Q30(z, z) >> 1);
    }

    return cosine;
}

#endif

#endif