A refined affine approximation method of multiplication for range analysis in word-length optimization

Affine arithmetic (AA) is widely used in range analysis in word-length optimization of hardware designs. To reduce the uncertainty in the AA and achieve efficient and accurate range analysis of multiplication, this paper presents a novel refined affine approximation method, Approximation Affine based on Space Extreme Estimation (AASEE). The affine form of multiplication is divided into two parts. The first part is the approximate affine form of the operation. In the second part, the equivalent affine form of the estimated range of the difference, which is introduced by the approximation, is represented by an extra noise symbol. In AASEE, it is proven that the proposed approximate affine form is the closest to the result of multiplication based on linear geometry. The proposed equivalent affine form of AASEE is more accurate since the extreme value theory of multivariable functions is used to minimize the difference between the result of multiplication and the approximate affine form. The computational complexity of AASEE is the same as that of trivial range estimation (AATRE) and lower than that of Chebyshev approximation (AACHA). The proposed affine form of multiplication is demonstrated with polynomial approximation, B-splines, and multivariate polynomial functions. In experiments, the average of the ranges derived by AASEE is 59% and 89% of that by AATRE and AACHA, respectively. The integer bits derived by AASEE are 2 and 1 b less than that by AATRE and AACHA at most, respectively.


Introduction
As a method of representing real numbers, floating point can support a wide dynamic range and high precision of values.It has been thus commonly used in signal processing, such as image processing, speech processing, and digital signals processing, to represent signals.When these applications are implemented on hardware for high speed and stability, the signals need to be represented in fixed point to optimize the performance of area, power, and speed of the hardware.Hence, the values in floatingpoint need to be converted to those in fixed point.This process is named as word-length optimization.Its goal is to achieve optimal system performance while satisfying the specification on the system output precision.Wordlength optimization involves range analysis and precision analysis.The former one is to find the minimum word length of the integer part of the value, while the latter one focuses on the optimization of the fractional part of the word length.
Word-length optimization has been proven to be an NPhard problem [1].It can be usually classified into dynamic analysis [2][3][4][5][6][7] and static analysis [8][9][10][11][12][13][14][15][16][17][18][19][20].By analyzing a large set of stimuli signals, dynamic analysis is applicable to all types of systems.However, it will take long time on simulation to provide sufficient confidence.Also, the precision for the signals without simulation cannot be guaranteed.Comparatively, the static analysis is an automated and efficient word-length optimization method and more applicable to large designs when compared to dynamic analysis.The static analysis mainly uses the characteristics of the input signals to estimate the word length conservatively, which can result in overestimation [12] to some extent.As a part of word-length optimization, the range analysis can also been classified in the same way.http://asp.eurasipjournals.com/content/2014/1/36Affine arithmetic (AA) [21] is often used for range analysis in static analysis.In AA, every signal must be represented in an affine form, which is a first-degree polynomial.As AA tracks the correlations among range intervals of signals, it can provide more accurate wordlength range.This makes it suitable for range analysis of the result of linear operations.It is noted that besides linear operations, nonlinear operations, such as multiplication, are also involved in hardware operations, typically in linear time invariant (LTI) systems.AA cannot provide an exact affine form for nonlinear operations.To solve this problem, Stolfi and de Figueiredo [22] proposed affine approximation methods for multiplication, which include trivial range estimation (AATRE) and Chebyshev approximation (AACHA).AATRE is efficient for computation, but the range produced by it can be four times of real range at most.The accumulation of the uncertainty of all signals in the computational chain may result in an error explosion, which is unacceptable in application.Such overestimation obviously cannot satisfy the accuracy requirement of the system, which limits the application of AATRE in large systems.The uncertainty of AACHA is less than AATRE, however, it is too complex to be used in large systems.Since LTI operations are accurately covered by AA, the proposed method is applied in the field of the range analysis of word-length optimization in this paper.
A novel affine approximation method, Approximation Affine based on Space Extreme Estimation (AASEE), is proposed to reduce the uncertainty of multiplication and achieve an accurate and efficient range analysis of multiplication in this paper.To analyze the uncertainty conveniently, we use two parts to divide the different parts of all the approximation methods for multiplication, which include AATRE, AACHA, and AASEE.The first part is named as approximate affine form, which is approximated to the nonlinear operation.The second part is named as equivalent affine form, which is the equivalent affine form of the estimated range of the difference between the result of multiplication and the approximate affine form.The more accurate the two parts are, the more accurate the approximation method is.Based on linear geometry [23], it is proven that the proposed approximate affine form is the closest to the result of multiplication.To derive the equivalent affine form, we use the extreme value theory of multivariable functions [24] to estimate the upper and lower bounds of the difference in space, and the difference is introduced by the approximation of the first part.The uncertainty of the proposed method is minimized.The accuracy of the resulting affine form by AASEE is higher than that by AATRE and averagely higher than that by AACHA.Meanwhile, the computational complexity of AASEE is equivalent to that of AATRE and lower than that of AACHA.The rest of this paper is organized as follows.Background of range analysis for multiplication is presented in Section 2. Section 3 presents the method of derivation of the two parts for multiplication.The refined affine form of multiplication, AASEE, is presented in next section.In Section 5, we compare the computational complexity and the accuracy among AASEE to AATRE and AACHA.The case studies and experimental results are demonstrated in Section 6. Section 7 concludes the paper.

