You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
 
 
 
 
 

71 line
1.7 KiB

  1. /* ----------------------------------------------------------------------
  2. * Project: CMSIS DSP Library
  3. * Title: arm_cosine_distance_f64.c
  4. * Description: Cosine distance between two vectors
  5. *
  6. * $Date: 13 September 2021
  7. * $Revision: V1.10.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. #include "dsp/distance_functions.h"
  29. #include <limits.h>
  30. #include <math.h>
  31. /**
  32. @addtogroup CosineDist
  33. @{
  34. */
  35. /**
  36. * @brief Cosine distance between two vectors
  37. *
  38. * @param[in] pA First vector
  39. * @param[in] pB Second vector
  40. * @param[in] blockSize vector length
  41. * @return distance
  42. *
  43. */
  44. float64_t arm_cosine_distance_f64(const float64_t *pA,const float64_t *pB, uint32_t blockSize)
  45. {
  46. float64_t pwra,pwrb,dot,tmp;
  47. arm_power_f64(pA, blockSize, &pwra);
  48. arm_power_f64(pB, blockSize, &pwrb);
  49. arm_dot_prod_f64(pA,pB,blockSize,&dot);
  50. tmp = sqrt(pwra * pwrb);
  51. return(1. - dot / tmp);
  52. }
  53. /**
  54. * @} end of CosineDist group
  55. */