forked from ARM-software/CMSIS_5
-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
CMSIS-DSP: Moved f16 bitreversal code into a new file
- Loading branch information
1 parent
63be1c0
commit 737e871
Showing
4 changed files
with
106 additions
and
75 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
100 changes: 100 additions & 0 deletions
100
CMSIS/DSP/Source/TransformFunctions/arm_bitreversal_f16.c
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,100 @@ | ||
| /* ---------------------------------------------------------------------- | ||
| * Project: CMSIS DSP Library | ||
| * Title: arm_bitreversal_f16.c | ||
| * Description: Bitreversal functions | ||
| * | ||
| * | ||
| * Target Processor: Cortex-M cores | ||
| * -------------------------------------------------------------------- */ | ||
| /* | ||
| * Copyright (C) 2010-2021 ARM Limited or its affiliates. All rights reserved. | ||
| * | ||
| * SPDX-License-Identifier: Apache-2.0 | ||
| * | ||
| * Licensed under the Apache License, Version 2.0 (the License); you may | ||
| * not use this file except in compliance with the License. | ||
| * You may obtain a copy of the License at | ||
| * | ||
| * www.apache.org/licenses/LICENSE-2.0 | ||
| * | ||
| * Unless required by applicable law or agreed to in writing, software | ||
| * distributed under the License is distributed on an AS IS BASIS, WITHOUT | ||
| * WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. | ||
| * See the License for the specific language governing permissions and | ||
| * limitations under the License. | ||
| */ | ||
|
|
||
| #include "dsp/transform_functions_f16.h" | ||
|
|
||
| /* | ||
| * @brief In-place bit reversal function. | ||
| * @param[in, out] *pSrc points to the in-place buffer of floating-point data type. | ||
| * @param[in] fftSize length of the FFT. | ||
| * @param[in] bitRevFactor bit reversal modifier that supports different size FFTs with the same bit reversal table. | ||
| * @param[in] *pBitRevTab points to the bit reversal table. | ||
| * @return none. | ||
| */ | ||
|
|
||
| #if defined(ARM_FLOAT16_SUPPORTED) | ||
|
|
||
| void arm_bitreversal_f16( | ||
| float16_t * pSrc, | ||
| uint16_t fftSize, | ||
| uint16_t bitRevFactor, | ||
| uint16_t * pBitRevTab) | ||
| { | ||
| uint16_t fftLenBy2, fftLenBy2p1; | ||
| uint16_t i, j; | ||
| float16_t in; | ||
|
|
||
| /* Initializations */ | ||
| j = 0U; | ||
| fftLenBy2 = fftSize >> 1U; | ||
| fftLenBy2p1 = (fftSize >> 1U) + 1U; | ||
|
|
||
| /* Bit Reversal Implementation */ | ||
| for (i = 0U; i <= (fftLenBy2 - 2U); i += 2U) | ||
| { | ||
| if (i < j) | ||
| { | ||
| /* pSrc[i] <-> pSrc[j]; */ | ||
| in = pSrc[2U * i]; | ||
| pSrc[2U * i] = pSrc[2U * j]; | ||
| pSrc[2U * j] = in; | ||
|
|
||
| /* pSrc[i+1U] <-> pSrc[j+1U] */ | ||
| in = pSrc[(2U * i) + 1U]; | ||
| pSrc[(2U * i) + 1U] = pSrc[(2U * j) + 1U]; | ||
| pSrc[(2U * j) + 1U] = in; | ||
|
|
||
| /* pSrc[i+fftLenBy2p1] <-> pSrc[j+fftLenBy2p1] */ | ||
| in = pSrc[2U * (i + fftLenBy2p1)]; | ||
| pSrc[2U * (i + fftLenBy2p1)] = pSrc[2U * (j + fftLenBy2p1)]; | ||
| pSrc[2U * (j + fftLenBy2p1)] = in; | ||
|
|
||
| /* pSrc[i+fftLenBy2p1+1U] <-> pSrc[j+fftLenBy2p1+1U] */ | ||
| in = pSrc[(2U * (i + fftLenBy2p1)) + 1U]; | ||
| pSrc[(2U * (i + fftLenBy2p1)) + 1U] = | ||
| pSrc[(2U * (j + fftLenBy2p1)) + 1U]; | ||
| pSrc[(2U * (j + fftLenBy2p1)) + 1U] = in; | ||
|
|
||
| } | ||
|
|
||
| /* pSrc[i+1U] <-> pSrc[j+1U] */ | ||
| in = pSrc[2U * (i + 1U)]; | ||
| pSrc[2U * (i + 1U)] = pSrc[2U * (j + fftLenBy2)]; | ||
| pSrc[2U * (j + fftLenBy2)] = in; | ||
|
|
||
| /* pSrc[i+2U] <-> pSrc[j+2U] */ | ||
| in = pSrc[(2U * (i + 1U)) + 1U]; | ||
| pSrc[(2U * (i + 1U)) + 1U] = pSrc[(2U * (j + fftLenBy2)) + 1U]; | ||
| pSrc[(2U * (j + fftLenBy2)) + 1U] = in; | ||
|
|
||
| /* Reading the index for the bit reversal */ | ||
| j = *pBitRevTab; | ||
|
|
||
| /* Updating the bit reversal index depending on the fft length */ | ||
| pBitRevTab += bitRevFactor; | ||
| } | ||
| } | ||
| #endif /* #if defined(ARM_FLOAT16_SUPPORTED) */ |