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 "Mixer_private.h"
23 #include "LVM_Macros.h"
24
25 /**********************************************************************************
26 FUNCTION CORE_MIXSOFT_1ST_D32C31_WRA
27 ***********************************************************************************/
28
29 #ifdef BUILD_FLOAT /* BUILD_FLOAT */
Core_MixInSoft_D32C31_SAT(Mix_1St_Cll_FLOAT_t * pInstance,const LVM_FLOAT * src,LVM_FLOAT * dst,LVM_INT16 n)30 void Core_MixInSoft_D32C31_SAT( Mix_1St_Cll_FLOAT_t *pInstance,
31 const LVM_FLOAT *src,
32 LVM_FLOAT *dst,
33 LVM_INT16 n)
34 {
35 LVM_FLOAT Temp1,Temp2,Temp3;
36 LVM_INT16 OutLoop;
37 LVM_INT16 InLoop;
38 LVM_FLOAT TargetTimesOneMinAlpha;
39 LVM_FLOAT CurrentTimesAlpha;
40 LVM_INT16 ii,jj;
41
42
43 InLoop = (LVM_INT16)(n >> 2); /* Process per 4 samples */
44 OutLoop = (LVM_INT16)(n - (InLoop << 2));
45
46 TargetTimesOneMinAlpha = ((1.0f -pInstance->Alpha) * pInstance->Target);
47 if (pInstance->Target >= pInstance->Current){
48 TargetTimesOneMinAlpha +=(LVM_FLOAT)(2.0f / 2147483647.0f); /* Ceil*/
49 }
50
51 if (OutLoop){
52
53 CurrentTimesAlpha = pInstance->Current * pInstance->Alpha;
54 pInstance->Current = TargetTimesOneMinAlpha + CurrentTimesAlpha;
55
56 for (ii = OutLoop; ii != 0; ii--){
57 Temp1 = *src++;
58 Temp2 = *dst;
59
60 Temp3 = Temp1 * (pInstance->Current);
61 Temp1 = Temp2 + Temp3;
62
63 if (Temp1 > 1.0f)
64 Temp1 = 1.0f;
65 else if (Temp1 < -1.0f)
66 Temp1 = -1.0f;
67
68 *dst++ = Temp1;
69 }
70 }
71
72 for (ii = InLoop; ii != 0; ii--){
73
74 CurrentTimesAlpha = pInstance->Current * pInstance->Alpha;
75 pInstance->Current = TargetTimesOneMinAlpha + CurrentTimesAlpha;
76
77 for (jj = 4; jj!=0 ; jj--){
78 Temp1 = *src++;
79 Temp2 = *dst;
80
81 Temp3 = Temp1 * (pInstance->Current);
82 Temp1 = Temp2 + Temp3;
83
84 if (Temp1 > 1.0f)
85 Temp1 = 1.0f;
86 else if (Temp1 < -1.0f)
87 Temp1 = -1.0f;
88 *dst++ = Temp1;
89 }
90 }
91 }
92 #else
Core_MixInSoft_D32C31_SAT(Mix_1St_Cll_t * pInstance,const LVM_INT32 * src,LVM_INT32 * dst,LVM_INT16 n)93 void Core_MixInSoft_D32C31_SAT( Mix_1St_Cll_t *pInstance,
94 const LVM_INT32 *src,
95 LVM_INT32 *dst,
96 LVM_INT16 n)
97 {
98 LVM_INT32 Temp1,Temp2,Temp3;
99 LVM_INT16 OutLoop;
100 LVM_INT16 InLoop;
101 LVM_INT32 TargetTimesOneMinAlpha;
102 LVM_INT32 CurrentTimesAlpha;
103 LVM_INT16 ii,jj;
104 LVM_INT16 CurrentShort;
105
106 InLoop = (LVM_INT16)(n >> 2); /* Process per 4 samples */
107 OutLoop = (LVM_INT16)(n - (InLoop << 2));
108
109 MUL32x32INTO32((0x7FFFFFFF-pInstance->Alpha),pInstance->Target,TargetTimesOneMinAlpha,31); /* Q31 * Q0 in Q0 */
110 if (pInstance->Target >= pInstance->Current){
111 TargetTimesOneMinAlpha +=2; /* Ceil*/
112 }
113
114 if (OutLoop){
115 MUL32x32INTO32(pInstance->Current,pInstance->Alpha,CurrentTimesAlpha,31); /* Q0 * Q31 in Q0 */
116 pInstance->Current = TargetTimesOneMinAlpha + CurrentTimesAlpha; /* Q0 + Q0 into Q0*/
117 CurrentShort = (LVM_INT16)(pInstance->Current>>16); /* From Q31 to Q15*/
118
119 for (ii = OutLoop; ii != 0; ii--){
120 Temp1=*src++;
121 Temp2=*dst;
122 MUL32x16INTO32(Temp1,CurrentShort,Temp3,15)
123 Temp1=(Temp2>>1)+(Temp3>>1);
124
125 if (Temp1 > 0x3FFFFFFF)
126 Temp1 = 0x7FFFFFFF;
127 else if (Temp1 < - 0x40000000)
128 Temp1 = 0x80000000;
129 else
130 Temp1=(Temp1<<1);
131 *dst++ = Temp1;
132 }
133 }
134
135 for (ii = InLoop; ii != 0; ii--){
136 MUL32x32INTO32(pInstance->Current,pInstance->Alpha,CurrentTimesAlpha,31); /* Q0 * Q31 in Q0 */
137 pInstance->Current = TargetTimesOneMinAlpha + CurrentTimesAlpha; /* Q0 + Q0 into Q0*/
138 CurrentShort = (LVM_INT16)(pInstance->Current>>16); /* From Q31 to Q15*/
139
140 for (jj = 4; jj!=0 ; jj--){
141 Temp1=*src++;
142 Temp2=*dst;
143 MUL32x16INTO32(Temp1,CurrentShort,Temp3,15)
144 Temp1=(Temp2>>1)+(Temp3>>1);
145
146 if (Temp1 > 0x3FFFFFFF)
147 Temp1 = 0x7FFFFFFF;
148 else if (Temp1 < - 0x40000000)
149 Temp1 = 0x80000000;
150 else
151 Temp1=(Temp1<<1);
152 *dst++ = Temp1;
153 }
154 }
155 }
156 #endif
157 /**********************************************************************************/
158