1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
|
/*
* Copyright (C) 2004-2010 NXP Software
* Copyright (C) 2010 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.
*/
/************************************************************************/
/* */
/* Project:: */
/* $Author: beq07716 $*/
/* $Revision: 1000 $*/
/* $Date: 2010-06-28 13:08:20 +0200 (Mon, 28 Jun 2010) $*/
/* */
/************************************************************************/
#include "LVM_Types.h"
#include "LVM_Macros.h"
#include "ScalarArithmetic.h"
/*-------------------------------------------------------------------------*/
/* FUNCTION: */
/* LVM_Polynomial */
/* */
/* DESCRIPTION: */
/* This function performs polynomial expansion */
/* Y = (A0 + A1*X + A2*X2 + A3*X3 + ….. + AN*xN) << AN+1 */
/* */
/* LVM_INT32 LVM_Polynomial(LVM_UINT16 N, */
/* LVM_INT32 *pCoefficients, */
/* LVM_INT32 X) */
/* */
/* PARAMETERS: */
/* */
/* N is the polynomial order */
/* pCoefficients is the ptr to polynomial coefficients A0,A1.. in Q.31 */
/* X is the input variable */
/* */
/* RETURNS: */
/* The result of the polynomial expansion in Q1.31 format */
/*-------------------------------------------------------------------------*/
LVM_INT32 LVM_Polynomial(LVM_UINT16 N,
LVM_INT32 *pCoefficients,
LVM_INT32 X)
{
LVM_INT32 i;
LVM_INT32 Y,A,XTemp,Temp,sign;
Y=*pCoefficients; /* Y=A0*/
pCoefficients++;
if((LVM_UINT32)X==0x80000000)
{
Temp=-1;
sign=Temp;
for(i=1;i<=N;i++)
{
Y+=((*pCoefficients)*sign);
pCoefficients++;
sign*=Temp;
}
}
else
{
XTemp=X;
for(i=N-1;i>=0;i--)
{
A=*pCoefficients;
pCoefficients++;
MUL32x32INTO32(A,XTemp,Temp,31)
Y+=Temp;
MUL32x32INTO32(XTemp,X,Temp,31)
XTemp=Temp;
}
}
A=*pCoefficients;
pCoefficients++;
if(A<0)
{
A=Abs_32(A);
Y=Y>>A;
}
else
{
Y = Y<<A;
}
return Y;
}
|