Related work
Interval arithmetic (IA) and affine arithmetic (AA) have been widely used in range analysis in word-length optimization.
IA [25] is a range arithmetic theory which is firstly presented by Moore in 1962.Cmar [2] employs it for range analysis of digital signal processing (DSP) systems.Carreras [20] presents a method based on IA.To reduce the oversized word length, the method provides the probability density functions that can be used when some truncation must be performed due to constraints in the specification.IA is not suitable for most real-world applications, since it could lead to drastic overestimation of the true range.
AA [21] is proposed to overcome the weakness of IA by Stolfi in 1993.In [8,9], Fang uses AA to analyze word-length optimization.Both range and precision are represented by the same affine form, which limits the optimization.Pu and Ha [10] also use AA for wordlength optimization.Simultaneously, they use two different affine forms for range analysis and precision analysis, respectively, and achieve more refined result of wordlength optimization.Similarly, Lee et al. [11] develop an automatic optimization approach, which is called MiniBit, to produce accuracy-guaranteed solutions, and area is minimized while meeting an error constraint.Osborne [12] uses both IA and AA for range analysis for different situations.Computation using either of the two methods in the design is time-consuming.The problem of overestimation is serious due to the approximation of the nonlinear operations.
Since AA cannot be used in the systems with infinite number of loops, an improved approach, quantized AA (QAA), has been proposed in [13] for linear time-invariant systems with feedback loops.This method can provide fast and tight estimation of the evolution of large sets of numerical inputs, using only an affine-based simulation, but it does not provide the exact bounds.
AATRE [22] is adopted for multiplication in most of the works for the low computational complexity.But the uncertainty of the range by AATRE is very large.To adjust the trade-off between the accuracy of approximation and computational complexity, Zhang [14] introduces a new http://asp.eurasipjournals.com/content/2014/1/36parameter N in the N-level simplified affine approximation (N-SAA).This method is faster than AACHA and more accurate than AATRE, but it is more complex than AATRE.Furthermore, it is troublesome to choose a suitable N. A method of range analysis is proposed by Pang [26].This method combines methods of IA, AATRE, and arithmetic transform (AT); and the result of the method is more accurate than AATRE, while the CPU implementation time is longer than AATRE.To deal with applications from the scientific computing domain, Kinsman [17,18] uses the computational methods based on Satisfiability Modulo Theory.Search efficiency of this method is improved leading to tighter bounds and thus smaller word length.
For all the existing methods, the accuracy of approximation is improved at the expense of the computational complexity.This paper presents an affine approximation method for multiplication, which achieves better trade-off between accuracy and computational complexity.

