1  /* ----------------------------------------------------------------------
  2   * Project:      CMSIS DSP Library
  3   * Title:        arm_vlog_f32.c
  4   * Description:  Fast vectorized log
  5   *
  6   * $Date:        23 April 2021
  7   * $Revision:    V1.9.0
  8   *
  9   * Target Processor: Cortex-M and Cortex-A cores
 10   * -------------------------------------------------------------------- */
 11  /*
 12   * Copyright (C) 2010-2021 ARM Limited or its affiliates. All rights reserved.
 13   *
 14   * SPDX-License-Identifier: Apache-2.0
 15   *
 16   * Licensed under the Apache License, Version 2.0 (the License); you may
 17   * not use this file except in compliance with the License.
 18   * You may obtain a copy of the License at
 19   *
 20   * www.apache.org/licenses/LICENSE-2.0
 21   *
 22   * Unless required by applicable law or agreed to in writing, software
 23   * distributed under the License is distributed on an AS IS BASIS, WITHOUT
 24   * WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 25   * See the License for the specific language governing permissions and
 26   * limitations under the License.
 27   */
 28  
 29  #include "dsp/fast_math_functions.h"
 30  #include "arm_common_tables.h"
 31  
 32  
 33  /**
 34    @ingroup groupFastMath
 35   */
 36  
 37  
 38  /**
 39    @defgroup vlog Vector Log
 40  
 41    Compute the log values of a vector of samples.
 42  
 43   */
 44  
 45  /**
 46    @addtogroup vlog
 47    @{
 48   */
 49  
 50  #if (defined(ARM_MATH_MVEF) || defined(ARM_MATH_HELIUM) || defined(ARM_MATH_NEON) || defined(ARM_MATH_NEON_EXPERIMENTAL)) && !defined(ARM_MATH_AUTOVECTORIZE)
 51  #include "arm_vec_math.h"
 52  #endif
 53  
 54  void arm_vlog_f32(
 55    const float32_t * pSrc,
 56          float32_t * pDst,
 57          uint32_t blockSize)
 58  {
 59     uint32_t blkCnt; 
 60  
 61  #if (defined(ARM_MATH_MVEF) || defined(ARM_MATH_HELIUM)) && !defined(ARM_MATH_AUTOVECTORIZE)
 62     f32x4_t src;
 63     f32x4_t dst;
 64  
 65     blkCnt = blockSize >> 2;
 66  
 67     while (blkCnt > 0U)
 68     {
 69        src = vld1q(pSrc);
 70        dst = vlogq_f32(src);
 71        vst1q(pDst, dst);
 72  
 73        pSrc += 4;
 74        pDst += 4;
 75        /* Decrement loop counter */
 76        blkCnt--;
 77     }
 78  
 79     blkCnt = blockSize & 3;
 80  #else
 81  #if (defined(ARM_MATH_NEON) || defined(ARM_MATH_NEON_EXPERIMENTAL)) && !defined(ARM_MATH_AUTOVECTORIZE)
 82     f32x4_t src;
 83     f32x4_t dst;
 84  
 85     blkCnt = blockSize >> 2;
 86  
 87     while (blkCnt > 0U)
 88     {
 89        src = vld1q_f32(pSrc);
 90        dst = vlogq_f32(src);
 91        vst1q_f32(pDst, dst);
 92  
 93        pSrc += 4;
 94        pDst += 4;
 95        /* Decrement loop counter */
 96        blkCnt--;
 97     }
 98  
 99     blkCnt = blockSize & 3;
100  #else
101     blkCnt = blockSize;
102  #endif
103  #endif
104  
105     while (blkCnt > 0U)
106     {
107        /* C = log(A) */
108    
109        /* Calculate log and store result in destination buffer. */
110        *pDst++ = logf(*pSrc++);
111    
112        /* Decrement loop counter */
113        blkCnt--;
114     }
115  }
116  
117  /**
118    @} end of vlog group
119   */