1 /* ----------------------------------------------------------------------
4  * Description:  Floating-point 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
36   @defgroup CmplxByCmplxMult Complex-by-Complex Multiplication
49       pDst[(2*n)+0] = pSrcA[(2*n)+0] * pSrcB[(2*n)+0] - pSrcA[(2*n)+1] * pSrcB[(2*n)+1];
54   There are separate functions for floating-point, Q15, and Q31 data types.
63   @brief         Floating-point complex-by-complex multiplication.
84     blkCnt -= 1;  in arm_cmplx_mult_cmplx_f32()
115             blkCnt--;  in arm_cmplx_mult_cmplx_f32()
149             blkCnt -= 4;  in arm_cmplx_mult_cmplx_f32()
173             blkCnt -= 4;  in arm_cmplx_mult_cmplx_f32()
188     float32_t a, b, c, d;  /* Temporary variables to store real and imaginary values */  in arm_cmplx_mult_cmplx_f32()  local
199         va = vld2q_f32(pSrcA);  // load & separate real/imag pSrcA (de-interleave 2)  in arm_cmplx_mult_cmplx_f32()
206 	/* Re{C} = Re{A}*Re{B} - Im{A}*Im{B} */  in arm_cmplx_mult_cmplx_f32()
220         blkCnt--;  in arm_cmplx_mult_cmplx_f32()
234     /* C[2 * i    ] = A[2 * i] * B[2 * i    ] - A[2 * i + 1] * B[2 * i + 1]. */  in arm_cmplx_mult_cmplx_f32()
240     d = *pSrcB++;  in arm_cmplx_mult_cmplx_f32()
242     *pDst++ = (a * c) - (b * d);  in arm_cmplx_mult_cmplx_f32()
243     *pDst++ = (a * d) + (b * c);  in arm_cmplx_mult_cmplx_f32()
248     d = *pSrcB++;  in arm_cmplx_mult_cmplx_f32()
249     *pDst++ = (a * c) - (b * d);  in arm_cmplx_mult_cmplx_f32()
250     *pDst++ = (a * d) + (b * c);  in arm_cmplx_mult_cmplx_f32()
255     d = *pSrcB++;  in arm_cmplx_mult_cmplx_f32()
256     *pDst++ = (a * c) - (b * d);  in arm_cmplx_mult_cmplx_f32()
257     *pDst++ = (a * d) + (b * c);  in arm_cmplx_mult_cmplx_f32()
262     d = *pSrcB++;  in arm_cmplx_mult_cmplx_f32()
263     *pDst++ = (a * c) - (b * d);  in arm_cmplx_mult_cmplx_f32()
264     *pDst++ = (a * d) + (b * c);  in arm_cmplx_mult_cmplx_f32()
267     blkCnt--;  in arm_cmplx_mult_cmplx_f32()
283     /* C[2 * i    ] = A[2 * i] * B[2 * i    ] - A[2 * i + 1] * B[2 * i + 1]. */  in arm_cmplx_mult_cmplx_f32()
289     d = *pSrcB++;  in arm_cmplx_mult_cmplx_f32()
292     *pDst++ = (a * c) - (b * d);  in arm_cmplx_mult_cmplx_f32()
293     *pDst++ = (a * d) + (b * c);  in arm_cmplx_mult_cmplx_f32()
296     blkCnt--;  in arm_cmplx_mult_cmplx_f32()