Range analysis
Range analysis involves studying the data range of every signal and minimizing the integer word lengths for signals on the premise that the signals in the design have enough bits to accommodate this range.The range of signal x is represented by x = [x min , x max ], where the two real numbers, x min and x max , denote the lower and upper bounds of x, respectively.The required integer part of the word length for signal x, which is represented as IWL x , can be derived by: where In ( 1), all the signals in the design are assumed to be expressed as signed numbers, and the sign bit is taken into account in IWL x .According to (1), once the range of a signal is decided, the integer part of word length of the signal can be derived.

Affine arithmetic
AA is widely applied for range analysis.In AA, an uncertain signal x is represented by an affine form as a firstdegree polynomial [22]: For the signal x, x 0 is the central value, and ε i is the ith noise symbol.ε i denotes an independent uncertainty source that contributes to the total uncertainty of the signal x, and x i is its coefficient.
The upper and lower bounds for the range of x can be represented as With x min and x max , the input interval x = [x min , x max ] can be converted into an equivalent affine form as (4), using only one independent noise symbol.
AA can keep correlations among the signals of the computational chain by contributing the sample noise symbol ε i to each signal [22].
For multiplication, AATRE and AACHA are typical approximation methods.
The affine form of AATRE is Suppose M 1 = max(n 1 , n 2 ), in which n 1 and n 2 denote the number of the noise symbol, whose coefficient is nonzero, of x and ŷ, respectively.The computational complexity of AATRE is O(M 1 ).
AACHA provides a better approximation result, but it is more complex.The affine form of AACHA is where a and b denote the minimum and the maximum of the range of The complexity of computing the both extremal values, a and b, is O(M 2 log M 2 ).As M 1 ≤ M 2 , the computational complexity of AATRE is lower than that of AACHA [22].

Extreme value theory
The proposed approximation is based on the extreme value theory of multivariable functions [24].
According to the extreme value theory of multivariable functions, the Hessian matrix of the function, H, and Jacobian matrix of the function, J, can be used to find http://asp.eurasipjournals.com/content/2014/1/36 the local maxima and the local minima.Hessian matrix of function Here we use The principal minor determinants are used to determine if a matrix is positive or negative definite or semidefinite.
It is necessary and sufficient for a positive semidefinite matrix that all the principal minor determinants of the matrix are nonnegative real numbers.
It is necessary and sufficient for a negative semidefinite matrix that all the odd order principal minor determinants of the matrix are non-positive real numbers and all the even order principal minor determinants of the matrix are nonnegative real numbers.

