geos/src/noding/Octant.cpp at main · Toucad/geos · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2006 Refractions Research Inc.
 *
 * This is free software; you can redistribute and/or modify it under
 * the terms of the GNU Lesser General Licence as published
 * by the Free Software Foundation.
 * See the COPYING file for more information.
 *
 **********************************************************************
 *
 * Last port: noding/Octant.java rev. 1.2 (JTS-1.7)
 *
 **********************************************************************/

#include <cmath>
#include <sstream>

#include <geos/util/IllegalArgumentException.h>
#include <geos/noding/Octant.h>
#include <geos/geom/Coordinate.h>

//
using namespace geos::geom;

namespace geos {
namespace noding { // geos.noding

/*public static*/
int
Octant::octant(double dx, double dy)
{
    if(dx == 0.0 && dy == 0.0) {
        std::ostringstream s;
        s << "Cannot compute the octant for point ( " << dx << ", " << dy << " )";
        throw util::IllegalArgumentException(s.str());
    }

    double adx = std::fabs(dx);
    double ady = std::fabs(dy);

    if(dx >= 0) {
        if(dy >= 0) {
            if(adx >= ady) {
                return 0;
            }
            else {
                return 1;
            }
        }
        else { // dy < 0
            if(adx >= ady) {
                return 7;
            }
            else {
                return 6;
            }
        }
    }
    else { // dx < 0
        if(dy >= 0) {
            if(adx >= ady) {
                return 3;
            }
            else {
                return 2;
            }
        }
        else { // dy < 0
            if(adx >= ady) {
                return 4;
            }
            else {
                return 5;
            }
        }
    }

}

/*public static*/
int
Octant::octant(const CoordinateXY& p0, const CoordinateXY& p1)
{
    double dx = p1.x - p0.x;
    double dy = p1.y - p0.y;

    return octant(dx, dy);
}

} // namespace geos.noding
} // namespace geos