diff options
| author | Jun Wako <wakojun@gmail.com> | 2015-04-24 16:26:14 +0900 | 
|---|---|---|
| committer | Jun Wako <wakojun@gmail.com> | 2015-04-24 16:26:14 +0900 | 
| commit | 1fe4406f374291ab2e86e95a97341fd9c475fcb8 (patch) | |
| tree | 1be0e16b4b07b5a31ea97ec50a9eb13a288c3d27 /tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions | |
| parent | a20ef7052c6e937d2f7672dd59456e55a5c08296 (diff) | |
Squashed 'tmk_core/' changes from 7967731..b9e0ea0
b9e0ea0 Merge commit '7fa9d8bdea3773d1195b04d98fcf27cf48ddd81d' as 'tool/mbed/mbed-sdk'
7fa9d8b Squashed 'tool/mbed/mbed-sdk/' content from commit 7c21ce5
git-subtree-dir: tmk_core
git-subtree-split: b9e0ea08cb940de20b3610ecdda18e9d8cd7c552
Diffstat (limited to 'tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions')
21 files changed, 3344 insertions, 0 deletions
| diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_copy_f32.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_copy_f32.c new file mode 100644 index 0000000000..f50cb532fd --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_copy_f32.c @@ -0,0 +1,135 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_copy_f32.c     +*     +* Description:	Copies the elements of a floating-point vector.   +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.   +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @defgroup copy Vector Copy     + *     + * Copies sample by sample from source vector to destination vector.     + *     + * <pre>     + * 	pDst[n] = pSrc[n];   0 <= n < blockSize.     + * </pre>     + *    + * There are separate functions for floating point, Q31, Q15, and Q7 data types.      + */ + +/**     + * @addtogroup copy     + * @{     + */ + +/**     + * @brief Copies the elements of a floating-point vector.      + * @param[in]       *pSrc points to input vector     + * @param[out]      *pDst points to output vector     + * @param[in]       blockSize length of the input vector    + * @return none.     + *     + */ + + +void arm_copy_f32( +  float32_t * pSrc, +  float32_t * pDst, +  uint32_t blockSize) +{ +  uint32_t blkCnt;                               /* loop counter */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ +  float32_t in1, in2, in3, in4; + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = A */ +    /* Copy and then store the results in the destination buffer */ +    in1 = *pSrc++; +    in2 = *pSrc++; +    in3 = *pSrc++; +    in4 = *pSrc++; + +    *pDst++ = in1; +    *pDst++ = in2; +    *pDst++ = in3; +    *pDst++ = in4; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = A */ +    /* Copy and then store the results in the destination buffer */ +    *pDst++ = *pSrc++; + +    /* Decrement the loop counter */ +    blkCnt--; +  } +} + +/**     + * @} end of BasicCopy group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_copy_q15.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_copy_q15.c new file mode 100644 index 0000000000..b60e68ac14 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_copy_q15.c @@ -0,0 +1,114 @@ +/* ----------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_copy_q15.c     +*     +* Description:	Copies the elements of a Q15 vector.    +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.    +* -------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup copy     + * @{     + */ +/**     + * @brief Copies the elements of a Q15 vector.      + * @param[in]       *pSrc points to input vector     + * @param[out]      *pDst points to output vector     + * @param[in]       blockSize length of the input vector    + * @return none.     + *     + */ + +void arm_copy_q15( +  q15_t * pSrc, +  q15_t * pDst, +  uint32_t blockSize) +{ +  uint32_t blkCnt;                               /* loop counter */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = A */ +    /* Read two inputs */ +    *__SIMD32(pDst)++ = *__SIMD32(pSrc)++; +    *__SIMD32(pDst)++ = *__SIMD32(pSrc)++; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = A */ +    /* Copy and then store the value in the destination buffer */ +    *pDst++ = *pSrc++; + +    /* Decrement the loop counter */ +    blkCnt--; +  } +} + +/**     + * @} end of BasicCopy group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_copy_q31.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_copy_q31.c new file mode 100644 index 0000000000..3654d3d30e --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_copy_q31.c @@ -0,0 +1,123 @@ +/* ----------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_copy_q31.c     +*     +* Description:	Copies the elements of a Q31 vector.    +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.    +* -------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup copy     + * @{     + */ + +/**     + * @brief Copies the elements of a Q31 vector.      + * @param[in]       *pSrc points to input vector     + * @param[out]      *pDst points to output vector     + * @param[in]       blockSize length of the input vector    + * @return none.     + *     + */ + +void arm_copy_q31( +  q31_t * pSrc, +  q31_t * pDst, +  uint32_t blockSize) +{ +  uint32_t blkCnt;                               /* loop counter */ + + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ +  q31_t in1, in2, in3, in4; + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = A */ +    /* Copy and then store the values in the destination buffer */ +    in1 = *pSrc++; +    in2 = *pSrc++; +    in3 = *pSrc++; +    in4 = *pSrc++; + +    *pDst++ = in1; +    *pDst++ = in2; +    *pDst++ = in3; +    *pDst++ = in4; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = A */ +    /* Copy and then store the value in the destination buffer */ +    *pDst++ = *pSrc++; + +    /* Decrement the loop counter */ +    blkCnt--; +  } +} + +/**     + * @} end of BasicCopy group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_copy_q7.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_copy_q7.c new file mode 100644 index 0000000000..303286fe5e --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_copy_q7.c @@ -0,0 +1,115 @@ +/* ----------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_copy_q7.c     +*     +* Description:	Copies the elements of a Q7 vector.     +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.     +* -------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup copy     + * @{     + */ + +/**     + * @brief Copies the elements of a Q7 vector.     + * @param[in]       *pSrc points to input vector     + * @param[out]      *pDst points to output vector     + * @param[in]       blockSize length of the input vector    + * @return none.     + *     + */ + +void arm_copy_q7( +  q7_t * pSrc, +  q7_t * pDst, +  uint32_t blockSize) +{ +  uint32_t blkCnt;                               /* loop counter */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = A */ +    /* Copy and then store the results in the destination buffer */ +    /* 4 samples are copied and stored at a time using SIMD */ +    *__SIMD32(pDst)++ = *__SIMD32(pSrc)++; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + + +  while(blkCnt > 0u) +  { +    /* C = A */ +    /* Copy and then store the results in the destination buffer */ +    *pDst++ = *pSrc++; + +    /* Decrement the loop counter */ +    blkCnt--; +  } +} + +/**     + * @} end of BasicCopy group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_fill_f32.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_fill_f32.c new file mode 100644 index 0000000000..3f5f86e0e5 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_fill_f32.c @@ -0,0 +1,134 @@ +/* ----------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_fill_f32.c     +*     +* Description:	Fills a constant value into a floating-point vector.    +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.     +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @defgroup Fill Vector Fill     + *     + * Fills the destination vector with a constant value.     + *     + * <pre>     + * 	pDst[n] = value;   0 <= n < blockSize.     + * </pre>     + *    + * There are separate functions for floating point, Q31, Q15, and Q7 data types.      + */ + +/**     + * @addtogroup Fill     + * @{     + */ + +/**     + * @brief Fills a constant value into a floating-point vector.      + * @param[in]       value input value to be filled    + * @param[out]      *pDst points to output vector     + * @param[in]       blockSize length of the output vector    + * @return none.     + *     + */ + + +void arm_fill_f32( +  float32_t value, +  float32_t * pDst, +  uint32_t blockSize) +{ +  uint32_t blkCnt;                               /* loop counter */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ +  float32_t in1 = value; +  float32_t in2 = value; +  float32_t in3 = value; +  float32_t in4 = value; + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = value */ +    /* Fill the value in the destination buffer */ +    *pDst++ = in1; +    *pDst++ = in2; +    *pDst++ = in3; +    *pDst++ = in4; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + + +  while(blkCnt > 0u) +  { +    /* C = value */ +    /* Fill the value in the destination buffer */ +    *pDst++ = value; + +    /* Decrement the loop counter */ +    blkCnt--; +  } +} + +/**     + * @} end of Fill group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_fill_q15.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_fill_q15.c new file mode 100644 index 0000000000..5c73cf627a --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_fill_q15.c @@ -0,0 +1,120 @@ +/* ----------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_fill_q15.c     +*     +* Description:	Fills a constant value into a Q15 vector.    +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.     +* -------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup Fill     + * @{     + */ + +/**     + * @brief Fills a constant value into a Q15 vector.     + * @param[in]       value input value to be filled    + * @param[out]      *pDst points to output vector     + * @param[in]       blockSize length of the output vector    + * @return none.     + *     + */ + +void arm_fill_q15( +  q15_t value, +  q15_t * pDst, +  uint32_t blockSize) +{ +  uint32_t blkCnt;                               /* loop counter */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ + +  q31_t packedValue;                             /* value packed to 32 bits */ + + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* Packing two 16 bit values to 32 bit value in order to use SIMD */ +  packedValue = __PKHBT(value, value, 16u); + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = value */ +    /* Fill the value in the destination buffer */ +    *__SIMD32(pDst)++ = packedValue; +    *__SIMD32(pDst)++ = packedValue; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = value */ +    /* Fill the value in the destination buffer */ +    *pDst++ = value; + +    /* Decrement the loop counter */ +    blkCnt--; +  } +} + +/**     + * @} end of Fill group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_fill_q31.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_fill_q31.c new file mode 100644 index 0000000000..2e8c133ae3 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_fill_q31.c @@ -0,0 +1,121 @@ +/* ----------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_fill_q31.c     +*     +* Description:	Fills a constant value into a Q31 vector.    +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.    +* -------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup Fill     + * @{     + */ + +/**     + * @brief Fills a constant value into a Q31 vector.     + * @param[in]       value input value to be filled    + * @param[out]      *pDst points to output vector     + * @param[in]       blockSize length of the output vector    + * @return none.     + *     + */ + +void arm_fill_q31( +  q31_t value, +  q31_t * pDst, +  uint32_t blockSize) +{ +  uint32_t blkCnt;                               /* loop counter */ + + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ +  q31_t in1 = value; +  q31_t in2 = value; +  q31_t in3 = value; +  q31_t in4 = value; + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = value */ +    /* Fill the value in the destination buffer */ +    *pDst++ = in1; +    *pDst++ = in2; +    *pDst++ = in3; +    *pDst++ = in4; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = value */ +    /* Fill the value in the destination buffer */ +    *pDst++ = value; + +    /* Decrement the loop counter */ +    blkCnt--; +  } +} + +/**     + * @} end of Fill group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_fill_q7.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_fill_q7.c new file mode 100644 index 0000000000..376b7a588b --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_fill_q7.c @@ -0,0 +1,118 @@ +/* ----------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_fill_q7.c     +*     +* Description:	Fills a constant value into a Q7 vector.    +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.     +* -------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup Fill     + * @{     + */ + +/**     + * @brief Fills a constant value into a Q7 vector.     + * @param[in]       value input value to be filled    + * @param[out]      *pDst points to output vector     + * @param[in]       blockSize length of the output vector    + * @return none.     + *     + */ + +void arm_fill_q7( +  q7_t value, +  q7_t * pDst, +  uint32_t blockSize) +{ +  uint32_t blkCnt;                               /* loop counter */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ + +  q31_t packedValue;                             /* value packed to 32 bits */ + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* Packing four 8 bit values to 32 bit value in order to use SIMD */ +  packedValue = __PACKq7(value, value, value, value); + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = value */ +    /* Fill the value in the destination buffer */ +    *__SIMD32(pDst)++ = packedValue; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = value */ +    /* Fill the value in the destination buffer */ +    *pDst++ = value; + +    /* Decrement the loop counter */ +    blkCnt--; +  } +} + +/**     + * @} end of Fill group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_float_to_q15.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_float_to_q15.c new file mode 100644 index 0000000000..cfa5ec6512 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_float_to_q15.c @@ -0,0 +1,204 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_float_to_q15.c     +*     +* Description:	Converts the elements of the floating-point vector to Q15 vector.     +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.     +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup float_to_x     + * @{     + */ + +/**     + * @brief Converts the elements of the floating-point vector to Q15 vector.     + * @param[in]       *pSrc points to the floating-point input vector     + * @param[out]      *pDst points to the Q15 output vector    + * @param[in]       blockSize length of the input vector     + * @return none.     + *     + * \par Description:     + * \par    + * The equation used for the conversion process is:     + * <pre>     + * 	pDst[n] = (q15_t)(pSrc[n] * 32768);   0 <= n < blockSize.     + * </pre>     + * \par Scaling and Overflow Behavior:     + * \par     + * The function uses saturating arithmetic.     + * Results outside of the allowable Q15 range [0x8000 0x7FFF] will be saturated.     + * \note    + * In order to apply rounding, the library should be rebuilt with the ROUNDING macro      + * defined in the preprocessor section of project options.      + *     + */ + + +void arm_float_to_q15( +  float32_t * pSrc, +  q15_t * pDst, +  uint32_t blockSize) +{ +  float32_t *pIn = pSrc;                         /* Src pointer */ +  uint32_t blkCnt;                               /* loop counter */ + +#ifdef ARM_MATH_ROUNDING + +  float32_t in; + +#endif /*      #ifdef ARM_MATH_ROUNDING        */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { + +#ifdef ARM_MATH_ROUNDING +    /* C = A * 32768 */ +    /* convert from float to q15 and then store the results in the destination buffer */ +    in = *pIn++; +    in = (in * 32768.0f); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = (q15_t) (__SSAT((q31_t) (in), 16)); + +    in = *pIn++; +    in = (in * 32768.0f); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = (q15_t) (__SSAT((q31_t) (in), 16)); + +    in = *pIn++; +    in = (in * 32768.0f); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = (q15_t) (__SSAT((q31_t) (in), 16)); + +    in = *pIn++; +    in = (in * 32768.0f); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = (q15_t) (__SSAT((q31_t) (in), 16)); + +#else + +    /* C = A * 32768 */ +    /* convert from float to q15 and then store the results in the destination buffer */ +    *pDst++ = (q15_t) __SSAT((q31_t) (*pIn++ * 32768.0f), 16); +    *pDst++ = (q15_t) __SSAT((q31_t) (*pIn++ * 32768.0f), 16); +    *pDst++ = (q15_t) __SSAT((q31_t) (*pIn++ * 32768.0f), 16); +    *pDst++ = (q15_t) __SSAT((q31_t) (*pIn++ * 32768.0f), 16); + +#endif /*      #ifdef ARM_MATH_ROUNDING        */ + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +  while(blkCnt > 0u) +  { + +#ifdef ARM_MATH_ROUNDING +    /* C = A * 32768 */ +    /* convert from float to q15 and then store the results in the destination buffer */ +    in = *pIn++; +    in = (in * 32768.0f); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = (q15_t) (__SSAT((q31_t) (in), 16)); + +#else + +    /* C = A * 32768 */ +    /* convert from float to q15 and then store the results in the destination buffer */ +    *pDst++ = (q15_t) __SSAT((q31_t) (*pIn++ * 32768.0f), 16); + +#endif /*      #ifdef ARM_MATH_ROUNDING        */ + +    /* Decrement the loop counter */ +    blkCnt--; +  } + + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +  while(blkCnt > 0u) +  { + +#ifdef ARM_MATH_ROUNDING +    /* C = A * 32768 */ +    /* convert from float to q15 and then store the results in the destination buffer */ +    in = *pIn++; +    in = (in * 32768.0f); +    in += in > 0 ? 0.5f : -0.5f; +    *pDst++ = (q15_t) (__SSAT((q31_t) (in), 16)); + +#else + +    /* C = A * 32768 */ +    /* convert from float to q15 and then store the results in the destination buffer */ +    *pDst++ = (q15_t) __SSAT((q31_t) (*pIn++ * 32768.0f), 16); + +#endif /*      #ifdef ARM_MATH_ROUNDING        */ + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +} + +/**     + * @} end of float_to_x group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_float_to_q31.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_float_to_q31.c new file mode 100644 index 0000000000..a39fbe7e39 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_float_to_q31.c @@ -0,0 +1,211 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_float_to_q31.c     +*     +* Description:	Converts the elements of the floating-point vector to Q31 vector.     +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.   +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @defgroup float_to_x  Convert 32-bit floating point value     + */ + +/**     + * @addtogroup float_to_x     + * @{     + */ + +/**     + * @brief Converts the elements of the floating-point vector to Q31 vector.     + * @param[in]       *pSrc points to the floating-point input vector     + * @param[out]      *pDst points to the Q31 output vector    + * @param[in]       blockSize length of the input vector     + * @return none.     + *     + *\par Description:     + * \par    + * The equation used for the conversion process is:     + *    + * <pre>     + * 	pDst[n] = (q31_t)(pSrc[n] * 2147483648);   0 <= n < blockSize.     + * </pre>     + * <b>Scaling and Overflow Behavior:</b>     + * \par     + * The function uses saturating arithmetic.     + * Results outside of the allowable Q31 range[0x80000000 0x7FFFFFFF] will be saturated.     + *    + * \note In order to apply rounding, the library should be rebuilt with the ROUNDING macro      + * defined in the preprocessor section of project options.      + */ + + +void arm_float_to_q31( +  float32_t * pSrc, +  q31_t * pDst, +  uint32_t blockSize) +{ +  float32_t *pIn = pSrc;                         /* Src pointer */ +  uint32_t blkCnt;                               /* loop counter */ + +#ifdef ARM_MATH_ROUNDING + +  float32_t in; + +#endif /*      #ifdef ARM_MATH_ROUNDING        */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { + +#ifdef ARM_MATH_ROUNDING + +    /* C = A * 32768 */ +    /* convert from float to Q31 and then store the results in the destination buffer */ +    in = *pIn++; +    in = (in * 2147483648.0f); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = clip_q63_to_q31((q63_t) (in)); + +    in = *pIn++; +    in = (in * 2147483648.0f); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = clip_q63_to_q31((q63_t) (in)); + +    in = *pIn++; +    in = (in * 2147483648.0f); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = clip_q63_to_q31((q63_t) (in)); + +    in = *pIn++; +    in = (in * 2147483648.0f); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = clip_q63_to_q31((q63_t) (in)); + +#else + +    /* C = A * 2147483648 */ +    /* convert from float to Q31 and then store the results in the destination buffer */ +    *pDst++ = clip_q63_to_q31((q63_t) (*pIn++ * 2147483648.0f)); +    *pDst++ = clip_q63_to_q31((q63_t) (*pIn++ * 2147483648.0f)); +    *pDst++ = clip_q63_to_q31((q63_t) (*pIn++ * 2147483648.0f)); +    *pDst++ = clip_q63_to_q31((q63_t) (*pIn++ * 2147483648.0f)); + +#endif /*      #ifdef ARM_MATH_ROUNDING        */ + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +  while(blkCnt > 0u) +  { + +#ifdef ARM_MATH_ROUNDING + +    /* C = A * 2147483648 */ +    /* convert from float to Q31 and then store the results in the destination buffer */ +    in = *pIn++; +    in = (in * 2147483648.0f); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = clip_q63_to_q31((q63_t) (in)); + +#else + +    /* C = A * 2147483648 */ +    /* convert from float to Q31 and then store the results in the destination buffer */ +    *pDst++ = clip_q63_to_q31((q63_t) (*pIn++ * 2147483648.0f)); + +#endif /*      #ifdef ARM_MATH_ROUNDING        */ + +    /* Decrement the loop counter */ +    blkCnt--; +  } + + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +  while(blkCnt > 0u) +  { + +#ifdef ARM_MATH_ROUNDING + +    /* C = A * 2147483648 */ +    /* convert from float to Q31 and then store the results in the destination buffer */ +    in = *pIn++; +    in = (in * 2147483648.0f); +    in += in > 0 ? 0.5f : -0.5f; +    *pDst++ = clip_q63_to_q31((q63_t) (in)); + +#else + +    /* C = A * 2147483648 */ +    /* convert from float to Q31 and then store the results in the destination buffer */ +    *pDst++ = clip_q63_to_q31((q63_t) (*pIn++ * 2147483648.0f)); + +#endif /*      #ifdef ARM_MATH_ROUNDING        */ + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +} + +/**     + * @} end of float_to_x group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_float_to_q7.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_float_to_q7.c new file mode 100644 index 0000000000..2820af7e66 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_float_to_q7.c @@ -0,0 +1,203 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_float_to_q7.c     +*     +* Description:	Converts the elements of the floating-point vector to Q7 vector.    +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.   +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup float_to_x     + * @{     + */ + +/**     + * @brief Converts the elements of the floating-point vector to Q7 vector.     + * @param[in]       *pSrc points to the floating-point input vector     + * @param[out]      *pDst points to the Q7 output vector    + * @param[in]       blockSize length of the input vector     + * @return none.     + *     + *\par Description:     + * \par    + * The equation used for the conversion process is:     + * <pre>     + * 	pDst[n] = (q7_t)(pSrc[n] * 128);   0 <= n < blockSize.     + * </pre>     + * \par Scaling and Overflow Behavior:     + * \par     + * The function uses saturating arithmetic.     + * Results outside of the allowable Q7 range [0x80 0x7F] will be saturated.     + * \note    + * In order to apply rounding, the library should be rebuilt with the ROUNDING macro      + * defined in the preprocessor section of project options.      + */ + + +void arm_float_to_q7( +  float32_t * pSrc, +  q7_t * pDst, +  uint32_t blockSize) +{ +  float32_t *pIn = pSrc;                         /* Src pointer */ +  uint32_t blkCnt;                               /* loop counter */ + +#ifdef ARM_MATH_ROUNDING + +  float32_t in; + +#endif /*      #ifdef ARM_MATH_ROUNDING        */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { + +#ifdef ARM_MATH_ROUNDING +    /* C = A * 128 */ +    /* convert from float to q7 and then store the results in the destination buffer */ +    in = *pIn++; +    in = (in * 128); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = (q7_t) (__SSAT((q15_t) (in), 8)); + +    in = *pIn++; +    in = (in * 128); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = (q7_t) (__SSAT((q15_t) (in), 8)); + +    in = *pIn++; +    in = (in * 128); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = (q7_t) (__SSAT((q15_t) (in), 8)); + +    in = *pIn++; +    in = (in * 128); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = (q7_t) (__SSAT((q15_t) (in), 8)); + +#else + +    /* C = A * 128 */ +    /* convert from float to q7 and then store the results in the destination buffer */ +    *pDst++ = __SSAT((q31_t) (*pIn++ * 128.0f), 8); +    *pDst++ = __SSAT((q31_t) (*pIn++ * 128.0f), 8); +    *pDst++ = __SSAT((q31_t) (*pIn++ * 128.0f), 8); +    *pDst++ = __SSAT((q31_t) (*pIn++ * 128.0f), 8); + +#endif /*      #ifdef ARM_MATH_ROUNDING        */ + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +  while(blkCnt > 0u) +  { + +#ifdef ARM_MATH_ROUNDING +    /* C = A * 128 */ +    /* convert from float to q7 and then store the results in the destination buffer */ +    in = *pIn++; +    in = (in * 128); +    in += in > 0 ? 0.5 : -0.5; +    *pDst++ = (q7_t) (__SSAT((q15_t) (in), 8)); + +#else + +    /* C = A * 128 */ +    /* convert from float to q7 and then store the results in the destination buffer */ +    *pDst++ = __SSAT((q31_t) (*pIn++ * 128.0f), 8); + +#endif /*      #ifdef ARM_MATH_ROUNDING        */ + +    /* Decrement the loop counter */ +    blkCnt--; +  } + + +#else + +  /* Run the below code for Cortex-M0 */ + + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +  while(blkCnt > 0u) +  { +#ifdef ARM_MATH_ROUNDING +    /* C = A * 128 */ +    /* convert from float to q7 and then store the results in the destination buffer */ +    in = *pIn++; +    in = (in * 128.0f); +    in += in > 0 ? 0.5f : -0.5f; +    *pDst++ = (q7_t) (__SSAT((q31_t) (in), 8)); + +#else + +    /* C = A * 128 */ +    /* convert from float to q7 and then store the results in the destination buffer */ +    *pDst++ = (q7_t) __SSAT((q31_t) (*pIn++ * 128.0f), 8); + +#endif /*      #ifdef ARM_MATH_ROUNDING        */ + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +} + +/**     + * @} end of float_to_x group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q15_to_float.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q15_to_float.c new file mode 100644 index 0000000000..2310b909dc --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q15_to_float.c @@ -0,0 +1,134 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_q15_to_float.c     +*     +* Description:	Converts the elements of the Q15 vector to floating-point vector.      +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.     +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @defgroup q15_to_x  Convert 16-bit Integer value     + */ + +/**     + * @addtogroup q15_to_x     + * @{     + */ + + + + +/**     + * @brief  Converts the elements of the Q15 vector to floating-point vector.      + * @param[in]       *pSrc points to the Q15 input vector     + * @param[out]      *pDst points to the floating-point output vector    + * @param[in]       blockSize length of the input vector     + * @return none.     + *     + * \par Description:     + *     + * The equation used for the conversion process is:     + *    + * <pre>     + * 	pDst[n] = (float32_t) pSrc[n] / 32768;   0 <= n < blockSize.     + * </pre>     + *    + */ + + +void arm_q15_to_float( +  q15_t * pSrc, +  float32_t * pDst, +  uint32_t blockSize) +{ +  q15_t *pIn = pSrc;                             /* Src pointer */ +  uint32_t blkCnt;                               /* loop counter */ + + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = (float32_t) A / 32768 */ +    /* convert from q15 to float and then store the results in the destination buffer */ +    *pDst++ = ((float32_t) * pIn++ / 32768.0f); +    *pDst++ = ((float32_t) * pIn++ / 32768.0f); +    *pDst++ = ((float32_t) * pIn++ / 32768.0f); +    *pDst++ = ((float32_t) * pIn++ / 32768.0f); + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = (float32_t) A / 32768 */ +    /* convert from q15 to float and then store the results in the destination buffer */ +    *pDst++ = ((float32_t) * pIn++ / 32768.0f); + +    /* Decrement the loop counter */ +    blkCnt--; +  } +} + +/**     + * @} end of q15_to_x group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q15_to_q31.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q15_to_q31.c new file mode 100644 index 0000000000..2d5c86e22a --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q15_to_q31.c @@ -0,0 +1,156 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_q15_to_q31.c     +*     +* Description:	Converts the elements of the Q15 vector to Q31 vector.   +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.      +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup q15_to_x     + * @{     + */ + +/**     + * @brief Converts the elements of the Q15 vector to Q31 vector.      + * @param[in]       *pSrc points to the Q15 input vector     + * @param[out]      *pDst points to the Q31 output vector    + * @param[in]       blockSize length of the input vector     + * @return none.     + *     + * \par Description:     + *     + * The equation used for the conversion process is:    + *    + * <pre>     + * 	pDst[n] = (q31_t) pSrc[n] << 16;   0 <= n < blockSize.     + * </pre>     + *    + */ + + +void arm_q15_to_q31( +  q15_t * pSrc, +  q31_t * pDst, +  uint32_t blockSize) +{ +  q15_t *pIn = pSrc;                             /* Src pointer */ +  uint32_t blkCnt;                               /* loop counter */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ +  q31_t in1, in2; +  q31_t out1, out2, out3, out4; + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = (q31_t)A << 16 */ +    /* convert from q15 to q31 and then store the results in the destination buffer */ +    in1 = *__SIMD32(pIn)++; +    in2 = *__SIMD32(pIn)++; + +#ifndef ARM_MATH_BIG_ENDIAN + +    /* extract lower 16 bits to 32 bit result */ +    out1 = in1 << 16u; +    /* extract upper 16 bits to 32 bit result */ +    out2 = in1 & 0xFFFF0000; +    /* extract lower 16 bits to 32 bit result */ +    out3 = in2 << 16u; +    /* extract upper 16 bits to 32 bit result */ +    out4 = in2 & 0xFFFF0000; + +#else + +    /* extract upper 16 bits to 32 bit result */ +    out1 = in1 & 0xFFFF0000; +    /* extract lower 16 bits to 32 bit result */ +    out2 = in1 << 16u; +    /* extract upper 16 bits to 32 bit result */ +    out3 = in2 & 0xFFFF0000; +    /* extract lower 16 bits to 32 bit result */ +    out4 = in2 << 16u; + +#endif //      #ifndef ARM_MATH_BIG_ENDIAN + +    *pDst++ = out1; +    *pDst++ = out2; +    *pDst++ = out3; +    *pDst++ = out4; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = (q31_t)A << 16 */ +    /* convert from q15 to q31 and then store the results in the destination buffer */ +    *pDst++ = (q31_t) * pIn++ << 16; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +} + +/**     + * @} end of q15_to_x group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q15_to_q7.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q15_to_q7.c new file mode 100644 index 0000000000..d261221505 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q15_to_q7.c @@ -0,0 +1,154 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_q15_to_q7.c     +*     +* Description:	Converts the elements of the Q15 vector to Q7 vector.   +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.     +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup q15_to_x     + * @{     + */ + + +/**     + * @brief Converts the elements of the Q15 vector to Q7 vector.      + * @param[in]       *pSrc points to the Q15 input vector     + * @param[out]      *pDst points to the Q7 output vector    + * @param[in]       blockSize length of the input vector     + * @return none.     + *     + * \par Description:     + *     + * The equation used for the conversion process is:     + *    + * <pre>     + * 	pDst[n] = (q7_t) pSrc[n] >> 8;   0 <= n < blockSize.     + * </pre>    + *    + */ + + +void arm_q15_to_q7( +  q15_t * pSrc, +  q7_t * pDst, +  uint32_t blockSize) +{ +  q15_t *pIn = pSrc;                             /* Src pointer */ +  uint32_t blkCnt;                               /* loop counter */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ +  q31_t in1, in2; +  q31_t out1, out2; + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = (q7_t) A >> 8 */ +    /* convert from q15 to q7 and then store the results in the destination buffer */ +    in1 = *__SIMD32(pIn)++; +    in2 = *__SIMD32(pIn)++; + +#ifndef ARM_MATH_BIG_ENDIAN + +    out1 = __PKHTB(in2, in1, 16); +    out2 = __PKHBT(in2, in1, 16); + +#else + +    out1 = __PKHTB(in1, in2, 16); +    out2 = __PKHBT(in1, in2, 16); + +#endif //      #ifndef ARM_MATH_BIG_ENDIAN + +    /* rotate packed value by 24 */ +    out2 = ((uint32_t) out2 << 8) | ((uint32_t) out2 >> 24); + +    /* anding with 0xff00ff00 to get two 8 bit values */ +    out1 = out1 & 0xFF00FF00; +    /* anding with 0x00ff00ff to get two 8 bit values */ +    out2 = out2 & 0x00FF00FF; + +    /* oring two values(contains two 8 bit values) to get four packed 8 bit values */ +    out1 = out1 | out2; + +    /* store 4 samples at a time to destiantion buffer */ +    *__SIMD32(pDst)++ = out1; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = (q7_t) A >> 8 */ +    /* convert from q15 to q7 and then store the results in the destination buffer */ +    *pDst++ = (q7_t) (*pIn++ >> 8); + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +} + +/**     + * @} end of q15_to_x group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q31_to_float.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q31_to_float.c new file mode 100644 index 0000000000..4f60511c22 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q31_to_float.c @@ -0,0 +1,131 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_q31_to_float.c     +*     +* Description:	Converts the elements of the Q31 vector to floating-point vector.       +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.   +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @defgroup q31_to_x  Convert 32-bit Integer value     + */ + +/**     + * @addtogroup q31_to_x     + * @{     + */ + +/**     + * @brief Converts the elements of the Q31 vector to floating-point vector.     + * @param[in]       *pSrc points to the Q31 input vector     + * @param[out]      *pDst points to the floating-point output vector    + * @param[in]       blockSize length of the input vector     + * @return none.     + *     + * \par Description:     + *     + * The equation used for the conversion process is:     + *    + * <pre>     + * 	pDst[n] = (float32_t) pSrc[n] / 2147483648;   0 <= n < blockSize.     + * </pre>     + *    + */ + + +void arm_q31_to_float( +  q31_t * pSrc, +  float32_t * pDst, +  uint32_t blockSize) +{ +  q31_t *pIn = pSrc;                             /* Src pointer */ +  uint32_t blkCnt;                               /* loop counter */ + + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = (float32_t) A / 2147483648 */ +    /* convert from q31 to float and then store the results in the destination buffer */ +    *pDst++ = ((float32_t) * pIn++ / 2147483648.0f); +    *pDst++ = ((float32_t) * pIn++ / 2147483648.0f); +    *pDst++ = ((float32_t) * pIn++ / 2147483648.0f); +    *pDst++ = ((float32_t) * pIn++ / 2147483648.0f); + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = (float32_t) A / 2147483648 */ +    /* convert from q31 to float and then store the results in the destination buffer */ +    *pDst++ = ((float32_t) * pIn++ / 2147483648.0f); + +    /* Decrement the loop counter */ +    blkCnt--; +  } +} + +/**     + * @} end of q31_to_x group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q31_to_q15.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q31_to_q15.c new file mode 100644 index 0000000000..a2b9fde747 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q31_to_q15.c @@ -0,0 +1,145 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_q31_to_q15.c     +*     +* Description:	Converts the elements of the Q31 vector to Q15 vector.     +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.     +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup q31_to_x     + * @{     + */ + +/**     + * @brief Converts the elements of the Q31 vector to Q15 vector.     + * @param[in]       *pSrc points to the Q31 input vector     + * @param[out]      *pDst points to the Q15 output vector    + * @param[in]       blockSize length of the input vector     + * @return none.     + *      + * \par Description:     + *     + * The equation used for the conversion process is:     + *    + * <pre>     + * 	pDst[n] = (q15_t) pSrc[n] >> 16;   0 <= n < blockSize.     + * </pre>     + *    + */ + + +void arm_q31_to_q15( +  q31_t * pSrc, +  q15_t * pDst, +  uint32_t blockSize) +{ +  q31_t *pIn = pSrc;                             /* Src pointer */ +  uint32_t blkCnt;                               /* loop counter */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ +  q31_t in1, in2, in3, in4; +  q31_t out1, out2; + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = (q15_t) A >> 16 */ +    /* convert from q31 to q15 and then store the results in the destination buffer */ +    in1 = *pIn++; +    in2 = *pIn++; +    in3 = *pIn++; +    in4 = *pIn++; + +    /* pack two higher 16-bit values from two 32-bit values */ +#ifndef ARM_MATH_BIG_ENDIAN + +    out1 = __PKHTB(in2, in1, 16); +    out2 = __PKHTB(in4, in3, 16); + +#else + +    out1 = __PKHTB(in1, in2, 16); +    out2 = __PKHTB(in3, in4, 16); + +#endif //      #ifdef ARM_MATH_BIG_ENDIAN + +    *__SIMD32(pDst)++ = out1; +    *__SIMD32(pDst)++ = out2; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = (q15_t) A >> 16 */ +    /* convert from q31 to q15 and then store the results in the destination buffer */ +    *pDst++ = (q15_t) (*pIn++ >> 16); + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +} + +/**     + * @} end of q31_to_x group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q31_to_q7.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q31_to_q7.c new file mode 100644 index 0000000000..c2f9b9a040 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q31_to_q7.c @@ -0,0 +1,136 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_q31_to_q7.c     +*     +* Description:	Converts the elements of the Q31 vector to Q7 vector.     +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.     +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup q31_to_x     + * @{     + */ + +/**     + * @brief Converts the elements of the Q31 vector to Q7 vector.     + * @param[in]       *pSrc points to the Q31 input vector     + * @param[out]      *pDst points to the Q7 output vector    + * @param[in]       blockSize length of the input vector     + * @return none.     + *     + * \par Description:     + *     + * The equation used for the conversion process is:     + *    + * <pre>     + * 	pDst[n] = (q7_t) pSrc[n] >> 24;   0 <= n < blockSize.      + * </pre>     + *    + */ + + +void arm_q31_to_q7( +  q31_t * pSrc, +  q7_t * pDst, +  uint32_t blockSize) +{ +  q31_t *pIn = pSrc;                             /* Src pointer */ +  uint32_t blkCnt;                               /* loop counter */ + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ +  q31_t in1, in2, in3, in4; +  q7_t out1, out2, out3, out4; + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = (q7_t) A >> 24 */ +    /* convert from q31 to q7 and then store the results in the destination buffer */ +    in1 = *pIn++; +    in2 = *pIn++; +    in3 = *pIn++; +    in4 = *pIn++; + +    out1 = (q7_t) (in1 >> 24); +    out2 = (q7_t) (in2 >> 24); +    out3 = (q7_t) (in3 >> 24); +    out4 = (q7_t) (in4 >> 24); + +    *__SIMD32(pDst)++ = __PACKq7(out1, out2, out3, out4); + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = (q7_t) A >> 24 */ +    /* convert from q31 to q7 and then store the results in the destination buffer */ +    *pDst++ = (q7_t) (*pIn++ >> 24); + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +} + +/**     + * @} end of q31_to_x group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q7_to_float.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q7_to_float.c new file mode 100644 index 0000000000..3b7f586a57 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q7_to_float.c @@ -0,0 +1,131 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_q7_to_float.c     +*     +* Description:	Converts the elements of the Q7 vector to floating-point vector.     +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.   +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @defgroup q7_to_x  Convert 8-bit Integer value     + */ + +/**     + * @addtogroup q7_to_x     + * @{     + */ + +/**     + * @brief Converts the elements of the Q7 vector to floating-point vector.     + * @param[in]       *pSrc points to the Q7 input vector     + * @param[out]      *pDst points to the floating-point output vector    + * @param[in]       blockSize length of the input vector     + * @return none.     + *		      + * \par Description:     + *     + * The equation used for the conversion process is:     + *    + * <pre>     + * 	pDst[n] = (float32_t) pSrc[n] / 128;   0 <= n < blockSize.     + * </pre>     + *    + */ + + +void arm_q7_to_float( +  q7_t * pSrc, +  float32_t * pDst, +  uint32_t blockSize) +{ +  q7_t *pIn = pSrc;                              /* Src pointer */ +  uint32_t blkCnt;                               /* loop counter */ + + +#ifndef ARM_MATH_CM0_FAMILY + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = (float32_t) A / 128 */ +    /* convert from q7 to float and then store the results in the destination buffer */ +    *pDst++ = ((float32_t) * pIn++ / 128.0f); +    *pDst++ = ((float32_t) * pIn++ / 128.0f); +    *pDst++ = ((float32_t) * pIn++ / 128.0f); +    *pDst++ = ((float32_t) * pIn++ / 128.0f); + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = (float32_t) A / 128 */ +    /* convert from q7 to float and then store the results in the destination buffer */ +    *pDst++ = ((float32_t) * pIn++ / 128.0f); + +    /* Decrement the loop counter */ +    blkCnt--; +  } +} + +/**     + * @} end of q7_to_x group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q7_to_q15.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q7_to_q15.c new file mode 100644 index 0000000000..444321c601 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q7_to_q15.c @@ -0,0 +1,157 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_q7_to_q15.c     +*     +* Description:	Converts the elements of the Q7 vector to Q15 vector.     +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.    +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup q7_to_x     + * @{     + */ + + + + +/**     + * @brief Converts the elements of the Q7 vector to Q15 vector.     + * @param[in]       *pSrc points to the Q7 input vector     + * @param[out]      *pDst points to the Q15 output vector    + * @param[in]       blockSize length of the input vector     + * @return none.     + *     + * \par Description:     + *     + * The equation used for the conversion process is:     + *    + * <pre>     + * 	pDst[n] = (q15_t) pSrc[n] << 8;   0 <= n < blockSize.     + * </pre>     + *    + */ + + +void arm_q7_to_q15( +  q7_t * pSrc, +  q15_t * pDst, +  uint32_t blockSize) +{ +  q7_t *pIn = pSrc;                              /* Src pointer */ +  uint32_t blkCnt;                               /* loop counter */ + +#ifndef ARM_MATH_CM0_FAMILY +  q31_t in; +  q31_t in1, in2; +  q31_t out1, out2; + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = (q15_t) A << 8 */ +    /* convert from q7 to q15 and then store the results in the destination buffer */ +    in = *__SIMD32(pIn)++; + +    /* rotatate in by 8 and extend two q7_t values to q15_t values */ +    in1 = __SXTB16(__ROR(in, 8)); + +    /* extend remainig two q7_t values to q15_t values */ +    in2 = __SXTB16(in); + +    in1 = in1 << 8u; +    in2 = in2 << 8u; + +    in1 = in1 & 0xFF00FF00; +    in2 = in2 & 0xFF00FF00; + +#ifndef ARM_MATH_BIG_ENDIAN + +    out2 = __PKHTB(in1, in2, 16); +    out1 = __PKHBT(in2, in1, 16); + +#else + +    out1 = __PKHTB(in1, in2, 16); +    out2 = __PKHBT(in2, in1, 16); + +#endif + +    *__SIMD32(pDst)++ = out1; +    *__SIMD32(pDst)++ = out2; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = (q15_t) A << 8 */ +    /* convert from q7 to q15 and then store the results in the destination buffer */ +    *pDst++ = (q15_t) * pIn++ << 8; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +} + +/**     + * @} end of q7_to_x group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q7_to_q31.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q7_to_q31.c new file mode 100644 index 0000000000..fefd78a017 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/arm_q7_to_q31.c @@ -0,0 +1,142 @@ +/* ----------------------------------------------------------------------------     +* Copyright (C) 2010-2013 ARM Limited. All rights reserved.     +*     +* $Date:        17. January 2013  +* $Revision: 	V1.4.1   +*     +* Project: 	    CMSIS DSP Library     +* Title:		arm_q7_to_q31.c     +*     +* Description:	Converts the elements of the Q7 vector to Q31 vector.   +*     +* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.    +* ---------------------------------------------------------------------------- */ + +#include "arm_math.h" + +/**     + * @ingroup groupSupport     + */ + +/**     + * @addtogroup q7_to_x     + * @{     + */ + +/**     + * @brief Converts the elements of the Q7 vector to Q31 vector.     + * @param[in]       *pSrc points to the Q7 input vector     + * @param[out]      *pDst points to the Q31 output vector    + * @param[in]       blockSize length of the input vector     + * @return none.     + *     + * \par Description:     + *     + * The equation used for the conversion process is:     + *    + * <pre>     + * 	pDst[n] = (q31_t) pSrc[n] << 24;   0 <= n < blockSize.    + * </pre>      + *    + */ + + +void arm_q7_to_q31( +  q7_t * pSrc, +  q31_t * pDst, +  uint32_t blockSize) +{ +  q7_t *pIn = pSrc;                              /* Src pointer */ +  uint32_t blkCnt;                               /* loop counter */ + +#ifndef ARM_MATH_CM0_FAMILY + +  q31_t in; + +  /* Run the below code for Cortex-M4 and Cortex-M3 */ + +  /*loop Unrolling */ +  blkCnt = blockSize >> 2u; + +  /* First part of the processing with loop unrolling.  Compute 4 outputs at a time.     +   ** a second loop below computes the remaining 1 to 3 samples. */ +  while(blkCnt > 0u) +  { +    /* C = (q31_t) A << 24 */ +    /* convert from q7 to q31 and then store the results in the destination buffer */ +    in = *__SIMD32(pIn)++; + +#ifndef ARM_MATH_BIG_ENDIAN + +    *pDst++ = (__ROR(in, 8)) & 0xFF000000; +    *pDst++ = (__ROR(in, 16)) & 0xFF000000; +    *pDst++ = (__ROR(in, 24)) & 0xFF000000; +    *pDst++ = (in & 0xFF000000); + +#else + +    *pDst++ = (in & 0xFF000000); +    *pDst++ = (__ROR(in, 24)) & 0xFF000000; +    *pDst++ = (__ROR(in, 16)) & 0xFF000000; +    *pDst++ = (__ROR(in, 8)) & 0xFF000000; + +#endif //              #ifndef ARM_MATH_BIG_ENDIAN + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +  /* If the blockSize is not a multiple of 4, compute any remaining output samples here.     +   ** No loop unrolling is used. */ +  blkCnt = blockSize % 0x4u; + +#else + +  /* Run the below code for Cortex-M0 */ + +  /* Loop over blockSize number of values */ +  blkCnt = blockSize; + +#endif /* #ifndef ARM_MATH_CM0_FAMILY */ + +  while(blkCnt > 0u) +  { +    /* C = (q31_t) A << 24 */ +    /* convert from q7 to q31 and then store the results in the destination buffer */ +    *pDst++ = (q31_t) * pIn++ << 24; + +    /* Decrement the loop counter */ +    blkCnt--; +  } + +} + +/**     + * @} end of q7_to_x group     + */ diff --git a/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/math_helper.c b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/math_helper.c new file mode 100644 index 0000000000..522f5a6762 --- /dev/null +++ b/tool/mbed/mbed-sdk/libraries/dsp/cmsis_dsp/SupportFunctions/math_helper.c @@ -0,0 +1,460 @@ +/* ----------------------------------------------------------------------    +* Copyright (C) 2010-2012 ARM Limited. All rights reserved.    +*    +* $Date:        17. January 2013   +* $Revision: 	V1.4.0     +*   +* Project: 	    CMSIS DSP Library  +* +* Title:	    math_helper.c +* +* Description:	Definition of all helper functions required.   +*   +* Target Processor: Cortex-M4/Cortex-M3 +*   +* Redistribution and use in source and binary forms, with or without  +* modification, are permitted provided that the following conditions +* are met: +*   - Redistributions of source code must retain the above copyright +*     notice, this list of conditions and the following disclaimer. +*   - Redistributions in binary form must reproduce the above copyright +*     notice, this list of conditions and the following disclaimer in +*     the documentation and/or other materials provided with the  +*     distribution. +*   - Neither the name of ARM LIMITED nor the names of its contributors +*     may be used to endorse or promote products derived from this +*     software without specific prior written permission. +* +* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS +* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT +* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS +* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE  +* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, +* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, +* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; +* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER +* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT +* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN +* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE +* POSSIBILITY OF SUCH DAMAGE.   +* -------------------------------------------------------------------- */ + +/* ---------------------------------------------------------------------- +*		Include standard header files   +* -------------------------------------------------------------------- */ +#include<math.h> + +/* ---------------------------------------------------------------------- +*		Include project header files   +* -------------------------------------------------------------------- */ +#include "math_helper.h" + +/**  + * @brief  Caluclation of SNR + * @param  float* 	Pointer to the reference buffer + * @param  float*	Pointer to the test buffer + * @param  uint32_t	total number of samples + * @return float	SNR + * The function Caluclates signal to noise ratio for the reference output  + * and test output  + */ + +float arm_snr_f32(float *pRef, float *pTest, uint32_t buffSize) +{ +  float EnergySignal = 0.0, EnergyError = 0.0; +  uint32_t i; +  float SNR; +  int temp; +  int *test; + +  for (i = 0; i < buffSize; i++) +    { + 	  /* Checking for a NAN value in pRef array */ +	  test =   (int *)(&pRef[i]); +      temp =  *test; + +	  if(temp == 0x7FC00000) +	  { +	  		return(0); +	  } + +	  /* Checking for a NAN value in pTest array */ +	  test =   (int *)(&pTest[i]); +      temp =  *test; + +	  if(temp == 0x7FC00000) +	  { +	  		return(0); +	  } +      EnergySignal += pRef[i] * pRef[i]; +      EnergyError += (pRef[i] - pTest[i]) * (pRef[i] - pTest[i]);  +    } + +	/* Checking for a NAN value in EnergyError */ +	test =   (int *)(&EnergyError); +    temp =  *test; + +    if(temp == 0x7FC00000) +    { +  		return(0); +    } +	 + +  SNR = 10 * log10 (EnergySignal / EnergyError); + +  return (SNR); + +} + + +/**  + * @brief  Provide guard bits for Input buffer + * @param  q15_t* 	    Pointer to input buffer + * @param  uint32_t 	blockSize + * @param  uint32_t 	guard_bits + * @return none + * The function Provides the guard bits for the buffer  + * to avoid overflow  + */ + +void arm_provide_guard_bits_q15 (q15_t * input_buf, uint32_t blockSize, +                            uint32_t guard_bits) +{ +  uint32_t i; + +  for (i = 0; i < blockSize; i++) +    { +      input_buf[i] = input_buf[i] >> guard_bits; +    } +} + +/**  + * @brief  Converts float to fixed in q12.20 format + * @param  uint32_t 	number of samples in the buffer + * @return none + * The function converts floating point values to fixed point(q12.20) values  + */ + +void arm_float_to_q12_20(float *pIn, q31_t * pOut, uint32_t numSamples) +{ +  uint32_t i; + +  for (i = 0; i < numSamples; i++) +    { +	  /* 1048576.0f corresponds to pow(2, 20) */ +      pOut[i] = (q31_t) (pIn[i] * 1048576.0f); + +      pOut[i] += pIn[i] > 0 ? 0.5 : -0.5; + +      if (pIn[i] == (float) 1.0) +        { +          pOut[i] = 0x000FFFFF; +        } +    } +} + +/**  + * @brief  Compare MATLAB Reference Output and ARM Test output + * @param  q15_t* 	Pointer to Ref buffer + * @param  q15_t* 	Pointer to Test buffer + * @param  uint32_t 	number of samples in the buffer + * @return none  + */ + +uint32_t arm_compare_fixed_q15(q15_t *pIn, q15_t * pOut, uint32_t numSamples) +{ +  uint32_t i;  +  int32_t diff; +  uint32_t diffCrnt = 0; +  uint32_t maxDiff = 0; + +  for (i = 0; i < numSamples; i++) +  { +  	diff = pIn[i] - pOut[i]; +  	diffCrnt = (diff > 0) ? diff : -diff; + +	if(diffCrnt > maxDiff) +	{ +		maxDiff = diffCrnt; +	}	 +  } + +  return(maxDiff); +} + +/**  + * @brief  Compare MATLAB Reference Output and ARM Test output + * @param  q31_t* 	Pointer to Ref buffer + * @param  q31_t* 	Pointer to Test buffer + * @param  uint32_t 	number of samples in the buffer + * @return none  + */ + +uint32_t arm_compare_fixed_q31(q31_t *pIn, q31_t * pOut, uint32_t numSamples) +{ +  uint32_t i;  +  int32_t diff; +  uint32_t diffCrnt = 0; +  uint32_t maxDiff = 0; + +  for (i = 0; i < numSamples; i++) +  { +  	diff = pIn[i] - pOut[i]; +  	diffCrnt = (diff > 0) ? diff : -diff; + +	if(diffCrnt > maxDiff) +	{ +		maxDiff = diffCrnt; +	} +  } + +  return(maxDiff); +} + +/**  + * @brief  Provide guard bits for Input buffer + * @param  q31_t* 	Pointer to input buffer + * @param  uint32_t 	blockSize + * @param  uint32_t 	guard_bits + * @return none + * The function Provides the guard bits for the buffer  + * to avoid overflow  + */ + +void arm_provide_guard_bits_q31 (q31_t * input_buf,  +								 uint32_t blockSize, +                                 uint32_t guard_bits) +{ +  uint32_t i; + +  for (i = 0; i < blockSize; i++) +    { +      input_buf[i] = input_buf[i] >> guard_bits; +    } +} + +/**  + * @brief  Provide guard bits for Input buffer + * @param  q31_t* 	Pointer to input buffer + * @param  uint32_t 	blockSize + * @param  uint32_t 	guard_bits + * @return none + * The function Provides the guard bits for the buffer  + * to avoid overflow  + */ + +void arm_provide_guard_bits_q7 (q7_t * input_buf,  +								uint32_t blockSize, +                                uint32_t guard_bits) +{ +  uint32_t i; + +  for (i = 0; i < blockSize; i++) +    { +      input_buf[i] = input_buf[i] >> guard_bits; +    } +} + + + +/**  + * @brief  Caluclates number of guard bits  + * @param  uint32_t 	number of additions + * @return none + * The function Caluclates the number of guard bits   + * depending on the numtaps  + */ + +uint32_t arm_calc_guard_bits (uint32_t num_adds) +{ +  uint32_t i = 1, j = 0; + +  if (num_adds == 1) +    { +      return (0); +    } + +  while (i < num_adds) +    { +      i = i * 2; +      j++; +    } + +  return (j); +} + +/**  + * @brief  Converts Q15 to floating-point + * @param  uint32_t 	number of samples in the buffer + * @return none + */ + +void arm_apply_guard_bits (float32_t * pIn,  +						   uint32_t numSamples,  +						   uint32_t guard_bits) +{ +  uint32_t i; + +  for (i = 0; i < numSamples; i++) +    { +      pIn[i] = pIn[i] * arm_calc_2pow(guard_bits); +    } +} + +/**  + * @brief  Calculates pow(2, numShifts) + * @param  uint32_t 	number of shifts + * @return pow(2, numShifts) + */ +uint32_t arm_calc_2pow(uint32_t numShifts) +{ + +  uint32_t i, val = 1; + +  for (i = 0; i < numShifts; i++) +    { +      val = val * 2; +    }	 + +  return(val); +} + + + +/**  + * @brief  Converts float to fixed q14  + * @param  uint32_t 	number of samples in the buffer + * @return none + * The function converts floating point values to fixed point values  + */ + +void arm_float_to_q14 (float *pIn, q15_t * pOut,  +                       uint32_t numSamples) +{ +  uint32_t i; + +  for (i = 0; i < numSamples; i++) +    { +	  /* 16384.0f corresponds to pow(2, 14) */ +      pOut[i] = (q15_t) (pIn[i] * 16384.0f); + +      pOut[i] += pIn[i] > 0 ? 0.5 : -0.5; + +      if (pIn[i] == (float) 2.0) +        { +          pOut[i] = 0x7FFF; +        } + +    } + +} + +  +/**  + * @brief  Converts float to fixed q30 format + * @param  uint32_t 	number of samples in the buffer + * @return none + * The function converts floating point values to fixed point values  + */ + +void arm_float_to_q30 (float *pIn, q31_t * pOut,  +					   uint32_t numSamples) +{ +  uint32_t i; + +  for (i = 0; i < numSamples; i++) +    { +	  /* 1073741824.0f corresponds to pow(2, 30) */ +      pOut[i] = (q31_t) (pIn[i] * 1073741824.0f); + +      pOut[i] += pIn[i] > 0 ? 0.5 : -0.5; + +      if (pIn[i] == (float) 2.0) +        { +          pOut[i] = 0x7FFFFFFF; +        } +    } +} + +/**  + * @brief  Converts float to fixed q30 format + * @param  uint32_t 	number of samples in the buffer + * @return none + * The function converts floating point values to fixed point values  + */ + +void arm_float_to_q29 (float *pIn, q31_t * pOut,  +					   uint32_t numSamples) +{ +  uint32_t i; + +  for (i = 0; i < numSamples; i++) +    { +	  /* 1073741824.0f corresponds to pow(2, 30) */ +      pOut[i] = (q31_t) (pIn[i] * 536870912.0f); + +      pOut[i] += pIn[i] > 0 ? 0.5 : -0.5; + +      if (pIn[i] == (float) 4.0) +        { +          pOut[i] = 0x7FFFFFFF; +        } +    } +} + + +/**  + * @brief  Converts float to fixed q28 format + * @param  uint32_t 	number of samples in the buffer + * @return none + * The function converts floating point values to fixed point values  + */ + +void arm_float_to_q28 (float *pIn, q31_t * pOut,  +                       uint32_t numSamples) +{ +  uint32_t i; + +  for (i = 0; i < numSamples; i++) +    { +	/* 268435456.0f corresponds to pow(2, 28) */ +      pOut[i] = (q31_t) (pIn[i] * 268435456.0f); + +      pOut[i] += pIn[i] > 0 ? 0.5 : -0.5; + +      if (pIn[i] == (float) 8.0) +        { +          pOut[i] = 0x7FFFFFFF; +        } +    } +} + +/**  + * @brief  Clip the float values to +/- 1  + * @param  pIn 	input buffer + * @param  numSamples 	number of samples in the buffer + * @return none + * The function converts floating point values to fixed point values  + */ + +void arm_clip_f32 (float *pIn, uint32_t numSamples) +{ +  uint32_t i; + +  for (i = 0; i < numSamples; i++) +    { +      if(pIn[i] > 1.0f) +	  { +	    pIn[i] = 1.0; +	  } +	  else if( pIn[i] < -1.0f) +	  { +	    pIn[i] = -1.0; +	  } +	        +    } +} + + + + | 
