1  
 2  /* ----------------------------------------------------------------------
 3   * Project:      CMSIS DSP Library
 4   * Title:        arm_sokalsneath_distance.c
 5   * Description:  Sokal-Sneath distance between two vectors
 6   *
 7   * $Date:        23 April 2021
 8   * $Revision:    V1.9.0
 9   *
10   * Target Processor: Cortex-M and Cortex-A cores
11   * -------------------------------------------------------------------- */
12  /*
13   * Copyright (C) 2010-2021 ARM Limited or its affiliates. All rights reserved.
14   *
15   * SPDX-License-Identifier: Apache-2.0
16   *
17   * Licensed under the Apache License, Version 2.0 (the License); you may
18   * not use this file except in compliance with the License.
19   * You may obtain a copy of the License at
20   *
21   * www.apache.org/licenses/LICENSE-2.0
22   *
23   * Unless required by applicable law or agreed to in writing, software
24   * distributed under the License is distributed on an AS IS BASIS, WITHOUT
25   * WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
26   * See the License for the specific language governing permissions and
27   * limitations under the License.
28   */
29  
30  #include "dsp/distance_functions.h"
31  #include <limits.h>
32  #include <math.h>
33  
34  
35  extern void arm_boolean_distance_TT_TF_FT(const uint32_t *pA
36         , const uint32_t *pB
37         , uint32_t numberOfBools
38         , uint32_t *cTT
39         , uint32_t *cTF
40         , uint32_t *cFT
41         );
42  
43  
44  /**
45    @addtogroup BoolDist
46    @{
47   */
48  
49  /**
50   * @brief        Sokal-Sneath distance between two vectors
51   *
52   * @param[in]    pA              First vector of packed booleans
53   * @param[in]    pB              Second vector of packed booleans
54   * @param[in]    numberOfBools   Number of booleans
55   * @return distance
56   *
57   */
58  
59  float32_t arm_sokalsneath_distance(const uint32_t *pA, const uint32_t *pB, uint32_t numberOfBools)
60  {
61      uint32_t ctt=0,cft=0,ctf=0;
62      float32_t r;
63  
64      arm_boolean_distance_TT_TF_FT(pA, pB, numberOfBools, &ctt, &ctf, &cft);
65  
66      r = 2.0*(ctf + cft);
67  
68      return(r / (r + ctt));
69  }
70  
71  
72  /**
73   * @} end of BoolDist group
74   */