/*
* Copyright (C) 2008 The Android Open Source Project
*
* 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.
*/
/* ---- includes ----------------------------------------------------------- */
#include "b_BasicEm/Phase.h"
#include "b_BasicEm/Math.h"
/* ------------------------------------------------------------------------- */
/* ========================================================================= */
/* */
/* ---- \ghd{ auxiliary functions } ---------------------------------------- */
/* */
/* ========================================================================= */
/* ------------------------------------------------------------------------- */
/* ========================================================================= */
/* */
/* ---- \ghd{ constructor / destructor } ----------------------------------- */
/* */
/* ========================================================================= */
/* ------------------------------------------------------------------------- */
/* ========================================================================= */
/* */
/* ---- \ghd{ operators } -------------------------------------------------- */
/* */
/* ========================================================================= */
/* ------------------------------------------------------------------------- */
/* ========================================================================= */
/* */
/* ---- \ghd{ query functions } -------------------------------------------- */
/* */
/* ========================================================================= */
/* ------------------------------------------------------------------------- */
/* ========================================================================= */
/* */
/* ---- \ghd{ modify functions } ------------------------------------------- */
/* */
/* ========================================================================= */
/* ------------------------------------------------------------------------- */
/* ========================================================================= */
/* */
/* ---- \ghd{ I/O } -------------------------------------------------------- */
/* */
/* ========================================================================= */
/* ------------------------------------------------------------------------- */
/* ========================================================================= */
/* */
/* ---- \ghd{ exec functions } --------------------------------------------- */
/* */
/* ========================================================================= */
/* ------------------------------------------------------------------------- */
/* ------------------------------------------------------------------------- */
#ifndef bbs_SIN_INTERPOLATION_METHOD_2
const int32 bbs_sin32_table1G[] =
{
0, 1608, 411648, 1607, 823040, 1606, 1234176, 1602,
1644288, 1599, 2053632, 1594, 2461696, 1588, 2868224, 1581,
3272960, 1574, 3675904, 1564, 4076288, 1556, 4474624, 1545,
4870144, 1533, 5262592, 1521, 5651968, 1508, 6038016, 1493,
6420224, 1478, 6798592, 1463, 7173120, 1445, 7543040, 1428,
7908608, 1409, 8269312, 1390, 8625152, 1369, 8975616, 1348,
9320704, 1327, 9660416, 1303, 9993984, 1280, 10321664, 1256,
10643200, 1231, 10958336, 1205, 11266816, 1178, 11568384, 1151,
11863040, 1124, 12150784, 1094, 12430848, 1066, 12703744, 1036,
12968960, 1005, 13226240, 974, 13475584, 942, 13716736, 910,
13949696, 877, 14174208, 844, 14390272, 810, 14597632, 775,
14796032, 741, 14985728, 705, 15166208, 670, 15337728, 634,
15500032, 597, 15652864, 561, 15796480, 523, 15930368, 486,
16054784, 448, 16169472, 409, 16274176, 372, 16369408, 333,
16454656, 295, 16530176, 255, 16595456, 217, 16651008, 177,
16696320, 138, 16731648, 99, 16756992, 59, 16772096, 20,
16777216, -20, 16772096, -59, 16756992, -99, 16731648, -138,
16696320, -177, 16651008, -217, 16595456, -255, 16530176, -295,
16454656, -333, 16369408, -372, 16274176, -409, 16169472, -448,
16054784, -486, 15930368, -523, 15796480, -561, 15652864, -597,
15500032, -634, 15337728, -670, 15166208, -705, 14985728, -741,
14796032, -775, 14597632, -810, 14390272, -844, 14174208, -877,
13949696, -910, 13716736, -942, 13475584, -974, 13226240, -1005,
12968960, -1036, 12703744, -1066, 12430848, -1094, 12150784, -1124,
11863040, -1151, 11568384, -1178, 11266816, -1205, 10958336, -1231,
10643200, -1256, 10321664, -1280, 9993984, -1303, 9660416, -1327,
9320704, -1348, 8975616, -1369, 8625152, -1390, 8269312, -1409,
7908608, -1428, 7543040, -1445, 7173120, -1463, 6798592, -1478,
6420224, -1493, 6038016, -1508, 5651968, -1521, 5262592, -1533,
4870144, -1545, 4474624, -1556, 4076288, -1564, 3675904, -1574,
3272960, -1581, 2868224, -1588, 2461696, -1594, 2053632, -1599,
1644288, -1602, 1234176, -1606, 823040, -1607, 411648, -1608
};
#else
const int32 bbs_sin32_table2G[] =
{
0, 12907, -122,
209469440, 12662, -368,
410894336, 11926, -596,
596525056, 10733, -802,
759234560, 9129, -978,
892780544, 7168, -1112,
992002048, 4939, -1210,
1053097984, 2516, -1256,
1073741824, -4, -1256,
1053097984, -2519, -1210,
992002048, -4944, -1112,
892780544, -7173, -978,
759234560, -9129, -802,
596525056, -10734, -596,
410894336, -11926, -368,
209469440, -12663, -122
};
#endif
int32 bbs_sin32( phase16 phaseA )
{
#ifndef bbs_SIN_INTERPOLATION_METHOD_2
int32 oL = ( phaseA & 0x00FF );
uint16 indexL = ( ( phaseA & 0x7F00 ) >> 8 ) << 1;
int32 sinL = bbs_sin32_table1G[ indexL ] + oL * bbs_sin32_table1G[ indexL + 1 ];
if( ( phaseA & 0x8000 ) != 0 )
{
return -sinL;
}
else
{
return sinL;
}
#else /*bbs_SIN_INTERPOLATION_METHOD_2*/
int32 o1L = ( phaseA & 0x07FF );
int32 o2L = ( o1L * o1L ) >> 8;
uint16 indexL = ( ( phaseA & 0x7800 ) >> 11 ) * 3;
int32 sinL = bbs_sin32_table2G[ indexL ] + ( ( o1L * bbs_sin32_table2G[ indexL + 1 ] ) << 3 ) + o2L * bbs_sin32_table2G[ indexL + 2 ];
if( ( phaseA & 0x8000 ) != 0 )
{
return -sinL >> 6;
}
else
{
return sinL >> 6;
}
#endif /*bbs_SIN_INTERPOLATION_METHOD_2*/
}
/** computation of sine tables (do not uncomment or remove)
void sin1Table()
{
long iL;
for( iL = 0; iL < 128; iL++ )
{
int32 phase1L = iL * 256;
int32 phase2L = phase1L + 256;
double angle1L = ( M_PI * phase1L ) / 32768;
double angle2L = ( M_PI * phase2L ) / 32768;
int32 sin1L = ( sin( angle1L ) * 65536 );
int32 sin2L = ( sin( angle2L ) * 65536 );
int32 diffL = sin2L - sin1L;
eout << iL << ": " << ( sin1L << 8 ) << " + oL * " << diffL << endl;
}
}
void sin2Table()
{
long iL;
for( iL = 0; iL < 16; iL++ )
{
int32 p0L = iL * ( 1 << 11 );
int32 p1L = p0L + ( 1 << 10 );
int32 p2L = p0L + ( 1 << 11 );
double a0L = ( M_PI * p0L ) / ( 1 << 15 );
double a1L = ( M_PI * p1L ) / ( 1 << 15 );
double a2L = ( M_PI * p2L ) / ( 1 << 15 );
int32 s0L = ( sin( a0L ) * ( 1 << 16 ) );
int32 s1L = ( sin( a1L ) * ( 1 << 16 ) );
int32 s2L = ( sin( a2L ) * ( 1 << 16 ) );
int32 aL = 4 * s1L - 3 * s0L - s2L;
int32 bL = 2 * s2L + 2 * s0L - 4 * s1L;
eout << iL << ": " << ( s0L << 14 ) << " + ( ( o1L * " << aL << " ) << 3 )"
<< " + o2L * " << bL << endl;
}
}
*/
/* ------------------------------------------------------------------------- */
int32 bbs_cos32( phase16 phaseA )
{
return bbs_sin32( ( phase16 )( phaseA + bbs_M_PI_2_16 ) );
}
/* ------------------------------------------------------------------------- */
int16 bbs_sin16( phase16 phaseA )
{
return bbs_sin32( phaseA ) >> 10;
}
/* ------------------------------------------------------------------------- */
int16 bbs_cos16( phase16 phaseA )
{
return bbs_sin32( ( phase16 )( phaseA + bbs_M_PI_2_16 ) ) >> 10;
}
/* ------------------------------------------------------------------------- */
const int32 bbs_atan16_tableG[] =
{
0, 325, 332800, 326, 666624, 326, 1000448, 325,
1333248, 324, 1665024, 323, 1995776, 323, 2326528, 322,
2656256, 320, 2983936, 319, 3310592, 317, 3635200, 316,
3958784, 314, 4280320, 312, 4599808, 310, 4917248, 308,
5232640, 306, 5545984, 303, 5856256, 301, 6164480, 298,
6469632, 296, 6772736, 292, 7071744, 291, 7369728, 287,
7663616, 284, 7954432, 281, 8242176, 279, 8527872, 275,
8809472, 272, 9088000, 269, 9363456, 265, 9634816, 263,
9904128, 259, 10169344, 256, 10431488, 252, 10689536, 249,
10944512, 246, 11196416, 243, 11445248, 239, 11689984, 236,
11931648, 233, 12170240, 230, 12405760, 226, 12637184, 223,
12865536, 219, 13089792, 217, 13312000, 213, 13530112, 210,
13745152, 207, 13957120, 204, 14166016, 201, 14371840, 198,
14574592, 195, 14774272, 192, 14970880, 189, 15164416, 186,
15354880, 183, 15542272, 180, 15726592, 178, 15908864, 175,
16088064, 172, 16264192, 169, 16437248, 167, 16608256, 165
};
phase16 bbs_atan16( uint32 valA )
{
uint32 oL = valA & 0x03FF;
uint16 indexL = ( valA >> 10 ) << 1;
uint32 phaseL = bbs_atan16_tableG[ indexL ] + oL * bbs_atan16_tableG[ indexL + 1 ];
return ( phase16 )( phaseL >> 11 );
}
/* ------------------------------------------------------------------------- */
phase16 bbs_phase16( int32 xA, int32 yA )
{
uint32 xL = ( xA > 0 ) ? xA : -xA;
uint32 yL = ( yA > 0 ) ? yA : -yA;
phase16 phaseL;
if( xL == 0 && yL == 0 ) return 0;
if( xL == yL )
{
phaseL = bbs_M_PI_4_16; /*PI/4*/
}
else if( xL > yL )
{
if( yL >= 65536 ) /* avoid overflow (1 << 16) */
{
uint32 shiftL = bbs_intLog2( yL ) - 15;
xL >>= shiftL;
yL >>= shiftL;
}
phaseL = bbs_atan16( ( yL << 16 ) / xL );
}
else
{
if( xL >= 65536 ) /* avoid overflow (1 << 16) */
{
uint32 shiftL = bbs_intLog2( xL ) - 15;
xL >>= shiftL;
yL >>= shiftL;
}
phaseL = bbs_M_PI_2_16 - bbs_atan16( ( xL << 16 ) / yL );
}
if( xA >= 0 )
{
if( yA >= 0 )
{
return phaseL;
}
else
{
return -phaseL;
}
}
else
{
if( yA >= 0 )
{
return bbs_M_PI_16 - phaseL;
}
else
{
return phaseL - bbs_M_PI_16;
}
}
}
/* ------------------------------------------------------------------------- */
/* ========================================================================= */