xref: /aosp_15_r20/external/libxaac/encoder/ixheaace_basic_ops.c (revision 15dc779a375ca8b5125643b829a8aa4b70d7f451)
1 /******************************************************************************
2  *                                                                            *
3  * Copyright (C) 2023 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  * Originally developed and contributed by Ittiam Systems Pvt. Ltd, Bangalore
19  */
20 
21 #include <float.h>
22 #include <math.h>
23 #include "ixheaac_type_def.h"
24 #include "ixheaac_constants.h"
25 #include "ixheaace_aac_constants.h"
26 #include "ixheaac_basic_ops32.h"
27 #include "ixheaac_basic_ops16.h"
28 #include "ixheaac_basic_ops40.h"
29 #include "ixheaac_basic_ops.h"
30 
ia_enhaacplus_enc_norm32_arr(const WORD32 * word32_arr,LOOPINDEX n)31 WORD ia_enhaacplus_enc_norm32_arr(const WORD32 *word32_arr, LOOPINDEX n) {
32   WORD32 max_bits = 0;
33   for (; n != 0; n--) {
34     max_bits = max_bits | ixheaac_abs32_sat(*word32_arr++);
35   }
36   return (ixheaac_pnorm32(max_bits));
37 }
38 
ixheaace_div32(FLOAT32 num,FLOAT32 den)39 FLOAT32 ixheaace_div32(FLOAT32 num, FLOAT32 den) {
40   if (fabs(den) < FLT_EPSILON) {
41     if (den < 0.0f) {
42       return -num;
43     }
44     else {
45       return num;
46     }
47   }
48   else {
49     return num / den;
50   }
51 }
52 
ixheaace_div64(FLOAT64 num,FLOAT64 den)53 FLOAT64 ixheaace_div64(FLOAT64 num, FLOAT64 den) {
54   if (fabs(den) < FLT_EPSILON) {
55     if (den < 0.0) {
56       return -num;
57     }
58     else {
59       return num;
60     }
61   }
62   else {
63     return num / den;
64   }
65 }
66