1 /* ----------------------------------------------------------------------
9  * Target Processor: Cortex-M and Cortex-A cores
10  * -------------------------------------------------------------------- */
12  * Copyright (C) 2010-2021 ARM Limited or its affiliates. All rights reserved.
14  * SPDX-License-Identifier: Apache-2.0
20  * www.apache.org/licenses/LICENSE-2.0
81     q15_t m;  in arm_mfcc_q15()  local
95     arm_absmax_q15(pSrc,S->fftLen,&m,&index);  in arm_mfcc_q15()
97     if ((m != 0) && (m != 0x7FFF))  in arm_mfcc_q15()
102        status = arm_divide_q15(0x7FFF,m,&quotient,&shift);  in arm_mfcc_q15()
108        arm_scale_q15(pSrc,quotient,shift,pSrc,S->fftLen);  in arm_mfcc_q15()
113     arm_mult_q15(pSrc,S->windowCoefs, pSrc, S->fftLen);  in arm_mfcc_q15()
118     fftShift = 31 - __CLZ(S->fftLen);  in arm_mfcc_q15()
120     /* some HW accelerator for CMSIS-DSP used in some boards  in arm_mfcc_q15()
128     for(i=0; i < S->fftLen ; i++)  in arm_mfcc_q15()
133     arm_cfft_q15(&(S->cfft),pTmp2,0,1);  in arm_mfcc_q15()
136     arm_rfft_q15(&(S->rfft),pSrc,pTmp2);  in arm_mfcc_q15()
138     filterLimit = 1 + (S->fftLen >> 1);  in arm_mfcc_q15()
141     // q15 - fftShift  in arm_mfcc_q15()
143     // q14 - fftShift  in arm_mfcc_q15()
147     for(i=0; i<S->nbMelFilters; i++)  in arm_mfcc_q15()
149       arm_dot_prod_q15(pSrc+S->filterPos[i],  in arm_mfcc_q15()
150         &(S->filterCoefs[coefsPos]),  in arm_mfcc_q15()
151         S->filterLengths[i],  in arm_mfcc_q15()
154       coefsPos += S->filterLengths[i];  in arm_mfcc_q15()
156       // q34.29 - fftShift  in arm_mfcc_q15()
159       // q34.29 - fftShift - satShift  in arm_mfcc_q15()
164     if ((m != 0) && (m != 0x7FFF))  in arm_mfcc_q15()
166       arm_scale_q31(pTmp,m<<16,0,pTmp,S->nbMelFilters);  in arm_mfcc_q15()
169     // q34.29 - fftShift - satShift  in arm_mfcc_q15()
171     arm_vlog_q31(pTmp,pTmp,S->nbMelFilters);  in arm_mfcc_q15()
181     arm_offset_q31(pTmp,logExponent,pTmp,S->nbMelFilters);  in arm_mfcc_q15()
182     arm_shift_q31(pTmp,-19,pTmp,S->nbMelFilters);  in arm_mfcc_q15()
183     for(i=0; i<S->nbMelFilters; i++)  in arm_mfcc_q15()
190     pDctMat.numRows=S->nbDctOutputs;  in arm_mfcc_q15()
191     pDctMat.numCols=S->nbMelFilters;  in arm_mfcc_q15()
192     pDctMat.pData=(q15_t*)S->dctCoefs;  in arm_mfcc_q15()