LCOV - code coverage report
Current view: top level - geospatial - latlong_metrics.cc (source / functions) Hit Total Coverage
Test: Test Coverage for xapian-core 954b5873a738 Lines: 61 67 91.0 %
Date: 2019-06-30 05:20:33 Functions: 10 10 100.0 %
Branches: 30 68 44.1 %

           Branch data     Line data    Source code
       1                 :            : /** @file latlong_metrics.cc
       2                 :            :  * @brief Geospatial distance metrics.
       3                 :            :  */
       4                 :            : /* Copyright 2008 Lemur Consulting Ltd
       5                 :            :  * Copyright 2011 Richard Boulton
       6                 :            :  *
       7                 :            :  * This program is free software; you can redistribute it and/or
       8                 :            :  * modify it under the terms of the GNU General Public License as
       9                 :            :  * published by the Free Software Foundation; either version 2 of the
      10                 :            :  * License, or (at your option) any later version.
      11                 :            :  *
      12                 :            :  * This program is distributed in the hope that it will be useful,
      13                 :            :  * but WITHOUT ANY WARRANTY; without even the implied warranty of
      14                 :            :  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
      15                 :            :  * GNU General Public License for more details.
      16                 :            :  *
      17                 :            :  * You should have received a copy of the GNU General Public License
      18                 :            :  * along with this program; if not, write to the Free Software
      19                 :            :  * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301
      20                 :            :  * USA
      21                 :            :  */
      22                 :            : 
      23                 :            : #include <config.h>
      24                 :            : 
      25                 :            : #include "xapian/geospatial.h"
      26                 :            : #include "xapian/error.h"
      27                 :            : #include "serialise-double.h"
      28                 :            : 
      29                 :            : #include <cmath>
      30                 :            : 
      31                 :            : using namespace Xapian;
      32                 :            : using namespace std;
      33                 :            : 
      34                 :            : /** Quadratic mean radius of the Earth in metres.
      35                 :            :  */
      36                 :            : #define QUAD_EARTH_RADIUS_METRES 6372797.6
      37                 :            : 
      38                 :       4717 : LatLongMetric::~LatLongMetric()
      39                 :            : {
      40         [ -  + ]:       4717 : }
      41                 :            : 
      42                 :            : double
      43                 :         15 : LatLongMetric::operator()(const LatLongCoords & a,
      44                 :            :                           const LatLongCoords &b) const
      45                 :            : {
      46 [ +  - ][ -  + ]:         15 :     if (a.empty() || b.empty()) {
                 [ -  + ]
      47 [ #  # ][ #  # ]:          0 :         throw InvalidArgumentError("Empty coordinate list supplied to LatLongMetric::operator()().");
                 [ #  # ]
      48                 :            :     }
      49                 :         15 :     double min_dist = 0.0;
      50                 :         15 :     bool have_min = false;
      51         [ +  + ]:         68 :     for (LatLongCoordsIterator a_iter = a.begin();
      52                 :         34 :          a_iter != a.end();
      53                 :            :          ++a_iter)
      54                 :            :     {
      55         [ +  + ]:         80 :         for (LatLongCoordsIterator b_iter = b.begin();
      56                 :         40 :              b_iter != b.end();
      57                 :            :              ++b_iter)
      58                 :            :         {
      59         [ +  - ]:         21 :             double dist = pointwise_distance(*a_iter, *b_iter);
      60         [ +  + ]:         21 :             if (!have_min) {
      61                 :         15 :                 min_dist = dist;
      62                 :         15 :                 have_min = true;
      63         [ +  + ]:          6 :             } else if (dist < min_dist) {
      64                 :          3 :                 min_dist = dist;
      65                 :            :             }
      66                 :            :         }
      67                 :            :     }
      68                 :         15 :     return min_dist;
      69                 :            : }
      70                 :            : 
      71                 :            : double
      72                 :         91 : LatLongMetric::operator()(const LatLongCoords & a,
      73                 :            :                           const char * b_ptr, size_t b_len) const
      74                 :            : {
      75 [ +  - ][ -  + ]:         91 :     if (a.empty() || b_len == 0) {
                 [ -  + ]
      76 [ #  # ][ #  # ]:          0 :         throw InvalidArgumentError("Empty coordinate list supplied to LatLongMetric::operator()().");
                 [ #  # ]
      77                 :            :     }
      78                 :         91 :     double min_dist = 0.0;
      79                 :         91 :     bool have_min = false;
      80                 :         91 :     LatLongCoord b;
      81                 :         91 :     const char * b_end = b_ptr + b_len;
      82         [ +  + ]:        182 :     while (b_ptr != b_end) {
      83         [ +  - ]:         91 :         b.unserialise(&b_ptr, b_end);
      84         [ +  + ]:        364 :         for (LatLongCoordsIterator a_iter = a.begin();
      85                 :        182 :              a_iter != a.end();
      86                 :            :              ++a_iter)
      87                 :            :         {
      88         [ +  - ]:         91 :             double dist = pointwise_distance(*a_iter, b);
      89         [ +  - ]:         91 :             if (!have_min) {
      90                 :         91 :                 min_dist = dist;
      91                 :         91 :                 have_min = true;
      92         [ #  # ]:          0 :             } else if (dist < min_dist) {
      93                 :          0 :                 min_dist = dist;
      94                 :            :             }
      95                 :            :         }
      96                 :            :     }
      97                 :         91 :     return min_dist;
      98                 :            : }
      99                 :            : 
     100                 :            : 
     101                 :       3138 : GreatCircleMetric::GreatCircleMetric()
     102                 :       3138 :         : radius(QUAD_EARTH_RADIUS_METRES)
     103                 :       3138 : {}
     104                 :            : 
     105                 :       1579 : GreatCircleMetric::GreatCircleMetric(double radius_)
     106                 :       1579 :         : radius(radius_)
     107                 :       1579 : {}
     108                 :            : 
     109                 :            : double
     110                 :        112 : GreatCircleMetric::pointwise_distance(const LatLongCoord & a,
     111                 :            :                                       const LatLongCoord & b) const
     112                 :            : {
     113                 :        112 :     double lata = a.latitude * (M_PI / 180.0);
     114                 :        112 :     double latb = b.latitude * (M_PI / 180.0);
     115                 :            : 
     116                 :        112 :     double latdiff = lata - latb;
     117                 :        112 :     double longdiff = (a.longitude - b.longitude) * (M_PI / 180.0);
     118                 :            : 
     119                 :        112 :     double sin_half_lat = sin(latdiff / 2);
     120                 :        112 :     double sin_half_long = sin(longdiff / 2);
     121                 :        112 :     double h = sin_half_lat * sin_half_lat +
     122                 :        112 :             sin_half_long * sin_half_long * cos(lata) * cos(latb);
     123         [ -  + ]:        112 :     if (rare(h > 1.0)) {
     124                 :            :         // Clamp to 1.0, asin(1.0) = M_PI / 2.0.
     125                 :          0 :         return radius * M_PI;
     126                 :            :     }
     127                 :        112 :     return 2 * radius * asin(sqrt(h));
     128                 :            : }
     129                 :            : 
     130                 :            : LatLongMetric *
     131                 :       1576 : GreatCircleMetric::clone() const
     132                 :            : {
     133         [ +  - ]:       1576 :     return new GreatCircleMetric(radius);
     134                 :            : }
     135                 :            : 
     136                 :            : string
     137                 :       1560 : GreatCircleMetric::name() const
     138                 :            : {
     139         [ +  - ]:       1560 :     return "Xapian::GreatCircleMetric";
     140                 :            : }
     141                 :            : 
     142                 :            : string
     143                 :          1 : GreatCircleMetric::serialise() const
     144                 :            : {
     145                 :          1 :     return serialise_double(radius);
     146                 :            : }
     147                 :            : 
     148                 :            : LatLongMetric *
     149                 :          1 : GreatCircleMetric::unserialise(const string & s) const
     150                 :            : {
     151                 :          1 :     const char * p = s.data();
     152                 :          1 :     const char * end = p + s.size();
     153                 :            : 
     154         [ +  - ]:          1 :     double new_radius = unserialise_double(&p, end);
     155         [ -  + ]:          1 :     if (p != end) {
     156 [ #  # ][ #  # ]:          0 :         throw Xapian::NetworkError("Bad serialised GreatCircleMetric - junk at end");
                 [ #  # ]
     157                 :            :     }
     158                 :            : 
     159 [ +  - ][ +  - ]:          1 :     return new GreatCircleMetric(new_radius);
     160                 :            : }

Generated by: LCOV version 1.11