CMSIS DSP Library from CMSIS 2.0. See http://www.onarm.com/cmsis/ for full details

Dependents:   K22F_DSP_Matrix_least_square BNO055-ELEC3810 1BNO055 ECE4180Project--Slave2 ... more

Committer:
simon
Date:
Thu Mar 10 15:07:50 2011 +0000
Revision:
0:1014af42efd9

        

Who changed what in which revision?

UserRevisionLine numberNew contents of line
simon 0:1014af42efd9 1 /* ----------------------------------------------------------------------
simon 0:1014af42efd9 2 * Copyright (C) 2010 ARM Limited. All rights reserved.
simon 0:1014af42efd9 3 *
simon 0:1014af42efd9 4 * $Date: 29. November 2010
simon 0:1014af42efd9 5 * $Revision: V1.0.3
simon 0:1014af42efd9 6 *
simon 0:1014af42efd9 7 * Project: CMSIS DSP Library
simon 0:1014af42efd9 8 * Title: arm_offset_q7.c
simon 0:1014af42efd9 9 *
simon 0:1014af42efd9 10 * Description: Q7 vector offset.
simon 0:1014af42efd9 11 *
simon 0:1014af42efd9 12 * Target Processor: Cortex-M4/Cortex-M3
simon 0:1014af42efd9 13 *
simon 0:1014af42efd9 14 * Version 1.0.3 2010/11/29
simon 0:1014af42efd9 15 * Re-organized the CMSIS folders and updated documentation.
simon 0:1014af42efd9 16 *
simon 0:1014af42efd9 17 * Version 1.0.2 2010/11/11
simon 0:1014af42efd9 18 * Documentation updated.
simon 0:1014af42efd9 19 *
simon 0:1014af42efd9 20 * Version 1.0.1 2010/10/05
simon 0:1014af42efd9 21 * Production release and review comments incorporated.
simon 0:1014af42efd9 22 *
simon 0:1014af42efd9 23 * Version 1.0.0 2010/09/20
simon 0:1014af42efd9 24 * Production release and review comments incorporated.
simon 0:1014af42efd9 25 *
simon 0:1014af42efd9 26 * Version 0.0.7 2010/06/10
simon 0:1014af42efd9 27 * Misra-C changes done
simon 0:1014af42efd9 28 * -------------------------------------------------------------------- */
simon 0:1014af42efd9 29
simon 0:1014af42efd9 30 #include "arm_math.h"
simon 0:1014af42efd9 31
simon 0:1014af42efd9 32 /**
simon 0:1014af42efd9 33 * @ingroup groupMath
simon 0:1014af42efd9 34 */
simon 0:1014af42efd9 35
simon 0:1014af42efd9 36 /**
simon 0:1014af42efd9 37 * @addtogroup offset
simon 0:1014af42efd9 38 * @{
simon 0:1014af42efd9 39 */
simon 0:1014af42efd9 40
simon 0:1014af42efd9 41 /**
simon 0:1014af42efd9 42 * @brief Adds a constant offset to a Q7 vector.
simon 0:1014af42efd9 43 * @param *pSrc points to the input vector
simon 0:1014af42efd9 44 * @param offset is the offset to be added
simon 0:1014af42efd9 45 * @param *pDst points to the output vector
simon 0:1014af42efd9 46 * @param blockSize number of samples in the vector
simon 0:1014af42efd9 47 * @return none.
simon 0:1014af42efd9 48 *
simon 0:1014af42efd9 49 * <b>Scaling and Overflow Behavior:</b>
simon 0:1014af42efd9 50 * \par
simon 0:1014af42efd9 51 * The function uses saturating arithmetic.
simon 0:1014af42efd9 52 * Results outside of the allowable Q7 range [0x80 0x7F] are saturated.
simon 0:1014af42efd9 53 */
simon 0:1014af42efd9 54
simon 0:1014af42efd9 55 void arm_offset_q7(
simon 0:1014af42efd9 56 q7_t * pSrc,
simon 0:1014af42efd9 57 q7_t offset,
simon 0:1014af42efd9 58 q7_t * pDst,
simon 0:1014af42efd9 59 uint32_t blockSize)
simon 0:1014af42efd9 60 {
simon 0:1014af42efd9 61 uint32_t blkCnt; /* loop counter */
simon 0:1014af42efd9 62 q31_t offset_packed; /* Offset packed to 32 bit */
simon 0:1014af42efd9 63
simon 0:1014af42efd9 64
simon 0:1014af42efd9 65 /*loop Unrolling */
simon 0:1014af42efd9 66 blkCnt = blockSize >> 2u;
simon 0:1014af42efd9 67
simon 0:1014af42efd9 68 /* Offset is packed to 32 bit in order to use SIMD32 for addition */
simon 0:1014af42efd9 69 offset_packed = __PACKq7(offset, offset, offset, offset);
simon 0:1014af42efd9 70
simon 0:1014af42efd9 71 /* First part of the processing with loop unrolling. Compute 4 outputs at a time.
simon 0:1014af42efd9 72 ** a second loop below computes the remaining 1 to 3 samples. */
simon 0:1014af42efd9 73 while(blkCnt > 0u)
simon 0:1014af42efd9 74 {
simon 0:1014af42efd9 75 /* C = A + offset */
simon 0:1014af42efd9 76 /* Add offset and then store the results in the destination bufferfor 4 samples at a time. */
simon 0:1014af42efd9 77 *__SIMD32(pDst)++ = __QADD8(*__SIMD32(pSrc)++, offset_packed);
simon 0:1014af42efd9 78
simon 0:1014af42efd9 79 /* Decrement the loop counter */
simon 0:1014af42efd9 80 blkCnt--;
simon 0:1014af42efd9 81 }
simon 0:1014af42efd9 82
simon 0:1014af42efd9 83 /* If the blockSize is not a multiple of 4, compute any remaining output samples here.
simon 0:1014af42efd9 84 ** No loop unrolling is used. */
simon 0:1014af42efd9 85 blkCnt = blockSize % 0x4u;
simon 0:1014af42efd9 86
simon 0:1014af42efd9 87 while(blkCnt > 0u)
simon 0:1014af42efd9 88 {
simon 0:1014af42efd9 89 /* C = A + offset */
simon 0:1014af42efd9 90 /* Add offset and then store the result in the destination buffer. */
simon 0:1014af42efd9 91 *pDst++ = (q7_t) __SSAT(*pSrc++ + offset, 8);
simon 0:1014af42efd9 92
simon 0:1014af42efd9 93 /* Decrement the loop counter */
simon 0:1014af42efd9 94 blkCnt--;
simon 0:1014af42efd9 95 }
simon 0:1014af42efd9 96 }
simon 0:1014af42efd9 97
simon 0:1014af42efd9 98 /**
simon 0:1014af42efd9 99 * @} end of offset group
simon 0:1014af42efd9 100 */