Derivation of the two parts for multiplication
A generic nonlinear operation z ← f (x, ŷ) proposed in [22] can be described by (8): Since the operation f is nonlinear, f * (ε 1 , . . ., ε n ) cannot be expressed exactly as an affine combination of the noise symbols, ε i .Under this case, an approximate affine form of the operation, which is represented as f z , must be used to approximate f * (ε 1 , . . ., ε n ).The difference introduced by this approximation, d f = f * − f z , can be expressed by an equivalent affine form of the estimated range of the difference, which is represented as d.Hence, the affine form of z can be expressed as In ( 9), f z is a first-degree function of ε i and can be expressed as (10) The computational complexity of computing the true range of d f is very high in a practical application.The estimated range of d f is utilized instead of the true range.Suppose d max and d min denote the upper and lower bounds of the estimated range of d f , respectively.According to (4), the d can be expressed as ( 11) With ( 10) and ( 11), the affine form of z can be represented as For multiplication, z can be expressed as The first three items of ( 13) form an affine form and the last term is a quadratic term.Its affine form can also be represented as (12).
According to the definition of f z in (10) and d in ( 11), AATRE and AACHA can also be represented by f z and d.For AATRE in (5), the f z and d are defined as For AACHA in (6), the f z and d are defined as In the existing affine approximation methods of AATRE and AACHA, d max and d min are estimated in the XY plane.In these methods, the same noise symbol of different variables is considered to be independent.Hence, the range of d is much larger than that of d f .The difference between d and d f will propagate to ẑ and result in uncertainty.
To describe the multiplication accurately, we use ε i as the input arguments and estimate the range of z in the (n + 1)-dimensional space E n+1 .The (n + 1)-dimensional space E n+1 is labeled as (ε 1 , . . ., ε n , z).In space E n+1 , a first-degree polynomial function can be expressed as a http://asp.eurasipjournals.com/content/2014/1/36(n + 1)-dimensional hyperplane and a nonlinear polynomial function denotes a (n + 1)-dimensional space curved surface.The approximate affine form in (10) denotes a (n + 1)-dimensional hyperplane in E n+1 .Each hyperplane in E n+1 can be viewed as a parallel translation of a tangent hyperplane at a certain point of (n + 1)-dimensional space curved surface.Hence, all possible approximate affine forms for z can be regarded as the (n + 1)-dimensional tangent hyperplanes at all points of (n + 1)-dimensional space curved surface in E n+1 .The translation amount is taken into account in d f , which is approximated by d.In space E n+1 , d f can be viewed as the function of the distance between the points of space curved surface and the tangent hyperplane.Figure 1 shows an example of x = 1 + ε 1 + 5ε 2 and ŷ = 3 − 6ε 1 + ε 2 .The space is labeled as (ε 1 , ε 2 , z).The red mesh surface represents the function z The blue plane represents the tangent plane f z , z = 3 − 3ε 1 + 16ε 2 , at the point z α = (0, 0, 3).All the possible approximate affine forms for z are the tangent planes of all the points.d f is a function of distance between z and f z .
Here we use f z α in (18) to represent the tangent hyperplane at the point Then, the possible approximate affine form can be represented as f z α , too.
In (18), z ε n are the partial derivatives of z with respect to the variables ε n at the point z α .
With the estimated range of d f , the maximum absolute error of d f can be expressed as To reduce the uncertainty, f z must be the most closed to the result of multiplication.Hence, f z is the tangent hyperplane whose maximum absolute error is minimum among that of all the possible affine form f z α , that is, The geometrical meaning of f z denotes the tangent hyperplane whose maximum absolute error is minimized.
f z is derived by the range of d f , while d is the equivalent affine form of d f .It is very complex to compute the true range of d f .With d in (11), the uncertainty in AA for nonlinear operations is generated due to the difference between the true range of d f and the estimated range of d f .
It is much tighter and easier to estimate range of d f in E n+1 space than in the XY plane.Based on the extreme value theory of multivariable functions, the estimated range of d f in AASEE is derived.
With more accurate d max and d min , f z and d can be calculated more precisely, and AASEE can achieve a refined affine approximation result.
In the next sections, the estimated range of d f will be derived firstly, and the two parts will be derived later.

Estimated range of the difference
For multiplication, which is expressed as (13), the value of z at the point z α is The partial derivatives of z with respect to the variable ε i at the point z α are Upon substitution for z α and z ε i , the tangent hyperplane f z α can be expressed as The difference between the tangent hyperplane f z α and (n + 1)-dimensional quadratic surface z is Suppose d emax and d emin denote the estimated maximum and minimum of the function value at the domain boundary respectively, and d fimax and d fimin denote the local maxima and the local minima, respectively.The estimated maximum and minimum of multivariable function d f , d max and d min , can be expressed as According to (24), the function value at the domain boundary, d fe , is represented by where To simplify, we observe the extreme case of ∀ε i = ±1.Under this case, for the first item, it is always positive when i = j.Hence, the estimated function value at the domain boundary, d e , is expressed as x i y j ε i ε α j where ∀ε i = ±1.
Hence, the maximum and minimum of d e , d emax and d emin are derived as To simply compare, d fimax and d fimin in ( 25) and ( 26) can be expressed as x i y j ε α i ε α j , where ε i , ε j = (−1, 1), and As the example in Section 3, Figure 2 shows the function of

