arm_chebyshev_distance_f64.c 1.8 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576
  1. /* ----------------------------------------------------------------------
  2. * Project: CMSIS DSP Library
  3. * Title: arm_chebyshev_distance_f64.c
  4. * Description: Chebyshev 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 Chebyshev
  33. @{
  34. */
  35. /**
  36. * @brief Chebyshev 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_chebyshev_distance_f64(const float64_t *pA,const float64_t *pB, uint32_t blockSize)
  44. {
  45. float64_t diff=0., maxVal,tmpA, tmpB;
  46. tmpA = *pA++;
  47. tmpB = *pB++;
  48. diff = fabs(tmpA - tmpB);
  49. maxVal = diff;
  50. blockSize--;
  51. while(blockSize > 0)
  52. {
  53. tmpA = *pA++;
  54. tmpB = *pB++;
  55. diff = fabs(tmpA - tmpB);
  56. if (diff > maxVal)
  57. {
  58. maxVal = diff;
  59. }
  60. blockSize --;
  61. }
  62. return(maxVal);
  63. }
  64. /**
  65. * @} end of Chebyshev group
  66. */