選択できるのは25トピックまでです。 トピックは、先頭が英数字で、英数字とダッシュ('-')を使用した35文字以内のものにしてください。
 
 
 
 
 

68 行
1.7 KiB

  1. /* ----------------------------------------------------------------------
  2. * Project: CMSIS DSP Library
  3. * Title: arm_cityblock_distance_f64.c
  4. * Description: Cityblock (Manhattan) 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 Manhattan
  33. @{
  34. */
  35. /**
  36. * @brief Cityblock (Manhattan) distance between two vectors
  37. * @param[in] pA First vector
  38. * @param[in] pB Second vector
  39. * @param[in] blockSize vector length
  40. * @return distance
  41. *
  42. */
  43. float64_t arm_cityblock_distance_f64(const float64_t *pA,const float64_t *pB, uint32_t blockSize)
  44. {
  45. float64_t accum,tmpA, tmpB;
  46. accum = 0.;
  47. while(blockSize > 0)
  48. {
  49. tmpA = *pA++;
  50. tmpB = *pB++;
  51. accum += fabs(tmpA - tmpB);
  52. blockSize --;
  53. }
  54. return(accum);
  55. }
  56. /**
  57. * @} end of Manhattan group
  58. */