Hessian matrix of function
From (33), we can see that H is independent of ε i .It is a expression of x i and y i .This means that H is same for all the points in the domain.
To determine if H is positive or negative definite or semidefinite, its principal minor determinants are derived as x i y j + x j y i x i y j + x j y i 2x j y j = −(x i y j − x j y i ) 2 (35) As introduced in Section 2.4, H is a positive semidefinite matrix, iff it satisfies H is a negative semidefinite matrix, iff it satisfies If it satisfies neither (37) nor (38), which means it satisfies (39), H is an indefinite matrix as According to (37), (38), and (39), we can compare d emax , d emin , d fimax , and d fimin , which are expressed as ( 29), (30), (31), and (32), respectively.Based on ( 25) and ( 26), d max and d min can be identified.

Lemma 1. The estimated maximum of function d f , d max equals to the estimated maximum of the function value at the domain boundary, and the estimated minimum of function d f , d min equals to the estimated minimum of the function value at the domain boundary. This can be expressed as
Proof.There are two cases to consider, as ∃x i y i < 0 and ∀x i y i ≥ 0. http://asp.eurasipjournals.com/content/2014/1/36 For ∃x i y i < 0, (39) is satisfied and H is indefinite.The stationary point is a saddle point, such as the point P in Figure 2. Neither d fimax nor d fimin exists in d f , that is, According to (41), Lemma 1 can be proven in this case.
For ∀x i y i ≥ 0, H may be positive semidefinite or negative semidefinite.d f may have local minima or local maxima under this condition. As , the following inequalities are established: If a local maximum lies at z α , the difference between d emax and d fimax is ∀x i y i ≥ 0, there exists According to (25) and (46), we can prove that Similarly, if a local minimum lies at z α , the difference between d emin and d fimin is As ∀x i y i ≥ 0 in (48), the inequality (49) can be proven: According to (26) and (49), we can prove that As ( 47) and (50) are established, Lemma 1 can be proven in the case of ∀x 1 y 1 ≥ 0. Combining these two cases, Lemma 1 is proven.
According to Lemma 1, d max and d min at a point z α can be computed as d emax and d emin in ( 29) and (30).
Proof.According to Lemma 1,(29), and (30), the maximum absolute error of d f is So the maximum absolute error between the tangent hyperplane f z 0 at the point z 0 = z 0 = (0, 0, . . ., 0) and Suppose that there is another point z α = z 0 , which is typically represented by , and ε i cannot be equal to 0 for all i, i = 1 . . .n.The maximum absolute error between the tangent hyperplane f z α at point z α and (n + 1)-dimensional quadratic surface xŷ is e a (z α ) and e a (z 0 ) can be compared by Because e a (z 0 ) ≤ e a (z α ), the tangent hyperplane f z 0 at the point z 0 = z 0 = (0, 0, . . ., 0) is the tangent hyperplane whose maximum absolute error is minimized.
It is proven that the chosen f z is a tangent hyperplane at the point z 0 = z 0 = (0, 0, . . ., 0).According to Lemma 2, f z of AASEE denotes the tangent hyperplane at the point z 0 = (0, 0, . . ., 0) and can be expressed as This f z is the same as the f z s in AATRE and AACHA.http://asp.eurasipjournals.com/content/2014/1/36

Expression of the equivalent affine form in AASEE
According to (55), the d f between the tangent hyperplane f z 0 and the quadratic surface is According to Lemma 1,(29), and (30), the estimated maximum and estimated minimum of d f , d max and d min can be expressed as (57) n = 1 is a special case and d max and d min can be optimized as By combining the two cases, d emax and d emin are rewritten as When n > 1, the range of d can be expressed as According to (11), the affine form of d can be expressed as When n = 1, the range of d can be expressed as The affine form of d can be expressed as (65)

