1 /*
2  * Copyright (C) 2004-2010 NXP Software
3  * Copyright (C) 2010 The Android Open Source Project
4  *
5  * Licensed under the Apache License, Version 2.0 (the "License");
6  * you may not use this file except in compliance with the License.
7  * You may obtain a copy of the License at
8  *
9  *      http://www.apache.org/licenses/LICENSE-2.0
10  *
11  * Unless required by applicable law or agreed to in writing, software
12  * distributed under the License is distributed on an "AS IS" BASIS,
13  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14  * See the License for the specific language governing permissions and
15  * limitations under the License.
16  */
17 
18 /**********************************************************************************
19    INCLUDE FILES
20 ***********************************************************************************/
21 
22 #include "VectorArithmetic.h"
23 
24 /**********************************************************************************
25    FUNCTION ADD2_SAT_32X32
26 ***********************************************************************************/
27 
Add2_Sat_32x32(const LVM_INT32 * src,LVM_INT32 * dst,LVM_INT16 n)28 void Add2_Sat_32x32( const LVM_INT32  *src,
29                            LVM_INT32  *dst,
30                            LVM_INT16  n )
31 {
32     LVM_INT32 a,b,c;
33     LVM_INT16 ii;
34     for (ii = n; ii != 0; ii--)
35     {
36         a=*src;
37         src++;
38 
39         b=*dst;
40         c=a+b;
41         if ((((c ^ a) & (c ^ b)) >> 31)!=0)     /* overflow / underflow */
42         {
43             if(a<0)
44             {
45                 c=0x80000000L;
46             }
47             else
48             {
49                 c=0x7FFFFFFFL;
50             }
51         }
52 
53         *dst = c;
54         dst++;
55     }
56     return;
57 }
58 
Add2_Sat_Float(const LVM_FLOAT * src,LVM_FLOAT * dst,LVM_INT16 n)59 void Add2_Sat_Float( const LVM_FLOAT  *src,
60                            LVM_FLOAT  *dst,
61                            LVM_INT16  n )
62 {
63     LVM_FLOAT Temp;
64     LVM_INT16 ii;
65     for (ii = n; ii != 0; ii--)
66     {
67         Temp = ((LVM_FLOAT) *src) + ((LVM_FLOAT) *dst);
68         src++;
69 
70         if (Temp > 1.000000f)
71         {
72             *dst = 1.000000f;
73         }
74         else if (Temp < -1.000000f)
75         {
76             *dst = -1.000000f;
77         }
78         else
79         {
80             *dst = Temp;
81         }
82         dst++;
83     }
84     return;
85 }
86 /**********************************************************************************/
87