123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177 |
- /* ----------------------------------------------------------------------
- * Copyright (C) 2010-2014 ARM Limited. All rights reserved.
- *
- * $Date: 12. March 2014
- * $Revision: V1.4.4
- *
- * Project: CMSIS DSP Library
- * Title: arm_min_q15.c
- *
- * Description: Minimum value 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 groupStats
- */
- /**
- * @addtogroup Min
- * @{
- */
- /**
- * @brief Minimum value of a Q15 vector.
- * @param[in] *pSrc points to the input vector
- * @param[in] blockSize length of the input vector
- * @param[out] *pResult minimum value returned here
- * @param[out] *pIndex index of minimum value returned here
- * @return none.
- *
- */
- void arm_min_q15(
- q15_t * pSrc,
- uint32_t blockSize,
- q15_t * pResult,
- uint32_t * pIndex)
- {
- #ifndef ARM_MATH_CM0_FAMILY
- /* Run the below code for Cortex-M4 and Cortex-M3 */
- q15_t minVal1, minVal2, out; /* Temporary variables to store the output value. */
- uint32_t blkCnt, outIndex, count; /* loop counter */
- /* Initialise the count value. */
- count = 0u;
- /* Initialise the index value to zero. */
- outIndex = 0u;
- /* Load first input value that act as reference value for comparision */
- out = *pSrc++;
- /* Loop unrolling */
- blkCnt = (blockSize - 1u) >> 2u;
- while(blkCnt > 0)
- {
- /* Initialize minVal to the next consecutive values one by one */
- minVal1 = *pSrc++;
- minVal2 = *pSrc++;
- /* compare for the minimum value */
- if(out > minVal1)
- {
- /* Update the minimum value and its index */
- out = minVal1;
- outIndex = count + 1u;
- }
- minVal1 = *pSrc++;
- /* compare for the minimum value */
- if(out > minVal2)
- {
- /* Update the minimum value and its index */
- out = minVal2;
- outIndex = count + 2u;
- }
- minVal2 = *pSrc++;
- /* compare for the minimum value */
- if(out > minVal1)
- {
- /* Update the minimum value and its index */
- out = minVal1;
- outIndex = count + 3u;
- }
- /* compare for the minimum value */
- if(out > minVal2)
- {
- /* Update the minimum value and its index */
- out = minVal2;
- outIndex = count + 4u;
- }
- count += 4u;
- blkCnt--;
- }
- /* if (blockSize - 1u ) is not multiple of 4 */
- blkCnt = (blockSize - 1u) % 4u;
- #else
- /* Run the below code for Cortex-M0 */
- q15_t minVal1, out; /* Temporary variables to store the output value. */
- uint32_t blkCnt, outIndex; /* loop counter */
- blkCnt = (blockSize - 1u);
- /* Initialise the index value to zero. */
- outIndex = 0u;
- /* Load first input value that act as reference value for comparision */
- out = *pSrc++;
- #endif // #ifndef ARM_MATH_CM0_FAMILY
- while(blkCnt > 0)
- {
- /* Initialize minVal to the next consecutive values one by one */
- minVal1 = *pSrc++;
- /* compare for the minimum value */
- if(out > minVal1)
- {
- /* Update the minimum value and it's index */
- out = minVal1;
- outIndex = blockSize - blkCnt;
- }
- blkCnt--;
- }
- /* Store the minimum value and its index into destination pointers */
- *pResult = out;
- *pIndex = outIndex;
- }
- /**
- * @} end of Min group
- */
|