Formulary of AASEE
According to (12), the affine form of AASEE for multiplication is It is impossible to obtain the exact affine form for multiplication in AA.The result of multiplication must be approximated to an affine form.Using ε i as the input arguments, the uncertainty of multiplication in AASEE is reduced.The proposed f z is the most closed to the result of multiplication among all the possible approximate affine forms, and the upper and lower bounds of d in AASEE are much closer to true bounds of d f .Hence, the uncertainty in AASEE is smaller than that in AATRE and AACHA.Formed by such f z and d, AASEE creates a refined affine form of multiplication.

Computational complexity
The computational complexity of an expression is determined by its most complex item.For n > 1, the most complex item is the coefficient of ε n+1 .To make the analysis convenient, we transform this coefficient: The computational complexity of the minuend is O(M 1 ), where M 1 is defined in Section 2.3, while the computational complexity of the subtrahend is less than O(M 1 ).
Hence, the computational complexity of AASEE is O(M 1 ).We can see that it is the same as that of AATRE and is lower than that of AACHA.

Accuracy
The accuracy of d is influential to the accuracy of the affine approximation methods of multiplication.The more accurate d will lead to a more accurate the affine approximation result.http://asp.eurasipjournals.com/content/2014/1/36 |y i |ε n+1 .In this method, the same noise symbol of different variables is considered to be independent.The range of this d is It is much larger than the range of d by AASEE, which is expressed in (62) and (64).
In AACHA, d = a+b 2 + b−a 2 ε n+1 , where a and b are represented the estimated range of d.In this method, a polygon in XY plane is used to find a and b.The domain of xŷ is bounded by the polygon.However, the polygon is larger than the true domain, and all the same noise symbols of different variables are not taken into account together.
All the same noise symbols of different variables are considered together by d of AASEE.It is more accurate than d of AATRE.In the most cases, it is more accurate than d of AACHA, too.

Case studies
The following nonlinear system cases are used to demonstrate the efficiency of the proposed refined affine form of multiplication.These cases are commonly used in signal processing.The first two cases are univariate cases and come from [11].The rest of cases are multivariate polynomial functions and come from [27][28][29].
Table 1 shows the variable ranges and the range intervals, (y max − y min ), of intermediates and output by the three methods.The true range of y lies in [0, 0.6931], and the range interval of output is 0.6931.Suppose R(T), R(C), and R(A) are represented as the ratios of range interval obtained by AATRE, AACHA, and AASEE to the true range interval, respectively.The closer this ratio converges to 1, the more accurate the method is.In this case, as R(T) = 1.33,R(C) = 1.15, and R(A) = 1.03, we can see the range by AASEE is closer to the true range than AATRE and AACHA.

