1 /* ----------------------------------------------------------------------
4  * Description:  Q31 complex-by-complex multiplication
9  * Target Processor: Cortex-M and Cortex-A cores
10  * -------------------------------------------------------------------- */
12  * Copyright (C) 2010-2021 ARM Limited or its affiliates. All rights reserved.
14  * SPDX-License-Identifier: Apache-2.0
20  * www.apache.org/licenses/LICENSE-2.0
41   @brief         Q31 complex-by-complex multiplication.
65     blkCnt -= 1;  in arm_cmplx_mult_cmplx_q31()
75             /* C[2 * i] = A[2 * i] * B[2 * i] - A[2 * i + 1] * B[2 * i + 1].  */  in arm_cmplx_mult_cmplx_q31()
102             blkCnt--;  in arm_cmplx_mult_cmplx_q31()
141             blkCnt -= 4;  in arm_cmplx_mult_cmplx_q31()
162             blkCnt -= 4;  in arm_cmplx_mult_cmplx_q31()
174         q31_t a, b, c, d;                              /* Temporary variables */  in arm_cmplx_mult_cmplx_q31()  local
183     /* C[2 * i    ] = A[2 * i] * B[2 * i    ] - A[2 * i + 1] * B[2 * i + 1]. */  in arm_cmplx_mult_cmplx_q31()
189     d = *pSrcB++;  in arm_cmplx_mult_cmplx_q31()
191     *pDst++ = (q31_t) ( (((q63_t) a * c) >> 33) - (((q63_t) b * d) >> 33) );  in arm_cmplx_mult_cmplx_q31()
192     *pDst++ = (q31_t) ( (((q63_t) a * d) >> 33) + (((q63_t) b * c) >> 33) );  in arm_cmplx_mult_cmplx_q31()
197     d = *pSrcB++;  in arm_cmplx_mult_cmplx_q31()
198     *pDst++ = (q31_t) ( (((q63_t) a * c) >> 33) - (((q63_t) b * d) >> 33) );  in arm_cmplx_mult_cmplx_q31()
199     *pDst++ = (q31_t) ( (((q63_t) a * d) >> 33) + (((q63_t) b * c) >> 33) );  in arm_cmplx_mult_cmplx_q31()
204     d = *pSrcB++;  in arm_cmplx_mult_cmplx_q31()
205     *pDst++ = (q31_t) ( (((q63_t) a * c) >> 33) - (((q63_t) b * d) >> 33) );  in arm_cmplx_mult_cmplx_q31()
206     *pDst++ = (q31_t) ( (((q63_t) a * d) >> 33) + (((q63_t) b * c) >> 33) );  in arm_cmplx_mult_cmplx_q31()
211     d = *pSrcB++;  in arm_cmplx_mult_cmplx_q31()
212     *pDst++ = (q31_t) ( (((q63_t) a * c) >> 33) - (((q63_t) b * d) >> 33) );  in arm_cmplx_mult_cmplx_q31()
213     *pDst++ = (q31_t) ( (((q63_t) a * d) >> 33) + (((q63_t) b * c) >> 33) );  in arm_cmplx_mult_cmplx_q31()
216     blkCnt--;  in arm_cmplx_mult_cmplx_q31()
231     /* C[2 * i    ] = A[2 * i] * B[2 * i    ] - A[2 * i + 1] * B[2 * i + 1]. */  in arm_cmplx_mult_cmplx_q31()
237     d = *pSrcB++;  in arm_cmplx_mult_cmplx_q31()
240     *pDst++ = (q31_t) ( (((q63_t) a * c) >> 33) - (((q63_t) b * d) >> 33) );  in arm_cmplx_mult_cmplx_q31()
241     *pDst++ = (q31_t) ( (((q63_t) a * d) >> 33) + (((q63_t) b * c) >> 33) );  in arm_cmplx_mult_cmplx_q31()
244     blkCnt--;  in arm_cmplx_mult_cmplx_q31()