Skip to content

/github/workspace/src/ComplexMathFunctions/kernels/plp_cmplx_mag_q16s_xpulpv2.c

Functions

Name
void plp_cmplx_mag_q16s_xpulpv2(const int16_t * pSrc, const uint32_t fracBits, int16_t * pRes, uint32_t numSamples)
complex magnitude for q16 on XPULPV2

Functions Documentation

function plp_cmplx_mag_q16s_xpulpv2

void plp_cmplx_mag_q16s_xpulpv2(
    const int16_t * pSrc,
    const uint32_t fracBits,
    int16_t * pRes,
    uint32_t numSamples
)

complex magnitude for q16 on XPULPV2

Parameters:

  • pSrc pointer to source
  • fracBits fractional bits -> Q(32-fracBits).fracBits
  • pRes pointer to result
  • numSamples The number of samples

Source code

/* =====================================================================
 * Project:      PULP DSP Library
 * Title:        plp_cmplx_mag_q16s_xpulpv2.c
 * Description:  q16 complex magnitude for XPULPV2
 *
 * $Date:        31. July 2020
 * $Revision:    V0
 *
 * Target Processor: PULP cores
 * ===================================================================== */
/*
 * Copyright (C) 2019 ETH Zurich and University of Bologna. 
 *
 * Author: Michael Rogenmoser, ETH Zurich
 *
 * SPDX-License-Identifier: Apache-2.0
 *
 * Licensed under the Apache License, Version 2.0 (the License); you may
 * not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 * www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an AS IS BASIS, WITHOUT
 * WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 *
 * Notice: project inspired by ARM CMSIS DSP and parts of source code
  ported and adopted for RISC-V PULP platform from ARM CMSIS DSP
 released under Copyright (C) 2010-2019 ARM Limited or its affiliates
  with Apache-2.0.
 */


#include "plp_math.h"


void plp_cmplx_mag_q16s_xpulpv2(const int16_t *pSrc,
                                const uint32_t fracBits,
                                int16_t *pRes,
                                uint32_t numSamples){

    int16_t real, cmplx, sqr;
    for (int i = 0; i < numSamples; i++) {

        real = (pSrc[2 * i] * pSrc[2 * i]) >> fracBits;
        cmplx = (pSrc[2 * i + 1] * pSrc[2 * i + 1]) >> fracBits;
        sqr = real + cmplx;
        plp_sqrt_q16(&sqr, fracBits, &pRes[i]);
    }
}

Updated on 2023-03-01 at 16:16:32 +0000