Comparison of range and computational complexity by the three cases
The output ranges by the three methods of case 2 and case 3 can be obtained according to the process of case 1. Table 2 demonstrates the ranges and the integer word lengths by AASEE and comparison among AATRE, AACHA and AASEE.Column c.fun shows the case study and the function of the row.The true output ranges, which are used as reference values, are obtained by numerical method or nonlinear programming technique, which are time-consuming and are not practical to solve the true bounds for large number of signals.From the table, we can see that the ranges, which are derived by AASEE, cover the true ranges and they are smaller than those by AATRE, for all the functions.For these thirteen functions, the ranges, which are derived by AASEE, are smaller than those by AACHA for nine functions, and equal to those by AACHA for two functions.According to (1), the integer word length can be decided by the range.The integer word-length, which is derived by AASEE, is 2 b less than that by AATRE and 1 b less than that by AACHA, at most.Comparing with AATRE, AASEE and AACHA can save 0.54 b on average.
To calculate the estimated range of d f , the values of ∃ε i = ±1, ∀i = 1, 2, . . ., n in (27) are substituted by ∀ε i = ±1 in AASEE.The difference between the estimated range and the true range of d f is introduced by this approximation.In most of the applications, the estimated ranges, which are computed by AASEE, are closer than those by AACHA.However, the estimated minimum and maximum of xŷ on the boundary of the polygon are independent of the value of ε i .In some applications such as functions f 2 and f 8 in Table 2, the results by AASEE are almost the same as those by AACHA.
In Table 3, ratios of range intervals and the computational complexity are compared among AATRE, AACHA, and AASEE.The computational complexity is calculated from the numbers of multiplications and additions.For AACHA, the extreme value of a quadratic function in one variable on a bounded interval needs to be calculated.N m , N a , and N e denote the numbers of multiplications, additions and the extreme value computations of each case, respectively.Table 3 shows that R(T) values are from 1.04 to 281.2, R(C) are from 1.03 to 233.7, and R(A) are from 1.03 to 192.9.The ratios of R(A) to R(T) and R(C) show the accuracy of AASEE compared to AATRE and AACHA, respectively.The average ratios can be used to evaluate the accuracy of the affine approximation methods.The ratios of R(A) to R(T) are from 0.18 to 0.99, and the average of these ratios is 0.59.The ratios of R(A) to R(C) are from 0.33 to 1.17, and the average of these ratios is 0.89.For these 13 cases, on average, the accuracy of AASEE is 1.69 times than that of AATRE and 1.12  times than that of AACHA.The extreme value computation, which is only necessary for AACHA, of the quadratic function is the most complex and time-consuming among the operations.Hence, the computational complexity of AACHA is much higher than that of AATRE and AASEE.The increase rate of the number of multiplications, N m , by AASEE to AATRE is from 0.091 to 1.75, and the average is 0.450.The increase rate of the number of multiplications, N m , by AASEE to AACHA is from 0.2 to 1.833, and the average is 0.567.The increase rate of the number of additions, N a , by AASEE to AATRE is from 0.05 to 3.4, and the average is 0.944.The increase rate of the number of additions, N a , by AASEE to AACHA is from 0 to 0.985, and the average is 0.157.The numbers of multipli- cations and additions of AASEE are increased a few.As shown in Table 3, AACHA is slightly more accurate for functions c 3 .f 2 and c 3 .f8 , but the computational complexity of AACHA is much higher than that of AASEE.

Comparison of the design cost by the three methods
To compare the design cost, the system area by the three methods, the fractional word lengths are obtained by the precise analysis in [11].Typically, we select the case of a random function of case 3, c 3 .f3 , for this section.The design of c 3 .f3 is synthesized on Xilinx Xc2vp30-7ff896 FPGA device (Xilinx, San Jose, CA, USA). Figure 3 shows the area variation for c 3 .f3 with increasing target precision.It can be seen that the area, which is calculated by AASEE, is less than that by AATRE and AACHA, and the area difference between them is increasing with the target precision.This difference is from 265 to 729 with the target precision increased.Such optimization of integer word length can save area.
Figure 4 shows the percentage area saving of AASEE over AATRE at different target precision for c 3 .f3 .The percentage area saving is from 14.34% to 5.62% with the target precision increased.Generally, we obtain increased relative saving for lower precision.

Conclusions
This paper presents a novel affine approximation method for multiplication, Approximation Affine based on Space Extreme Estimation.In this method, an extra noise symbol is added to an approximated affine form.
To reduce the uncertainty in AA, we derive this method in the (n + 1)-dimensional space E n+1 .In space E n+1 , approximate affine form can be regarded as the tangent hyperplane at a certain point of (n + 1)-dimensional space

Figure 4
Figure 4 Percentage area saving of AASEE over AATRE at different target precision for c 3 .f3 .

Table 2 Comparison of analytical ranges and bits by the three methods
c.fun, case study and the function of the row.

Table 3 Comparison of range ratios and computational complexity by the three methods
.fun, case study and the function of the row.http://asp.eurasipjournals.com/content/2014/1/36 c