ExtensionFunctionsGeo.hpp File Reference

#include "Geospatial/Compression.h"
#include "Geospatial/CompressionRuntime.h"
#include "Geospatial/Types.h"
#include "Geospatial/Utm.h"
#include <iomanip>
#include <string>
#include "Shared/likely.h"

Include dependency graph for ExtensionFunctionsGeo.hpp:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
struct	anonymous_namespace{ExtensionFunctionsGeo.hpp}::CoordData
struct	anonymous_namespace{ExtensionFunctionsGeo.hpp}::Point2D

Namespaces
	anonymous_namespace{ExtensionFunctionsGeo.hpp}

Macros
#define	DEBUG_STMT(x)

Enumerations
enum	anonymous_namespace{ExtensionFunctionsGeo.hpp}::Coords { anonymous_namespace{ExtensionFunctionsGeo.hpp}::None = 0, anonymous_namespace{ExtensionFunctionsGeo.hpp}::Y, anonymous_namespace{ExtensionFunctionsGeo.hpp}::X, anonymous_namespace{ExtensionFunctionsGeo.hpp}::XY }
enum	EdgeBehavior { EdgeBehavior::kIncludePointOnEdge, EdgeBehavior::kExcludePointOnEdge }

Functions
DEVICE ALWAYS_INLINE double	tol (int32_t ic)
DEVICE ALWAYS_INLINE double	tol (const int32_t ic1, const int32_t ic2)
DEVICE ALWAYS_INLINE bool	tol_zero (const double x, const double tolerance=TOLERANCE_DEFAULT)
DEVICE ALWAYS_INLINE bool	tol_eq (const double x, const double y, const double tolerance=TOLERANCE_DEFAULT)
DEVICE ALWAYS_INLINE bool	tol_le (const double x, const double y, const double tolerance=TOLERANCE_DEFAULT)
DEVICE ALWAYS_INLINE bool	tol_ge (const double x, const double y, const double tolerance=TOLERANCE_DEFAULT)
template<Coords C>
DEVICE ALWAYS_INLINE double	decompress_coord (const int8_t *data, const int32_t x_index, const int32_t ic)
DEVICE ALWAYS_INLINE int32_t	compression_unit_size (const int32_t ic)
template<Coords C>
DEVICE ALWAYS_INLINE Point2D	conv_4326_900913 (const Point2D point)
template<Coords C>
DEVICE ALWAYS_INLINE Point2D	conv_4326_utm (const Point2D point, const int32_t utm_srid)
template<Coords C>
DEVICE ALWAYS_INLINE Point2D	conv_utm_4326 (const Point2D point, const int32_t utm_srid)
template<Coords C>
DEVICE ALWAYS_INLINE Point2D	transform_point (const Point2D point, const int32_t isr, const int32_t osr)
DEVICE ALWAYS_INLINE bool	x_and_y_are_dependent (const int32_t isr, const int32_t osr)
DEVICE ALWAYS_INLINE Point2D	get_point (const int8_t *data, const int32_t x_index, const int32_t ic, const int32_t isr, const int32_t osr)
template<Coords C>
DEVICE ALWAYS_INLINE Point2D	coord (const int8_t *data, const int32_t x_index, const int32_t ic, const int32_t isr, const int32_t osr)
DEVICE ALWAYS_INLINE int32_t	compressed_coord (const int8_t *data, const int32_t index)
DEVICE ALWAYS_INLINE double	distance_point_point_squared (double p1x, double p1y, double p2x, double p2y)
DEVICE ALWAYS_INLINE double	distance_point_point (double p1x, double p1y, double p2x, double p2y)
DEVICE double	distance_point_line (double px, double py, double l1x, double l1y, double l2x, double l2y)
DEVICE double	distance_point_line_squared (double px, double py, double l1x, double l1y, double l2x, double l2y)
DEVICE ALWAYS_INLINE bool	on_segment (double px, double py, double qx, double qy, double rx, double ry)
DEVICE ALWAYS_INLINE int16_t	orientation (double px, double py, double qx, double qy, double rx, double ry)
DEVICE bool	line_intersects_line (double l11x, double l11y, double l12x, double l12y, double l21x, double l21y, double l22x, double l22y)
DEVICE bool	linestring_intersects_line (int8_t *l, int32_t lnum_coords, double l1x, double l1y, double l2x, double l2y, int32_t ic1, int32_t isr1, int32_t osr)
DEVICE bool	ring_intersects_line (int8_t *ring, int32_t ring_num_coords, double l1x, double l1y, double l2x, double l2y, int32_t ic1, int32_t isr1, int32_t osr)
DEVICE bool	linestring_intersects_linestring (int8_t *l, int32_t lnum_coords, double l1x, double l1y, double l2x, double l2y, int32_t ic1, int32_t isr1, int32_t osr)
DEVICE double	distance_line_line (double l11x, double l11y, double l12x, double l12y, double l21x, double l21y, double l22x, double l22y)
DEVICE double	distance_line_line_squared (double l11x, double l11y, double l12x, double l12y, double l21x, double l21y, double l22x, double l22y)
DEVICE double	distance_ring_linestring (int8_t *ring, int32_t ring_num_coords, int8_t *l, int32_t lnum_coords, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
DEVICE double	distance_ring_ring (int8_t *ring1, int32_t ring1_num_coords, int8_t *ring2, int32_t ring2_num_coords, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
template<typename T >
DEVICE ALWAYS_INLINE T	is_left (const T lx0, const T ly0, const T lx1, const T ly1, const T px, const T py)
template<typename T >
DEVICE ALWAYS_INLINE bool	tol_zero_template (const T x, const T tolerance=TOLERANCE_DEFAULT)
template<typename T , EdgeBehavior TEdgeBehavior>
DEVICE ALWAYS_INLINE bool	point_in_polygon_winding_number (const int8_t *poly, const int32_t poly_num_coords, const T px, const T py, const int32_t ic1, const int32_t isr1, const int32_t osr)
DEVICE ALWAYS_INLINE bool	polygon_contains_point (const int8_t *poly, const int32_t poly_num_coords, const double px, const double py, const int32_t ic1, const int32_t isr1, const int32_t osr)
DEVICE bool	polygon_contains_linestring (int8_t *poly, int32_t poly_num_coords, int8_t *l, int64_t lnum_coords, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
DEVICE ALWAYS_INLINE bool	box_contains_point (const double *bounds, const int64_t bounds_size, const double px, const double py)
EXTENSION_NOINLINE bool	Point_Overlaps_Box (double *bounds, int64_t bounds_size, double px, double py)
DEVICE ALWAYS_INLINE bool	box_contains_box (double *bounds1, int64_t bounds1_size, double *bounds2, int64_t bounds2_size)
DEVICE ALWAYS_INLINE bool	box_contains_box_vertex (double *bounds1, int64_t bounds1_size, double *bounds2, int64_t bounds2_size)
DEVICE ALWAYS_INLINE bool	box_overlaps_box (double *bounds1, int64_t bounds1_size, double *bounds2, int64_t bounds2_size)
DEVICE ALWAYS_INLINE bool	point_dwithin_box (int8_t *p1, int64_t p1size, int32_t ic1, int32_t isr1, double *bounds2, int64_t bounds2_size, int32_t isr2, int32_t osr, double distance)
DEVICE ALWAYS_INLINE bool	box_dwithin_box (double *bounds1, int64_t bounds1_size, int32_t isr1, double *bounds2, int64_t bounds2_size, int32_t isr2, int32_t osr, double distance)
DEVICE ALWAYS_INLINE CoordData	trim_linestring_to_buffered_box (int8_t *l1, int64_t l1size, int32_t ic1, int32_t isr1, double *bounds2, int64_t bounds2_size, int32_t isr2, int32_t osr, double distance)
EXTENSION_NOINLINE double	ST_X_Point (int8_t *p, int64_t psize, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE double	ST_Y_Point (int8_t *p, int64_t psize, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE double	ST_XMin (int8_t *coords, int64_t size, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE double	ST_YMin (int8_t *coords, int64_t size, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE double	ST_XMax (int8_t *coords, int64_t size, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE double	ST_YMax (int8_t *coords, int64_t size, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_INLINE double	ST_XMin_Bounds (double *bounds, int64_t size, int32_t isr, int32_t osr)
EXTENSION_INLINE double	ST_YMin_Bounds (double *bounds, int64_t size, int32_t isr, int32_t osr)
EXTENSION_INLINE double	ST_XMax_Bounds (double *bounds, int64_t size, int32_t isr, int32_t osr)
EXTENSION_INLINE double	ST_YMax_Bounds (double *bounds, int64_t size, int32_t isr, int32_t osr)
DEVICE ALWAYS_INLINE double	length_linestring (int8_t *l, int32_t lsize, int32_t ic, int32_t isr, int32_t osr, bool geodesic, bool check_closed)
EXTENSION_NOINLINE double	ST_Length_LineString (int8_t *coords, int64_t coords_sz, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE double	ST_Length_LineString_Geodesic (int8_t *coords, int64_t coords_sz, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE double	ST_Length_MultiLineString (int8_t *coords, int64_t coords_sz, int8_t *linestring_sizes_in, int64_t linestring_sizes_sz, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE double	ST_Perimeter_Polygon (int8_t *poly, int32_t polysize, int8_t *poly_ring_sizes, int32_t poly_num_rings, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE double	ST_Perimeter_Polygon_Geodesic (int8_t *poly, int32_t polysize, int8_t *poly_ring_sizes_in, int32_t poly_num_rings, int32_t ic, int32_t isr, int32_t osr)
DEVICE ALWAYS_INLINE double	perimeter_multipolygon (int8_t *mpoly_coords, int32_t mpoly_coords_size, int8_t *mpoly_ring_sizes_in, int32_t mpoly_num_rings, int8_t *mpoly_poly_sizes, int32_t mpoly_num_polys, int32_t ic, int32_t isr, int32_t osr, bool geodesic)
EXTENSION_NOINLINE double	ST_Perimeter_MultiPolygon (int8_t *mpoly_coords, int32_t mpoly_coords_size, int8_t *mpoly_ring_sizes, int32_t mpoly_num_rings, int8_t *mpoly_poly_sizes, int32_t mpoly_num_polys, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE double	ST_Perimeter_MultiPolygon_Geodesic (int8_t *mpoly_coords, int32_t mpoly_coords_size, int8_t *mpoly_ring_sizes, int32_t mpoly_num_rings, int8_t *mpoly_poly_sizes, int32_t mpoly_num_polys, int32_t ic, int32_t isr, int32_t osr)
DEVICE ALWAYS_INLINE double	area_triangle (double x1, double y1, double x2, double y2, double x3, double y3)
DEVICE ALWAYS_INLINE double	area_ring (int8_t *ring, int64_t ringsize, int32_t ic, int32_t isr, int32_t osr)
DEVICE ALWAYS_INLINE double	area_polygon (int8_t *poly_coords, int32_t poly_coords_size, int8_t *poly_ring_sizes_in, int32_t poly_num_rings, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE double	ST_Area_Polygon (int8_t *poly_coords, int32_t poly_coords_size, int8_t *poly_ring_sizes, int32_t poly_num_rings, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE double	ST_Area_MultiPolygon (int8_t *mpoly_coords, int32_t mpoly_coords_size, int8_t *mpoly_ring_sizes, int32_t mpoly_num_rings, int8_t *mpoly_poly_sizes_in, int32_t mpoly_num_polys, int32_t ic, int32_t isr, int32_t osr)
EXTENSION_NOINLINE void	ST_Centroid_Point (int8_t *p, int32_t psize, int32_t ic, int32_t isr, int32_t osr, double *point_centroid)
EXTENSION_NOINLINE void	ST_Centroid_MultiPoint (int8_t *mp, int32_t mpsize, int32_t ic, int32_t isr, int32_t osr, double *multipoint_centroid)
DEVICE ALWAYS_INLINE bool	centroid_add_segment (double x1, double y1, double x2, double y2, double *length, double *linestring_centroid_sum)
DEVICE ALWAYS_INLINE bool	centroid_add_linestring (int8_t *l, int64_t lsize, int32_t ic, int32_t isr, int32_t osr, bool closed, double *total_length, double *linestring_centroid_sum, int64_t *num_points, double *point_centroid_sum)
EXTENSION_NOINLINE void	ST_Centroid_LineString (int8_t *coords, int32_t coords_sz, int32_t ic, int32_t isr, int32_t osr, double *linestring_centroid)
DEVICE ALWAYS_INLINE bool	centroid_add_triangle (double x1, double y1, double x2, double y2, double x3, double y3, double sign, double *total_area2, double *cg3)
DEVICE ALWAYS_INLINE bool	centroid_add_ring (int8_t *ring, int64_t ringsize, int32_t ic, int32_t isr, int32_t osr, double sign, double *total_area2, double *cg3, double *total_length, double *linestring_centroid_sum, int64_t *num_points, double *point_centroid_sum)
DEVICE ALWAYS_INLINE bool	centroid_add_polygon (int8_t *poly_coords, int64_t poly_coords_size, int32_t *poly_ring_sizes, int64_t poly_num_rings, int32_t ic, int32_t isr, int32_t osr, double *total_area2, double *cg3, double *total_length, double *linestring_centroid_sum, int64_t *num_points, double *point_centroid_sum)
EXTENSION_NOINLINE void	ST_Centroid_Polygon (int8_t *poly_coords, int32_t poly_coords_size, int32_t *poly_ring_sizes, int32_t poly_num_rings, int32_t ic, int32_t isr, int32_t osr, double *poly_centroid)
EXTENSION_NOINLINE void	ST_Centroid_MultiPolygon (int8_t *mpoly_coords, int32_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int32_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int32_t mpoly_num_polys, int32_t ic, int32_t isr, int32_t osr, double *mpoly_centroid)
EXTENSION_NOINLINE double	ST_Distance_Point_Point (int8_t *p1, int64_t p1size, int8_t *p2, int64_t p2size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE double	ST_Distance_Point_Point_Squared (int8_t *p1, int64_t p1size, int8_t *p2, int64_t p2size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE double	ST_Distance_Point_Point_Geodesic (int8_t *p1, int64_t p1size, int8_t *p2, int64_t p2size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE double	ST_Distance_Point_LineString_Geodesic (int8_t *p, int64_t psize, int8_t *l, int64_t lsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE double	ST_Distance_LineString_Point_Geodesic (int8_t *l, int64_t lsize, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE double	ST_Distance_Point_MultiPoint (int8_t *p, int64_t psize, int8_t *mp, int64_t mpsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_Point_MultiPoint_Squared (int8_t *p, int64_t psize, int8_t *mp, int64_t mpsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
DEVICE ALWAYS_INLINE double	distance_point_linestring (int8_t *p, int64_t psize, int8_t *l, int64_t lsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, bool check_closed, double threshold)
EXTENSION_NOINLINE double	ST_Distance_Point_ClosedLineString (int8_t *p, int64_t psize, int8_t *l, int64_t lsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_Point_LineString (int8_t *p, int64_t psize, int8_t *l, int64_t lsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
DEVICE ALWAYS_INLINE double	distance_point_multilinestring (int8_t *p, int64_t psize, int8_t *mls, int64_t mls_size, int32_t *mls_ls_sizes, int64_t mls_ls_num, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_Point_MultiLineString (int8_t *p, int64_t psize, int8_t *mls, int64_t mls_size, int32_t *mls_ls_sizes, int64_t mls_ls_num, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
DEVICE ALWAYS_INLINE double	distance_point_polygon (int8_t *p, int64_t psize, int8_t *poly, int64_t polysize, int32_t *poly_ring_sizes, int64_t poly_num_rings, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
DEVICE ALWAYS_INLINE double	distance_point_multipolygon (int8_t *p, int64_t psize, int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_Point_Polygon (int8_t *p, int64_t psize, int8_t *poly, int64_t polysize, int32_t *poly_ring_sizes, int64_t poly_num_rings, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_Point_MultiPolygon (int8_t *p, int64_t psize, int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_MultiPoint_Polygon (int8_t *mp, int64_t mpsize, int8_t *poly, int64_t polysize, int32_t *poly_ring_sizes, int64_t poly_num_rings, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_MultiPoint_MultiPolygon (int8_t *mp, int64_t mpsize, int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_Polygon_MultiPoint (int8_t *poly, int64_t polysize, int32_t *poly_ring_sizes, int64_t poly_num_rings, int8_t *mp, int64_t mpsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_MultiPolygon_MultiPoint (int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, int8_t *mp, int64_t mpsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_MultiPoint_Point (int8_t *mp, int64_t mpsize, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_MultiPoint_Point_Squared (int8_t *mp, int64_t mpsize, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_MultiPoint_MultiPoint (int8_t *mp1, int64_t mp1size, int8_t *mp2, int64_t mp2size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_MultiPoint_MultiPoint_Squared (int8_t *mp1, int64_t mp1size, int8_t *mp2, int64_t mp2size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_MultiPoint_LineString (int8_t *mp, int64_t mpsize, int8_t *l, int64_t lsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_MultiPoint_MultiLineString (int8_t *mp, int64_t mpsize, int8_t *mls, int64_t mls_size, int32_t *mls_ls_sizes, int64_t mls_ls_num, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_LineString_Point (int8_t *l, int64_t lsize, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_LineString_MultiPoint (int8_t *l, int64_t lsize, int8_t *mp, int64_t mpsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_MultiLineString_MultiPoint (int8_t *mls, int64_t mls_size, int32_t *mls_ls_sizes, int64_t mls_ls_num, int8_t *mp, int64_t mpsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_LineString_LineString (int8_t *l1, int64_t l1size, int8_t *l2, int64_t l2size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_LineString_Polygon (int8_t *l, int64_t lsize, int8_t *poly_coords, int64_t poly_coords_size, int32_t *poly_ring_sizes, int64_t poly_num_rings, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_LineString_MultiPolygon (int8_t *l, int64_t lsize, int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_MultiLineString_Point (int8_t *mls, int64_t mls_size, int32_t *mls_ls_sizes, int64_t mls_ls_num, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_Polygon_Point (int8_t *poly_coords, int64_t poly_coords_size, int32_t *poly_ring_sizes, int64_t poly_num_rings, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_Polygon_LineString (int8_t *poly_coords, int64_t poly_coords_size, int32_t *poly_ring_sizes, int64_t poly_num_rings, int8_t *l, int64_t lsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_Polygon_Polygon (int8_t *poly1_coords, int64_t poly1_coords_size, int32_t *poly1_ring_sizes, int64_t poly1_num_rings, int8_t *poly2_coords, int64_t poly2_coords_size, int32_t *poly2_ring_sizes, int64_t poly2_num_rings, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_Polygon_MultiPolygon (int8_t *poly1_coords, int64_t poly1_coords_size, int32_t *poly1_ring_sizes, int64_t poly1_num_rings, int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_MultiPolygon_Point (int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_MultiPolygon_LineString (int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, int8_t *l, int64_t lsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE double	ST_Distance_MultiPolygon_Polygon (int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, int8_t *poly1_coords, int64_t poly1_coords_size, int32_t *poly1_ring_sizes, int64_t poly1_num_rings, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_MultiPolygon_MultiPolygon (int8_t *mpoly1_coords, int64_t mpoly1_coords_size, int32_t *mpoly1_ring_sizes, int64_t mpoly1_num_rings, int32_t *mpoly1_poly_sizes, int64_t mpoly1_num_polys, int8_t *mpoly2_coords, int64_t mpoly2_coords_size, int32_t *mpoly2_ring_sizes, int64_t mpoly2_num_rings, int32_t *mpoly2_poly_sizes, int64_t mpoly2_num_polys, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_NOINLINE double	ST_Distance_MultiLineString_MultiLineString (int8_t *mls1, int64_t mls1size, int32_t *mls1_ls_sizes, int64_t mls1_ls_num, int8_t *mls2, int64_t mls2size, int32_t *mls2_ls_sizes, int64_t mls2_ls_num, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double threshold)
EXTENSION_INLINE bool	ST_DWithin_Point_Point (int8_t *p1, int64_t p1size, int8_t *p2, int64_t p2size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double distance_within)
EXTENSION_INLINE bool	ST_DWithin_Point_Point_Geodesic (int8_t *p1, int64_t p1size, int8_t *p2, int64_t p2size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double distance_within)
EXTENSION_INLINE bool	ST_DWithin_Point_LineString (int8_t *p1, int64_t p1size, int8_t *l2, int64_t l2size, double *l2bounds, int64_t l2bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double distance_within)
EXTENSION_INLINE bool	ST_DWithin_Point_Polygon (int8_t *p, int64_t psize, int8_t *poly_coords, int64_t poly_coords_size, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double distance_within)
EXTENSION_INLINE bool	ST_DWithin_Point_MultiPolygon (int8_t *p, int64_t psize, int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double distance_within)
EXTENSION_INLINE bool	ST_DWithin_LineString_LineString (int8_t *l1, int64_t l1size, double *l1bounds, int64_t l1bounds_size, int8_t *l2, int64_t l2size, double *l2bounds, int64_t l2bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double distance_within)
EXTENSION_INLINE bool	ST_DWithin_LineString_Polygon (int8_t *l1, int64_t l1size, double *l1bounds, int64_t l1bounds_size, int8_t *poly_coords, int64_t poly_coords_size, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double distance_within)
EXTENSION_INLINE bool	ST_DWithin_LineString_MultiPolygon (int8_t *l1, int64_t l1size, double *l1bounds, int64_t l1bounds_size, int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double distance_within)
EXTENSION_INLINE bool	ST_DWithin_Polygon_Polygon (int8_t *poly1_coords, int64_t poly1_coords_size, int32_t *poly1_ring_sizes, int64_t poly1_num_rings, double *poly1_bounds, int64_t poly1_bounds_size, int8_t *poly2_coords, int64_t poly2_coords_size, int32_t *poly2_ring_sizes, int64_t poly2_num_rings, double *poly2_bounds, int64_t poly2_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double distance_within)
EXTENSION_INLINE bool	ST_DWithin_Polygon_MultiPolygon (int8_t *poly_coords, int64_t poly_coords_size, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double distance_within)
EXTENSION_INLINE bool	ST_DWithin_MultiPolygon_MultiPolygon (int8_t *mpoly1_coords, int64_t mpoly1_coords_size, int32_t *mpoly1_ring_sizes, int64_t mpoly1_num_rings, int32_t *mpoly1_poly_sizes, int64_t mpoly1_num_polys, double *mpoly1_bounds, int64_t mpoly1_bounds_size, int8_t *mpoly2_coords, int64_t mpoly2_coords_size, int32_t *mpoly2_ring_sizes, int64_t mpoly2_num_rings, int32_t *mpoly2_poly_sizes, int64_t mpoly2_num_polys, double *mpoly2_bounds, int64_t mpoly2_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, double distance_within)
DEVICE double	max_distance_point_line (double px, double py, double l1x, double l1y, double l2x, double l2y)
DEVICE ALWAYS_INLINE double	max_distance_point_linestring (int8_t *p, int64_t psize, int8_t *l, int64_t lsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr, bool check_closed)
EXTENSION_NOINLINE double	ST_MaxDistance_Point_LineString (int8_t *p, int64_t psize, int8_t *l, int64_t lsize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE double	ST_MaxDistance_LineString_Point (int8_t *l, int64_t lsize, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_Contains_Point_Point (int8_t *p1, int64_t p1size, int8_t *p2, int64_t p2size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_Contains_Point_LineString (int8_t *p, int64_t psize, int8_t *l, int64_t lsize, double *lbounds, int64_t lbounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_Contains_Point_Polygon (int8_t *p, int64_t psize, int8_t *poly_coords, int64_t poly_coords_size, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_Contains_LineString_Point (int8_t *l, int64_t lsize, double *lbounds, int64_t lbounds_size, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_Contains_LineString_Polygon (int8_t *l, int64_t lsize, double *lbounds, int64_t lbounds_size, int8_t *poly_coords, int64_t poly_coords_size, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
template<typename T , EdgeBehavior TEdgeBehavior>
DEVICE ALWAYS_INLINE bool	Contains_Polygon_Point_Impl (const int8_t *poly_coords, const int64_t poly_coords_size, const int32_t *poly_ring_sizes, const int64_t poly_num_rings, const double *poly_bounds, const int64_t poly_bounds_size, const int8_t *p, const int64_t psize, const int32_t ic1, const int32_t isr1, const int32_t ic2, const int32_t isr2, const int32_t osr)
EXTENSION_NOINLINE bool	ST_Contains_Polygon_Point (const int8_t *poly_coords, const int64_t poly_coords_size, const int32_t *poly_ring_sizes, const int64_t poly_num_rings, const double *poly_bounds, const int64_t poly_bounds_size, const int8_t *p, const int64_t psize, const int32_t ic1, const int32_t isr1, const int32_t ic2, const int32_t isr2, const int32_t osr)
EXTENSION_NOINLINE bool	ST_cContains_Polygon_Point (const int8_t *poly_coords, const int64_t poly_coords_size, const int32_t *poly_ring_sizes, const int64_t poly_num_rings, const double *poly_bounds, const int64_t poly_bounds_size, const int8_t *p, const int64_t psize, const int32_t ic1, const int32_t isr1, const int32_t ic2, const int32_t isr2, const int32_t osr)
EXTENSION_NOINLINE bool	ST_Contains_Polygon_LineString (int8_t *poly_coords, int64_t poly_coords_size, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int8_t *l, int64_t lsize, double *lbounds, int64_t lbounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_Contains_Polygon_Polygon (int8_t *poly1_coords, int64_t poly1_coords_size, int32_t *poly1_ring_sizes, int64_t poly1_num_rings, double *poly1_bounds, int64_t poly1_bounds_size, int8_t *poly2_coords, int64_t poly2_coords_size, int32_t *poly2_ring_sizes, int64_t poly2_num_rings, double *poly2_bounds, int64_t poly2_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
template<typename T , EdgeBehavior TEdgeBehavior>
DEVICE ALWAYS_INLINE bool	Contains_MultiPolygon_Point_Impl (const int8_t *mpoly_coords, const int64_t mpoly_coords_size, const int32_t *mpoly_ring_sizes, const int64_t mpoly_num_rings, const int32_t *mpoly_poly_sizes, const int64_t mpoly_num_polys, const double *mpoly_bounds, const int64_t mpoly_bounds_size, const int8_t *p, const int64_t psize, const int32_t ic1, const int32_t isr1, const int32_t ic2, const int32_t isr2, const int32_t osr)
EXTENSION_NOINLINE bool	ST_Contains_MultiPolygon_Point (int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_cContains_MultiPolygon_Point (int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_Contains_MultiPolygon_LineString (int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int8_t *l, int64_t lsize, double *lbounds, int64_t lbounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_Intersects_Point_Point (int8_t *p1, int64_t p1size, int8_t *p2, int64_t p2size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_Intersects_Point_LineString (int8_t *p, int64_t psize, int8_t *l, int64_t lsize, double *lbounds, int64_t lbounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_Intersects_Point_Polygon (int8_t *p, int64_t psize, int8_t *poly, int64_t polysize, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_Intersects_Point_MultiPolygon (int8_t *p, int64_t psize, int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_Intersects_LineString_Point (int8_t *l, int64_t lsize, double *lbounds, int64_t lbounds_size, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_Intersects_LineString_Linestring (int8_t *l1, int64_t l1size, double *l1bounds, int64_t l1bounds_size, int8_t *l2, int64_t l2size, double *l2bounds, int64_t l2bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_Intersects_LineString_Polygon (int8_t *l, int64_t lsize, double *lbounds, int64_t lbounds_size, int8_t *poly, int64_t polysize, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_Intersects_LineString_MultiPolygon (int8_t *l, int64_t lsize, double *lbounds, int64_t lbounds_size, int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_Intersects_Polygon_Point (int8_t *poly, int64_t polysize, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_cIntersects_Polygon_Point (int8_t *poly, int64_t polysize, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_Intersects_Polygon_LineString (int8_t *poly, int64_t polysize, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int8_t *l, int64_t lsize, double *lbounds, int64_t lbounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_Intersects_Polygon_Polygon (int8_t *poly1_coords, int64_t poly1_coords_size, int32_t *poly1_ring_sizes, int64_t poly1_num_rings, double *poly1_bounds, int64_t poly1_bounds_size, int8_t *poly2_coords, int64_t poly2_coords_size, int32_t *poly2_ring_sizes, int64_t poly2_num_rings, double *poly2_bounds, int64_t poly2_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_NOINLINE bool	ST_Intersects_Polygon_MultiPolygon (int8_t *poly_coords, int64_t poly_coords_size, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_Intersects_MultiPolygon_Point (int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_cIntersects_MultiPolygon_Point (int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int8_t *p, int64_t psize, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_Intersects_MultiPolygon_LineString (int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int8_t *l, int64_t lsize, double *lbounds, int64_t lbounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_Intersects_MultiPolygon_Polygon (int8_t *mpoly_coords, int64_t mpoly_coords_size, int32_t *mpoly_ring_sizes, int64_t mpoly_num_rings, int32_t *mpoly_poly_sizes, int64_t mpoly_num_polys, double *mpoly_bounds, int64_t mpoly_bounds_size, int8_t *poly_coords, int64_t poly_coords_size, int32_t *poly_ring_sizes, int64_t poly_num_rings, double *poly_bounds, int64_t poly_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE bool	ST_Intersects_MultiPolygon_MultiPolygon (int8_t *mpoly1_coords, int64_t mpoly1_coords_size, int32_t *mpoly1_ring_sizes, int64_t mpoly1_num_rings, int32_t *mpoly1_poly_sizes, int64_t mpoly1_num_polys, double *mpoly1_bounds, int64_t mpoly1_bounds_size, int8_t *mpoly2_coords, int64_t mpoly2_coords_size, int32_t *mpoly2_ring_sizes, int64_t mpoly2_num_rings, int32_t *mpoly2_poly_sizes, int64_t mpoly2_num_polys, double *mpoly2_bounds, int64_t mpoly2_bounds_size, int32_t ic1, int32_t isr1, int32_t ic2, int32_t isr2, int32_t osr)
EXTENSION_INLINE int64_t	HeavyDB_Geo_PolyBoundsPtr (double *bounds, int64_t size)
EXTENSION_NOINLINE double	convert_meters_to_pixel_width (const double meters, int8_t *p, const int64_t psize, const int32_t ic, const int32_t isr, const int32_t osr, const double min_lon, const double max_lon, const int32_t img_width, const double min_width)
EXTENSION_NOINLINE double	convert_meters_to_pixel_height (const double meters, int8_t *p, const int64_t psize, const int32_t ic, const int32_t isr, const int32_t osr, const double min_lat, const double max_lat, const int32_t img_height, const double min_height)
EXTENSION_NOINLINE bool	is_point_in_view (int8_t *p, const int64_t psize, const int32_t ic, const double min_lon, const double max_lon, const double min_lat, const double max_lat)
EXTENSION_NOINLINE bool	is_point_size_in_view (int8_t *p, const int64_t psize, const int32_t ic, const double meters, const double min_lon, const double max_lon, const double min_lat, const double max_lat)
Geospatial::GeoPoint	to_Geospatial_GeoPoint (GeoPoint &p)
Geospatial::GeoMultiPoint	to_Geospatial_GeoMultiPoint (GeoMultiPoint &mp)
Geospatial::GeoLineString	to_Geospatial_GeoLineString (GeoLineString &l)
Geospatial::GeoMultiLineString	to_Geospatial_GeoMultiLineString (GeoMultiLineString &ml)
Geospatial::GeoPolygon	to_Geospatial_GeoPolygon (GeoPolygon &p)
Geospatial::GeoMultiPolygon	to_Geospatial_GeoMultiPolygon (GeoMultiPolygon &mp)
EXTENSION_NOINLINE TextEncodingNone	ST_AsText__GeoPoint__cpu_ (RowFunctionManager &mgr, GeoPoint &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsText__GeoMultiPoint__cpu_ (RowFunctionManager &mgr, GeoMultiPoint &mp)
EXTENSION_NOINLINE TextEncodingNone	ST_AsText__GeoLineString__cpu_ (RowFunctionManager &mgr, GeoLineString &l)
EXTENSION_NOINLINE TextEncodingNone	ST_AsText__GeoMultiLineString__cpu_ (RowFunctionManager &mgr, GeoMultiLineString &ml)
EXTENSION_NOINLINE TextEncodingNone	ST_AsText__GeoPolygon__cpu_ (RowFunctionManager &mgr, GeoPolygon &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsText__GeoMultiPolygon__cpu_ (RowFunctionManager &mgr, GeoMultiPolygon &mp)
EXTENSION_NOINLINE TextEncodingNone	ST_AsWkt__GeoPoint__cpu_ (RowFunctionManager &mgr, GeoPoint &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsWkt__GeoMultiPoint__cpu_ (RowFunctionManager &mgr, GeoMultiPoint &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsWkt__GeoLineString__cpu_ (RowFunctionManager &mgr, GeoLineString &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsWkt__GeoMultiLineString__cpu_ (RowFunctionManager &mgr, GeoMultiLineString &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsWkt__GeoPolygon__cpu_ (RowFunctionManager &mgr, GeoPolygon &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsWkt__GeoMultiPolygon__cpu_ (RowFunctionManager &mgr, GeoMultiPolygon &p)
template<typename T >
ALWAYS_INLINE std::string	__wkb_to_str (T &geometry)
EXTENSION_NOINLINE TextEncodingNone	ST_AsBinary__GeoPoint__cpu_ (RowFunctionManager &mgr, GeoPoint &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsBinary__GeoMultiPoint__cpu_ (RowFunctionManager &mgr, GeoMultiPoint &mp)
EXTENSION_NOINLINE TextEncodingNone	ST_AsBinary__GeoLineString__cpu_ (RowFunctionManager &mgr, GeoLineString &l)
EXTENSION_NOINLINE TextEncodingNone	ST_AsBinary__GeoMultiLineString__cpu_ (RowFunctionManager &mgr, GeoMultiLineString &ml)
EXTENSION_NOINLINE TextEncodingNone	ST_AsBinary__GeoPolygon__cpu_ (RowFunctionManager &mgr, GeoPolygon &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsBinary__GeoMultiPolygon__cpu_ (RowFunctionManager &mgr, GeoMultiPolygon &mp)
EXTENSION_NOINLINE TextEncodingNone	ST_AsWkb__GeoPoint__cpu_ (RowFunctionManager &mgr, GeoPoint &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsWkb__GeoMultiPoint__cpu_ (RowFunctionManager &mgr, GeoMultiPoint &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsWkb__GeoLineString__cpu_ (RowFunctionManager &mgr, GeoLineString &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsWkb__GeoMultiLineString__cpu_ (RowFunctionManager &mgr, GeoMultiLineString &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsWkb__GeoPolygon__cpu_ (RowFunctionManager &mgr, GeoPolygon &p)
EXTENSION_NOINLINE TextEncodingNone	ST_AsWkb__GeoMultiPolygon__cpu_ (RowFunctionManager &mgr, GeoMultiPolygon &p)

Macro Definition Documentation

#define DEBUG_STMT

(

x

)

Definition at line 21 of file ExtensionFunctionsGeo.hpp.

Referenced by Contains_Polygon_Point_Impl(), and point_in_polygon_winding_number().

Enumeration Type Documentation

enum EdgeBehavior

strong

Enumerator
kIncludePointOnEdge
kExcludePointOnEdge

Definition at line 571 of file ExtensionFunctionsGeo.hpp.

{ kIncludePointOnEdge, kExcludePointOnEdge };

Function Documentation

template<typename T >

ALWAYS_INLINE std::string __wkb_to_str

(

T &

geometry

)

Definition at line 5296 of file ExtensionFunctionsGeo.hpp.

Referenced by ST_AsBinary__GeoLineString__cpu_(), ST_AsBinary__GeoMultiLineString__cpu_(), ST_AsBinary__GeoMultiPoint__cpu_(), ST_AsBinary__GeoMultiPolygon__cpu_(), ST_AsBinary__GeoPoint__cpu_(), and ST_AsBinary__GeoPolygon__cpu_().

                                                  {
  std::vector<uint8_t> bytes;

  bool ok = geometry.getWkb(bytes);
  if (!ok) {
    return std::string();
  }

  std::ostringstream ss;
  for (auto b : bytes) {
    ss << std::setfill('0') << std::setw(2) << std::hex << static_cast<int>(b);
  }
  return ss.str();
}

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double area_polygon	(	int8_t *	poly_coords,
		int32_t	poly_coords_size,
		int8_t *	poly_ring_sizes_in,
		int32_t	poly_num_rings,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1381 of file ExtensionFunctionsGeo.hpp.

References area_ring(), and compression_unit_size().

Referenced by ST_Area_MultiPolygon(), and ST_Area_Polygon().

                                                      {
  if (poly_num_rings <= 0) {
    return 0.0;
  }

  double area = 0.0;
  auto ring_coords = poly_coords;
  auto poly_ring_sizes = reinterpret_cast<int32_t*>(poly_ring_sizes_in);

  // Add up the areas of all rings.
  // External ring is CCW, open - positive area.
  // Internal rings (holes) are CW, open - negative areas.
  for (auto r = 0; r < poly_num_rings; r++) {
    auto ring_coords_size = poly_ring_sizes[r] * 2 * compression_unit_size(ic);
    area += area_ring(ring_coords, ring_coords_size, ic, isr, osr);
    // Advance to the next ring.
    ring_coords += ring_coords_size;
  }
  return area;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double area_ring	(	int8_t *	ring,
		int64_t	ringsize,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1358 of file ExtensionFunctionsGeo.hpp.

References area_triangle(), compression_unit_size(), and get_point().

Referenced by area_polygon().

                                                   {
  auto ring_num_coords = ringsize / compression_unit_size(ic);

  if (ring_num_coords < 6) {
    return 0.0;
  }

  double area = 0.0;

  Point2D p1 = get_point(ring, 0, ic, isr, osr);
  Point2D p2 = get_point(ring, 2, ic, isr, osr);
  for (int32_t i = 4; i < ring_num_coords; i += 2) {
    Point2D p3 = get_point(ring, i, ic, isr, osr);
    area += area_triangle(p1.x, p1.y, p2.x, p2.y, p3.x, p3.y);
    p2 = p3;
  }
  return area;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double area_triangle	(	double	x1,
		double	y1,
		double	x2,
		double	y2,
		double	x3,
		double	y3
	)

Definition at line 1349 of file ExtensionFunctionsGeo.hpp.

Referenced by area_ring().

                                                     {
  return (x1 * y2 - x2 * y1 + x3 * y1 - x1 * y3 + x2 * y3 - x3 * y2) / 2.0;
}

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool box_contains_box	(	double *	bounds1,
		int64_t	bounds1_size,
		double *	bounds2,
		int64_t	bounds2_size
	)

Definition at line 878 of file ExtensionFunctionsGeo.hpp.

References box_contains_point().

Referenced by ST_Contains_MultiPolygon_LineString(), ST_Contains_Polygon_LineString(), and ST_Contains_Polygon_Polygon().

                                                                 {
  // TODO: mixed spatial references
  return (
      box_contains_point(
          bounds1, bounds1_size, bounds2[0], bounds2[1]) &&  // box1 <- box2: xmin, ymin
      box_contains_point(
          bounds1, bounds1_size, bounds2[2], bounds2[3]));  // box1 <- box2: xmax, ymax
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool box_contains_box_vertex	(	double *	bounds1,
		int64_t	bounds1_size,
		double *	bounds2,
		int64_t	bounds2_size
	)

Definition at line 890 of file ExtensionFunctionsGeo.hpp.

References box_contains_point().

                                                                        {
  // TODO: mixed spatial references
  return (
      box_contains_point(
          bounds1, bounds1_size, bounds2[0], bounds2[1]) ||  // box1 <- box2: xmin, ymin
      box_contains_point(
          bounds1, bounds1_size, bounds2[2], bounds2[3]) ||  // box1 <- box2: xmax, ymax
      box_contains_point(
          bounds1, bounds1_size, bounds2[0], bounds2[3]) ||  // box1 <- box2: xmin, ymax
      box_contains_point(
          bounds1, bounds1_size, bounds2[2], bounds2[1]));  // box1 <- box2: xmax, ymin
}

Here is the call graph for this function:

DEVICE ALWAYS_INLINE bool box_contains_point	(	const double *	bounds,
		const int64_t	bounds_size,
		const double	px,
		const double	py
	)

Definition at line 861 of file ExtensionFunctionsGeo.hpp.

References tol_ge(), and tol_le().

Referenced by box_contains_box(), box_contains_box_vertex(), Contains_MultiPolygon_Point_Impl(), Contains_Polygon_Point_Impl(), Point_Overlaps_Box(), and ST_Intersects_Point_LineString().

                                                              {
  // TODO: mixed spatial references
  return (tol_ge(px, bounds[0]) && tol_ge(py, bounds[1]) && tol_le(px, bounds[2]) &&
          tol_le(py, bounds[3]));
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool box_dwithin_box	(	double *	bounds1,
		int64_t	bounds1_size,
		int32_t	isr1,
		double *	bounds2,
		int64_t	bounds2_size,
		int32_t	isr2,
		int32_t	osr,
		double	distance
	)

Definition at line 946 of file ExtensionFunctionsGeo.hpp.

Referenced by ST_DWithin_LineString_LineString(), ST_DWithin_LineString_MultiPolygon(), ST_DWithin_LineString_Polygon(), ST_DWithin_MultiPolygon_MultiPolygon(), ST_DWithin_Polygon_MultiPolygon(), and ST_DWithin_Polygon_Polygon().

                                                           {
  // TODO: tolerance

  // Bounds come in corresponding input spatial references.
  // For example, here the first bounding box may come in 4326, second in 900913:
  //   ST_DWithin(ST_Transform(linestring4326, 900913), linestring900913, 10)
  // Distance is given in osr spatial reference units, in this case 10 meters.
  // Bounds should be converted to osr before calculations, if isr != osr.

  // TODO: revise all other functions that process bounds

  Point2D const lb1 = transform_point<XY>({bounds1[0], bounds1[1]}, isr1, osr);
  Point2D const rt1 = transform_point<XY>({bounds1[2], bounds1[3]}, isr1, osr);
  Point2D const lb2 = transform_point<XY>({bounds2[0], bounds2[1]}, isr2, osr);
  Point2D const rt2 = transform_point<XY>({bounds2[2], bounds2[3]}, isr2, osr);
  if (
      // box1 is left of box2:  box1.xmax + distance < box2.xmin
      rt1.x + distance < lb2.x ||
      // box1 is right of box2: box1.xmin - distance > box2.xmax
      lb1.x - distance > rt2.x ||
      // box1 is below box2:    box1.ymax + distance < box2.ymin
      rt1.y + distance < lb2.y ||
      // box1 is above box2:    box1.ymin - distance > box1.ymax
      lb1.y - distance > rt2.y) {
    return false;
  }
  return true;
}

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool box_overlaps_box	(	double *	bounds1,
		int64_t	bounds1_size,
		double *	bounds2,
		int64_t	bounds2_size
	)

Definition at line 906 of file ExtensionFunctionsGeo.hpp.

Referenced by ST_Intersects_LineString_Linestring(), ST_Intersects_LineString_MultiPolygon(), ST_Intersects_LineString_Polygon(), ST_Intersects_MultiPolygon_MultiPolygon(), ST_Intersects_Polygon_MultiPolygon(), and ST_Intersects_Polygon_Polygon().

                                                                 {
  // TODO: tolerance
  // TODO: mixed spatial references
  if (bounds1[2] < bounds2[0] ||  // box1 is left of box2:  box1.xmax < box2.xmin
      bounds1[0] > bounds2[2] ||  // box1 is right of box2: box1.xmin > box2.xmax
      bounds1[3] < bounds2[1] ||  // box1 is below box2:    box1.ymax < box2.ymin
      bounds1[1] > bounds2[3]) {  // box1 is above box2:    box1.ymin > box1.ymax
    return false;
  }
  return true;
}

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool centroid_add_linestring	(	int8_t *	l,
		int64_t	lsize,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr,
		bool	closed,
		double *	total_length,
		double *	linestring_centroid_sum,
		int64_t *	num_points,
		double *	point_centroid_sum
	)

Definition at line 1527 of file ExtensionFunctionsGeo.hpp.

References centroid_add_segment(), compression_unit_size(), and get_point().

Referenced by centroid_add_ring(), and ST_Centroid_LineString().

                                                                              {
  auto l_num_coords = lsize / compression_unit_size(ic);
  double length = 0.0;
  Point2D const l0 = get_point(l, 0, ic, isr, osr);
  Point2D l2 = l0;
  for (int32_t i = 2; i < l_num_coords; i += 2) {
    Point2D const l1 = l2;
    l2 = get_point(l, i, ic, isr, osr);
    centroid_add_segment(l1.x, l1.y, l2.x, l2.y, &length, linestring_centroid_sum);
  }
  if (l_num_coords > 4 && closed) {
    // Also add the closing segment between the last and the first points
    centroid_add_segment(l2.x, l2.y, l0.x, l0.y, &length, linestring_centroid_sum);
  }
  *total_length += length;
  if (length == 0.0 && l_num_coords > 0) {
    *num_points += 1;
    point_centroid_sum[0] += l0.x;
    point_centroid_sum[1] += l0.y;
  }
  return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool centroid_add_polygon	(	int8_t *	poly_coords,
		int64_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr,
		double *	total_area2,
		double *	cg3,
		double *	total_length,
		double *	linestring_centroid_sum,
		int64_t *	num_points,
		double *	point_centroid_sum
	)

Definition at line 1646 of file ExtensionFunctionsGeo.hpp.

References centroid_add_ring(), and compression_unit_size().

Referenced by ST_Centroid_MultiPolygon(), and ST_Centroid_Polygon().

                                                                           {
  if (poly_num_rings <= 0) {
    return false;
  }

  auto ring_coords = poly_coords;

  for (auto r = 0; r < poly_num_rings; r++) {
    auto ring_coords_size = poly_ring_sizes[r] * 2 * compression_unit_size(ic);
    // Shape - positive area, holes - negative areas
    double sign = (r == 0) ? 1.0 : -1.0;
    centroid_add_ring(ring_coords,
                      ring_coords_size,
                      ic,
                      isr,
                      osr,
                      sign,
                      total_area2,
                      cg3,
                      total_length,
                      linestring_centroid_sum,
                      num_points,
                      point_centroid_sum);
    // Advance to the next ring.
    ring_coords += ring_coords_size;
  }
  return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool centroid_add_ring	(	int8_t *	ring,
		int64_t	ringsize,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr,
		double	sign,
		double *	total_area2,
		double *	cg3,
		double *	total_length,
		double *	linestring_centroid_sum,
		int64_t *	num_points,
		double *	point_centroid_sum
	)

Definition at line 1607 of file ExtensionFunctionsGeo.hpp.

References centroid_add_linestring(), centroid_add_triangle(), compression_unit_size(), and get_point().

Referenced by centroid_add_polygon().

                                                                        {
  auto ring_num_coords = ringsize / compression_unit_size(ic);

  if (ring_num_coords < 6) {
    return false;
  }

  Point2D p1 = get_point(ring, 0, ic, isr, osr);
  Point2D p2 = get_point(ring, 2, ic, isr, osr);
  for (int32_t i = 4; i < ring_num_coords; i += 2) {
    Point2D p3 = get_point(ring, i, ic, isr, osr);
    centroid_add_triangle(p1.x, p1.y, p2.x, p2.y, p3.x, p3.y, sign, total_area2, cg3);
    p2 = p3;
  }

  centroid_add_linestring(ring,
                          ringsize,
                          ic,
                          isr,
                          osr,
                          true,  // closed
                          total_length,
                          linestring_centroid_sum,
                          num_points,
                          point_centroid_sum);
  return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool centroid_add_segment	(	double	x1,
		double	y1,
		double	x2,
		double	y2,
		double *	length,
		double *	linestring_centroid_sum
	)

Definition at line 1512 of file ExtensionFunctionsGeo.hpp.

References distance_point_point().

Referenced by centroid_add_linestring().

                                                                                {
  double ldist = distance_point_point(x1, y1, x2, y2);
  *length += ldist;
  double segment_midpoint_x = (x1 + x2) / 2.0;
  double segment_midpoint_y = (y1 + y2) / 2.0;
  linestring_centroid_sum[0] += ldist * segment_midpoint_x;
  linestring_centroid_sum[1] += ldist * segment_midpoint_y;
  return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool centroid_add_triangle	(	double	x1,
		double	y1,
		double	x2,
		double	y2,
		double	x3,
		double	y3,
		double	sign,
		double *	total_area2,
		double *	cg3
	)

Definition at line 1589 of file ExtensionFunctionsGeo.hpp.

Referenced by centroid_add_ring().

                                                             {
  double cx = x1 + x2 + x3;
  double cy = y1 + y2 + y3;
  double area2 = x1 * y2 - x2 * y1 + x3 * y1 - x1 * y3 + x2 * y3 - x3 * y2;
  cg3[0] += sign * area2 * cx;
  cg3[1] += sign * area2 * cy;
  *total_area2 += sign * area2;
  return true;
}

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE int32_t compressed_coord	(	const int8_t *	data,
		const int32_t	index
	)

Definition at line 226 of file ExtensionFunctionsGeo.hpp.

Referenced by Geospatial::compress_latitude_coord_geoint32(), Geospatial::compress_longitude_coord_geoint32(), Contains_Polygon_Point_Impl(), decompress_coord(), and point_in_polygon_winding_number().

                                                                                       {
  const auto compressed_coords = reinterpret_cast<const int32_t*>(data);
  return compressed_coords[index];
}

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE int32_t compression_unit_size

(

const int32_t

ic

)

Definition at line 105 of file ExtensionFunctionsGeo.hpp.

References COMPRESSION_GEOINT32.

Referenced by area_polygon(), area_ring(), centroid_add_linestring(), centroid_add_polygon(), centroid_add_ring(), Contains_MultiPolygon_Point_Impl(), Contains_Polygon_Point_Impl(), distance_point_linestring(), distance_point_multilinestring(), distance_point_multipolygon(), distance_point_polygon(), length_linestring(), max_distance_point_linestring(), perimeter_multipolygon(), ST_Area_MultiPolygon(), ST_Centroid_MultiPoint(), ST_Centroid_MultiPolygon(), ST_Contains_MultiPolygon_LineString(), ST_Contains_Point_LineString(), ST_Contains_Point_Polygon(), ST_Contains_Polygon_LineString(), ST_Contains_Polygon_Polygon(), ST_Distance_LineString_LineString(), ST_Distance_LineString_MultiPolygon(), ST_Distance_LineString_Polygon(), ST_Distance_MultiPoint_LineString(), ST_Distance_MultiPoint_MultiLineString(), ST_Distance_MultiPoint_MultiPoint(), ST_Distance_MultiPoint_MultiPoint_Squared(), ST_Distance_MultiPoint_MultiPolygon(), ST_Distance_MultiPoint_Polygon(), ST_Distance_MultiPolygon_MultiPolygon(), ST_Distance_Point_LineString_Geodesic(), ST_Distance_Point_MultiPoint(), ST_Distance_Point_MultiPoint_Squared(), ST_Distance_Polygon_MultiPolygon(), ST_Distance_Polygon_Polygon(), ST_DWithin_LineString_MultiPolygon(), ST_DWithin_LineString_Polygon(), ST_Length_MultiLineString(), ST_Perimeter_Polygon(), ST_Perimeter_Polygon_Geodesic(), ST_XMax(), ST_XMin(), ST_YMax(), ST_YMin(), and trim_linestring_to_buffered_box().

                                                                     {
  if (ic == COMPRESSION_GEOINT32) {
    return 4;
  }
  return 8;
}

Here is the caller graph for this function:

template<typename T , EdgeBehavior TEdgeBehavior>

DEVICE ALWAYS_INLINE bool Contains_MultiPolygon_Point_Impl	(	const int8_t *	mpoly_coords,
		const int64_t	mpoly_coords_size,
		const int32_t *	mpoly_ring_sizes,
		const int64_t	mpoly_num_rings,
		const int32_t *	mpoly_poly_sizes,
		const int64_t	mpoly_num_polys,
		const double *	mpoly_bounds,
		const int64_t	mpoly_bounds_size,
		const int8_t *	p,
		const int64_t	psize,
		const int32_t	ic1,
		const int32_t	isr1,
		const int32_t	ic2,
		const int32_t	isr2,
		const int32_t	osr
	)

Definition at line 4187 of file ExtensionFunctionsGeo.hpp.

References box_contains_point(), compression_unit_size(), and get_point().

                       {
  if (mpoly_num_polys <= 0) {
    return false;
  }

  // TODO: mpoly_bounds could contain individual bounding boxes too:
  // first two points - box for the entire multipolygon, then a pair for each polygon
  if (mpoly_bounds) {
    Point2D const pt = get_point(p, 0, ic2, isr2, osr);
    if (!box_contains_point(mpoly_bounds, mpoly_bounds_size, pt.x, pt.y)) {
      return false;
    }
  }

  // Set specific poly pointers as we move through the coords/ringsizes/polyrings
  // arrays.
  auto next_poly_coords = mpoly_coords;
  auto next_poly_ring_sizes = mpoly_ring_sizes;

  for (auto poly = 0; poly < mpoly_num_polys; poly++) {
    const auto poly_coords = next_poly_coords;
    const auto poly_ring_sizes = next_poly_ring_sizes;
    const auto poly_num_rings = mpoly_poly_sizes[poly];
    // Count number of coords in all of poly's rings, advance ring size pointer.
    int32_t poly_num_coords = 0;
    for (auto ring = 0; ring < poly_num_rings; ring++) {
      poly_num_coords += 2 * *next_poly_ring_sizes++;
    }
    const auto poly_coords_size = poly_num_coords * compression_unit_size(ic1);
    next_poly_coords += poly_coords_size;
    // TODO: pass individual bounding boxes for each polygon
    if (Contains_Polygon_Point_Impl<T, TEdgeBehavior>(poly_coords,
                                                      poly_coords_size,
                                                      poly_ring_sizes,
                                                      poly_num_rings,
                                                      nullptr,
                                                      0,
                                                      p,
                                                      psize,
                                                      ic1,
                                                      isr1,
                                                      ic2,
                                                      isr2,
                                                      osr)) {
      return true;
    }
  }

  return false;
}

Here is the call graph for this function:

template<typename T , EdgeBehavior TEdgeBehavior>

DEVICE ALWAYS_INLINE bool Contains_Polygon_Point_Impl	(	const int8_t *	poly_coords,
		const int64_t	poly_coords_size,
		const int32_t *	poly_ring_sizes,
		const int64_t	poly_num_rings,
		const double *	poly_bounds,
		const int64_t	poly_bounds_size,
		const int8_t *	p,
		const int64_t	psize,
		const int32_t	ic1,
		const int32_t	isr1,
		const int32_t	ic2,
		const int32_t	isr2,
		const int32_t	osr
	)

Definition at line 3974 of file ExtensionFunctionsGeo.hpp.

References box_contains_point(), compressed_coord(), compression_unit_size(), DEBUG_STMT, get_point(), and heavydb.dtypes::T.

                                                                         {
  if (poly_bounds) {
    // TODO: codegen
    Point2D const pt = get_point(p, 0, ic2, isr2, osr);
    if (!box_contains_point(poly_bounds, poly_bounds_size, pt.x, pt.y)) {
      DEBUG_STMT(printf("Bounding box does not contain point, exiting.\n"));
      return false;
    }
  }

  auto get_x_coord = [=]() -> T {
    if constexpr (std::is_floating_point<T>::value) {
      return get_point(p, 0, ic2, isr2, osr).x;
    } else {
      return compressed_coord(p, 0);
    }
    return T{};  // https://stackoverflow.com/a/64561686/2700898
  };

  auto get_y_coord = [=]() -> T {
    if constexpr (std::is_floating_point<T>::value) {
      return get_point(p, 0, ic2, isr2, osr).y;
    } else {
      return compressed_coord(p, 1);
    }
    return T{};  // https://stackoverflow.com/a/64561686/2700898
  };

  const T px = get_x_coord();
  const T py = get_y_coord();
  DEBUG_STMT(printf("pt: %f, %f\n", (double)px, (double)py));

  const auto poly_num_coords = poly_coords_size / compression_unit_size(ic1);
  auto exterior_ring_num_coords =
      poly_num_rings > 0 ? poly_ring_sizes[0] * 2 : poly_num_coords;

  auto poly = poly_coords;
  if (point_in_polygon_winding_number<T, TEdgeBehavior>(
          poly, exterior_ring_num_coords, px, py, ic1, isr1, osr)) {
    // Inside exterior ring
    poly += exterior_ring_num_coords * compression_unit_size(ic1);
    // Check that none of the polygon's holes contain that point
    for (auto r = 1; r < poly_num_rings; r++) {
      const int64_t interior_ring_num_coords = poly_ring_sizes[r] * 2;
      if (point_in_polygon_winding_number<T, TEdgeBehavior>(
              poly, interior_ring_num_coords, px, py, ic1, isr1, osr)) {
        return false;
      }
      poly += interior_ring_num_coords * compression_unit_size(ic1);
    }
    return true;
  }
  return false;
}

Here is the call graph for this function:

template<Coords C>

DEVICE ALWAYS_INLINE Point2D conv_4326_900913

(

const Point2D

point

)

Definition at line 113 of file ExtensionFunctionsGeo.hpp.

References conv_4326_900913_x(), conv_4326_900913_y(), anonymous_namespace{ExtensionFunctionsGeo.hpp}::X, and anonymous_namespace{ExtensionFunctionsGeo.hpp}::Y.

                                                                   {
  Point2D retval;
  if constexpr (static_cast<bool>(C & Coords::X)) {
    retval.x = conv_4326_900913_x(point.x);
  }
  if constexpr (static_cast<bool>(C & Coords::Y)) {
    retval.y = conv_4326_900913_y(point.y);
  }
  return retval;
}

Here is the call graph for this function:

template<Coords C>

DEVICE ALWAYS_INLINE Point2D conv_4326_utm	(	const Point2D	point,
		const int32_t	utm_srid
	)

Definition at line 125 of file ExtensionFunctionsGeo.hpp.

References Transform4326ToUTM::calculateX(), Transform4326ToUTM::calculateY(), anonymous_namespace{ExtensionFunctionsGeo.hpp}::X, and anonymous_namespace{ExtensionFunctionsGeo.hpp}::Y.

                                                                                        {
  Point2D retval;
  Transform4326ToUTM const utm(utm_srid, point.x, point.y);
  if constexpr (static_cast<bool>(C & Coords::X)) {
    retval.x = utm.calculateX();
  }
  if constexpr (static_cast<bool>(C & Coords::Y)) {
    retval.y = utm.calculateY();
  }
  return retval;
}

Here is the call graph for this function:

template<Coords C>

DEVICE ALWAYS_INLINE Point2D conv_utm_4326	(	const Point2D	point,
		const int32_t	utm_srid
	)

Definition at line 138 of file ExtensionFunctionsGeo.hpp.

References TransformUTMTo4326::calculateX(), TransformUTMTo4326::calculateY(), anonymous_namespace{ExtensionFunctionsGeo.hpp}::X, and anonymous_namespace{ExtensionFunctionsGeo.hpp}::Y.

                                                                                        {
  Point2D retval;
  TransformUTMTo4326 const wgs84(utm_srid, point.x, point.y);
  if constexpr (static_cast<bool>(C & Coords::X)) {
    retval.x = wgs84.calculateX();
  }
  if constexpr (static_cast<bool>(C & Coords::Y)) {
    retval.y = wgs84.calculateY();
  }
  return retval;
}

Here is the call graph for this function:

EXTENSION_NOINLINE double convert_meters_to_pixel_height	(	const double	meters,
		int8_t *	p,
		const int64_t	psize,
		const int32_t	ic,
		const int32_t	isr,
		const int32_t	osr,
		const double	min_lat,
		const double	max_lat,
		const int32_t	img_height,
		const double	min_height
	)

Definition at line 5033 of file ExtensionFunctionsGeo.hpp.

Referenced by com.mapd.calcite.parser.HeavyDBSqlOperatorTable::addUDF().

                                                               {
  const double const1 = 0.017453292519943295769236907684886;
  const double const2 = 6372797.560856;
  const auto lat = decompress_coord<Y>(p, 0, ic);
  const double latdiff = meters / (const1 * const2);
  const double newlat =
      (lat < 0) ? lat + latdiff : lat - latdiff;  // assumes a lat range of [-90, 90]
  double t1 = transform_point<Y>({{}, lat}, isr, osr).y;
  double t2 = transform_point<Y>({{}, newlat}, isr, osr).y;
  const double min_domain_y = transform_point<Y>({{}, min_lat}, isr, osr).y;
  const double max_domain_y = transform_point<Y>({{}, max_lat}, isr, osr).y;
  const double domain_diff = max_domain_y - min_domain_y;
  t1 = ((t1 - min_domain_y) / domain_diff) * static_cast<double>(img_height);
  t2 = ((t2 - min_domain_y) / domain_diff) * static_cast<double>(img_height);

  // TODO(croot): need to account for edge cases, such as getting close to the poles.
  const double sz = fabs(t1 - t2);
  return (sz < min_height ? min_height : sz);
}

Here is the caller graph for this function:

EXTENSION_NOINLINE double convert_meters_to_pixel_width	(	const double	meters,
		int8_t *	p,
		const int64_t	psize,
		const int32_t	ic,
		const int32_t	isr,
		const int32_t	osr,
		const double	min_lon,
		const double	max_lon,
		const int32_t	img_width,
		const double	min_width
	)

Definition at line 5001 of file ExtensionFunctionsGeo.hpp.

Referenced by com.mapd.calcite.parser.HeavyDBSqlOperatorTable::addUDF().

                                                             {
  const double const1 = 0.017453292519943295769236907684886;
  const double const2 = 6372797.560856;
  const auto lon = decompress_coord<X>(p, 0, ic);
  const auto lat = decompress_coord<Y>(p, 0, ic);
  double t1 = sinf(meters / (2.0 * const2));
  double t2 = cosf(const1 * lat);
  const double newlon = lon - (2.0 * asinf(t1 / t2)) / const1;
  t1 = transform_point<X>({lon, {}}, isr, osr).x;
  t2 = transform_point<X>({newlon, {}}, isr, osr).x;
  const double min_domain_x = transform_point<X>({min_lon, {}}, isr, osr).x;
  const double max_domain_x = transform_point<X>({max_lon, {}}, isr, osr).x;
  const double domain_diff = max_domain_x - min_domain_x;
  t1 = ((t1 - min_domain_x) / domain_diff) * static_cast<double>(img_width);
  t2 = ((t2 - min_domain_x) / domain_diff) * static_cast<double>(img_width);

  // TODO(croot): need to account for edge cases, such as getting close to the poles.
  const double sz = fabs(t1 - t2);
  return (sz < min_width ? min_width : sz);
}

Here is the caller graph for this function:

template<Coords C>

DEVICE ALWAYS_INLINE Point2D coord	(	const int8_t *	data,
		const int32_t	x_index,
		const int32_t	ic,
		const int32_t	isr,
		const int32_t	osr
	)

Definition at line 205 of file ExtensionFunctionsGeo.hpp.

References anonymous_namespace{ExtensionFunctionsGeo.hpp}::X, x_and_y_are_dependent(), and anonymous_namespace{ExtensionFunctionsGeo.hpp}::Y.

Referenced by Geospatial::compress_coords().

                                                      {
  Point2D decompressed;
  static_assert(C == Coords::X || C == Coords::Y, "Use get_point() instead of XY.");
  if (x_and_y_are_dependent(isr, osr)) {
    decompressed = {decompress_coord<X>(data, x_index, ic),
                    decompress_coord<Y>(data, x_index, ic)};
  } else {
    if constexpr (C == Coords::X) {
      decompressed.x = decompress_coord<X>(data, x_index, ic);
    } else {
      decompressed.y = decompress_coord<Y>(data, x_index, ic);
    }
  }
  return transform_point<C>(decompressed, isr, osr);
}

Here is the call graph for this function:

Here is the caller graph for this function:

template<Coords C>

DEVICE ALWAYS_INLINE double decompress_coord	(	const int8_t *	data,
		const int32_t	x_index,
		const int32_t	ic
	)

Definition at line 87 of file ExtensionFunctionsGeo.hpp.

References compressed_coord(), COMPRESSION_GEOINT32, Geospatial::decompress_latitude_coord_geoint32(), Geospatial::decompress_longitude_coord_geoint32(), anonymous_namespace{ExtensionFunctionsGeo.hpp}::X, and anonymous_namespace{ExtensionFunctionsGeo.hpp}::Y.

                                                               {
  static_assert(C == Coords::X || C == Coords::Y);
  const int32_t adjusted_index = x_index + (C == Coords::Y);
  if (ic == COMPRESSION_GEOINT32) {
    auto compressed_coords = reinterpret_cast<const int32_t*>(data);
    auto compressed_coord = compressed_coords[adjusted_index];
    if constexpr (C == Coords::X) {
      return Geospatial::decompress_longitude_coord_geoint32(compressed_coord);
    } else {
      return Geospatial::decompress_latitude_coord_geoint32(compressed_coord);
    }
  }
  auto double_coords = reinterpret_cast<const double*>(data);
  return double_coords[adjusted_index];
}

Here is the call graph for this function:

DEVICE double distance_line_line	(	double	l11x,
		double	l11y,
		double	l12x,
		double	l12y,
		double	l21x,
		double	l21y,
		double	l22x,
		double	l22y
	)

Definition at line 441 of file ExtensionFunctionsGeo.hpp.

References distance_point_line(), and line_intersects_line().

Referenced by distance_ring_linestring(), distance_ring_ring(), and ST_Distance_LineString_LineString().

                                       {
  if (line_intersects_line(l11x, l11y, l12x, l12y, l21x, l21y, l22x, l22y)) {
    return 0.0;
  }
  double dist12 = fmin(distance_point_line(l11x, l11y, l21x, l21y, l22x, l22y),
                       distance_point_line(l12x, l12y, l21x, l21y, l22x, l22y));
  double dist21 = fmin(distance_point_line(l21x, l21y, l11x, l11y, l12x, l12y),
                       distance_point_line(l22x, l22y, l11x, l11y, l12x, l12y));
  return fmin(dist12, dist21);
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE double distance_line_line_squared	(	double	l11x,
		double	l11y,
		double	l12x,
		double	l12y,
		double	l21x,
		double	l21y,
		double	l22x,
		double	l22y
	)

Definition at line 461 of file ExtensionFunctionsGeo.hpp.

References distance_point_line_squared(), and line_intersects_line().

                                               {
  if (line_intersects_line(l11x, l11y, l12x, l12y, l21x, l21y, l22x, l22y)) {
    return 0.0;
  }
  double dist12 = fmin(distance_point_line_squared(l11x, l11y, l21x, l21y, l22x, l22y),
                       distance_point_line_squared(l12x, l12y, l21x, l21y, l22x, l22y));
  double dist21 = fmin(distance_point_line_squared(l21x, l21y, l11x, l11y, l12x, l12y),
                       distance_point_line_squared(l22x, l22y, l11x, l11y, l12x, l12y));
  return fmin(dist12, dist21);
}

Here is the call graph for this function:

DEVICE double distance_point_line	(	double	px,
		double	py,
		double	l1x,
		double	l1y,
		double	l2x,
		double	l2y
	)

Definition at line 258 of file ExtensionFunctionsGeo.hpp.

References distance_point_point(), and tol_zero().

Referenced by distance_line_line(), distance_point_linestring(), and polygon_contains_point().

                                       {
  double length = distance_point_point(l1x, l1y, l2x, l2y);
  if (tol_zero(length)) {
    return distance_point_point(px, py, l1x, l1y);
  }

  // Find projection of point P onto the line segment AB:
  // Line containing that segment: A + k * (B - A)
  // Projection of point P onto the line touches it at
  //   k = dot(P-A,B-A) / length^2
  // AB segment is represented by k = [0,1]
  // Clamping k to [0,1] will give the shortest distance from P to AB segment
  double dotprod = (px - l1x) * (l2x - l1x) + (py - l1y) * (l2y - l1y);
  double k = dotprod / (length * length);
  k = fmax(0.0, fmin(1.0, k));
  double projx = l1x + k * (l2x - l1x);
  double projy = l1y + k * (l2y - l1y);
  return distance_point_point(px, py, projx, projy);
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE double distance_point_line_squared	(	double	px,
		double	py,
		double	l1x,
		double	l1y,
		double	l2x,
		double	l2y
	)

Definition at line 285 of file ExtensionFunctionsGeo.hpp.

References distance_point_point_squared(), tol_zero(), and TOLERANCE_DEFAULT_SQUARED.

Referenced by distance_line_line_squared().

                                               {
  double length = distance_point_point_squared(l1x, l1y, l2x, l2y);
  if (tol_zero(length, TOLERANCE_DEFAULT_SQUARED)) {
    return distance_point_point_squared(px, py, l1x, l1y);
  }

  // Find projection of point P onto the line segment AB:
  // Line containing that segment: A + k * (B - A)
  // Projection of point P onto the line touches it at
  //   k = dot(P-A,B-A) / length^2
  // AB segment is represented by k = [0,1]
  // Clamping k to [0,1] will give the shortest distance from P to AB segment
  double dotprod = (px - l1x) * (l2x - l1x) + (py - l1y) * (l2y - l1y);
  double k = dotprod / (length * length);
  k = fmax(0.0, fmin(1.0, k));
  double projx = l1x + k * (l2x - l1x);
  double projy = l1y + k * (l2y - l1y);
  return distance_point_point_squared(px, py, projx, projy);
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double distance_point_linestring	(	int8_t *	p,
		int64_t	psize,
		int8_t *	l,
		int64_t	lsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		bool	check_closed,
		double	threshold
	)

Definition at line 1947 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), distance_point_line(), and get_point().

Referenced by distance_point_multilinestring(), ST_Distance_MultiPoint_LineString(), ST_Distance_Point_ClosedLineString(), and ST_Distance_Point_LineString().

                                                                        {
  Point2D pt = get_point(p, 0, ic1, isr1, osr);

  auto l_num_coords = lsize / compression_unit_size(ic2);

  Point2D l1 = get_point(l, 0, ic2, isr2, osr);
  Point2D l2 = get_point(l, 2, ic2, isr2, osr);

  double dist = distance_point_line(pt.x, pt.y, l1.x, l1.y, l2.x, l2.y);
  for (int32_t i = 4; i < l_num_coords; i += 2) {
    l1 = l2;  // advance one point
    l2 = get_point(l, i, ic2, isr2, osr);
    double ldist = distance_point_line(pt.x, pt.y, l1.x, l1.y, l2.x, l2.y);
    if (dist > ldist) {
      dist = ldist;
    }
    if (dist <= threshold) {
      return dist;
    }
  }
  if (l_num_coords > 4 && check_closed) {
    // Also check distance to the closing edge between the first and the last points
    l1 = get_point(l, 0, ic2, isr2, osr);
    double ldist = distance_point_line(pt.x, pt.y, l1.x, l1.y, l2.x, l2.y);
    if (dist > ldist) {
      dist = ldist;
    }
  }
  return dist;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double distance_point_multilinestring	(	int8_t *	p,
		int64_t	psize,
		int8_t *	mls,
		int64_t	mls_size,
		int32_t *	mls_ls_sizes,
		int64_t	mls_ls_num,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2018 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), distance_point_linestring(), and tol_zero().

Referenced by ST_Distance_MultiPoint_MultiLineString(), and ST_Distance_Point_MultiLineString().

                                                                             {
  if (mls_ls_num <= 0) {
    return 0.0;
  }
  double min_distance = 0.0;

  auto next_linestring_coords = mls;

  for (auto l = 0; l < mls_ls_num; l++) {
    auto linestring_coords = next_linestring_coords;
    auto linestring_num_points = mls_ls_sizes[l];
    auto linestring_num_coords = 2 * linestring_num_points;
    auto linestring_coords_size = linestring_num_coords * compression_unit_size(ic2);
    next_linestring_coords += linestring_coords_size;
    double distance = distance_point_linestring(p,
                                                psize,
                                                linestring_coords,
                                                linestring_coords_size,
                                                ic1,
                                                isr1,
                                                ic2,
                                                isr2,
                                                osr,
                                                false,
                                                threshold);
    if (l == 0 || min_distance > distance) {
      min_distance = distance;
      if (tol_zero(min_distance)) {
        min_distance = 0.0;
        break;
      }
      if (min_distance <= threshold) {
        break;
      }
    }
  }

  return min_distance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double distance_point_multipolygon	(	int8_t *	p,
		int64_t	psize,
		int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2148 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), distance_point_polygon(), and tol_zero().

Referenced by ST_Distance_MultiPoint_MultiPolygon(), and ST_Distance_Point_MultiPolygon().

                                                                          {
  if (mpoly_num_polys <= 0) {
    return 0.0;
  }
  double min_distance = 0.0;

  // Set specific poly pointers as we move through the coords/ringsizes/polyrings arrays.
  auto next_poly_coords = mpoly_coords;
  auto next_poly_ring_sizes = mpoly_ring_sizes;

  for (auto poly = 0; poly < mpoly_num_polys; poly++) {
    auto poly_coords = next_poly_coords;
    auto poly_ring_sizes = next_poly_ring_sizes;
    auto poly_num_rings = mpoly_poly_sizes[poly];
    // Count number of coords in all of poly's rings, advance ring size pointer.
    int32_t poly_num_coords = 0;
    for (auto ring = 0; ring < poly_num_rings; ring++) {
      poly_num_coords += 2 * *next_poly_ring_sizes++;
    }
    auto poly_coords_size = poly_num_coords * compression_unit_size(ic2);
    next_poly_coords += poly_coords_size;
    double distance = distance_point_polygon(p,
                                             psize,
                                             poly_coords,
                                             poly_coords_size,
                                             poly_ring_sizes,
                                             poly_num_rings,
                                             ic1,
                                             isr1,
                                             ic2,
                                             isr2,
                                             osr,
                                             threshold);
    if (poly == 0 || min_distance > distance) {
      min_distance = distance;
      if (tol_zero(min_distance)) {
        min_distance = 0.0;
        break;
      }
      if (min_distance <= threshold) {
        break;
      }
    }
  }

  return min_distance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double distance_point_point	(	double	p1x,
		double	p1y,
		double	p2x,
		double	p2y
	)

Definition at line 248 of file ExtensionFunctionsGeo.hpp.

References distance_point_point_squared().

Referenced by centroid_add_segment(), distance_point_line(), length_linestring(), max_distance_point_line(), ST_Distance_MultiPoint_MultiPoint(), ST_Distance_Point_MultiPoint(), and ST_Distance_Point_Point().

                                                             {
  auto d2 = distance_point_point_squared(p1x, p1y, p2x, p2y);
  return sqrt(d2);
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double distance_point_point_squared	(	double	p1x,
		double	p1y,
		double	p2x,
		double	p2y
	)

Definition at line 232 of file ExtensionFunctionsGeo.hpp.

References tol_zero(), and TOLERANCE_DEFAULT_SQUARED.

Referenced by distance_point_line_squared(), distance_point_point(), ST_Distance_MultiPoint_MultiPoint_Squared(), ST_Distance_Point_MultiPoint_Squared(), and ST_Distance_Point_Point_Squared().

                                                                     {
  auto x = p1x - p2x;
  auto y = p1y - p2y;
  auto x2 = x * x;
  auto y2 = y * y;
  auto d2 = x2 + y2;
  if (tol_zero(d2, TOLERANCE_DEFAULT_SQUARED)) {
    return 0.0;
  }
  return d2;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double distance_point_polygon	(	int8_t *	p,
		int64_t	psize,
		int8_t *	poly,
		int64_t	polysize,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2096 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), get_point(), polygon_contains_point(), and ST_Distance_Point_ClosedLineString().

Referenced by distance_point_multipolygon(), ST_Distance_MultiPoint_Polygon(), and ST_Distance_Point_Polygon().

                                                                     {
  auto exterior_ring_num_coords = polysize / compression_unit_size(ic2);
  if (poly_num_rings > 0) {
    exterior_ring_num_coords = poly_ring_sizes[0] * 2;
  }
  auto exterior_ring_coords_size = exterior_ring_num_coords * compression_unit_size(ic2);

  Point2D pt = get_point(p, 0, ic1, isr1, osr);
  if (!polygon_contains_point(
          poly, exterior_ring_num_coords, pt.x, pt.y, ic2, isr2, osr)) {
    // Outside the exterior ring
    return ST_Distance_Point_ClosedLineString(
        p, psize, poly, exterior_ring_coords_size, ic1, isr1, ic2, isr2, osr, threshold);
  }
  // Inside exterior ring
  // Advance to first interior ring
  poly += exterior_ring_coords_size;
  // Check if one of the polygon's holes contains that point
  for (auto r = 1; r < poly_num_rings; r++) {
    auto interior_ring_num_coords = poly_ring_sizes[r] * 2;
    auto interior_ring_coords_size =
        interior_ring_num_coords * compression_unit_size(ic2);
    if (polygon_contains_point(
            poly, interior_ring_num_coords, pt.x, pt.y, ic2, isr2, osr)) {
      // Inside an interior ring
      return ST_Distance_Point_ClosedLineString(p,
                                                psize,
                                                poly,
                                                interior_ring_coords_size,
                                                ic1,
                                                isr1,
                                                ic2,
                                                isr2,
                                                osr,
                                                threshold);
    }
    poly += interior_ring_coords_size;
  }
  return 0.0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE double distance_ring_linestring	(	int8_t *	ring,
		int32_t	ring_num_coords,
		int8_t *	l,
		int32_t	lnum_coords,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 480 of file ExtensionFunctionsGeo.hpp.

References distance_line_line(), get_point(), and tol_zero().

Referenced by ST_Distance_LineString_Polygon().

                                                  {
  double min_distance = 0.0;

  Point2D re1 = get_point(ring, ring_num_coords - 2, ic1, isr1, osr);
  for (auto i = 0; i < ring_num_coords; i += 2) {
    Point2D re2 = get_point(ring, i, ic1, isr1, osr);
    Point2D le1 = get_point(l, 0, ic2, isr2, osr);
    for (auto j = 2; j < lnum_coords; j += 2) {
      Point2D le2 = get_point(l, j, ic2, isr2, osr);
      auto distance =
          distance_line_line(re1.x, re1.y, re2.x, re2.y, le1.x, le1.y, le2.x, le2.y);
      if ((i == 0 && j == 2) || min_distance > distance) {
        min_distance = distance;
        if (tol_zero(min_distance)) {
          return 0.0;
        }
        if (min_distance <= threshold) {
          return min_distance;
        }
      }
      le1 = le2;
    }
    re1 = re2;
  }
  return min_distance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE double distance_ring_ring	(	int8_t *	ring1,
		int32_t	ring1_num_coords,
		int8_t *	ring2,
		int32_t	ring2_num_coords,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 517 of file ExtensionFunctionsGeo.hpp.

References distance_line_line(), get_point(), and tol_zero().

Referenced by ST_Distance_Polygon_Polygon().

                                            {
  double min_distance = 0.0;
  Point2D e11 = get_point(ring1, ring1_num_coords - 2, ic1, isr1, osr);
  for (auto i = 0; i < ring1_num_coords; i += 2) {
    Point2D e12 = get_point(ring1, i, ic1, isr1, osr);
    Point2D e21 = get_point(ring2, ring2_num_coords - 2, ic2, isr2, osr);
    for (auto j = 0; j < ring2_num_coords; j += 2) {
      Point2D e22 = get_point(ring2, j, ic2, isr2, osr);
      auto distance =
          distance_line_line(e11.x, e11.y, e12.x, e12.y, e21.x, e21.y, e22.x, e22.y);
      if ((i == 0 && j == 0) || min_distance > distance) {
        min_distance = distance;
        if (tol_zero(min_distance)) {
          return 0.0;
        }
        if (min_distance <= threshold) {
          return min_distance;
        }
      }
      e21 = e22;
    }
    e11 = e12;
  }
  return min_distance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE Point2D get_point	(	const int8_t *	data,
		const int32_t	x_index,
		const int32_t	ic,
		const int32_t	isr,
		const int32_t	osr
	)

Definition at line 192 of file ExtensionFunctionsGeo.hpp.

Referenced by area_ring(), centroid_add_linestring(), centroid_add_ring(), Contains_MultiPolygon_Point_Impl(), Contains_Polygon_Point_Impl(), distance_point_linestring(), distance_point_polygon(), distance_ring_linestring(), distance_ring_ring(), length_linestring(), linestring_intersects_line(), linestring_intersects_linestring(), max_distance_point_linestring(), point_dwithin_box(), point_in_polygon_winding_number(), polygon_contains_linestring(), polygon_contains_point(), ring_intersects_line(), ST_Centroid_MultiPoint(), ST_Centroid_MultiPolygon(), ST_Centroid_Point(), ST_Centroid_Polygon(), ST_Contains_Point_LineString(), ST_Contains_Point_Point(), ST_Distance_LineString_LineString(), ST_Distance_MultiPoint_MultiPoint(), ST_Distance_MultiPoint_MultiPoint_Squared(), ST_Distance_Point_LineString_Geodesic(), ST_Distance_Point_MultiPoint(), ST_Distance_Point_MultiPoint_Squared(), ST_Distance_Point_Point(), ST_Distance_Point_Point_Geodesic(), ST_Distance_Point_Point_Squared(), ST_Intersects_Point_LineString(), and trim_linestring_to_buffered_box().

                                                          {
  Point2D const decompressed{decompress_coord<X>(data, x_index, ic),
                             decompress_coord<Y>(data, x_index, ic)};
  return transform_point<XY>(decompressed, isr, osr);
}

Here is the caller graph for this function:

EXTENSION_INLINE int64_t HeavyDB_Geo_PolyBoundsPtr	(	double *	bounds,
		int64_t	size
	)

Definition at line 4995 of file ExtensionFunctionsGeo.hpp.

Referenced by com.mapd.calcite.parser.HeavyDBSqlOperatorTable::addUDF().

                                                                {
  return reinterpret_cast<int64_t>(bounds);
}

Here is the caller graph for this function:

template<typename T >

DEVICE ALWAYS_INLINE T is_left	(	const T	lx0,
		const T	ly0,
		const T	lx1,
		const T	ly1,
		const T	px,
		const T	py
	)

Tests if a point is left, on, or to the right of a 2D line

Returns

>0 if p is left of l, 0 if p is on l, and < 0 if p is right of l

Definition at line 558 of file ExtensionFunctionsGeo.hpp.

Referenced by point_in_polygon_winding_number().

                                                                                    {
  return ((lx1 - lx0) * (py - ly0) - (px - lx0) * (ly1 - ly0));
}

Here is the caller graph for this function:

EXTENSION_NOINLINE bool is_point_in_view	(	int8_t *	p,
		const int64_t	psize,
		const int32_t	ic,
		const double	min_lon,
		const double	max_lon,
		const double	min_lat,
		const double	max_lat
	)

Definition at line 5062 of file ExtensionFunctionsGeo.hpp.

Referenced by com.mapd.calcite.parser.HeavyDBSqlOperatorTable::addUDF().

                                                               {
  const auto lon = decompress_coord<X>(p, 0, ic);
  const auto lat = decompress_coord<Y>(p, 0, ic);
  return !(lon < min_lon || lon > max_lon || lat < min_lat || lat > max_lat);
}

Here is the caller graph for this function:

EXTENSION_NOINLINE bool is_point_size_in_view	(	int8_t *	p,
		const int64_t	psize,
		const int32_t	ic,
		const double	meters,
		const double	min_lon,
		const double	max_lon,
		const double	min_lat,
		const double	max_lat
	)

Definition at line 5074 of file ExtensionFunctionsGeo.hpp.

Referenced by com.mapd.calcite.parser.HeavyDBSqlOperatorTable::addUDF().

                                                                    {
  const double const1 = 0.017453292519943295769236907684886;
  const double const2 = 6372797.560856;
  const auto lon = decompress_coord<X>(p, 0, ic);
  const auto lat = decompress_coord<Y>(p, 0, ic);
  const double latdiff = meters / (const1 * const2);
  const double t1 = sinf(meters / (2.0 * const2));
  const double t2 = cosf(const1 * lat);
  const double londiff = (2.0 * asinf(t1 / t2)) / const1;
  return !(lon + londiff < min_lon || lon - londiff > max_lon ||
           lat + latdiff < min_lat || lat - latdiff > max_lat);
}

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double length_linestring	(	int8_t *	l,
		int32_t	lsize,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr,
		bool	geodesic,
		bool	check_closed
	)

Definition at line 1140 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), distance_in_meters(), distance_point_point(), and get_point().

Referenced by perimeter_multipolygon(), ST_Length_LineString(), ST_Length_LineString_Geodesic(), ST_Length_MultiLineString(), ST_Perimeter_Polygon(), and ST_Perimeter_Polygon_Geodesic().

                                                                 {
  auto l_num_coords = lsize / compression_unit_size(ic);

  double length = 0.0;

  Point2D l0 = get_point(l, 0, ic, isr, osr);
  Point2D l2 = l0;
  for (int32_t i = 2; i < l_num_coords; i += 2) {
    Point2D l1 = l2;
    l2 = get_point(l, i, ic, isr, osr);
    double ldist = geodesic ? distance_in_meters(l1.x, l1.y, l2.x, l2.y)
                            : distance_point_point(l1.x, l1.y, l2.x, l2.y);
    length += ldist;
  }
  if (check_closed) {
    double ldist = geodesic ? distance_in_meters(l2.x, l2.y, l0.x, l0.y)
                            : distance_point_point(l2.x, l2.y, l0.x, l0.y);
    length += ldist;
  }
  return length;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE bool line_intersects_line	(	double	l11x,
		double	l11y,
		double	l12x,
		double	l12y,
		double	l21x,
		double	l21y,
		double	l22x,
		double	l22y
	)

Definition at line 336 of file ExtensionFunctionsGeo.hpp.

References on_segment(), and orientation().

Referenced by distance_line_line(), distance_line_line_squared(), linestring_intersects_line(), linestring_intersects_linestring(), polygon_contains_point(), and ring_intersects_line().

                                       {
  auto o1 = orientation(l11x, l11y, l12x, l12y, l21x, l21y);
  auto o2 = orientation(l11x, l11y, l12x, l12y, l22x, l22y);
  auto o3 = orientation(l21x, l21y, l22x, l22y, l11x, l11y);
  auto o4 = orientation(l21x, l21y, l22x, l22y, l12x, l12y);

  // General case
  if (o1 != o2 && o3 != o4) {
    return true;
  }

  // Special Cases
  // l11, l12 and l21 are colinear and l21 lies on segment l11-l12
  if (o1 == 0 && on_segment(l11x, l11y, l21x, l21y, l12x, l12y)) {
    return true;
  }

  // l11, l12 and l21 are colinear and l22 lies on segment l11-l12
  if (o2 == 0 && on_segment(l11x, l11y, l22x, l22y, l12x, l12y)) {
    return true;
  }

  // l21, l22 and l11 are colinear and l11 lies on segment l21-l22
  if (o3 == 0 && on_segment(l21x, l21y, l11x, l11y, l22x, l22y)) {
    return true;
  }

  // l21, l22 and l12 are colinear and l12 lies on segment l21-l22
  if (o4 == 0 && on_segment(l21x, l21y, l12x, l12y, l22x, l22y)) {
    return true;
  }

  return false;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE bool linestring_intersects_line	(	int8_t *	l,
		int32_t	lnum_coords,
		double	l1x,
		double	l1y,
		double	l2x,
		double	l2y,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	osr
	)

Definition at line 379 of file ExtensionFunctionsGeo.hpp.

References get_point(), and line_intersects_line().

Referenced by ring_intersects_line().

                                             {
  Point2D e1 = get_point(l, 0, ic1, isr1, osr);
  for (int64_t i = 2; i < lnum_coords; i += 2) {
    Point2D e2 = get_point(l, i, ic1, isr1, osr);
    if (line_intersects_line(e1.x, e1.y, e2.x, e2.y, l1x, l1y, l2x, l2y)) {
      return true;
    }
    e1 = e2;
  }
  return false;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE bool linestring_intersects_linestring	(	int8_t *	l,
		int32_t	lnum_coords,
		double	l1x,
		double	l1y,
		double	l2x,
		double	l2y,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	osr
	)

Definition at line 419 of file ExtensionFunctionsGeo.hpp.

References get_point(), and line_intersects_line().

                                                   {
  Point2D e1 = get_point(l, 0, ic1, isr1, osr);
  for (int64_t i = 2; i < lnum_coords; i += 2) {
    Point2D e2 = get_point(l, i, ic1, isr1, osr);
    if (line_intersects_line(e1.x, e1.y, e2.x, e2.y, l1x, l1y, l2x, l2y)) {
      return true;
    }
    e1 = e2;
  }
  return false;
}

Here is the call graph for this function:

DEVICE double max_distance_point_line	(	double	px,
		double	py,
		double	l1x,
		double	l1y,
		double	l2x,
		double	l2y
	)

Definition at line 3769 of file ExtensionFunctionsGeo.hpp.

References distance_point_point().

Referenced by max_distance_point_linestring().

                                           {
  // TODO: switch to squared distances
  double length1 = distance_point_point(px, py, l1x, l1y);
  double length2 = distance_point_point(px, py, l2x, l2y);
  if (length1 > length2) {
    return length1;
  }
  return length2;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double max_distance_point_linestring	(	int8_t *	p,
		int64_t	psize,
		int8_t *	l,
		int64_t	lsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		bool	check_closed
	)

Definition at line 3784 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), get_point(), and max_distance_point_line().

Referenced by ST_MaxDistance_LineString_Point(), and ST_MaxDistance_Point_LineString().

                                                                             {
  // TODO: switch to squared distances
  Point2D pt = get_point(p, 0, ic1, isr1, osr);

  auto l_num_coords = lsize / compression_unit_size(ic2);

  Point2D l1 = get_point(l, 0, ic2, isr2, osr);
  Point2D l2 = get_point(l, 2, ic2, isr2, osr);

  double max_dist = max_distance_point_line(pt.x, pt.y, l1.x, l1.y, l2.x, l2.y);
  for (int32_t i = 4; i < l_num_coords; i += 2) {
    l1 = l2;  // advance one point
    l2 = get_point(l, i, ic2, isr2, osr);
    double ldist = max_distance_point_line(pt.x, pt.y, l1.x, l1.y, l2.x, l2.y);
    if (max_dist < ldist) {
      max_dist = ldist;
    }
  }
  if (l_num_coords > 4 && check_closed) {
    // Also check distance to the closing edge between the first and the last points
    l1 = get_point(l, 0, ic2, isr2, osr);
    double ldist = max_distance_point_line(pt.x, pt.y, l1.x, l1.y, l2.x, l2.y);
    if (max_dist < ldist) {
      max_dist = ldist;
    }
  }
  return max_dist;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool on_segment	(	double	px,
		double	py,
		double	qx,
		double	qy,
		double	rx,
		double	ry
	)

Definition at line 312 of file ExtensionFunctionsGeo.hpp.

References tol_ge(), and tol_le().

Referenced by line_intersects_line().

                                                {
  return (tol_le(qx, fmax(px, rx)) && tol_ge(qx, fmin(px, rx)) &&
          tol_le(qy, fmax(py, ry)) && tol_ge(qy, fmin(py, ry)));
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE int16_t orientation	(	double	px,
		double	py,
		double	qx,
		double	qy,
		double	rx,
		double	ry
	)

Definition at line 323 of file ExtensionFunctionsGeo.hpp.

References tol_zero().

Referenced by line_intersects_line().

                                                                              {
  auto val = ((qy - py) * (rx - qx) - (qx - px) * (ry - qy));
  if (tol_zero(val)) {
    return 0;  // Points p, q and r are colinear
  }
  if (val > 0.0) {
    return 1;  // Clockwise point orientation
  }
  return 2;  // Counterclockwise point orientation
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double perimeter_multipolygon	(	int8_t *	mpoly_coords,
		int32_t	mpoly_coords_size,
		int8_t *	mpoly_ring_sizes_in,
		int32_t	mpoly_num_rings,
		int8_t *	mpoly_poly_sizes,
		int32_t	mpoly_num_polys,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr,
		bool	geodesic
	)

Definition at line 1257 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), and length_linestring().

Referenced by ST_Perimeter_MultiPolygon(), and ST_Perimeter_MultiPolygon_Geodesic().

                                                                  {
  if (mpoly_num_polys <= 0 || mpoly_num_rings <= 0) {
    return 0.0;
  }

  double perimeter = 0.0;

  auto mpoly_ring_sizes = reinterpret_cast<int32_t*>(mpoly_ring_sizes_in);

  // Set specific poly pointers as we move through the coords/ringsizes/polyrings arrays.
  auto next_poly_coords = mpoly_coords;
  auto next_poly_ring_sizes = mpoly_ring_sizes;

  for (auto poly = 0; poly < mpoly_num_polys; poly++) {
    auto poly_coords = next_poly_coords;
    auto poly_ring_sizes = next_poly_ring_sizes;
    auto poly_num_rings = reinterpret_cast<int32_t*>(mpoly_poly_sizes)[poly];
    // Count number of coords in all of poly's rings, advance ring size pointer.
    int32_t poly_num_coords = 0;
    for (auto ring = 0; ring < poly_num_rings; ring++) {
      poly_num_coords += 2 * *next_poly_ring_sizes++;
    }
    auto poly_coords_size = poly_num_coords * compression_unit_size(ic);
    next_poly_coords += poly_coords_size;

    auto exterior_ring_num_coords = poly_ring_sizes[0] * 2;
    auto exterior_ring_coords_size = exterior_ring_num_coords * compression_unit_size(ic);

    perimeter += length_linestring(
        poly_coords, exterior_ring_coords_size, ic, isr, osr, geodesic, true);
  }

  return perimeter;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool point_dwithin_box	(	int8_t *	p1,
		int64_t	p1size,
		int32_t	ic1,
		int32_t	isr1,
		double *	bounds2,
		int64_t	bounds2_size,
		int32_t	isr2,
		int32_t	osr,
		double	distance
	)

Definition at line 921 of file ExtensionFunctionsGeo.hpp.

References get_point().

Referenced by ST_DWithin_Point_LineString(), ST_DWithin_Point_MultiPolygon(), and ST_DWithin_Point_Polygon().

                                                             {
  // TODO: tolerance

  // Point has to be uncompressed and transformed to output SR.
  // Bounding box has to be transformed to output SR.
  Point2D p = get_point(p1, 0, ic1, isr1, osr);
  Point2D const lb = transform_point<XY>({bounds2[0], bounds2[1]}, isr2, osr);
  Point2D const rt = transform_point<XY>({bounds2[2], bounds2[3]}, isr2, osr);
  if (p.x < lb.x - distance ||  // point is left of box2:  px < box2.xmin - distance
      p.x > rt.x + distance ||  // point is right of box2: px > box2.xmax + distance
      p.y < lb.y - distance ||  // point is below box2:    py < box2.ymin - distance
      p.y > rt.y + distance) {  // point is above box2:    py > box2.ymax + distance
    return false;
  }
  return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

template<typename T , EdgeBehavior TEdgeBehavior>

DEVICE ALWAYS_INLINE bool point_in_polygon_winding_number	(	const int8_t *	poly,
		const int32_t	poly_num_coords,
		const T	px,
		const T	py,
		const int32_t	ic1,
		const int32_t	isr1,
		const int32_t	osr
	)

Computes whether the point p is inside the polygon poly using the winding number algorithm. The winding number algorithm efficiently computes the winding number, which is the number of revolutions made around a point P while traveling around the polygon. If the winding number is non-zero, then P must be inside the polygon. The actual algorithm loops over all edges in the polygon, and determines the location of the crossing between a ray from the point P and each edge E of the polygon. In the implementation below, this crossing calculation is performed by assigning a fixed orientation to each edge, then computing the location of the point with respect to that edge. The parameter include_point_on_edge allows the algorithm to work for both point-inside-polygon, or point-on-polygon calculations. The original source and additional detail for this implementation is here: http://geomalgorithms.com/a03-_inclusion.html

Definition at line 587 of file ExtensionFunctionsGeo.hpp.

References compressed_coord(), DEBUG_STMT, get_point(), is_left(), kIncludePointOnEdge, heavydb.dtypes::T, tol_zero_template(), and UNLIKELY.

                                                                             {
  int32_t wn = 0;

  constexpr bool include_point_on_edge =
      TEdgeBehavior == EdgeBehavior::kIncludePointOnEdge;

  auto get_x_coord = [=](const auto data, const auto index) -> T {
    if constexpr (std::is_floating_point<T>::value) {
      return get_point(data, index, ic1, isr1, osr).x;
    } else {
      return compressed_coord(data, index);
    }
    return T{};  // https://stackoverflow.com/a/64561686/2700898
  };

  auto get_y_coord = [=](const auto data, const auto index) -> T {
    if constexpr (std::is_floating_point<T>::value) {
      return get_point(data, index, ic1, isr1, osr).y;
    } else {
      return compressed_coord(data, index + 1);
    }
    return T{};  // https://stackoverflow.com/a/64561686/2700898
  };

  DEBUG_STMT(printf("Poly num coords: %d\n", poly_num_coords));
  const int64_t num_edges = poly_num_coords / 2;

  int64_t e0_index = 0;
  T e0x = get_x_coord(poly, e0_index * 2);
  T e0y = get_y_coord(poly, e0_index * 2);

  int64_t e1_index;
  T e1x, e1y;

  for (int64_t edge = 1; edge <= num_edges; edge++) {
    DEBUG_STMT(printf("Edge %ld\n", edge));
    // build the edge
    e1_index = edge % num_edges;
    e1x = get_x_coord(poly, e1_index * 2);
    e1y = get_y_coord(poly, e1_index * 2);

    DEBUG_STMT(printf("edge 0: %ld : %f, %f\n", e0_index, (double)e0x, (double)e0y));
    DEBUG_STMT(printf("edge 1: %ld : %f, %f\n", e1_index, (double)e1x, (double)e1y));

    auto get_epsilon = [=]() -> T {
      if constexpr (std::is_floating_point<T>::value) {
        const T edge_vec_magnitude =
            (e1x - e0x) * (e1x - e0x) + (e1y - e0y) * (e1y - e0y);
        return 0.003 * edge_vec_magnitude;
      } else {
        return T(0);
      }
      return T{};  // https://stackoverflow.com/a/64561686/2700898
    };

    const T epsilon = get_epsilon();
    DEBUG_STMT(printf("using epsilon value %e\n", (double)epsilon));

    if constexpr (include_point_on_edge) {
      const T xp = (px - e1x) * (e0y - e1y) - (py - e1y) * (e0x - e1x);
      const T pt_vec_magnitude = (px - e0x) * (px - e0x) + (py - e0y) * (py - e0y);
      const T edge_vec_magnitude = (e1x - e0x) * (e1x - e0x) + (e1y - e0y) * (e1y - e0y);
      if (tol_zero_template(xp, epsilon) && (pt_vec_magnitude <= edge_vec_magnitude)) {
        DEBUG_STMT(printf("point is on edge: %e && %e <= %e\n",
                          (double)xp,
                          (double)pt_vec_magnitude,
                          (double)edge_vec_magnitude));
        return include_point_on_edge;
      }
    }

    if (e0y <= py) {
      if (e1y > py) {
        // upward crossing
        DEBUG_STMT(printf("upward crossing\n"));
        const auto is_left_val = is_left(e0x, e0y, e1x, e1y, px, py);
        DEBUG_STMT(printf("Left val %f and tol zero left val %d\n",
                          (double)is_left_val,
                          tol_zero_template(is_left_val, epsilon)));
        if (UNLIKELY(tol_zero_template(is_left_val, epsilon))) {
          // p is on the edge
          return include_point_on_edge;
        } else if (is_left_val > T(0)) {  // p left of edge
          // valid up intersection
          DEBUG_STMT(printf("++wn\n"));
          ++wn;
        }
      }
    } else if (e1y <= py) {
      // downward crossing
      DEBUG_STMT(printf("downward crossing\n"));
      const auto is_left_val = is_left(e0x, e0y, e1x, e1y, px, py);
      DEBUG_STMT(printf("Left val %f and tol zero left val %d\n",
                        (double)is_left_val,
                        tol_zero_template(is_left_val, epsilon)));
      if (UNLIKELY(tol_zero_template(is_left_val, epsilon))) {
        // p is on the edge
        return include_point_on_edge;
      } else if (is_left_val < T(0)) {  // p right of edge
        // valid down intersection
        DEBUG_STMT(printf("--wn\n"));
        --wn;
      }
    }

    e0_index = e1_index;
    e0x = e1x;
    e0y = e1y;
  }
  DEBUG_STMT(printf("wn: %d\n", wn));
  // wn == 0 --> P is outside the polygon
  return !(wn == 0);
}

Here is the call graph for this function:

EXTENSION_NOINLINE bool Point_Overlaps_Box	(	double *	bounds,
		int64_t	bounds_size,
		double	px,
		double	py
	)

Definition at line 870 of file ExtensionFunctionsGeo.hpp.

References box_contains_point().

                                                      {
  // TODO: mixed spatial references
  return box_contains_point(bounds, bounds_size, px, py);
}

Here is the call graph for this function:

DEVICE bool polygon_contains_linestring	(	int8_t *	poly,
		int32_t	poly_num_coords,
		int8_t *	l,
		int64_t	lnum_coords,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 833 of file ExtensionFunctionsGeo.hpp.

References get_point(), polygon_contains_point(), and ring_intersects_line().

Referenced by ST_Contains_Polygon_LineString().

                                              {
  // Check that the first point is in the polygon
  Point2D l1 = get_point(l, 0, ic2, isr2, osr);
  if (!polygon_contains_point(poly, poly_num_coords, l1.x, l1.y, ic1, isr1, osr)) {
    return false;
  }

  // Go through line segments and check if there are no intersections with poly edges,
  // i.e. linestring doesn't escape
  for (int32_t i = 2; i < lnum_coords; i += 2) {
    Point2D l2 = get_point(l, i, ic2, isr2, osr);
    if (ring_intersects_line(
            poly, poly_num_coords, l1.x, l1.y, l2.x, l2.y, ic1, isr1, osr)) {
      return false;
    }
    l1 = l2;
  }
  return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool polygon_contains_point	(	const int8_t *	poly,
		const int32_t	poly_num_coords,
		const double	px,
		const double	py,
		const int32_t	ic1,
		const int32_t	isr1,
		const int32_t	osr
	)

Definition at line 721 of file ExtensionFunctionsGeo.hpp.

References distance_point_line(), get_point(), line_intersects_line(), run_benchmark_import::result, tol_eq(), and tol_zero().

Referenced by distance_point_polygon(), and polygon_contains_linestring().

                                                                    {
  bool result = false;
  int xray_touch = 0;
  bool horizontal_edge = false;
  bool yray_intersects = false;

  Point2D e1 = get_point(poly, poly_num_coords - 2, ic1, isr1, osr);
  for (int64_t i = 0; i < poly_num_coords; i += 2) {
    Point2D e2 = get_point(poly, i, ic1, isr1, osr);

    // Check if point sits on an edge.
    if (tol_zero(distance_point_line(px, py, e1.x, e1.y, e2.x, e2.y))) {
      return true;
    }

    // Before flipping the switch, check if xray hit a horizontal edge
    // - If an edge lays on the xray, one of the previous edges touched it
    //   so while moving horizontally we're in 'xray_touch' state
    // - Last edge that touched xray at (e2.x,e2.y) didn't register intersection
    // - Next edge that diverges from xray at (e1,e1y) will register intersection
    // - Can have several horizontal edges, one after the other, keep moving though
    //   in 'xray_touch' state without flipping the switch
    horizontal_edge = (xray_touch != 0) && tol_eq(py, e1.y) && tol_eq(py, e2.y);

    // Main probe: xray
    // Overshoot the xray to detect an intersection if there is one.
    double xray = fmax(e2.x, e1.x) + 1.0;
    if (px <= xray &&        // Only check for intersection if the edge is on the right
        !horizontal_edge &&  // Keep moving through horizontal edges
        line_intersects_line(px,  // xray shooting from point p to the right
                             py,
                             xray,
                             py,
                             e1.x,  // polygon edge
                             e1.y,
                             e2.x,
                             e2.y)) {
      // Register intersection
      result = !result;

      // Adjust for special cases
      if (xray_touch == 0) {
        if (tol_zero(distance_point_line(e2.x, e2.y, px, py, xray + 1.0, py))) {
          // Xray goes through the edge's second vertex, unregister intersection -
          // that vertex will be crossed again when we look at the following edge(s)
          result = !result;
          // Enter the xray-touch state:
          // (1) - xray was touched by the edge from above, (-1) from below
          xray_touch = (e1.y > py) ? 1 : -1;
        }
      } else {
        // Previous edge touched the xray, intersection hasn't been registered,
        // it has to be registered now if this edge continues across the xray.
        if (xray_touch > 0) {
          // Previous edge touched the xray from above
          if (e2.y <= py) {
            // Current edge crosses under xray: intersection is already registered
          } else {
            // Current edge just touched the xray and pulled up: unregister intersection
            result = !result;
          }
        } else {
          // Previous edge touched the xray from below
          if (e2.y > py) {
            // Current edge crosses over xray: intersection is already registered
          } else {
            // Current edge just touched the xray and pulled down: unregister intersection
            result = !result;
          }
        }
        // Exit the xray-touch state
        xray_touch = 0;
      }
    }

    // Redundancy: vertical yray down
    // Main probe xray may hit multiple complex fragments which increases a chance of
    // error. Perform a simple secondary check for edge intersections to see if point is
    // outside.
    if (!yray_intersects) {  // Continue checking on yray until intersection is found
      double yray = fmin(e2.y, e1.y) - 1.0;
      if (yray <= py) {  // Only check for yray intersection if point P is above the edge
        yray_intersects = line_intersects_line(px,  // yray shooting from point P down
                                               py,
                                               px,
                                               yray,
                                               e1.x,  // polygon edge
                                               e1.y,
                                               e2.x,
                                               e2.y);
      }
    }

    // Advance to the next vertex
    e1 = e2;
  }
  if (!yray_intersects) {
    // yray has zero intersections - point is outside the polygon
    return false;
  }
  // Otherwise rely on the main probe
  return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE bool ring_intersects_line	(	int8_t *	ring,
		int32_t	ring_num_coords,
		double	l1x,
		double	l1y,
		double	l2x,
		double	l2y,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	osr
	)

Definition at line 400 of file ExtensionFunctionsGeo.hpp.

References get_point(), line_intersects_line(), and linestring_intersects_line().

Referenced by polygon_contains_linestring().

                                       {
  Point2D e1 = get_point(ring, ring_num_coords - 2, ic1, isr1, osr);
  Point2D e2 = get_point(ring, 0, ic1, isr1, osr);
  if (line_intersects_line(e1.x, e1.y, e2.x, e2.y, l1x, l1y, l2x, l2y)) {
    return true;
  }
  return linestring_intersects_line(
      ring, ring_num_coords, l1x, l1y, l2x, l2y, ic1, isr1, osr);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Area_MultiPolygon	(	int8_t *	mpoly_coords,
		int32_t	mpoly_coords_size,
		int8_t *	mpoly_ring_sizes,
		int32_t	mpoly_num_rings,
		int8_t *	mpoly_poly_sizes_in,
		int32_t	mpoly_num_polys,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1421 of file ExtensionFunctionsGeo.hpp.

References area_polygon(), and compression_unit_size().

                                         {
  if (mpoly_num_rings <= 0 || mpoly_num_polys <= 0) {
    return 0.0;
  }

  double area = 0.0;

  auto mpoly_poly_sizes = reinterpret_cast<int32_t*>(mpoly_poly_sizes_in);

  // Set specific poly pointers as we move through the coords/ringsizes/polyrings arrays.
  auto next_poly_coords = mpoly_coords;
  auto next_poly_ring_sizes = reinterpret_cast<int32_t*>(mpoly_ring_sizes);

  for (auto poly = 0; poly < mpoly_num_polys; poly++) {
    auto poly_coords = next_poly_coords;
    auto poly_ring_sizes = next_poly_ring_sizes;
    auto poly_num_rings = mpoly_poly_sizes[poly];
    // Count number of coords in all of poly's rings, advance ring size pointer.
    int32_t poly_num_coords = 0;
    for (auto ring = 0; ring < poly_num_rings; ring++) {
      poly_num_coords += 2 * *next_poly_ring_sizes++;
    }
    auto poly_coords_size = poly_num_coords * compression_unit_size(ic);
    next_poly_coords += poly_coords_size;

    area += area_polygon(poly_coords,
                         poly_coords_size,
                         reinterpret_cast<int8_t*>(poly_ring_sizes),
                         poly_num_rings,
                         ic,
                         isr,
                         osr);
  }
  return area;
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Area_Polygon	(	int8_t *	poly_coords,
		int32_t	poly_coords_size,
		int8_t *	poly_ring_sizes,
		int32_t	poly_num_rings,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1409 of file ExtensionFunctionsGeo.hpp.

References area_polygon().

                                    {
  return area_polygon(
      poly_coords, poly_coords_size, poly_ring_sizes, poly_num_rings, ic, isr, osr);
}

Here is the call graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsBinary__GeoLineString__cpu_	(	RowFunctionManager &	mgr,
		GeoLineString &	l
	)

Definition at line 5326 of file ExtensionFunctionsGeo.hpp.

References __wkb_to_str(), and to_Geospatial_GeoLineString().

Referenced by ST_AsWkb__GeoLineString__cpu_().

                                                                            {
  Geospatial::GeoLineString linestring = to_Geospatial_GeoLineString(l);
  std::string wkb = __wkb_to_str(linestring);
  return TextEncodingNone(mgr, wkb);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsBinary__GeoMultiLineString__cpu_	(	RowFunctionManager &	mgr,
		GeoMultiLineString &	ml
	)

Definition at line 5333 of file ExtensionFunctionsGeo.hpp.

References __wkb_to_str(), and to_Geospatial_GeoMultiLineString().

Referenced by ST_AsWkb__GeoMultiLineString__cpu_().

                                                                                       {
  Geospatial::GeoMultiLineString multilinestring = to_Geospatial_GeoMultiLineString(ml);
  std::string wkb = __wkb_to_str(multilinestring);
  return TextEncodingNone(mgr, wkb);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsBinary__GeoMultiPoint__cpu_	(	RowFunctionManager &	mgr,
		GeoMultiPoint &	mp
	)

Definition at line 5319 of file ExtensionFunctionsGeo.hpp.

References __wkb_to_str(), and to_Geospatial_GeoMultiPoint().

Referenced by ST_AsWkb__GeoMultiPoint__cpu_().

                                                                             {
  Geospatial::GeoMultiPoint multipoint = to_Geospatial_GeoMultiPoint(mp);
  std::string wkb = __wkb_to_str(multipoint);
  return TextEncodingNone(mgr, wkb);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsBinary__GeoMultiPolygon__cpu_	(	RowFunctionManager &	mgr,
		GeoMultiPolygon &	mp
	)

Definition at line 5347 of file ExtensionFunctionsGeo.hpp.

References __wkb_to_str(), and to_Geospatial_GeoMultiPolygon().

Referenced by ST_AsWkb__GeoMultiPolygon__cpu_().

                                                                                 {
  Geospatial::GeoMultiPolygon multipolygon = to_Geospatial_GeoMultiPolygon(mp);
  std::string wkb = __wkb_to_str(multipolygon);
  return TextEncodingNone(mgr, wkb);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsBinary__GeoPoint__cpu_	(	RowFunctionManager &	mgr,
		GeoPoint &	p
	)

Definition at line 5311 of file ExtensionFunctionsGeo.hpp.

References __wkb_to_str(), and to_Geospatial_GeoPoint().

Referenced by ST_AsWkb__GeoPoint__cpu_().

                                                                             {
  Geospatial::GeoPoint point = to_Geospatial_GeoPoint(p);
  std::string wkb = __wkb_to_str(point);
  return TextEncodingNone(mgr, wkb);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsBinary__GeoPolygon__cpu_	(	RowFunctionManager &	mgr,
		GeoPolygon &	p
	)

Definition at line 5339 of file ExtensionFunctionsGeo.hpp.

References __wkb_to_str(), and to_Geospatial_GeoPolygon().

Referenced by ST_AsWkb__GeoPolygon__cpu_().

                                                                                 {
  Geospatial::GeoPolygon polygon = to_Geospatial_GeoPolygon(p);
  std::string wkb = __wkb_to_str(polygon);
  return TextEncodingNone(mgr, wkb);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsText__GeoLineString__cpu_	(	RowFunctionManager &	mgr,
		GeoLineString &	l
	)

Definition at line 5220 of file ExtensionFunctionsGeo.hpp.

References Geospatial::GeoBase::getWktString(), and to_Geospatial_GeoLineString().

Referenced by ST_AsWkt__GeoLineString__cpu_().

                                                                          {
  Geospatial::GeoLineString linestring = to_Geospatial_GeoLineString(l);
  std::string wkt = linestring.getWktString();
  return TextEncodingNone(mgr, wkt);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsText__GeoMultiLineString__cpu_	(	RowFunctionManager &	mgr,
		GeoMultiLineString &	ml
	)

Definition at line 5227 of file ExtensionFunctionsGeo.hpp.

References Geospatial::GeoBase::getWktString(), and to_Geospatial_GeoMultiLineString().

Referenced by ST_AsWkt__GeoMultiLineString__cpu_().

                                                                                     {
  Geospatial::GeoMultiLineString multilinestring = to_Geospatial_GeoMultiLineString(ml);
  std::string wkt = multilinestring.getWktString();
  return TextEncodingNone(mgr, wkt);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsText__GeoMultiPoint__cpu_	(	RowFunctionManager &	mgr,
		GeoMultiPoint &	mp
	)

Definition at line 5213 of file ExtensionFunctionsGeo.hpp.

References Geospatial::GeoBase::getWktString(), and to_Geospatial_GeoMultiPoint().

Referenced by ST_AsWkt__GeoMultiPoint__cpu_().

                                                                           {
  Geospatial::GeoMultiPoint multipoint = to_Geospatial_GeoMultiPoint(mp);
  std::string wkt = multipoint.getWktString();
  return TextEncodingNone(mgr, wkt);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsText__GeoMultiPolygon__cpu_	(	RowFunctionManager &	mgr,
		GeoMultiPolygon &	mp
	)

Definition at line 5241 of file ExtensionFunctionsGeo.hpp.

References Geospatial::GeoBase::getWktString(), and to_Geospatial_GeoMultiPolygon().

Referenced by ST_AsWkt__GeoMultiPolygon__cpu_().

                                                                               {
  Geospatial::GeoMultiPolygon multipolygon = to_Geospatial_GeoMultiPolygon(mp);
  std::string wkt = multipolygon.getWktString();
  return TextEncodingNone(mgr, wkt);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsText__GeoPoint__cpu_	(	RowFunctionManager &	mgr,
		GeoPoint &	p
	)

Definition at line 5205 of file ExtensionFunctionsGeo.hpp.

References Geospatial::GeoBase::getWktString(), and to_Geospatial_GeoPoint().

Referenced by ST_AsWkt__GeoPoint__cpu_().

                                                                           {
  Geospatial::GeoPoint point = to_Geospatial_GeoPoint(p);
  std::string wkt = point.getWktString();
  return TextEncodingNone(mgr, wkt);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsText__GeoPolygon__cpu_	(	RowFunctionManager &	mgr,
		GeoPolygon &	p
	)

Definition at line 5233 of file ExtensionFunctionsGeo.hpp.

References Geospatial::GeoBase::getWktString(), and to_Geospatial_GeoPolygon().

Referenced by ST_AsWkt__GeoPolygon__cpu_().

                                                                               {
  Geospatial::GeoPolygon polygon = to_Geospatial_GeoPolygon(p);
  std::string wkt = polygon.getWktString();
  return TextEncodingNone(mgr, wkt);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsWkb__GeoLineString__cpu_	(	RowFunctionManager &	mgr,
		GeoLineString &	p
	)

Definition at line 5371 of file ExtensionFunctionsGeo.hpp.

References ST_AsBinary__GeoLineString__cpu_().

                                                                                    {
  return ST_AsBinary__GeoLineString__cpu_(mgr, p);
}

Here is the call graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsWkb__GeoMultiLineString__cpu_	(	RowFunctionManager &	mgr,
		GeoMultiLineString &	p
	)

Definition at line 5377 of file ExtensionFunctionsGeo.hpp.

References ST_AsBinary__GeoMultiLineString__cpu_().

                                                                                   {
  return ST_AsBinary__GeoMultiLineString__cpu_(mgr, p);
}

Here is the call graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsWkb__GeoMultiPoint__cpu_	(	RowFunctionManager &	mgr,
		GeoMultiPoint &	p
	)

Definition at line 5366 of file ExtensionFunctionsGeo.hpp.

References ST_AsBinary__GeoMultiPoint__cpu_().

                                                                                    {
  return ST_AsBinary__GeoMultiPoint__cpu_(mgr, p);
}

Here is the call graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsWkb__GeoMultiPolygon__cpu_	(	RowFunctionManager &	mgr,
		GeoMultiPolygon &	p
	)

Definition at line 5387 of file ExtensionFunctionsGeo.hpp.

References ST_AsBinary__GeoMultiPolygon__cpu_().

                                                                             {
  return ST_AsBinary__GeoMultiPolygon__cpu_(mgr, p);
}

Here is the call graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsWkb__GeoPoint__cpu_	(	RowFunctionManager &	mgr,
		GeoPoint &	p
	)

Definition at line 5361 of file ExtensionFunctionsGeo.hpp.

References ST_AsBinary__GeoPoint__cpu_().

                                                                          {
  return ST_AsBinary__GeoPoint__cpu_(mgr, p);
}

Here is the call graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsWkb__GeoPolygon__cpu_	(	RowFunctionManager &	mgr,
		GeoPolygon &	p
	)

Definition at line 5381 of file ExtensionFunctionsGeo.hpp.

References ST_AsBinary__GeoPolygon__cpu_().

                                                                              {
  return ST_AsBinary__GeoPolygon__cpu_(mgr, p);
}

Here is the call graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsWkt__GeoLineString__cpu_	(	RowFunctionManager &	mgr,
		GeoLineString &	p
	)

Definition at line 5265 of file ExtensionFunctionsGeo.hpp.

References ST_AsText__GeoLineString__cpu_().

                                                                                    {
  return ST_AsText__GeoLineString__cpu_(mgr, p);
}

Here is the call graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsWkt__GeoMultiLineString__cpu_	(	RowFunctionManager &	mgr,
		GeoMultiLineString &	p
	)

Definition at line 5271 of file ExtensionFunctionsGeo.hpp.

References ST_AsText__GeoMultiLineString__cpu_().

                                                                                   {
  return ST_AsText__GeoMultiLineString__cpu_(mgr, p);
}

Here is the call graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsWkt__GeoMultiPoint__cpu_	(	RowFunctionManager &	mgr,
		GeoMultiPoint &	p
	)

Definition at line 5260 of file ExtensionFunctionsGeo.hpp.

References ST_AsText__GeoMultiPoint__cpu_().

                                                                                    {
  return ST_AsText__GeoMultiPoint__cpu_(mgr, p);
}

Here is the call graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsWkt__GeoMultiPolygon__cpu_	(	RowFunctionManager &	mgr,
		GeoMultiPolygon &	p
	)

Definition at line 5281 of file ExtensionFunctionsGeo.hpp.

References ST_AsText__GeoMultiPolygon__cpu_().

                                                                             {
  return ST_AsText__GeoMultiPolygon__cpu_(mgr, p);
}

Here is the call graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsWkt__GeoPoint__cpu_	(	RowFunctionManager &	mgr,
		GeoPoint &	p
	)

Definition at line 5255 of file ExtensionFunctionsGeo.hpp.

References ST_AsText__GeoPoint__cpu_().

                                                                          {
  return ST_AsText__GeoPoint__cpu_(mgr, p);
}

Here is the call graph for this function:

EXTENSION_NOINLINE TextEncodingNone ST_AsWkt__GeoPolygon__cpu_	(	RowFunctionManager &	mgr,
		GeoPolygon &	p
	)

Definition at line 5275 of file ExtensionFunctionsGeo.hpp.

References ST_AsText__GeoPolygon__cpu_().

                                                                              {
  return ST_AsText__GeoPolygon__cpu_(mgr, p);
}

Here is the call graph for this function:

EXTENSION_NOINLINE bool ST_cContains_MultiPolygon_Point	(	int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4287 of file ExtensionFunctionsGeo.hpp.

                                                                     {
  return Contains_MultiPolygon_Point_Impl<int64_t, EdgeBehavior::kExcludePointOnEdge>(
      mpoly_coords,
      mpoly_coords_size,
      mpoly_ring_sizes,
      mpoly_num_rings,
      mpoly_poly_sizes,
      mpoly_num_polys,
      mpoly_bounds,
      mpoly_bounds_size,
      p,
      psize,
      ic1,
      isr1,
      ic2,
      isr2,
      osr);
}

EXTENSION_NOINLINE bool ST_cContains_Polygon_Point	(	const int8_t *	poly_coords,
		const int64_t	poly_coords_size,
		const int32_t *	poly_ring_sizes,
		const int64_t	poly_num_rings,
		const double *	poly_bounds,
		const int64_t	poly_bounds_size,
		const int8_t *	p,
		const int64_t	psize,
		const int32_t	ic1,
		const int32_t	isr1,
		const int32_t	ic2,
		const int32_t	isr2,
		const int32_t	osr
	)

Definition at line 4070 of file ExtensionFunctionsGeo.hpp.

                                                                      {
  return Contains_Polygon_Point_Impl<int64_t, EdgeBehavior::kExcludePointOnEdge>(
      poly_coords,
      poly_coords_size,
      poly_ring_sizes,
      poly_num_rings,
      poly_bounds,
      poly_bounds_size,
      p,
      psize,
      ic1,
      isr1,
      ic2,
      isr2,
      osr);
}

EXTENSION_NOINLINE void ST_Centroid_LineString	(	int8_t *	coords,
		int32_t	coords_sz,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr,
		double *	linestring_centroid
	)

Definition at line 1560 of file ExtensionFunctionsGeo.hpp.

References centroid_add_linestring().

                                                         {
  double length = 0.0;
  double linestring_centroid_sum[2] = {0.0, 0.0};
  int64_t num_points = 0;
  double point_centroid_sum[2] = {0.0, 0.0};
  centroid_add_linestring(coords,
                          coords_sz,
                          ic,
                          isr,
                          osr,
                          false,  // not closed
                          &length,
                          &linestring_centroid_sum[0],
                          &num_points,
                          &point_centroid_sum[0]);
  if (length > 0) {
    linestring_centroid[0] = linestring_centroid_sum[0] / length;
    linestring_centroid[1] = linestring_centroid_sum[1] / length;
  } else if (num_points > 0) {
    linestring_centroid[0] = point_centroid_sum[0] / num_points;
    linestring_centroid[1] = point_centroid_sum[1] / num_points;
  }
}

Here is the call graph for this function:

EXTENSION_NOINLINE void ST_Centroid_MultiPoint	(	int8_t *	mp,
		int32_t	mpsize,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr,
		double *	multipoint_centroid
	)

Definition at line 1493 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), and get_point().

                                                         {
  auto mp_num_coords = mpsize / compression_unit_size(ic);
  double x = 0.0;
  double y = 0.0;
  for (int32_t i = 0; i < mp_num_coords; i += 2) {
    Point2D mpp = get_point(mp, i, ic, isr, osr);
    x += mpp.x;  // In case of overflows on large geometries,
    y += mpp.y;  // move division into the loop
  }
  auto mp_num_points = mp_num_coords >> 1;
  multipoint_centroid[0] = x / mp_num_points;
  multipoint_centroid[1] = y / mp_num_points;
}

Here is the call graph for this function:

EXTENSION_NOINLINE void ST_Centroid_MultiPolygon	(	int8_t *	mpoly_coords,
		int32_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int32_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int32_t	mpoly_num_polys,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr,
		double *	mpoly_centroid
	)

Definition at line 1739 of file ExtensionFunctionsGeo.hpp.

References centroid_add_polygon(), compression_unit_size(), and get_point().

                                                      {
  if (mpoly_num_rings <= 0 || mpoly_num_polys <= 0) {
    mpoly_centroid[0] = 0.0;
    mpoly_centroid[1] = 0.0;
  }

  double total_area2 = 0.0;
  double cg3[2] = {0.0, 0.0};
  double total_length = 0.0;
  double linestring_centroid_sum[2] = {0.0, 0.0};
  int64_t num_points = 0;
  double point_centroid_sum[2] = {0.0, 0.0};

  // Set specific poly pointers as we move through the coords/ringsizes/polyrings arrays.
  auto next_poly_coords = mpoly_coords;
  auto next_poly_ring_sizes = mpoly_ring_sizes;

  for (auto poly = 0; poly < mpoly_num_polys; poly++) {
    auto poly_coords = next_poly_coords;
    auto poly_ring_sizes = next_poly_ring_sizes;
    auto poly_num_rings = mpoly_poly_sizes[poly];
    // Count number of coords in all of poly's rings, advance ring size pointer.
    int32_t poly_num_coords = 0;
    for (auto ring = 0; ring < poly_num_rings; ring++) {
      poly_num_coords += 2 * *next_poly_ring_sizes++;
    }
    auto poly_coords_size = poly_num_coords * compression_unit_size(ic);
    next_poly_coords += poly_coords_size;

    centroid_add_polygon(poly_coords,
                         poly_coords_size,
                         poly_ring_sizes,
                         poly_num_rings,
                         ic,
                         isr,
                         osr,
                         &total_area2,
                         &cg3[0],
                         &total_length,
                         &linestring_centroid_sum[0],
                         &num_points,
                         &point_centroid_sum[0]);
  }

  if (total_area2 != 0.0) {
    mpoly_centroid[0] = cg3[0] / 3 / total_area2;
    mpoly_centroid[1] = cg3[1] / 3 / total_area2;
  } else if (total_length > 0.0) {
    // zero-area multipolygon, fall back to linear centroid
    mpoly_centroid[0] = linestring_centroid_sum[0] / total_length;
    mpoly_centroid[1] = linestring_centroid_sum[1] / total_length;
  } else if (num_points > 0) {
    // zero-area zero-length multipolygon, fall further back to point centroid
    mpoly_centroid[0] = point_centroid_sum[0] / num_points;
    mpoly_centroid[1] = point_centroid_sum[1] / num_points;
  } else {
    Point2D centroid = get_point(mpoly_coords, 0, ic, isr, osr);
    mpoly_centroid[0] = centroid.x;
    mpoly_centroid[1] = centroid.y;
  }
}

Here is the call graph for this function:

EXTENSION_NOINLINE void ST_Centroid_Point	(	int8_t *	p,
		int32_t	psize,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr,
		double *	point_centroid
	)

Definition at line 1481 of file ExtensionFunctionsGeo.hpp.

References get_point().

                                               {
  Point2D const centroid = get_point(p, 0, ic, isr, osr);
  point_centroid[0] = centroid.x;
  point_centroid[1] = centroid.y;
}

Here is the call graph for this function:

EXTENSION_NOINLINE void ST_Centroid_Polygon	(	int8_t *	poly_coords,
		int32_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int32_t	poly_num_rings,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr,
		double *	poly_centroid
	)

Definition at line 1688 of file ExtensionFunctionsGeo.hpp.

References centroid_add_polygon(), and get_point().

                                                {
  if (poly_num_rings <= 0) {
    poly_centroid[0] = 0.0;
    poly_centroid[1] = 0.0;
  }
  double total_area2 = 0.0;
  double cg3[2] = {0.0, 0.0};
  double total_length = 0.0;
  double linestring_centroid_sum[2] = {0.0, 0.0};
  int64_t num_points = 0;
  double point_centroid_sum[2] = {0.0, 0.0};
  centroid_add_polygon(poly_coords,
                       poly_coords_size,
                       poly_ring_sizes,
                       poly_num_rings,
                       ic,
                       isr,
                       osr,
                       &total_area2,
                       &cg3[0],
                       &total_length,
                       &linestring_centroid_sum[0],
                       &num_points,
                       &point_centroid_sum[0]);

  if (total_area2 != 0.0) {
    poly_centroid[0] = cg3[0] / 3 / total_area2;
    poly_centroid[1] = cg3[1] / 3 / total_area2;
  } else if (total_length > 0.0) {
    // zero-area polygon, fall back to linear centroid
    poly_centroid[0] = linestring_centroid_sum[0] / total_length;
    poly_centroid[1] = linestring_centroid_sum[1] / total_length;
  } else if (num_points > 0) {
    // zero-area zero-length polygon, fall further back to point centroid
    poly_centroid[0] = point_centroid_sum[0] / num_points;
    poly_centroid[1] = point_centroid_sum[1] / num_points;
  } else {
    Point2D centroid = get_point(poly_coords, 0, ic, isr, osr);
    poly_centroid[0] = centroid.x;
    poly_centroid[1] = centroid.y;
  }
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_cIntersects_MultiPolygon_Point	(	int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4834 of file ExtensionFunctionsGeo.hpp.

                                                    {
  return Contains_MultiPolygon_Point_Impl<int64_t, EdgeBehavior::kIncludePointOnEdge>(
      mpoly_coords,
      mpoly_coords_size,
      mpoly_ring_sizes,
      mpoly_num_rings,
      mpoly_poly_sizes,
      mpoly_num_polys,
      mpoly_bounds,
      mpoly_bounds_size,
      p,
      psize,
      ic1,
      isr1,
      ic2,
      isr2,
      osr);
}

EXTENSION_INLINE bool ST_cIntersects_Polygon_Point	(	int8_t *	poly,
		int64_t	polysize,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4651 of file ExtensionFunctionsGeo.hpp.

                                               {
  return Contains_Polygon_Point_Impl<int64_t, EdgeBehavior::kIncludePointOnEdge>(
      poly,
      polysize,
      poly_ring_sizes,
      poly_num_rings,
      poly_bounds,
      poly_bounds_size,
      p,
      psize,
      ic1,
      isr1,
      ic2,
      isr2,
      osr);
}

EXTENSION_INLINE bool ST_Contains_LineString_Point	(	int8_t *	l,
		int64_t	lsize,
		double *	lbounds,
		int64_t	lbounds_size,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 3938 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_Point_LineString(), and tol_zero().

                                               {
  return tol_zero(
      ST_Distance_Point_LineString(p, psize, l, lsize, ic2, isr2, ic1, isr1, osr, 0.0));
}

Here is the call graph for this function:

EXTENSION_NOINLINE bool ST_Contains_LineString_Polygon	(	int8_t *	l,
		int64_t	lsize,
		double *	lbounds,
		int64_t	lbounds_size,
		int8_t *	poly_coords,
		int64_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 3954 of file ExtensionFunctionsGeo.hpp.

                                                 {
  // TODO
  return false;
}

EXTENSION_NOINLINE bool ST_Contains_MultiPolygon_LineString	(	int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int8_t *	l,
		int64_t	lsize,
		double *	lbounds,
		int64_t	lbounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4321 of file ExtensionFunctionsGeo.hpp.

References box_contains_box(), compression_unit_size(), and ST_Contains_Polygon_LineString().

                                                      {
  if (mpoly_num_polys <= 0) {
    return false;
  }

  if (mpoly_bounds && lbounds) {
    if (!box_contains_box(mpoly_bounds, mpoly_bounds_size, lbounds, lbounds_size)) {
      return false;
    }
  }

  // Set specific poly pointers as we move through the coords/ringsizes/polyrings arrays.
  auto next_poly_coords = mpoly_coords;
  auto next_poly_ring_sizes = mpoly_ring_sizes;

  for (auto poly = 0; poly < mpoly_num_polys; poly++) {
    auto poly_coords = next_poly_coords;
    auto poly_ring_sizes = next_poly_ring_sizes;
    auto poly_num_rings = mpoly_poly_sizes[poly];
    // Count number of coords in all of poly's rings, advance ring size pointer.
    int32_t poly_num_coords = 0;
    for (auto ring = 0; ring < poly_num_rings; ring++) {
      poly_num_coords += 2 * *next_poly_ring_sizes++;
    }
    auto poly_coords_size = poly_num_coords * compression_unit_size(ic1);
    next_poly_coords += poly_coords_size;

    if (ST_Contains_Polygon_LineString(poly_coords,
                                       poly_coords_size,
                                       poly_ring_sizes,
                                       poly_num_rings,
                                       nullptr,
                                       0,
                                       l,
                                       lsize,
                                       nullptr,
                                       0,
                                       ic1,
                                       isr1,
                                       ic2,
                                       isr2,
                                       osr)) {
      return true;
    }
  }

  return false;
}

Here is the call graph for this function:

EXTENSION_NOINLINE bool ST_Contains_MultiPolygon_Point	(	int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4254 of file ExtensionFunctionsGeo.hpp.

                                                 {
  return Contains_MultiPolygon_Point_Impl<double, EdgeBehavior::kExcludePointOnEdge>(
      mpoly_coords,
      mpoly_coords_size,
      mpoly_ring_sizes,
      mpoly_num_rings,
      mpoly_poly_sizes,
      mpoly_num_polys,
      mpoly_bounds,
      mpoly_bounds_size,
      p,
      psize,
      ic1,
      isr1,
      ic2,
      isr2,
      osr);
}

EXTENSION_NOINLINE bool ST_Contains_Point_LineString	(	int8_t *	p,
		int64_t	psize,
		int8_t *	l,
		int64_t	lsize,
		double *	lbounds,
		int64_t	lbounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 3873 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), get_point(), and tol_eq().

Referenced by ST_Contains_Point_Polygon().

                                               {
  Point2D const pt = get_point(p, 0, ic1, isr1, osr);

  if (lbounds) {
    if (tol_eq(pt.x, lbounds[0]) && tol_eq(pt.y, lbounds[1]) &&
        tol_eq(pt.x, lbounds[2]) && tol_eq(pt.y, lbounds[3])) {
      return true;
    }
  }

  auto l_num_coords = lsize / compression_unit_size(ic2);
  for (int i = 0; i < l_num_coords; i += 2) {
    Point2D pl = get_point(l, i, ic2, isr2, osr);
    if (tol_eq(pt.x, pl.x) && tol_eq(pt.y, pl.y)) {
      continue;
    }
    return false;
  }
  return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE bool ST_Contains_Point_Point	(	int8_t *	p1,
		int64_t	p1size,
		int8_t *	p2,
		int64_t	p2size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 3857 of file ExtensionFunctionsGeo.hpp.

References get_point(), tol(), and tol_eq().

                                          {
  Point2D const pt1 = get_point(p1, 0, ic1, isr1, osr);
  Point2D const pt2 = get_point(p2, 0, ic2, isr2, osr);
  double tolerance = tol(ic1, ic2);
  return tol_eq(pt1.x, pt2.x, tolerance) && tol_eq(pt1.y, pt2.y, tolerance);
}

Here is the call graph for this function:

EXTENSION_NOINLINE bool ST_Contains_Point_Polygon	(	int8_t *	p,
		int64_t	psize,
		int8_t *	poly_coords,
		int64_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 3905 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), and ST_Contains_Point_LineString().

                                            {
  auto exterior_ring_num_coords = poly_coords_size / compression_unit_size(ic2);
  if (poly_num_rings > 0) {
    exterior_ring_num_coords = poly_ring_sizes[0] * 2;
  }
  auto exterior_ring_coords_size = exterior_ring_num_coords * compression_unit_size(ic2);

  return ST_Contains_Point_LineString(p,
                                      psize,
                                      poly_coords,
                                      exterior_ring_coords_size,
                                      poly_bounds,
                                      poly_bounds_size,
                                      ic1,
                                      isr1,
                                      ic2,
                                      isr2,
                                      osr);
}

Here is the call graph for this function:

EXTENSION_NOINLINE bool ST_Contains_Polygon_LineString	(	int8_t *	poly_coords,
		int64_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int8_t *	l,
		int64_t	lsize,
		double *	lbounds,
		int64_t	lbounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4100 of file ExtensionFunctionsGeo.hpp.

References box_contains_box(), compression_unit_size(), and polygon_contains_linestring().

Referenced by ST_Contains_MultiPolygon_LineString(), and ST_Contains_Polygon_Polygon().

                                                 {
  if (poly_num_rings > 1) {
    return false;  // TODO: support polygons with interior rings
  }

  // Bail out if poly bounding box doesn't contain linestring bounding box
  if (poly_bounds && lbounds) {
    if (!box_contains_box(poly_bounds, poly_bounds_size, lbounds, lbounds_size)) {
      return false;
    }
  }

  const auto poly_num_coords = poly_coords_size / compression_unit_size(ic1);
  const auto lnum_coords = lsize / compression_unit_size(ic2);

  return polygon_contains_linestring(
      poly_coords, poly_num_coords, l, lnum_coords, ic1, isr1, ic2, isr2, osr);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE bool ST_Contains_Polygon_Point	(	const int8_t *	poly_coords,
		const int64_t	poly_coords_size,
		const int32_t *	poly_ring_sizes,
		const int64_t	poly_num_rings,
		const double *	poly_bounds,
		const int64_t	poly_bounds_size,
		const int8_t *	p,
		const int64_t	psize,
		const int32_t	ic1,
		const int32_t	isr1,
		const int32_t	ic2,
		const int32_t	isr2,
		const int32_t	osr
	)

Definition at line 4041 of file ExtensionFunctionsGeo.hpp.

                                                                     {
  return Contains_Polygon_Point_Impl<double, EdgeBehavior::kExcludePointOnEdge>(
      poly_coords,
      poly_coords_size,
      poly_ring_sizes,
      poly_num_rings,
      poly_bounds,
      poly_bounds_size,
      p,
      psize,
      ic1,
      isr1,
      ic2,
      isr2,
      osr);
}

EXTENSION_NOINLINE bool ST_Contains_Polygon_Polygon	(	int8_t *	poly1_coords,
		int64_t	poly1_coords_size,
		int32_t *	poly1_ring_sizes,
		int64_t	poly1_num_rings,
		double *	poly1_bounds,
		int64_t	poly1_bounds_size,
		int8_t *	poly2_coords,
		int64_t	poly2_coords_size,
		int32_t *	poly2_ring_sizes,
		int64_t	poly2_num_rings,
		double *	poly2_bounds,
		int64_t	poly2_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4134 of file ExtensionFunctionsGeo.hpp.

References box_contains_box(), compression_unit_size(), and ST_Contains_Polygon_LineString().

                                              {
  // TODO: needs to be extended, cover more cases
  // Right now only checking if simple poly1 (no holes) contains poly2's exterior shape
  if (poly1_num_rings > 1) {
    return false;  // TODO: support polygons with interior rings
  }

  if (poly1_bounds && poly2_bounds) {
    if (!box_contains_box(
            poly1_bounds, poly1_bounds_size, poly2_bounds, poly2_bounds_size)) {
      return false;
    }
  }

  int64_t poly2_exterior_ring_coords_size = poly2_coords_size;
  if (poly2_num_rings > 0) {
    poly2_exterior_ring_coords_size =
        2 * poly2_ring_sizes[0] * compression_unit_size(ic2);
  }
  return ST_Contains_Polygon_LineString(poly1_coords,
                                        poly1_coords_size,
                                        poly1_ring_sizes,
                                        poly1_num_rings,
                                        poly1_bounds,
                                        poly1_bounds_size,
                                        poly2_coords,
                                        poly2_exterior_ring_coords_size,
                                        poly2_bounds,
                                        poly2_bounds_size,
                                        ic1,
                                        isr1,
                                        ic2,
                                        isr2,
                                        osr);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Distance_LineString_LineString	(	int8_t *	l1,
		int64_t	l1size,
		int8_t *	l2,
		int64_t	l2size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2669 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), distance_line_line(), get_point(), tol_zero(), and TOLERANCE_DEFAULT_SQUARED.

Referenced by ST_DWithin_LineString_LineString(), and ST_Intersects_LineString_Linestring().

                                                           {
  auto l1_num_coords = l1size / compression_unit_size(ic1);
  auto l2_num_coords = l2size / compression_unit_size(ic2);

  double threshold_squared = threshold * threshold;
  double dist_squared = 0.0;
  Point2D l11 = get_point(l1, 0, ic1, isr1, osr);
  for (int32_t i1 = 2; i1 < l1_num_coords; i1 += 2) {
    Point2D l12 = get_point(l1, i1, ic1, isr1, osr);
    Point2D l21 = get_point(l2, 0, ic2, isr2, osr);
    for (int32_t i2 = 2; i2 < l2_num_coords; i2 += 2) {
      Point2D l22 = get_point(l2, i2, ic2, isr2, osr);

      // double ldist_squared =
      //    distance_line_line_squared(l11x, l11y, l12x, l12y, l21x, l21y, l22x, l22y);
      // TODO: fix distance_line_line_squared
      double ldist =
          distance_line_line(l11.x, l11.y, l12.x, l12.y, l21.x, l21.y, l22.x, l22.y);
      double ldist_squared = ldist * ldist;

      if (i1 == 2 && i2 == 2) {
        dist_squared = ldist_squared;  // initialize dist with distance between the first
                                       // two segments
      } else if (dist_squared > ldist_squared) {
        dist_squared = ldist_squared;
      }
      if (tol_zero(dist_squared, TOLERANCE_DEFAULT_SQUARED)) {
        return 0.0;  // Segments touch, short-circuit
      }
      if (dist_squared <= threshold_squared) {
        // If threashold_squared is defined and the calculated dist_squared dips under it,
        // short-circuit and return it
        return sqrt(dist_squared);
      }

      l21 = l22;  // advance to the next point on l2
    }

    l11 = l12;  // advance to the next point on l1
  }
  return sqrt(dist_squared);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE double ST_Distance_LineString_MultiPoint	(	int8_t *	l,
		int64_t	lsize,
		int8_t *	mp,
		int64_t	mpsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2627 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_MultiPoint_LineString().

                                                           {
  return ST_Distance_MultiPoint_LineString(
      mp, mpsize, l, lsize, ic2, isr2, ic1, isr1, osr, threshold);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Distance_LineString_MultiPolygon	(	int8_t *	l,
		int64_t	lsize,
		int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2791 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), ST_Distance_LineString_Polygon(), and tol_zero().

Referenced by ST_Distance_MultiPolygon_LineString(), ST_DWithin_LineString_MultiPolygon(), and ST_Intersects_LineString_MultiPolygon().

                                                             {
  // TODO: revisit implementation, cover all cases
  double min_distance = 0.0;

  // Set specific poly pointers as we move through the coords/ringsizes/polyrings arrays.
  auto next_poly_coords = mpoly_coords;
  auto next_poly_ring_sizes = mpoly_ring_sizes;

  for (auto poly = 0; poly < mpoly_num_polys; poly++) {
    auto poly_coords = next_poly_coords;
    auto poly_ring_sizes = next_poly_ring_sizes;
    auto poly_num_rings = mpoly_poly_sizes[poly];
    // Count number of coords in all of poly's rings, advance ring size pointer.
    int32_t poly_num_coords = 0;
    for (auto ring = 0; ring < poly_num_rings; ring++) {
      poly_num_coords += 2 * *next_poly_ring_sizes++;
    }
    auto poly_coords_size = poly_num_coords * compression_unit_size(ic2);
    next_poly_coords += poly_coords_size;
    double distance = ST_Distance_LineString_Polygon(l,
                                                     lsize,
                                                     poly_coords,
                                                     poly_coords_size,
                                                     poly_ring_sizes,
                                                     poly_num_rings,
                                                     ic1,
                                                     isr1,
                                                     ic2,
                                                     isr2,
                                                     osr,
                                                     threshold);
    if (poly == 0 || min_distance > distance) {
      min_distance = distance;
      if (tol_zero(min_distance)) {
        min_distance = 0.0;
        break;
      }
      if (min_distance <= threshold) {
        return min_distance;
      }
    }
  }

  return min_distance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE double ST_Distance_LineString_Point	(	int8_t *	l,
		int64_t	lsize,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2612 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_Point_LineString().

                                                      {
  return ST_Distance_Point_LineString(
      p, psize, l, lsize, ic2, isr2, ic1, isr1, osr, threshold);
}

Here is the call graph for this function:

EXTENSION_INLINE double ST_Distance_LineString_Point_Geodesic	(	int8_t *	l,
		int64_t	lsize,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 1877 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_Point_LineString_Geodesic().

                                                          {
  // Currently only statically indexed LineString is supported
  return ST_Distance_Point_LineString_Geodesic(
      p, psize, l, lsize, ic2, isr2, ic1, isr1, osr);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Distance_LineString_Polygon	(	int8_t *	l,
		int64_t	lsize,
		int8_t *	poly_coords,
		int64_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2722 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), distance_ring_linestring(), ST_Distance_Point_Polygon(), and tol_zero().

Referenced by ST_Distance_LineString_MultiPolygon(), ST_Distance_Polygon_LineString(), ST_DWithin_LineString_Polygon(), and ST_Intersects_LineString_Polygon().

                                                        {
  auto lnum_coords = lsize / compression_unit_size(ic1);
  auto psize = 2 * compression_unit_size(ic1);
  auto min_distance =
      ST_Distance_Point_Polygon(l,  // pointer to start of linestring for first point
                                psize,
                                poly_coords,
                                poly_coords_size,
                                poly_ring_sizes,
                                poly_num_rings,
                                ic1,
                                isr1,
                                ic2,
                                isr2,
                                osr,
                                threshold);
  if (tol_zero(min_distance)) {
    // Linestring's first point is inside the poly
    return 0.0;
  }
  if (min_distance <= threshold) {
    return min_distance;
  }

  // Otherwise, linestring's first point is outside the external ring or inside
  // an internal ring. Measure minimum distance between linestring segments and
  // poly rings. Crossing a ring zeroes the distance and causes an early return.
  auto poly_ring_coords = poly_coords;
  for (auto r = 0; r < poly_num_rings; r++) {
    int64_t poly_ring_num_coords = poly_ring_sizes[r] * 2;

    auto distance = distance_ring_linestring(poly_ring_coords,
                                             poly_ring_num_coords,
                                             l,
                                             lnum_coords,
                                             ic2,
                                             isr2,
                                             ic1,
                                             isr1,
                                             osr,
                                             threshold);
    if (min_distance > distance) {
      min_distance = distance;
      if (tol_zero(min_distance)) {
        return 0.0;
      }
      if (min_distance <= threshold) {
        return min_distance;
      }
    }

    poly_ring_coords += poly_ring_num_coords * compression_unit_size(ic2);
  }

  return min_distance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_MultiLineString_MultiLineString	(	int8_t *	mls1,
		int64_t	mls1size,
		int32_t *	mls1_ls_sizes,
		int64_t	mls1_ls_num,
		int8_t *	mls2,
		int64_t	mls2size,
		int32_t *	mls2_ls_sizes,
		int64_t	mls2_ls_num,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 3245 of file ExtensionFunctionsGeo.hpp.

                                                                     {
  return 0.0;
}

EXTENSION_INLINE double ST_Distance_MultiLineString_MultiPoint	(	int8_t *	mls,
		int64_t	mls_size,
		int32_t *	mls_ls_sizes,
		int64_t	mls_ls_num,
		int8_t *	mp,
		int64_t	mpsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2642 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_MultiPoint_MultiLineString().

                                                                {
  return ST_Distance_MultiPoint_MultiLineString(mp,
                                                mpsize,
                                                mls,
                                                mls_size,
                                                mls_ls_sizes,
                                                mls_ls_num,
                                                ic2,
                                                isr2,
                                                ic1,
                                                isr1,
                                                osr,
                                                threshold);
}

Here is the call graph for this function:

EXTENSION_INLINE double ST_Distance_MultiLineString_Point	(	int8_t *	mls,
		int64_t	mls_size,
		int32_t *	mls_ls_sizes,
		int64_t	mls_ls_num,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2851 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_Point_MultiLineString().

                                                           {
  return ST_Distance_Point_MultiLineString(p,
                                           psize,
                                           mls,
                                           mls_size,
                                           mls_ls_sizes,
                                           mls_ls_num,
                                           ic2,
                                           isr2,
                                           ic1,
                                           isr1,
                                           osr,
                                           threshold);
}

Here is the call graph for this function:

EXTENSION_INLINE double ST_Distance_MultiPoint_LineString	(	int8_t *	mp,
		int64_t	mpsize,
		int8_t *	l,
		int64_t	lsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2528 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), and distance_point_linestring().

Referenced by ST_Distance_LineString_MultiPoint().

                                                           {
  auto mp_num_coords = mpsize / compression_unit_size(ic1);
  auto* p = mp;
  auto const psize = compression_unit_size(ic1) * 2;
  auto dist = distance_point_linestring(
      p, psize, l, lsize, ic1, isr1, ic2, isr2, osr, true, threshold);
  p += psize;
  for (int32_t i = 2; i < mp_num_coords; i += 2) {
    auto ldist = distance_point_linestring(
        p, psize, l, lsize, ic1, isr1, ic2, isr2, osr, true, threshold);
    p += psize;
    if (dist > ldist) {
      dist = ldist;
    }
    if (dist <= threshold) {
      return dist;
    }
  }
  return dist;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE double ST_Distance_MultiPoint_MultiLineString	(	int8_t *	mp,
		int64_t	mpsize,
		int8_t *	mls,
		int64_t	mls_size,
		int32_t *	mls_ls_sizes,
		int64_t	mls_ls_num,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2559 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), and distance_point_multilinestring().

Referenced by ST_Distance_MultiLineString_MultiPoint().

                                                                {
  auto mp_num_coords = mpsize / compression_unit_size(ic1);
  auto* p = mp;
  auto const psize = compression_unit_size(ic1) * 2;
  auto dist = distance_point_multilinestring(p,
                                             psize,
                                             mls,
                                             mls_size,
                                             mls_ls_sizes,
                                             mls_ls_num,
                                             ic1,
                                             isr1,
                                             ic2,
                                             isr2,
                                             osr,
                                             threshold);
  p += psize;
  for (int32_t i = 2; i < mp_num_coords; i += 2) {
    auto ldist = distance_point_multilinestring(p,
                                                psize,
                                                mls,
                                                mls_size,
                                                mls_ls_sizes,
                                                mls_ls_num,
                                                ic1,
                                                isr1,
                                                ic2,
                                                isr2,
                                                osr,
                                                threshold);
    p += psize;
    if (dist > ldist) {
      dist = ldist;
    }
    if (dist <= threshold) {
      return dist;
    }
  }
  return dist;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_MultiPoint_MultiPoint	(	int8_t *	mp1,
		int64_t	mp1size,
		int8_t *	mp2,
		int64_t	mp2size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2468 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), distance_point_point(), and get_point().

                                                           {
  auto mp1_num_coords = mp1size / compression_unit_size(ic1);
  auto mp2_num_coords = mp2size / compression_unit_size(ic2);
  double dist = -1.0;
  for (int32_t i = 0; i < mp1_num_coords; i += 2) {
    Point2D mp1p = get_point(mp1, 0, ic1, isr1, osr);
    for (int32_t j = 0; j < mp2_num_coords; j += 2) {
      Point2D mp2p = get_point(mp2, 0, ic2, isr2, osr);
      double ldist = distance_point_point(mp1p.x, mp1p.y, mp2p.x, mp2p.y);
      if (dist > ldist || dist < 0.0) {
        dist = ldist;
      }
      if (dist <= threshold) {
        return dist;
      }
    }
  }
  return dist;
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Distance_MultiPoint_MultiPoint_Squared	(	int8_t *	mp1,
		int64_t	mp1size,
		int8_t *	mp2,
		int64_t	mp2size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2498 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), distance_point_point_squared(), and get_point().

                                                                   {
  auto mp1_num_coords = mp1size / compression_unit_size(ic1);
  auto mp2_num_coords = mp2size / compression_unit_size(ic2);
  double dist = -1.0;
  for (int32_t i = 0; i < mp1_num_coords; i += 2) {
    Point2D mp1p = get_point(mp1, i, ic1, isr1, osr);
    for (int32_t j = 0; j < mp2_num_coords; j += 2) {
      Point2D mp2p = get_point(mp2, j, ic2, isr2, osr);
      double ldist = distance_point_point_squared(mp1p.x, mp1p.y, mp2p.x, mp2p.y);
      if (dist > ldist || dist < 0.0) {
        dist = ldist;
      }
      if (dist <= threshold) {
        return dist;
      }
    }
  }
  return dist;
}

Here is the call graph for this function:

EXTENSION_INLINE double ST_Distance_MultiPoint_MultiPolygon	(	int8_t *	mp,
		int64_t	mpsize,
		int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2321 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), and distance_point_multipolygon().

Referenced by ST_Distance_MultiPolygon_MultiPoint().

                                                             {
  auto mp_num_coords = mpsize / compression_unit_size(ic1);
  auto* p = mp;
  auto const psize = compression_unit_size(ic1) * 2;
  auto dist = distance_point_multipolygon(p,
                                          psize,
                                          mpoly_coords,
                                          mpoly_coords_size,
                                          mpoly_ring_sizes,
                                          mpoly_num_rings,
                                          mpoly_poly_sizes,
                                          mpoly_num_polys,
                                          ic1,
                                          isr1,
                                          ic2,
                                          isr2,
                                          osr,
                                          threshold);
  p += psize;
  for (int32_t i = 2; i < mp_num_coords; i += 2) {
    auto ldist = distance_point_multipolygon(p,
                                             psize,
                                             mpoly_coords,
                                             mpoly_coords_size,
                                             mpoly_ring_sizes,
                                             mpoly_num_rings,
                                             mpoly_poly_sizes,
                                             mpoly_num_polys,
                                             ic1,
                                             isr1,
                                             ic2,
                                             isr2,
                                             osr,
                                             threshold);
    p += psize;
    if (dist > ldist) {
      dist = ldist;
    }
    if (dist <= threshold) {
      return dist;
    }
  }
  return dist;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE double ST_Distance_MultiPoint_Point	(	int8_t *	mp,
		int64_t	mpsize,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2438 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_Point_MultiPoint().

                                                      {
  return ST_Distance_Point_MultiPoint(
      p, psize, mp, mpsize, ic2, isr2, ic1, isr1, osr, threshold);
}

Here is the call graph for this function:

EXTENSION_INLINE double ST_Distance_MultiPoint_Point_Squared	(	int8_t *	mp,
		int64_t	mpsize,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2453 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_Point_MultiPoint_Squared().

                                                              {
  return ST_Distance_Point_MultiPoint_Squared(
      p, psize, mp, mpsize, ic2, isr2, ic1, isr1, osr, threshold);
}

Here is the call graph for this function:

EXTENSION_INLINE double ST_Distance_MultiPoint_Polygon	(	int8_t *	mp,
		int64_t	mpsize,
		int8_t *	poly,
		int64_t	polysize,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2268 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), and distance_point_polygon().

Referenced by ST_Distance_Polygon_MultiPoint().

                                                        {
  auto mp_num_coords = mpsize / compression_unit_size(ic1);
  auto* p = mp;
  auto const psize = compression_unit_size(ic1) * 2;
  auto dist = distance_point_polygon(p,
                                     psize,
                                     poly,
                                     polysize,
                                     poly_ring_sizes,
                                     poly_num_rings,
                                     ic1,
                                     isr1,
                                     ic2,
                                     isr2,
                                     osr,
                                     threshold);
  p += psize;
  for (int32_t i = 2; i < mp_num_coords; i += 2) {
    auto ldist = distance_point_polygon(p,
                                        psize,
                                        poly,
                                        polysize,
                                        poly_ring_sizes,
                                        poly_num_rings,
                                        ic1,
                                        isr1,
                                        ic2,
                                        isr2,
                                        osr,
                                        threshold);
    p += psize;
    if (dist > ldist) {
      dist = ldist;
    }
    if (dist <= threshold) {
      return dist;
    }
  }
  return dist;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE double ST_Distance_MultiPolygon_LineString	(	int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		int8_t *	l,
		int64_t	lsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 3111 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_LineString_MultiPolygon().

                                                             {
  return ST_Distance_LineString_MultiPolygon(l,
                                             lsize,
                                             mpoly_coords,
                                             mpoly_coords_size,
                                             mpoly_ring_sizes,
                                             mpoly_num_rings,
                                             mpoly_poly_sizes,
                                             mpoly_num_polys,
                                             ic2,
                                             isr2,
                                             ic1,
                                             isr1,
                                             osr,
                                             threshold);
}

Here is the call graph for this function:

EXTENSION_INLINE double ST_Distance_MultiPolygon_MultiPoint	(	int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		int8_t *	mp,
		int64_t	mpsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2407 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_MultiPoint_MultiPolygon().

                                                             {
  return ST_Distance_MultiPoint_MultiPolygon(mp,
                                             mpsize,
                                             mpoly_coords,
                                             mpoly_coords_size,
                                             mpoly_ring_sizes,
                                             mpoly_num_rings,
                                             mpoly_poly_sizes,
                                             mpoly_num_polys,
                                             ic2,
                                             isr2,
                                             ic1,
                                             isr1,
                                             osr,
                                             threshold);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Distance_MultiPolygon_MultiPolygon	(	int8_t *	mpoly1_coords,
		int64_t	mpoly1_coords_size,
		int32_t *	mpoly1_ring_sizes,
		int64_t	mpoly1_num_rings,
		int32_t *	mpoly1_poly_sizes,
		int64_t	mpoly1_num_polys,
		int8_t *	mpoly2_coords,
		int64_t	mpoly2_coords_size,
		int32_t *	mpoly2_ring_sizes,
		int64_t	mpoly2_num_rings,
		int32_t *	mpoly2_poly_sizes,
		int64_t	mpoly2_num_polys,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 3177 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), ST_Distance_Polygon_MultiPolygon(), and tol_zero().

Referenced by ST_DWithin_MultiPolygon_MultiPolygon(), and ST_Intersects_MultiPolygon_MultiPolygon().

                                                               {
  double min_distance = 0.0;

  // Set specific poly pointers as we move through mpoly1's coords/ringsizes/polyrings
  // arrays.
  auto next_poly_coords = mpoly1_coords;
  auto next_poly_ring_sizes = mpoly1_ring_sizes;

  for (auto poly = 0; poly < mpoly1_num_polys; poly++) {
    auto poly_coords = next_poly_coords;
    auto poly_ring_sizes = next_poly_ring_sizes;
    auto poly_num_rings = mpoly1_poly_sizes[poly];
    // Count number of coords in all of poly's rings, advance ring size pointer.
    int32_t poly_num_coords = 0;
    for (auto ring = 0; ring < poly_num_rings; ring++) {
      poly_num_coords += 2 * *next_poly_ring_sizes++;
    }
    auto poly_coords_size = poly_num_coords * compression_unit_size(ic1);
    next_poly_coords += poly_coords_size;
    double distance = ST_Distance_Polygon_MultiPolygon(poly_coords,
                                                       poly_coords_size,
                                                       poly_ring_sizes,
                                                       poly_num_rings,
                                                       mpoly2_coords,
                                                       mpoly2_coords_size,
                                                       mpoly2_ring_sizes,
                                                       mpoly2_num_rings,
                                                       mpoly2_poly_sizes,
                                                       mpoly2_num_polys,
                                                       ic1,
                                                       isr1,
                                                       ic2,
                                                       isr2,
                                                       osr,
                                                       threshold);
    if (poly == 0 || min_distance > distance) {
      min_distance = distance;
      if (tol_zero(min_distance)) {
        min_distance = 0.0;
        break;
      }
      if (min_distance <= threshold) {
        break;
      }
    }
  }

  return min_distance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE double ST_Distance_MultiPolygon_Point	(	int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 3080 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_Point_MultiPolygon().

                                                        {
  return ST_Distance_Point_MultiPolygon(p,
                                        psize,
                                        mpoly_coords,
                                        mpoly_coords_size,
                                        mpoly_ring_sizes,
                                        mpoly_num_rings,
                                        mpoly_poly_sizes,
                                        mpoly_num_polys,
                                        ic2,
                                        isr2,
                                        ic1,
                                        isr1,
                                        osr,
                                        threshold);
}

Here is the call graph for this function:

EXTENSION_INLINE double ST_Distance_MultiPolygon_Polygon	(	int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		int8_t *	poly1_coords,
		int64_t	poly1_coords_size,
		int32_t *	poly1_ring_sizes,
		int64_t	poly1_num_rings,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 3142 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_Polygon_MultiPolygon().

                                                          {
  return ST_Distance_Polygon_MultiPolygon(poly1_coords,
                                          poly1_coords_size,
                                          poly1_ring_sizes,
                                          poly1_num_rings,
                                          mpoly_coords,
                                          mpoly_coords_size,
                                          mpoly_ring_sizes,
                                          mpoly_num_rings,
                                          mpoly_poly_sizes,
                                          mpoly_num_polys,
                                          ic2,
                                          isr2,
                                          ic1,
                                          isr1,
                                          osr,
                                          threshold);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Distance_Point_ClosedLineString	(	int8_t *	p,
		int64_t	psize,
		int8_t *	l,
		int64_t	lsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 1989 of file ExtensionFunctionsGeo.hpp.

References distance_point_linestring().

Referenced by distance_point_polygon().

                                                            {
  return distance_point_linestring(
      p, psize, l, lsize, ic1, isr1, ic2, isr2, osr, true, threshold);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_Point_LineString	(	int8_t *	p,
		int64_t	psize,
		int8_t *	l,
		int64_t	lsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2004 of file ExtensionFunctionsGeo.hpp.

References distance_point_linestring().

Referenced by ST_Contains_LineString_Point(), ST_Distance_LineString_Point(), ST_DWithin_Point_LineString(), and ST_Intersects_Point_LineString().

                                                      {
  return distance_point_linestring(
      p, psize, l, lsize, ic1, isr1, ic2, isr2, osr, false, threshold);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_Point_LineString_Geodesic	(	int8_t *	p,
		int64_t	psize,
		int8_t *	l,
		int64_t	lsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 1860 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), distance_in_meters(), and get_point().

Referenced by ST_Distance_LineString_Point_Geodesic().

                                                          {
  // Currently only statically indexed LineString is supported
  Point2D const pt = get_point(p, 0, ic1, 4326, 4326);
  const auto lpoints = lsize / (2 * compression_unit_size(ic2));
  Point2D const pl = get_point(l, 2 * (lpoints - 1), ic2, 4326, 4326);
  return distance_in_meters(pt.x, pt.y, pl.x, pl.y);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_Point_MultiLineString	(	int8_t *	p,
		int64_t	psize,
		int8_t *	mls,
		int64_t	mls_size,
		int32_t *	mls_ls_sizes,
		int64_t	mls_ls_num,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2070 of file ExtensionFunctionsGeo.hpp.

References distance_point_multilinestring().

Referenced by ST_Distance_MultiLineString_Point().

                                                           {
  return distance_point_multilinestring(p,
                                        psize,
                                        mls,
                                        mls_size,
                                        mls_ls_sizes,
                                        mls_ls_num,
                                        ic1,
                                        isr1,
                                        ic2,
                                        isr2,
                                        osr,
                                        threshold);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_Point_MultiPoint	(	int8_t *	p,
		int64_t	psize,
		int8_t *	mp,
		int64_t	mpsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 1892 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), distance_point_point(), and get_point().

Referenced by ST_Distance_MultiPoint_Point().

                                                      {
  Point2D pt = get_point(p, 0, ic1, isr1, osr);
  auto mp_num_coords = mpsize / compression_unit_size(ic2);
  Point2D mpp = get_point(mp, 0, ic2, isr2, osr);
  double dist = distance_point_point(pt.x, pt.y, mpp.x, mpp.y);
  for (int32_t i = 2; i < mp_num_coords; i += 2) {
    mpp = get_point(mp, i, ic2, isr2, osr);
    double ldist = distance_point_point(pt.x, pt.y, mpp.x, mpp.y);
    if (dist > ldist) {
      dist = ldist;
    }
    if (dist <= threshold) {
      return dist;
    }
  }
  return dist;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_Point_MultiPoint_Squared	(	int8_t *	p,
		int64_t	psize,
		int8_t *	mp,
		int64_t	mpsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 1920 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), distance_point_point_squared(), and get_point().

Referenced by ST_Distance_MultiPoint_Point_Squared().

                                                              {
  Point2D pt = get_point(p, 0, ic1, isr1, osr);
  auto mp_num_coords = mpsize / compression_unit_size(ic2);
  Point2D mpp = get_point(mp, 0, ic2, isr2, osr);
  double dist = distance_point_point_squared(pt.x, pt.y, mpp.x, mpp.y);
  for (int32_t i = 2; i < mp_num_coords; i += 2) {
    mpp = get_point(mp, i, ic2, isr2, osr);
    double ldist = distance_point_point_squared(pt.x, pt.y, mpp.x, mpp.y);
    if (dist > ldist) {
      dist = ldist;
    }
    if (dist <= threshold) {
      return dist;
    }
  }
  return dist;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_Point_MultiPolygon	(	int8_t *	p,
		int64_t	psize,
		int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2237 of file ExtensionFunctionsGeo.hpp.

References distance_point_multipolygon().

Referenced by ST_Distance_MultiPolygon_Point(), and ST_DWithin_Point_MultiPolygon().

                                                        {
  return distance_point_multipolygon(p,
                                     psize,
                                     mpoly_coords,
                                     mpoly_coords_size,
                                     mpoly_ring_sizes,
                                     mpoly_num_rings,
                                     mpoly_poly_sizes,
                                     mpoly_num_polys,
                                     ic1,
                                     isr1,
                                     ic2,
                                     isr2,
                                     osr,
                                     threshold);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_Point_Point	(	int8_t *	p1,
		int64_t	p1size,
		int8_t *	p2,
		int64_t	p2size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 1815 of file ExtensionFunctionsGeo.hpp.

References distance_point_point(), and get_point().

Referenced by ST_Intersects_Point_Point().

                                            {
  Point2D const pt1 = get_point(p1, 0, ic1, isr1, osr);
  Point2D const pt2 = get_point(p2, 0, ic2, isr2, osr);
  return distance_point_point(pt1.x, pt1.y, pt2.x, pt2.y);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_Point_Point_Geodesic	(	int8_t *	p1,
		int64_t	p1size,
		int8_t *	p2,
		int64_t	p2size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 1845 of file ExtensionFunctionsGeo.hpp.

References distance_in_meters(), and get_point().

Referenced by ST_DWithin_Point_Point_Geodesic().

                                                     {
  Point2D const pt1 = get_point(p1, 0, ic1, 4326, 4326);
  Point2D const pt2 = get_point(p2, 0, ic2, 4326, 4326);
  return distance_in_meters(pt1.x, pt1.y, pt2.x, pt2.y);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_Point_Point_Squared	(	int8_t *	p1,
		int64_t	p1size,
		int8_t *	p2,
		int64_t	p2size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 1830 of file ExtensionFunctionsGeo.hpp.

References distance_point_point_squared(), and get_point().

Referenced by ST_DWithin_Point_Point().

                                                    {
  Point2D const pt1 = get_point(p1, 0, ic1, isr1, osr);
  Point2D const pt2 = get_point(p2, 0, ic2, isr2, osr);
  return distance_point_point_squared(pt1.x, pt1.y, pt2.x, pt2.y);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_Point_Polygon	(	int8_t *	p,
		int64_t	psize,
		int8_t *	poly,
		int64_t	polysize,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2210 of file ExtensionFunctionsGeo.hpp.

References distance_point_polygon().

Referenced by ST_Distance_LineString_Polygon(), ST_Distance_Polygon_Point(), and ST_DWithin_Point_Polygon().

                                                   {
  return distance_point_polygon(p,
                                psize,
                                poly,
                                polysize,
                                poly_ring_sizes,
                                poly_num_rings,
                                ic1,
                                isr1,
                                ic2,
                                isr2,
                                osr,
                                threshold);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE double ST_Distance_Polygon_LineString	(	int8_t *	poly_coords,
		int64_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		int8_t *	l,
		int64_t	lsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2905 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_LineString_Polygon().

                                                        {
  return ST_Distance_LineString_Polygon(l,
                                        lsize,
                                        poly_coords,
                                        poly_coords_size,
                                        poly_ring_sizes,
                                        poly_num_rings,
                                        ic2,
                                        isr2,
                                        ic1,
                                        isr2,
                                        osr,
                                        threshold);
}

Here is the call graph for this function:

EXTENSION_INLINE double ST_Distance_Polygon_MultiPoint	(	int8_t *	poly,
		int64_t	polysize,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		int8_t *	mp,
		int64_t	mpsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2380 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_MultiPoint_Polygon().

                                                        {
  return ST_Distance_MultiPoint_Polygon(mp,
                                        mpsize,
                                        poly,
                                        polysize,
                                        poly_ring_sizes,
                                        poly_num_rings,
                                        ic2,
                                        isr2,
                                        ic1,
                                        isr1,
                                        osr,
                                        threshold);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Distance_Polygon_MultiPolygon	(	int8_t *	poly1_coords,
		int64_t	poly1_coords_size,
		int32_t *	poly1_ring_sizes,
		int64_t	poly1_num_rings,
		int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 3017 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), ST_Distance_Polygon_Polygon(), and tol_zero().

Referenced by ST_Distance_MultiPolygon_MultiPolygon(), ST_Distance_MultiPolygon_Polygon(), ST_DWithin_Polygon_MultiPolygon(), and ST_Intersects_Polygon_MultiPolygon().

                                                          {
  double min_distance = 0.0;

  // Set specific poly pointers as we move through the coords/ringsizes/polyrings arrays.
  auto next_poly_coords = mpoly_coords;
  auto next_poly_ring_sizes = mpoly_ring_sizes;

  for (auto poly = 0; poly < mpoly_num_polys; poly++) {
    auto poly_coords = next_poly_coords;
    auto poly_ring_sizes = next_poly_ring_sizes;
    auto poly_num_rings = mpoly_poly_sizes[poly];
    // Count number of coords in all of poly's rings, advance ring size pointer.
    int32_t poly_num_coords = 0;
    for (auto ring = 0; ring < poly_num_rings; ring++) {
      poly_num_coords += 2 * *next_poly_ring_sizes++;
    }
    auto poly_coords_size = poly_num_coords * compression_unit_size(ic2);
    next_poly_coords += poly_coords_size;
    double distance = ST_Distance_Polygon_Polygon(poly1_coords,
                                                  poly1_coords_size,
                                                  poly1_ring_sizes,
                                                  poly1_num_rings,
                                                  poly_coords,
                                                  poly_coords_size,
                                                  poly_ring_sizes,
                                                  poly_num_rings,
                                                  ic1,
                                                  isr1,
                                                  ic2,
                                                  isr2,
                                                  osr,
                                                  threshold);
    if (poly == 0 || min_distance > distance) {
      min_distance = distance;
      if (tol_zero(min_distance)) {
        min_distance = 0.0;
        break;
      }
      if (min_distance <= threshold) {
        break;
      }
    }
  }

  return min_distance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE double ST_Distance_Polygon_Point	(	int8_t *	poly_coords,
		int64_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2878 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_Point_Polygon().

Referenced by ST_Distance_Polygon_Polygon().

                                                   {
  return ST_Distance_Point_Polygon(p,
                                   psize,
                                   poly_coords,
                                   poly_coords_size,
                                   poly_ring_sizes,
                                   poly_num_rings,
                                   ic2,
                                   isr2,
                                   ic1,
                                   isr1,
                                   osr,
                                   threshold);
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_Distance_Polygon_Polygon	(	int8_t *	poly1_coords,
		int64_t	poly1_coords_size,
		int32_t *	poly1_ring_sizes,
		int64_t	poly1_num_rings,
		int8_t *	poly2_coords,
		int64_t	poly2_coords_size,
		int32_t *	poly2_ring_sizes,
		int64_t	poly2_num_rings,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	threshold
	)

Definition at line 2932 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), distance_ring_ring(), ST_Distance_Polygon_Point(), and tol_zero().

Referenced by ST_Distance_Polygon_MultiPolygon(), ST_DWithin_Polygon_Polygon(), and ST_Intersects_Polygon_Polygon().

                                                     {
  // Check if poly1 contains the first point of poly2's shape, i.e. the external ring
  auto poly2_first_point_coords = poly2_coords;
  auto poly2_first_point_coords_size = compression_unit_size(ic2) * 2;
  auto min_distance = ST_Distance_Polygon_Point(poly1_coords,
                                                poly1_coords_size,
                                                poly1_ring_sizes,
                                                poly1_num_rings,
                                                poly2_first_point_coords,
                                                poly2_first_point_coords_size,
                                                ic1,
                                                isr1,
                                                ic2,
                                                isr2,
                                                osr,
                                                threshold);
  if (tol_zero(min_distance)) {
    // Polygons overlap
    return 0.0;
  }
  if (min_distance <= threshold) {
    return min_distance;
  }

  // Poly2's first point is either outside poly1's external ring or inside one of the
  // internal rings. Measure the smallest distance between a poly1 ring (external or
  // internal) and a poly2 ring (external or internal). If poly2 is completely outside
  // poly1, then the min distance would be between poly1's and poly2's external rings. If
  // poly2 is completely inside one of poly1 internal rings then the min distance would be
  // between that poly1 internal ring and poly2's external ring. If poly1 is completely
  // inside one of poly2 internal rings, min distance is between that internal ring and
  // poly1's external ring. In each case other rings don't get in the way. Any ring
  // intersection means zero distance - short-circuit and return.

  auto poly1_ring_coords = poly1_coords;
  for (auto r1 = 0; r1 < poly1_num_rings; r1++) {
    int64_t poly1_ring_num_coords = poly1_ring_sizes[r1] * 2;

    auto poly2_ring_coords = poly2_coords;
    for (auto r2 = 0; r2 < poly2_num_rings; r2++) {
      int64_t poly2_ring_num_coords = poly2_ring_sizes[r2] * 2;

      auto distance = distance_ring_ring(poly1_ring_coords,
                                         poly1_ring_num_coords,
                                         poly2_ring_coords,
                                         poly2_ring_num_coords,
                                         ic1,
                                         isr1,
                                         ic2,
                                         isr2,
                                         osr,
                                         threshold);
      if (min_distance > distance) {
        min_distance = distance;
        if (tol_zero(min_distance)) {
          return 0.0;
        }
        if (min_distance <= threshold) {
          return min_distance;
        }
      }

      poly2_ring_coords += poly2_ring_num_coords * compression_unit_size(ic2);
    }

    poly1_ring_coords += poly1_ring_num_coords * compression_unit_size(ic1);
  }

  return min_distance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE bool ST_DWithin_LineString_LineString	(	int8_t *	l1,
		int64_t	l1size,
		double *	l1bounds,
		int64_t	l1bounds_size,
		int8_t *	l2,
		int64_t	l2size,
		double *	l2bounds,
		int64_t	l2bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	distance_within
	)

Definition at line 3421 of file ExtensionFunctionsGeo.hpp.

References box_dwithin_box(), and ST_Distance_LineString_LineString().

                                                              {
  if (l1bounds && l2bounds) {
    // Bounding boxes need to be transformed to output SR before proximity check
    if (!box_dwithin_box(l1bounds,
                         l1bounds_size,
                         isr1,
                         l2bounds,
                         l2bounds_size,
                         isr2,
                         osr,
                         distance_within)) {
      return false;
    }
  }

  // May need to adjust the threshold by TOLERANCE_DEFAULT
  const double threshold = distance_within;
  return ST_Distance_LineString_LineString(
             l1, l1size, l2, l2size, ic1, isr1, ic2, isr2, osr, threshold) <=
         distance_within;
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_DWithin_LineString_MultiPolygon	(	int8_t *	l1,
		int64_t	l1size,
		double *	l1bounds,
		int64_t	l1bounds_size,
		int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	distance_within
	)

Definition at line 3522 of file ExtensionFunctionsGeo.hpp.

References box_dwithin_box(), compression_unit_size(), ST_Distance_LineString_MultiPolygon(), and trim_linestring_to_buffered_box().

                                                                {
  if (l1bounds && mpoly_bounds) {
    // Bounding boxes need to be transformed to output SR before proximity check
    if (!box_dwithin_box(l1bounds,
                         l1bounds_size,
                         isr1,
                         mpoly_bounds,
                         mpoly_bounds_size,
                         isr2,
                         osr,
                         distance_within)) {
      return false;
    }
  }

  // First, assume the entire linestring is relevant.
  CoordData l1_relevant_section{l1, l1size};
  // Mitigation for very long linestrings (5+ segments: think highways, powerlines)
  if (mpoly_bounds && l1size > 12 * compression_unit_size(ic1)) {
    // Before diving into distance calculations, try to trim the linestring just to
    // segments whose bounding boxes overlap the threshold-buffered mpoly bounding box.
    l1_relevant_section = trim_linestring_to_buffered_box(l1,
                                                          l1size,
                                                          ic1,
                                                          isr1,
                                                          mpoly_bounds,
                                                          mpoly_bounds_size,
                                                          isr2,
                                                          osr,
                                                          distance_within);
    if (!l1_relevant_section.ptr) {
      // The big short-circuit: if not a single segment's bounding box overlaped
      // with buffered mpoly bounding box, then it means that the entire linestring
      // doesn't come close enough to the multipolygon to be within distance.
      // The global bounding boxes for the input geometries may overlap but the linestring
      // never gets close enough to the mpoly to require an actual distance calculation.
      return false;
    }
  }

  // Run distance calculation but only on the linestring section that is relevant.
  // May need to adjust the threshold by TOLERANCE_DEFAULT
  const double threshold = distance_within;
  return ST_Distance_LineString_MultiPolygon(l1_relevant_section.ptr,
                                             l1_relevant_section.size,
                                             mpoly_coords,
                                             mpoly_coords_size,
                                             mpoly_ring_sizes,
                                             mpoly_num_rings,
                                             mpoly_poly_sizes,
                                             mpoly_num_polys,
                                             ic1,
                                             isr1,
                                             ic2,
                                             isr2,
                                             osr,
                                             threshold) <= distance_within;
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_DWithin_LineString_Polygon	(	int8_t *	l1,
		int64_t	l1size,
		double *	l1bounds,
		int64_t	l1bounds_size,
		int8_t *	poly_coords,
		int64_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	distance_within
	)

Definition at line 3457 of file ExtensionFunctionsGeo.hpp.

References box_dwithin_box(), compression_unit_size(), ST_Distance_LineString_Polygon(), and trim_linestring_to_buffered_box().

                                                           {
  if (l1bounds && poly_bounds) {
    // Bounding boxes need to be transformed to output SR before proximity check
    if (!box_dwithin_box(l1bounds,
                         l1bounds_size,
                         isr1,
                         poly_bounds,
                         poly_bounds_size,
                         isr2,
                         osr,
                         distance_within)) {
      return false;
    }
  }

  // First, assume the entire linestring is relevant.
  CoordData l1_relevant_section{l1, l1size};
  // Mitigation for very long linestrings (5+ segments: think highways, powerlines)
  if (poly_bounds && l1size > 12 * compression_unit_size(ic1)) {
    // Before diving into distance calculations, try to trim the linestring just to
    // segments whose bounding boxes overlap the threshold-buffered poly bounding box.
    l1_relevant_section = trim_linestring_to_buffered_box(
        l1, l1size, ic1, isr1, poly_bounds, poly_bounds_size, isr2, osr, distance_within);
    if (!l1_relevant_section.ptr) {
      // The big short-circuit: if not a single segment's bounding box overlaped
      // with buffered poly bounding box, then it means that the entire linestring
      // doesn't come close enough to the polygon to be within distance.
      // The global bounding boxes for the input geometries may overlap but the linestring
      // never gets close enough to the poly to require an actual distance calculation.
      return false;
    }
  }

  // May need to adjust the threshold by TOLERANCE_DEFAULT
  const double threshold = distance_within;
  return ST_Distance_LineString_Polygon(l1_relevant_section.ptr,
                                        l1_relevant_section.size,
                                        poly_coords,
                                        poly_coords_size,
                                        poly_ring_sizes,
                                        poly_num_rings,
                                        ic1,
                                        isr1,
                                        ic2,
                                        isr2,
                                        osr,
                                        threshold) <= distance_within;
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_DWithin_MultiPolygon_MultiPolygon	(	int8_t *	mpoly1_coords,
		int64_t	mpoly1_coords_size,
		int32_t *	mpoly1_ring_sizes,
		int64_t	mpoly1_num_rings,
		int32_t *	mpoly1_poly_sizes,
		int64_t	mpoly1_num_polys,
		double *	mpoly1_bounds,
		int64_t	mpoly1_bounds_size,
		int8_t *	mpoly2_coords,
		int64_t	mpoly2_coords_size,
		int32_t *	mpoly2_ring_sizes,
		int64_t	mpoly2_num_rings,
		int32_t *	mpoly2_poly_sizes,
		int64_t	mpoly2_num_polys,
		double *	mpoly2_bounds,
		int64_t	mpoly2_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	distance_within
	)

Definition at line 3705 of file ExtensionFunctionsGeo.hpp.

References box_dwithin_box(), and ST_Distance_MultiPolygon_MultiPolygon().

                                                                  {
  if (mpoly1_bounds && mpoly2_bounds) {
    // Bounding boxes need to be transformed to output SR before proximity check
    if (!box_dwithin_box(mpoly1_bounds,
                         mpoly1_bounds_size,
                         isr1,
                         mpoly2_bounds,
                         mpoly2_bounds_size,
                         isr2,
                         osr,
                         distance_within)) {
      return false;
    }
  }

  // May need to adjust the threshold by TOLERANCE_DEFAULT
  const double threshold = distance_within;
  return ST_Distance_MultiPolygon_MultiPolygon(mpoly1_coords,
                                               mpoly1_coords_size,
                                               mpoly1_ring_sizes,
                                               mpoly1_num_rings,
                                               mpoly1_poly_sizes,
                                               mpoly1_num_polys,
                                               mpoly2_coords,
                                               mpoly2_coords_size,
                                               mpoly2_ring_sizes,
                                               mpoly2_num_rings,
                                               mpoly2_poly_sizes,
                                               mpoly2_num_polys,
                                               ic1,
                                               isr1,
                                               ic2,
                                               isr2,
                                               osr,
                                               threshold) <= distance_within;
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_DWithin_Point_LineString	(	int8_t *	p1,
		int64_t	p1size,
		int8_t *	l2,
		int64_t	l2size,
		double *	l2bounds,
		int64_t	l2bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	distance_within
	)

Definition at line 3299 of file ExtensionFunctionsGeo.hpp.

References point_dwithin_box(), and ST_Distance_Point_LineString().

                                                         {
  if (l2bounds) {
    // Bounding boxes need to be transformed to output SR before proximity check
    if (!point_dwithin_box(
            p1, p1size, ic1, isr1, l2bounds, l2bounds_size, isr2, osr, distance_within)) {
      return false;
    }
  }

  // May need to adjust the threshold by TOLERANCE_DEFAULT
  const double threshold = distance_within;
  return ST_Distance_Point_LineString(
             p1, p1size, l2, l2size, ic1, isr1, ic2, isr2, osr, threshold) <=
         distance_within;
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_DWithin_Point_MultiPolygon	(	int8_t *	p,
		int64_t	psize,
		int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	distance_within
	)

Definition at line 3372 of file ExtensionFunctionsGeo.hpp.

References point_dwithin_box(), and ST_Distance_Point_MultiPolygon().

                                                           {
  if (mpoly_bounds) {
    if (!point_dwithin_box(p,
                           psize,
                           ic1,
                           isr1,
                           mpoly_bounds,
                           mpoly_bounds_size,
                           isr2,
                           osr,
                           distance_within)) {
      return false;
    }
  }

  // May need to adjust the threshold by TOLERANCE_DEFAULT
  const double threshold = distance_within;
  return ST_Distance_Point_MultiPolygon(p,
                                        psize,
                                        mpoly_coords,
                                        mpoly_coords_size,
                                        mpoly_ring_sizes,
                                        mpoly_num_rings,
                                        mpoly_poly_sizes,
                                        mpoly_num_polys,
                                        ic1,
                                        isr1,
                                        ic2,
                                        isr2,
                                        osr,
                                        threshold) <= distance_within;
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_DWithin_Point_Point	(	int8_t *	p1,
		int64_t	p1size,
		int8_t *	p2,
		int64_t	p2size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	distance_within
	)

Definition at line 3267 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_Point_Point_Squared().

                                                    {
  return ST_Distance_Point_Point_Squared(
             p1, p1size, p2, p2size, ic1, isr1, ic2, isr2, osr) <=
         distance_within * distance_within;
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_DWithin_Point_Point_Geodesic	(	int8_t *	p1,
		int64_t	p1size,
		int8_t *	p2,
		int64_t	p2size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	distance_within
	)

Definition at line 3283 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_Point_Point_Geodesic().

                                                             {
  auto dist_meters =
      ST_Distance_Point_Point_Geodesic(p1, p1size, p2, p2size, ic1, isr1, ic2, isr2, osr);
  return (dist_meters <= distance_within);
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_DWithin_Point_Polygon	(	int8_t *	p,
		int64_t	psize,
		int8_t *	poly_coords,
		int64_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	distance_within
	)

Definition at line 3327 of file ExtensionFunctionsGeo.hpp.

References point_dwithin_box(), and ST_Distance_Point_Polygon().

                                                      {
  if (poly_bounds) {
    if (!point_dwithin_box(p,
                           psize,
                           ic1,
                           isr1,
                           poly_bounds,
                           poly_bounds_size,
                           isr2,
                           osr,
                           distance_within)) {
      return false;
    }
  }

  // May need to adjust the threshold by TOLERANCE_DEFAULT
  const double threshold = distance_within;
  return ST_Distance_Point_Polygon(p,
                                   psize,
                                   poly_coords,
                                   poly_coords_size,
                                   poly_ring_sizes,
                                   poly_num_rings,
                                   ic1,
                                   isr1,
                                   ic2,
                                   isr2,
                                   osr,
                                   threshold) <= distance_within;
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_DWithin_Polygon_MultiPolygon	(	int8_t *	poly_coords,
		int64_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	distance_within
	)

Definition at line 3650 of file ExtensionFunctionsGeo.hpp.

References box_dwithin_box(), and ST_Distance_Polygon_MultiPolygon().

                                                             {
  if (poly_bounds && mpoly_bounds) {
    // Bounding boxes need to be transformed to output SR before proximity check
    if (!box_dwithin_box(poly_bounds,
                         poly_bounds_size,
                         isr1,
                         mpoly_bounds,
                         mpoly_bounds_size,
                         isr2,
                         osr,
                         distance_within)) {
      return false;
    }
  }

  // May need to adjust the threshold by TOLERANCE_DEFAULT
  const double threshold = distance_within;
  return ST_Distance_Polygon_MultiPolygon(poly_coords,
                                          poly_coords_size,
                                          poly_ring_sizes,
                                          poly_num_rings,
                                          mpoly_coords,
                                          mpoly_coords_size,
                                          mpoly_ring_sizes,
                                          mpoly_num_rings,
                                          mpoly_poly_sizes,
                                          mpoly_num_polys,
                                          ic1,
                                          isr1,
                                          ic2,
                                          isr2,
                                          osr,
                                          threshold) <= distance_within;
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_DWithin_Polygon_Polygon	(	int8_t *	poly1_coords,
		int64_t	poly1_coords_size,
		int32_t *	poly1_ring_sizes,
		int64_t	poly1_num_rings,
		double *	poly1_bounds,
		int64_t	poly1_bounds_size,
		int8_t *	poly2_coords,
		int64_t	poly2_coords_size,
		int32_t *	poly2_ring_sizes,
		int64_t	poly2_num_rings,
		double *	poly2_bounds,
		int64_t	poly2_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr,
		double	distance_within
	)

Definition at line 3599 of file ExtensionFunctionsGeo.hpp.

References box_dwithin_box(), and ST_Distance_Polygon_Polygon().

                                                        {
  if (poly1_bounds && poly2_bounds) {
    // Bounding boxes need to be transformed to output SR before proximity check
    if (!box_dwithin_box(poly1_bounds,
                         poly1_bounds_size,
                         isr1,
                         poly2_bounds,
                         poly2_bounds_size,
                         isr2,
                         osr,
                         distance_within)) {
      return false;
    }
  }

  // May need to adjust the threshold by TOLERANCE_DEFAULT
  const double threshold = distance_within;
  return ST_Distance_Polygon_Polygon(poly1_coords,
                                     poly1_coords_size,
                                     poly1_ring_sizes,
                                     poly1_num_rings,
                                     poly2_coords,
                                     poly2_coords_size,
                                     poly2_ring_sizes,
                                     poly2_num_rings,
                                     ic1,
                                     isr1,
                                     ic2,
                                     isr2,
                                     osr,
                                     threshold) <= distance_within;
}

Here is the call graph for this function:

EXTENSION_NOINLINE bool ST_Intersects_LineString_Linestring	(	int8_t *	l1,
		int64_t	l1size,
		double *	l1bounds,
		int64_t	l1bounds_size,
		int8_t *	l2,
		int64_t	l2size,
		double *	l2bounds,
		int64_t	l2bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4507 of file ExtensionFunctionsGeo.hpp.

References box_overlaps_box(), ST_Distance_LineString_LineString(), and tol_zero().

                                                      {
  if (l1bounds && l2bounds) {
    if (!box_overlaps_box(l1bounds, l1bounds_size, l2bounds, l2bounds_size)) {
      return false;
    }
  }

  return tol_zero(ST_Distance_LineString_LineString(
      l1, l1size, l2, l2size, ic1, isr1, ic2, isr2, osr, 0.0));
}

Here is the call graph for this function:

EXTENSION_NOINLINE bool ST_Intersects_LineString_MultiPolygon	(	int8_t *	l,
		int64_t	lsize,
		double *	lbounds,
		int64_t	lbounds_size,
		int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4574 of file ExtensionFunctionsGeo.hpp.

References box_overlaps_box(), ST_Distance_LineString_MultiPolygon(), and tol_zero().

Referenced by ST_Intersects_MultiPolygon_LineString().

                                                        {
  if (lbounds && mpoly_bounds) {
    if (!box_overlaps_box(lbounds, lbounds_size, mpoly_bounds, mpoly_bounds_size)) {
      return false;
    }
  }

  // Check for spatial intersection.
  // One way to do that would be to start with linestring's first point, if it's inside
  // any of the polygons - it means intersection. Otherwise follow the linestring, segment
  // by segment, checking for intersections with polygon shapes/holes, bail as soon as we
  // cross into a polygon.

  // Or, alternatively, just measure the distance:
  return tol_zero(ST_Distance_LineString_MultiPolygon(l,
                                                      lsize,
                                                      mpoly_coords,
                                                      mpoly_coords_size,
                                                      mpoly_ring_sizes,
                                                      mpoly_num_rings,
                                                      mpoly_poly_sizes,
                                                      mpoly_num_polys,
                                                      ic1,
                                                      isr1,
                                                      ic2,
                                                      isr2,
                                                      osr,
                                                      0.0));
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE bool ST_Intersects_LineString_Point	(	int8_t *	l,
		int64_t	lsize,
		double *	lbounds,
		int64_t	lbounds_size,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4491 of file ExtensionFunctionsGeo.hpp.

References ST_Intersects_Point_LineString().

                                                 {
  return ST_Intersects_Point_LineString(
      p, psize, l, lsize, lbounds, lbounds_size, ic2, isr2, ic1, isr1, osr);
}

Here is the call graph for this function:

EXTENSION_NOINLINE bool ST_Intersects_LineString_Polygon	(	int8_t *	l,
		int64_t	lsize,
		double *	lbounds,
		int64_t	lbounds_size,
		int8_t *	poly,
		int64_t	polysize,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4531 of file ExtensionFunctionsGeo.hpp.

References box_overlaps_box(), ST_Distance_LineString_Polygon(), and tol_zero().

Referenced by ST_Intersects_Polygon_LineString().

                                                   {
  if (lbounds && poly_bounds) {
    if (!box_overlaps_box(lbounds, lbounds_size, poly_bounds, poly_bounds_size)) {
      return false;
    }
  }

  // Check for spatial intersection.
  // One way to do that would be to start with linestring's first point, if it's inside
  // the polygon - it means intersection. Otherwise follow the linestring, segment by
  // segment, checking for intersections with polygon rings, bail as soon as we cross into
  // the polygon.

  // Or, alternatively, just measure the distance:
  return tol_zero(ST_Distance_LineString_Polygon(l,
                                                 lsize,
                                                 poly,
                                                 polysize,
                                                 poly_ring_sizes,
                                                 poly_num_rings,
                                                 ic1,
                                                 isr1,
                                                 ic2,
                                                 isr2,
                                                 osr,
                                                 0.0));
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE bool ST_Intersects_MultiPolygon_LineString	(	int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int8_t *	l,
		int64_t	lsize,
		double *	lbounds,
		int64_t	lbounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4868 of file ExtensionFunctionsGeo.hpp.

References ST_Intersects_LineString_MultiPolygon().

                                                        {
  return ST_Intersects_LineString_MultiPolygon(l,
                                               lsize,
                                               lbounds,
                                               lbounds_size,
                                               mpoly_coords,
                                               mpoly_coords_size,
                                               mpoly_ring_sizes,
                                               mpoly_num_rings,
                                               mpoly_poly_sizes,
                                               mpoly_num_polys,
                                               mpoly_bounds,
                                               mpoly_bounds_size,
                                               ic2,
                                               isr2,
                                               ic1,
                                               isr1,
                                               osr);
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_Intersects_MultiPolygon_MultiPolygon	(	int8_t *	mpoly1_coords,
		int64_t	mpoly1_coords_size,
		int32_t *	mpoly1_ring_sizes,
		int64_t	mpoly1_num_rings,
		int32_t *	mpoly1_poly_sizes,
		int64_t	mpoly1_num_polys,
		double *	mpoly1_bounds,
		int64_t	mpoly1_bounds_size,
		int8_t *	mpoly2_coords,
		int64_t	mpoly2_coords_size,
		int32_t *	mpoly2_ring_sizes,
		int64_t	mpoly2_num_rings,
		int32_t *	mpoly2_poly_sizes,
		int64_t	mpoly2_num_polys,
		double *	mpoly2_bounds,
		int64_t	mpoly2_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4946 of file ExtensionFunctionsGeo.hpp.

References box_overlaps_box(), ST_Distance_MultiPolygon_MultiPolygon(), and tol_zero().

                                                          {
  if (mpoly1_bounds && mpoly2_bounds) {
    if (!box_overlaps_box(
            mpoly1_bounds, mpoly1_bounds_size, mpoly2_bounds, mpoly2_bounds_size)) {
      return false;
    }
  }

  return tol_zero(ST_Distance_MultiPolygon_MultiPolygon(mpoly1_coords,
                                                        mpoly1_coords_size,
                                                        mpoly1_ring_sizes,
                                                        mpoly1_num_rings,
                                                        mpoly1_poly_sizes,
                                                        mpoly1_num_polys,
                                                        mpoly2_coords,
                                                        mpoly2_coords_size,
                                                        mpoly2_ring_sizes,
                                                        mpoly2_num_rings,
                                                        mpoly2_poly_sizes,
                                                        mpoly2_num_polys,
                                                        ic1,
                                                        isr1,
                                                        ic2,
                                                        isr2,
                                                        osr,
                                                        0.0));
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_Intersects_MultiPolygon_Point	(	int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4800 of file ExtensionFunctionsGeo.hpp.

                                                   {
  return Contains_MultiPolygon_Point_Impl<double, EdgeBehavior::kIncludePointOnEdge>(
      mpoly_coords,
      mpoly_coords_size,
      mpoly_ring_sizes,
      mpoly_num_rings,
      mpoly_poly_sizes,
      mpoly_num_polys,
      mpoly_bounds,
      mpoly_bounds_size,
      p,
      psize,
      ic1,
      isr1,
      ic2,
      isr2,
      osr);
}

EXTENSION_INLINE bool ST_Intersects_MultiPolygon_Polygon	(	int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int8_t *	poly_coords,
		int64_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4905 of file ExtensionFunctionsGeo.hpp.

References ST_Intersects_Polygon_MultiPolygon().

                                                     {
  return ST_Intersects_Polygon_MultiPolygon(poly_coords,
                                            poly_coords_size,
                                            poly_ring_sizes,
                                            poly_num_rings,
                                            poly_bounds,
                                            poly_bounds_size,
                                            mpoly_coords,
                                            mpoly_coords_size,
                                            mpoly_ring_sizes,
                                            mpoly_num_rings,
                                            mpoly_poly_sizes,
                                            mpoly_num_polys,
                                            mpoly_bounds,
                                            mpoly_bounds_size,
                                            ic2,
                                            isr2,
                                            ic1,
                                            isr1,
                                            osr);
}

Here is the call graph for this function:

EXTENSION_NOINLINE bool ST_Intersects_Point_LineString	(	int8_t *	p,
		int64_t	psize,
		int8_t *	l,
		int64_t	lsize,
		double *	lbounds,
		int64_t	lbounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4405 of file ExtensionFunctionsGeo.hpp.

References box_contains_point(), get_point(), ST_Distance_Point_LineString(), and tol_zero().

Referenced by ST_Intersects_LineString_Point().

                                                 {
  if (lbounds) {
    Point2D const pt = get_point(p, 0, ic1, isr1, osr);
    if (!box_contains_point(lbounds, lbounds_size, pt.x, pt.y)) {
      return false;
    }
  }
  return tol_zero(
      ST_Distance_Point_LineString(p, psize, l, lsize, ic1, isr1, ic2, isr2, osr, 0.0));
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE bool ST_Intersects_Point_MultiPolygon	(	int8_t *	p,
		int64_t	psize,
		int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4457 of file ExtensionFunctionsGeo.hpp.

                                                   {
  return Contains_MultiPolygon_Point_Impl<double, EdgeBehavior::kIncludePointOnEdge>(
      mpoly_coords,
      mpoly_coords_size,
      mpoly_ring_sizes,
      mpoly_num_rings,
      mpoly_poly_sizes,
      mpoly_num_polys,
      mpoly_bounds,
      mpoly_bounds_size,
      p,
      psize,
      ic1,
      isr1,
      ic2,
      isr2,
      osr);
}

EXTENSION_INLINE bool ST_Intersects_Point_Point	(	int8_t *	p1,
		int64_t	p1size,
		int8_t *	p2,
		int64_t	p2size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4391 of file ExtensionFunctionsGeo.hpp.

References ST_Distance_Point_Point(), and tol_zero().

                                            {
  return tol_zero(
      ST_Distance_Point_Point(p1, p1size, p2, p2size, ic1, isr1, ic2, isr2, osr));
}

Here is the call graph for this function:

EXTENSION_INLINE bool ST_Intersects_Point_Polygon	(	int8_t *	p,
		int64_t	psize,
		int8_t *	poly,
		int64_t	polysize,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4427 of file ExtensionFunctionsGeo.hpp.

                                              {
  return Contains_Polygon_Point_Impl<double, EdgeBehavior::kIncludePointOnEdge>(
      poly,
      polysize,
      poly_ring_sizes,
      poly_num_rings,
      poly_bounds,
      poly_bounds_size,
      p,
      psize,
      ic2,
      isr2,
      ic1,
      isr1,
      osr);
}

EXTENSION_INLINE bool ST_Intersects_Polygon_LineString	(	int8_t *	poly,
		int64_t	polysize,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int8_t *	l,
		int64_t	lsize,
		double *	lbounds,
		int64_t	lbounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4681 of file ExtensionFunctionsGeo.hpp.

References ST_Intersects_LineString_Polygon().

                                                   {
  return ST_Intersects_LineString_Polygon(l,
                                          lsize,
                                          lbounds,
                                          lbounds_size,
                                          poly,
                                          polysize,
                                          poly_ring_sizes,
                                          poly_num_rings,
                                          poly_bounds,
                                          poly_bounds_size,
                                          ic2,
                                          isr2,
                                          ic1,
                                          isr1,
                                          osr);
}

Here is the call graph for this function:

EXTENSION_NOINLINE bool ST_Intersects_Polygon_MultiPolygon	(	int8_t *	poly_coords,
		int64_t	poly_coords_size,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int8_t *	mpoly_coords,
		int64_t	mpoly_coords_size,
		int32_t *	mpoly_ring_sizes,
		int64_t	mpoly_num_rings,
		int32_t *	mpoly_poly_sizes,
		int64_t	mpoly_num_polys,
		double *	mpoly_bounds,
		int64_t	mpoly_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4755 of file ExtensionFunctionsGeo.hpp.

References box_overlaps_box(), ST_Distance_Polygon_MultiPolygon(), and tol_zero().

Referenced by ST_Intersects_MultiPolygon_Polygon().

                                                     {
  if (poly_bounds && mpoly_bounds) {
    if (!box_overlaps_box(
            poly_bounds, poly_bounds_size, mpoly_bounds, mpoly_bounds_size)) {
      return false;
    }
  }

  return tol_zero(ST_Distance_Polygon_MultiPolygon(poly_coords,
                                                   poly_coords_size,
                                                   poly_ring_sizes,
                                                   poly_num_rings,
                                                   mpoly_coords,
                                                   mpoly_coords_size,
                                                   mpoly_ring_sizes,
                                                   mpoly_num_rings,
                                                   mpoly_poly_sizes,
                                                   mpoly_num_polys,
                                                   ic1,
                                                   isr1,
                                                   ic2,
                                                   isr2,
                                                   osr,
                                                   0.0));
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE bool ST_Intersects_Polygon_Point	(	int8_t *	poly,
		int64_t	polysize,
		int32_t *	poly_ring_sizes,
		int64_t	poly_num_rings,
		double *	poly_bounds,
		int64_t	poly_bounds_size,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4621 of file ExtensionFunctionsGeo.hpp.

                                              {
  return Contains_Polygon_Point_Impl<double, EdgeBehavior::kIncludePointOnEdge>(
      poly,
      polysize,
      poly_ring_sizes,
      poly_num_rings,
      poly_bounds,
      poly_bounds_size,
      p,
      psize,
      ic1,
      isr1,
      ic2,
      isr2,
      osr);
}

EXTENSION_NOINLINE bool ST_Intersects_Polygon_Polygon	(	int8_t *	poly1_coords,
		int64_t	poly1_coords_size,
		int32_t *	poly1_ring_sizes,
		int64_t	poly1_num_rings,
		double *	poly1_bounds,
		int64_t	poly1_bounds_size,
		int8_t *	poly2_coords,
		int64_t	poly2_coords_size,
		int32_t *	poly2_ring_sizes,
		int64_t	poly2_num_rings,
		double *	poly2_bounds,
		int64_t	poly2_bounds_size,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 4714 of file ExtensionFunctionsGeo.hpp.

References box_overlaps_box(), ST_Distance_Polygon_Polygon(), and tol_zero().

                                                {
  if (poly1_bounds && poly2_bounds) {
    if (!box_overlaps_box(
            poly1_bounds, poly1_bounds_size, poly2_bounds, poly2_bounds_size)) {
      return false;
    }
  }

  return tol_zero(ST_Distance_Polygon_Polygon(poly1_coords,
                                              poly1_coords_size,
                                              poly1_ring_sizes,
                                              poly1_num_rings,
                                              poly2_coords,
                                              poly2_coords_size,
                                              poly2_ring_sizes,
                                              poly2_num_rings,
                                              ic1,
                                              isr1,
                                              ic2,
                                              isr2,
                                              osr,
                                              0.0));
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Length_LineString	(	int8_t *	coords,
		int64_t	coords_sz,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1169 of file ExtensionFunctionsGeo.hpp.

References length_linestring().

                                         {
  return length_linestring(coords, coords_sz, ic, isr, osr, false, false);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Length_LineString_Geodesic	(	int8_t *	coords,
		int64_t	coords_sz,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1178 of file ExtensionFunctionsGeo.hpp.

References length_linestring().

                                                  {
  return length_linestring(coords, coords_sz, ic, isr, osr, true, false);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Length_MultiLineString	(	int8_t *	coords,
		int64_t	coords_sz,
		int8_t *	linestring_sizes_in,
		int64_t	linestring_sizes_sz,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1187 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), and length_linestring().

                                              {
  if (linestring_sizes_sz <= 0) {
    return 0.0;
  }
  double mls_length = 0.0;

  auto next_linestring_coords = coords;

  for (auto l = 0; l < linestring_sizes_sz; l++) {
    auto linestring_coords = next_linestring_coords;
    auto linestring_sizes = reinterpret_cast<int32_t*>(linestring_sizes_in);
    auto linestring_num_points = linestring_sizes[l];
    auto linestring_num_coords = 2 * linestring_num_points;
    auto linestring_coords_size = linestring_num_coords * compression_unit_size(ic);
    next_linestring_coords += linestring_coords_size;
    double ls_length = length_linestring(
        linestring_coords, linestring_coords_size, ic, isr, osr, false, false);
    mls_length += ls_length;
  }

  return mls_length;
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_MaxDistance_LineString_Point	(	int8_t *	l,
		int64_t	lsize,
		int8_t *	p,
		int64_t	psize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 3837 of file ExtensionFunctionsGeo.hpp.

References max_distance_point_linestring().

                                                    {
  return max_distance_point_linestring(
      p, psize, l, lsize, ic2, isr2, ic1, isr1, osr, false);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_MaxDistance_Point_LineString	(	int8_t *	p,
		int64_t	psize,
		int8_t *	l,
		int64_t	lsize,
		int32_t	ic1,
		int32_t	isr1,
		int32_t	ic2,
		int32_t	isr2,
		int32_t	osr
	)

Definition at line 3823 of file ExtensionFunctionsGeo.hpp.

References max_distance_point_linestring().

                                                    {
  return max_distance_point_linestring(
      p, psize, l, lsize, ic1, isr1, ic2, isr2, osr, false);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Perimeter_MultiPolygon	(	int8_t *	mpoly_coords,
		int32_t	mpoly_coords_size,
		int8_t *	mpoly_ring_sizes,
		int32_t	mpoly_num_rings,
		int8_t *	mpoly_poly_sizes,
		int32_t	mpoly_num_polys,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1302 of file ExtensionFunctionsGeo.hpp.

References perimeter_multipolygon().

                                              {
  return perimeter_multipolygon(mpoly_coords,
                                mpoly_coords_size,
                                mpoly_ring_sizes,
                                mpoly_num_rings,
                                mpoly_poly_sizes,
                                mpoly_num_polys,
                                ic,
                                isr,
                                osr,
                                false);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Perimeter_MultiPolygon_Geodesic	(	int8_t *	mpoly_coords,
		int32_t	mpoly_coords_size,
		int8_t *	mpoly_ring_sizes,
		int32_t	mpoly_num_rings,
		int8_t *	mpoly_poly_sizes,
		int32_t	mpoly_num_polys,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1324 of file ExtensionFunctionsGeo.hpp.

References perimeter_multipolygon().

                                                       {
  return perimeter_multipolygon(mpoly_coords,
                                mpoly_coords_size,
                                mpoly_ring_sizes,
                                mpoly_num_rings,
                                mpoly_poly_sizes,
                                mpoly_num_polys,
                                ic,
                                isr,
                                osr,
                                true);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Perimeter_Polygon	(	int8_t *	poly,
		int32_t	polysize,
		int8_t *	poly_ring_sizes,
		int32_t	poly_num_rings,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1221 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), and length_linestring().

                                         {
  if (poly_num_rings <= 0) {
    return 0.0;
  }

  auto exterior_ring_num_coords = poly_ring_sizes[0] * 2;
  auto exterior_ring_coords_size = exterior_ring_num_coords * compression_unit_size(ic);

  return length_linestring(poly, exterior_ring_coords_size, ic, isr, osr, false, true);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_Perimeter_Polygon_Geodesic	(	int8_t *	poly,
		int32_t	polysize,
		int8_t *	poly_ring_sizes_in,
		int32_t	poly_num_rings,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1239 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), and length_linestring().

                                                  {
  if (poly_num_rings <= 0) {
    return 0.0;
  }

  auto poly_ring_sizes = reinterpret_cast<int32_t*>(poly_ring_sizes_in);
  auto exterior_ring_num_coords = poly_ring_sizes[0] * 2;
  auto exterior_ring_coords_size = exterior_ring_num_coords * compression_unit_size(ic);

  return length_linestring(poly, exterior_ring_coords_size, ic, isr, osr, true, true);
}

Here is the call graph for this function:

EXTENSION_NOINLINE double ST_X_Point	(	int8_t *	p,
		int64_t	psize,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1055 of file ExtensionFunctionsGeo.hpp.

Referenced by to_Geospatial_GeoPoint().

                                                                                  {
  return coord<X>(p, 0, ic, isr, osr).x;
}

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_XMax	(	int8_t *	coords,
		int64_t	size,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1091 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size().

Referenced by com.mapd.calcite.parser.HeavyDBSqlOperatorTable::addUDF().

                                                                                   {
  auto num_coords = size / compression_unit_size(ic);
  double xmax = 0.0;
  for (int32_t i = 0; i < num_coords; i += 2) {
    double x = coord<X>(coords, i, ic, isr, osr).x;
    if (i == 0 || x > xmax) {
      xmax = x;
    }
  }
  return xmax;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE double ST_XMax_Bounds	(	double *	bounds,
		int64_t	size,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1127 of file ExtensionFunctionsGeo.hpp.

                                                                              {
  return transform_point<X>({bounds[2], bounds[3]}, isr, osr).x;
}

EXTENSION_NOINLINE double ST_XMin	(	int8_t *	coords,
		int64_t	size,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1065 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size().

Referenced by com.mapd.calcite.parser.HeavyDBSqlOperatorTable::addUDF().

                                                                                   {
  auto num_coords = size / compression_unit_size(ic);
  double xmin = 0.0;
  for (int32_t i = 0; i < num_coords; i += 2) {
    double x = coord<X>(coords, i, ic, isr, osr).x;
    if (i == 0 || x < xmin) {
      xmin = x;
    }
  }
  return xmin;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE double ST_XMin_Bounds	(	double *	bounds,
		int64_t	size,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1117 of file ExtensionFunctionsGeo.hpp.

                                                                              {
  return transform_point<X>({bounds[0], bounds[1]}, isr, osr).x;
}

EXTENSION_NOINLINE double ST_Y_Point	(	int8_t *	p,
		int64_t	psize,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1060 of file ExtensionFunctionsGeo.hpp.

Referenced by to_Geospatial_GeoPoint().

                                                                                  {
  return coord<Y>(p, 0, ic, isr, osr).y;
}

Here is the caller graph for this function:

EXTENSION_NOINLINE double ST_YMax	(	int8_t *	coords,
		int64_t	size,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1104 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size().

Referenced by com.mapd.calcite.parser.HeavyDBSqlOperatorTable::addUDF().

                                                                                   {
  auto num_coords = size / compression_unit_size(ic);
  double ymax = 0.0;
  for (int32_t i = 0; i < num_coords; i += 2) {
    double y = coord<Y>(coords, i, ic, isr, osr).y;
    if (i == 0 || y > ymax) {
      ymax = y;
    }
  }
  return ymax;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE double ST_YMax_Bounds	(	double *	bounds,
		int64_t	size,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1132 of file ExtensionFunctionsGeo.hpp.

                                                                              {
  return transform_point<Y>({bounds[2], bounds[3]}, isr, osr).y;
}

EXTENSION_NOINLINE double ST_YMin	(	int8_t *	coords,
		int64_t	size,
		int32_t	ic,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1078 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size().

Referenced by com.mapd.calcite.parser.HeavyDBSqlOperatorTable::addUDF().

                                                                                   {
  auto num_coords = size / compression_unit_size(ic);
  double ymin = 0.0;
  for (int32_t i = 0; i < num_coords; i += 2) {
    double y = coord<Y>(coords, i, ic, isr, osr).y;
    if (i == 0 || y < ymin) {
      ymin = y;
    }
  }
  return ymin;
}

Here is the call graph for this function:

Here is the caller graph for this function:

EXTENSION_INLINE double ST_YMin_Bounds	(	double *	bounds,
		int64_t	size,
		int32_t	isr,
		int32_t	osr
	)

Definition at line 1122 of file ExtensionFunctionsGeo.hpp.

                                                                              {
  return transform_point<Y>({bounds[0], bounds[1]}, isr, osr).y;
}

Geospatial::GeoLineString to_Geospatial_GeoLineString

(

GeoLineString &

l

)

Definition at line 5123 of file ExtensionFunctionsGeo.hpp.

References Geospatial::decompress_coords< double, int32_t >(), GeoLineStringStruct::getCompression(), GeoLineStringStruct::getInputSrid(), GeoLineStringStruct::getOutputSrid(), GeoLineStringStruct::getSize(), and GeoLineStringStruct::ptr.

Referenced by ST_AsBinary__GeoLineString__cpu_(), and ST_AsText__GeoLineString__cpu_().

                                                                      {
  int8_t* ptr = l.ptr;
  int32_t lsize = l.getSize();
  int32_t ic = l.getCompression();
  int32_t isr = l.getInputSrid();
  int32_t osr = l.getOutputSrid();

  std::vector<double> coords =
      *Geospatial::decompress_coords<double, int32_t>(ic, ptr, lsize);
  auto linestring = Geospatial::GeoLineString(coords);
  linestring.transform(isr, osr);
  return linestring;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Geospatial::GeoMultiLineString to_Geospatial_GeoMultiLineString

(

GeoMultiLineString &

ml

)

Definition at line 5137 of file ExtensionFunctionsGeo.hpp.

References Geospatial::decompress_coords< double, int32_t >(), GeoMultiLineStringStruct::getCompression(), GeoMultiLineStringStruct::getCoordsSize(), GeoMultiLineStringStruct::getInputSrid(), GeoMultiLineStringStruct::getLineStringSizes(), GeoMultiLineStringStruct::getNumLineStrings(), GeoMultiLineStringStruct::getOutputSrid(), and GeoMultiLineStringStruct::ptr.

Referenced by ST_AsBinary__GeoMultiLineString__cpu_(), and ST_AsText__GeoMultiLineString__cpu_().

                                                                                      {
  int8_t* ptr = ml.ptr;
  int32_t ml_coords_size = ml.getCoordsSize();
  int32_t* linestring_sizes = reinterpret_cast<int32_t*>(ml.getLineStringSizes());
  int32_t num_linestrings = ml.getNumLineStrings();
  int32_t ic = ml.getCompression();
  int32_t isr = ml.getInputSrid();
  int32_t osr = ml.getOutputSrid();

  std::vector<double> coords =
      *Geospatial::decompress_coords<double, int32_t>(ic, ptr, ml_coords_size);

  std::vector<int32_t> v_linestring_sizes(linestring_sizes,
                                          linestring_sizes + num_linestrings);

  auto multilinestring = Geospatial::GeoMultiLineString(coords, v_linestring_sizes);
  multilinestring.transform(isr, osr);
  return multilinestring;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Geospatial::GeoMultiPoint to_Geospatial_GeoMultiPoint

(

GeoMultiPoint &

mp

)

Definition at line 5109 of file ExtensionFunctionsGeo.hpp.

References Geospatial::decompress_coords< double, int32_t >(), GeoMultiPointStruct::getCompression(), GeoMultiPointStruct::getInputSrid(), GeoMultiPointStruct::getOutputSrid(), GeoMultiPointStruct::getSize(), and GeoMultiPointStruct::ptr.

Referenced by ST_AsBinary__GeoMultiPoint__cpu_(), and ST_AsText__GeoMultiPoint__cpu_().

                                                                       {
  int8_t* ptr = mp.ptr;
  int32_t lsize = mp.getSize();
  int32_t ic = mp.getCompression();
  int32_t isr = mp.getInputSrid();
  int32_t osr = mp.getOutputSrid();

  std::vector<double> coords =
      *Geospatial::decompress_coords<double, int32_t>(ic, ptr, lsize);
  auto multipoint = Geospatial::GeoMultiPoint(coords);
  multipoint.transform(isr, osr);
  return multipoint;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Geospatial::GeoMultiPolygon to_Geospatial_GeoMultiPolygon

(

GeoMultiPolygon &

mp

)

Definition at line 5175 of file ExtensionFunctionsGeo.hpp.

References Geospatial::decompress_coords< double, int32_t >().

Referenced by ST_AsBinary__GeoMultiPolygon__cpu_(), and ST_AsText__GeoMultiPolygon__cpu_().

                                                                             {
  int8_t* poly_coords = mp.ptr_coords;
  int32_t poly_coords_size = mp.getCoordsSize();
  int32_t* poly_ring_sizes = reinterpret_cast<int32_t*>(mp.getRingSizes());
  int32_t poly_num_rings = mp.getNumRings();
  int32_t* poly_sizes = reinterpret_cast<int32_t*>(mp.getPolygonSizes());
  int32_t num_polys = mp.getNumPolygons();
  int32_t ic = mp.getCompression();
  int32_t isr = mp.getInputSrid();
  int32_t osr = mp.getOutputSrid();

  std::vector<double> coords =
      *Geospatial::decompress_coords<double, int32_t>(ic, poly_coords, poly_coords_size);

  std::vector<int32_t> ring_sizes(poly_ring_sizes, poly_ring_sizes + poly_num_rings);
  std::vector<int32_t> poly_rings(poly_sizes, poly_sizes + num_polys);

  auto multipolygon = Geospatial::GeoMultiPolygon(coords, ring_sizes, poly_rings);
  multipolygon.transform(isr, osr);
  return multipolygon;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Geospatial::GeoPoint to_Geospatial_GeoPoint

(

GeoPoint &

p

)

Definition at line 5096 of file ExtensionFunctionsGeo.hpp.

References GeoPointStruct::getCompression(), GeoPointStruct::getInputSrid(), GeoPointStruct::getOutputSrid(), GeoPointStruct::getSize(), GeoPointStruct::ptr, ST_X_Point(), and ST_Y_Point().

Referenced by ST_AsBinary__GeoPoint__cpu_(), and ST_AsText__GeoPoint__cpu_().

                                                       {
  int8_t* ptr = p.ptr;
  int32_t lsize = p.getSize();
  int32_t ic = p.getCompression();
  int32_t isr = p.getInputSrid();
  int32_t osr = p.getOutputSrid();

  double x = ST_X_Point(ptr, lsize, ic, isr, osr);
  double y = ST_Y_Point(ptr, lsize, ic, isr, osr);
  auto point = Geospatial::GeoPoint(std::vector<double>{x, y});
  return point;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Geospatial::GeoPolygon to_Geospatial_GeoPolygon

(

GeoPolygon &

p

)

Definition at line 5157 of file ExtensionFunctionsGeo.hpp.

References Geospatial::decompress_coords< double, int32_t >().

Referenced by ST_AsBinary__GeoPolygon__cpu_(), and ST_AsText__GeoPolygon__cpu_().

                                                             {
  int8_t* poly_coords = p.ptr_coords;
  int32_t poly_coords_size = p.getCoordsSize();
  int32_t* poly_ring_sizes = reinterpret_cast<int32_t*>(p.getRingSizes());
  int32_t poly_num_rings = p.getNumRings();
  int32_t ic = p.getCompression();
  int32_t isr = p.getInputSrid();
  int32_t osr = p.getOutputSrid();

  std::vector<double> coords =
      *Geospatial::decompress_coords<double, int32_t>(ic, poly_coords, poly_coords_size);

  std::vector<int32_t> ring_sizes(poly_ring_sizes, poly_ring_sizes + poly_num_rings);
  auto polygon = Geospatial::GeoPolygon(coords, ring_sizes);
  polygon.transform(isr, osr);
  return polygon;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double tol

(

int32_t

ic

)

Definition at line 35 of file ExtensionFunctionsGeo.hpp.

References COMPRESSION_GEOINT32, TOLERANCE_DEFAULT, and TOLERANCE_GEOINT32.

Referenced by ST_Contains_Point_Point(), and tol().

                                            {
  if (ic == COMPRESSION_GEOINT32) {
    return TOLERANCE_GEOINT32;
  }
  return TOLERANCE_DEFAULT;
}

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE double tol	(	const int32_t	ic1,
		const int32_t	ic2
	)

Definition at line 43 of file ExtensionFunctionsGeo.hpp.

References tol().

                                                                      {
  return fmax(tol(ic1), tol(ic2));
}

Here is the call graph for this function:

DEVICE ALWAYS_INLINE bool tol_eq	(	const double	x,
		const double	y,
		const double	tolerance = TOLERANCE_DEFAULT
	)

Definition at line 52 of file ExtensionFunctionsGeo.hpp.

Referenced by polygon_contains_point(), ST_Contains_Point_LineString(), and ST_Contains_Point_Point().

                                                                             {
  auto diff = x - y;
  return (-tolerance <= diff) && (diff <= tolerance);
}

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool tol_ge	(	const double	x,
		const double	y,
		const double	tolerance = TOLERANCE_DEFAULT
	)

Definition at line 65 of file ExtensionFunctionsGeo.hpp.

Referenced by box_contains_point(), and on_segment().

                                                                             {
  return (x + tolerance) >= y;
}

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool tol_le	(	const double	x,
		const double	y,
		const double	tolerance = TOLERANCE_DEFAULT
	)

Definition at line 59 of file ExtensionFunctionsGeo.hpp.

Referenced by box_contains_point(), and on_segment().

                                                                             {
  return x <= (y + tolerance);
}

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool tol_zero	(	const double	x,
		const double	tolerance = TOLERANCE_DEFAULT
	)

Definition at line 47 of file ExtensionFunctionsGeo.hpp.

Referenced by distance_point_line(), distance_point_line_squared(), distance_point_multilinestring(), distance_point_multipolygon(), distance_point_point_squared(), distance_ring_linestring(), distance_ring_ring(), orientation(), polygon_contains_point(), ST_Contains_LineString_Point(), ST_Distance_LineString_LineString(), ST_Distance_LineString_MultiPolygon(), ST_Distance_LineString_Polygon(), ST_Distance_MultiPolygon_MultiPolygon(), ST_Distance_Polygon_MultiPolygon(), ST_Distance_Polygon_Polygon(), ST_Intersects_LineString_Linestring(), ST_Intersects_LineString_MultiPolygon(), ST_Intersects_LineString_Polygon(), ST_Intersects_MultiPolygon_MultiPolygon(), ST_Intersects_Point_LineString(), ST_Intersects_Point_Point(), ST_Intersects_Polygon_MultiPolygon(), and ST_Intersects_Polygon_Polygon().

                                                                               {
  return (-tolerance <= x) && (x <= tolerance);
}

Here is the caller graph for this function:

template<typename T >

DEVICE ALWAYS_INLINE bool tol_zero_template	(	const T	x,
		const T	tolerance = TOLERANCE_DEFAULT
	)

Definition at line 563 of file ExtensionFunctionsGeo.hpp.

Referenced by point_in_polygon_winding_number().

                                                                                   {
  if constexpr (!std::is_floating_point<T>::value) {
    return x == 0;  // assume integer 0s are always 0
  }
  return (-tolerance <= x) && (x <= tolerance);
}

Here is the caller graph for this function:

template<Coords C>

DEVICE ALWAYS_INLINE Point2D transform_point	(	const Point2D	point,
		const int32_t	isr,
		const int32_t	osr
	)

Definition at line 155 of file ExtensionFunctionsGeo.hpp.

References is_utm_srid(), and to_string().

                                                                {
  if (isr == osr || osr == 0) {
    return point;
  } else if (isr == 4326) {
    if (osr == 900913) {
      return conv_4326_900913<C>(point);  // WGS 84 --> Web Mercator
#ifdef ENABLE_UTM_TRANSFORM
    } else if (is_utm_srid(osr)) {
      return conv_4326_utm<C>(point, osr);  // WGS 84 --> UTM
    }
  } else if (is_utm_srid(isr)) {
    if (osr == 4326) {
      return conv_utm_4326<C>(point, osr);  // UTM --> WGS 84
#endif
    }
  }
#ifdef __CUDACC__
  return {};  // (NaN,NaN)
#else
  throw std::runtime_error("Unhandled geo transformation from " + std::to_string(isr) +
                           " to " + std::to_string(osr) + '.');
#endif
}

Here is the call graph for this function:

DEVICE ALWAYS_INLINE CoordData trim_linestring_to_buffered_box	(	int8_t *	l1,
		int64_t	l1size,
		int32_t	ic1,
		int32_t	isr1,
		double *	bounds2,
		int64_t	bounds2_size,
		int32_t	isr2,
		int32_t	osr,
		double	distance
	)

Definition at line 982 of file ExtensionFunctionsGeo.hpp.

References compression_unit_size(), and get_point().

Referenced by ST_DWithin_LineString_MultiPolygon(), and ST_DWithin_LineString_Polygon().

                                                                                {
  // TODO: tolerance

  // Bounds come in corresponding input spatial references.
  // For example, here the first bounding box may come in 4326, second in 900913:
  //   ST_DWithin(ST_Transform(linestring4326, 900913), linestring900913, 10)
  // Distance is given in osr spatial reference units, in this case 10 meters.
  // Bounds should be converted to osr before calculations, if isr != osr.

  Point2D const lb = transform_point<XY>({bounds2[0], bounds2[1]}, isr2, osr);
  Point2D const rt = transform_point<XY>({bounds2[2], bounds2[3]}, isr2, osr);
  CoordData l1_relevant_section;
  // Interate through linestring segments.
  // Mark the start of the relevant section FIRST TIME WE ENTER the buffered box.
  // Set the size of the relevant section the LAST TIME WE GET OUT of the box.
  auto l1_num_coords = l1size / compression_unit_size(ic1);
  Point2D l11 = get_point(l1, 0, ic1, isr1, osr);
  for (int32_t i1 = 2; i1 < l1_num_coords; i1 += 2) {
    Point2D l12 = get_point(l1, i1, ic1, isr1, osr);
    // Check for overlap between the line box and the buffered box
    if (
        // line box is left of buffered box
        fmax(l11.x, l12.x) + distance < lb.x ||
        // line box is right of buffered box
        fmin(l11.x, l12.x) - distance > rt.x ||
        // line box is below buffered box
        fmax(l11.y, l12.y) + distance < lb.y ||
        // line box is above buffered box
        fmin(l11.y, l12.y) - distance > rt.y) {
      // No overlap
      if (l1_relevant_section.ptr && l1_relevant_section.size == 0) {
        // Can finalize relevant section size
        l1_relevant_section.size =
            l1 + i1 * compression_unit_size(ic1) - l1_relevant_section.ptr;
      }
    } else {
      // Overlap
      if (!l1_relevant_section.ptr) {
        // Start tracking relevant section
        l1_relevant_section.ptr = l1 + (i1 - 2) * compression_unit_size(ic1);
      }
      // Keep relevant section size unfinalized while there is overlap,
      // linestring may reenter the buffered box.
      l1_relevant_section.size = 0;
    }
    // Advance to the next line segment.
    l11.x = l12.x;
    l11.y = l12.y;
  }
  if (l1_relevant_section.ptr && l1_relevant_section.size == 0) {
    // If we did start tracking the relevant section (i.e. entered the buffered bbox)
    // and haven't finalized the size (i.e. stayed in the bbox until the end)
    // finalize now based on the last point.
    l1_relevant_section.size = l1size - (l1_relevant_section.ptr - l1);
  }
  if (l1_relevant_section.ptr &&
      (l1_relevant_section.size < 4 * compression_unit_size(ic1))) {
    // Haven't found a section with at least 1 segment: zoom out
    l1_relevant_section = {l1, l1size};
  }

  return l1_relevant_section;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DEVICE ALWAYS_INLINE bool x_and_y_are_dependent	(	const int32_t	isr,
		const int32_t	osr
	)

Definition at line 182 of file ExtensionFunctionsGeo.hpp.

References is_utm_srid().

Referenced by coord().

                                                                                      {
#ifdef ENABLE_UTM_TRANSFORM
  return isr != osr && (is_utm_srid(isr) || is_utm_srid(osr));
#else
  return false;
#endif
}

Here is the call graph for this function:

Here is the caller graph for this function: