postgraph-postgis.sql.in

/*
 * Copyright (C) 2023 PostGraphDB
 *  
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU Affero General Public License as
 * published by the Free Software Foundation, either version 3 of the
 * License, or (at your option) any later version.
 *  
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU Affero General Public License for more details.
 *
 * You should have received a copy of the GNU Affero General Public License
 * along with this program.  If not, see <https://www.gnu.org/licenses/>.
 */

-----------------------------------------------------------------------------
-- GEOMETRY Operators
-----------------------------------------------------------------------------


-- Availability: 2.2.0
CREATE OR REPLACE FUNCTION ST_ShiftLongitude(gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_longitude_shift';

-- TODO: ST_Wrapx




-----------------------------------------------------------------------------
--  BOX3D FUNCTIONS
-----------------------------------------------------------------------------

CREATE FUNCTION ST_XMin (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_xmin';

CREATE FUNCTION ST_XMax (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_xmax';

CREATE FUNCTION ST_ZMin (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_zmin';

CREATE FUNCTION ST_YMin (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_ymin';

CREATE FUNCTION ST_YMax (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_ymax';

CREATE FUNCTION ST_ZMax (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_zmax';



-----------------------------------------------------------------------------
--  Overloaded ST_Expand FUNCTIONS
-----------------------------------------------------------------------------

-- Availability: 1.2.2
CREATE OR REPLACE FUNCTION ST_Expand(gtype,gtype)
RETURNS gtype
AS 'MODULE_PATHNAME', 'gtype_ST_Expand'
LANGUAGE 'c' IMMUTABLE STRICT PARALLEL SAFE;

CREATE OR REPLACE FUNCTION ST_Expand(gtype, dx float8, dy float8)
RETURNS gtype
AS 'MODULE_PATHNAME', 'gtype_ST_Expand'
LANGUAGE 'c' IMMUTABLE STRICT PARALLEL SAFE;

-- Availability: 2.3.0
CREATE OR REPLACE FUNCTION ST_Expand(box gtype, dx float8, dy float8, dz float8 DEFAULT 0)
RETURNS gtype
AS 'MODULE_PATHNAME', 'gtype_ST_Expand'
LANGUAGE 'c' IMMUTABLE STRICT PARALLEL SAFE;

-- Availability: 2.3.0
CREATE OR REPLACE FUNCTION ST_Expand(geom gtype, dx gtype, dy gtype, dz gtype DEFAULT '0.0'::gtype, dm gtype DEFAULT '0.0'::gtype)
RETURNS gtype
AS 'MODULE_PATHNAME', 'gtype_ST_Expand'
LANGUAGE 'c' IMMUTABLE STRICT PARALLEL SAFE;


-----------------------------------------------------------------------------
--  BOX2D FUNCTIONS
-----------------------------------------------------------------------------

-- Availability: 1.5.0
CREATE OR REPLACE FUNCTION postgis_getbbox(gtype)
RETURNS gtype
AS 'MODULE_PATHNAME','gtype_to_BOX2DF'
LANGUAGE 'c' IMMUTABLE STRICT PARALLEL SAFE;


-- Availability: 1.2.2
CREATE OR REPLACE FUNCTION ST_MakeBox2d(geom1 gtype, geom2 gtype)
RETURNS gtype
AS 'MODULE_PATHNAME', 'gtype_construct'
LANGUAGE 'c' IMMUTABLE STRICT PARALLEL SAFE;

--
-- PostGIS 2D Geometry Operators
--
CREATE OR REPLACE FUNCTION gtype_geometry_same(geom1 gtype, geom2 gtype)
RETURNS boolean
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_same_2d';

CREATE OPERATOR ~= (
    LEFTARG = gtype, 
    RIGHTARG = gtype,
    PROCEDURE = gtype_geometry_same,
    RESTRICT = contsel,
    JOIN = contjoinsel
);

CREATE OR REPLACE FUNCTION gtype_geometry_within(geom1 gtype, geom2 gtype)
RETURNS bool
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_within_2d';

CREATE OPERATOR @ (
LEFTARG = gtype,
RIGHTARG = gtype,
PROCEDURE = gtype_geometry_within--,
-- Updated: 3.4.0 to use selectivity estimator
--RESTRICT = gtype_gserialized_gist_sel_2d,
-- Updated: 3.4.0 to use join selectivity estimator
--JOIN = gtype_gserialized_gist_joinsel_2d
);

CREATE OR REPLACE FUNCTION gtype_geometry_overleft(geom1 gtype, geom2 gtype)
RETURNS bool
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_overleft_2d';

CREATE OPERATOR &< (
LEFTARG = gtype,
RIGHTARG = gtype,
PROCEDURE = gtype_geometry_overleft,
RESTRICT = positionsel,
JOIN = positionjoinsel
);

CREATE OR REPLACE FUNCTION gtype_geometry_below(geom1 gtype, geom2 gtype)
RETURNS bool
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_below_2d';


CREATE OPERATOR <<| (
LEFTARG = gtype,
RIGHTARG = gtype,
PROCEDURE = gtype_geometry_below,
--COMMuTAOR = '<<|',
RESTRICT = positionsel,
JOIN = positionjoinsel
);


CREATE OR REPLACE FUNCTION gtype_geometry_overbelow(geom1 gtype, geom2 gtype)
RETURNS bool
LANGUAGE c
IMMUTABLE 
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_overbelow_2d';

CREATE OPERATOR &<| (
LEFTARG = gtype,
RIGHTARG = gtype,
PROCEDURE = gtype_geometry_overbelow,
RESTRICT = positionsel,
JOIN = positionjoinsel
);


CREATE OR REPLACE FUNCTION gtype_geometry_overright(geom1 gtype, geom2 gtype)
RETURNS bool
LANGUAGE c
IMMUTABLE 
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_overright_2d';

CREATE OPERATOR &> (
LEFTARG = gtype,
RIGHTARG = gtype,
PROCEDURE = gtype_geometry_overright,
RESTRICT = positionsel,
JOIN = positionjoinsel
);

CREATE OR REPLACE FUNCTION gtype_geometry_overabove(geom1 gtype, geom2 gtype)
RETURNS bool
LANGUAGE c
IMMUTABLE 
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_overabove_2d';

CREATE OPERATOR |&> (
LEFTARG = gtype,
RIGHTARG = gtype,
PROCEDURE = gtype_geometry_overabove,
RESTRICT = positionsel,
JOIN = positionjoinsel
);

CREATE OR REPLACE FUNCTION gtype_geometry_above(geom1 gtype, geom2 gtype)
RETURNS bool
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_above_2d';

CREATE OPERATOR |>> (
LEFTARG = gtype,
RIGHTARG = gtype,
PROCEDURE = gtype_geometry_above,
--COMMuTAOR = '<<|',
RESTRICT = positionsel,
JOIN = positionjoinsel
);

--
-- 2D Gist Index
--
CREATE OR REPLACE FUNCTION geometry_gist_distance_2d(internal, gtype, integer)
RETURNS float8
AS 'MODULE_PATHNAME' ,'gtype_gserialized_gist_distance_2d'
LANGUAGE 'c' PARALLEL SAFE;

CREATE OR REPLACE FUNCTION geometry_gist_consistent_2d(internal,gtype,integer)
RETURNS bool
AS 'MODULE_PATHNAME' ,'gtype_gserialized_gist_consistent_2d'
LANGUAGE 'c' PARALLEL SAFE;

CREATE OR REPLACE FUNCTION geometry_gist_compress_2d(internal)
RETURNS internal
AS 'MODULE_PATHNAME','gtype_gserialized_gist_compress_2d'
LANGUAGE 'c' PARALLEL SAFE;

CREATE OR REPLACE FUNCTION geometry_gist_penalty_2d(internal,internal,internal)
RETURNS internal
AS 'MODULE_PATHNAME' ,'gtype_gserialized_gist_penalty_2d'
LANGUAGE 'c' PARALLEL SAFE;

CREATE OR REPLACE FUNCTION geometry_gist_picksplit_2d(internal, internal)
RETURNS internal
AS 'MODULE_PATHNAME' ,'gtype_gserialized_gist_picksplit_2d'
LANGUAGE 'c' PARALLEL SAFE;

CREATE OR REPLACE FUNCTION geometry_gist_union_2d(bytea, internal)
RETURNS internal
AS 'MODULE_PATHNAME' ,'gtype_gserialized_gist_union_2d'
LANGUAGE 'c' PARALLEL SAFE;

CREATE OR REPLACE FUNCTION geometry_gist_same_2d(geom1 gtype, geom2 gtype, internal)
RETURNS internal
AS 'MODULE_PATHNAME' ,'gtype_gserialized_gist_same_2d'
LANGUAGE 'c' PARALLEL SAFE;

CREATE OR REPLACE FUNCTION geometry_gist_decompress_2d(internal)
RETURNS internal
AS 'MODULE_PATHNAME' ,'gtype_gserialized_gist_decompress_2d'
LANGUAGE 'c' PARALLEL SAFE;

CREATE OR REPLACE FUNCTION geometry_gist_sortsupport_2d(internal)
RETURNS void
AS 'MODULE_PATHNAME', 'gtype_gserialized_gist_sortsupport_2d'
LANGUAGE 'c' STRICT;

CREATE OPERATOR CLASS gist_geometry_ops_2d
        DEFAULT FOR TYPE gtype USING GIST AS
        STORAGE public.box2df,
        OPERATOR        1        <<  ,
        OPERATOR        2        &<      ,
        OPERATOR        3        &&  ,
        OPERATOR        4        &>      ,
        OPERATOR        5        >>      ,
        OPERATOR        6        ~=      ,
        OPERATOR        7        ~       ,
        OPERATOR        8        @       ,
        OPERATOR        9        &<| ,
        OPERATOR        10       <<| ,
        OPERATOR        11       |>> ,
        OPERATOR        12       |&> ,
        OPERATOR        13       <-> FOR ORDER BY gtype_ops_btree,
        OPERATOR        14       <#> FOR ORDER BY gtype_ops_btree,
--
-- Sort support in bulk indexing not included in the default
-- opclass for PostgreSQL versions <15 due to query performance
-- degradation caused by GiST index page overlap.
-- Since PostgreSQL 15 sorting build uses picksplit function to
-- find better partitioning for index records. This allows to
-- build indices performing similarly to those produced by default
-- method while still reducing index build time significantly.
--
-- To enable sortsupport:
--   alter operator family gist_geometry_ops_2d using gist
--     add function 11 (geometry)
--     geometry_gist_sortsupport_2d (internal);
--
-- To remove sortsupport:
--   alter operator family gist_geometry_ops_2d using gist
--     drop function 11 (geometry);
--
        FUNCTION        8        geometry_gist_distance_2d (internal, gtype, integer),
        FUNCTION        1        geometry_gist_consistent_2d (internal, gtype, integer),
        FUNCTION        2        geometry_gist_union_2d (bytea, internal),
        FUNCTION        3        geometry_gist_compress_2d (internal),
        FUNCTION        4        geometry_gist_decompress_2d (internal),
        FUNCTION        5        geometry_gist_penalty_2d (internal, internal, internal),
        FUNCTION        6        geometry_gist_picksplit_2d (internal, internal),
        FUNCTION        7        geometry_gist_same_2d (geom1 gtype, geom2 gtype, internal)/*,
        FUNCTION        11       geometry_gist_sortsupport_2d(gtype)*/;

--
-- PostGIS N-D Geometry Operators
--
CREATE OR REPLACE FUNCTION geometry_overlaps_nd(gtype, gtype)
RETURNS boolean
AS 'MODULE_PATHNAME' ,'gtype_gserialized_overlaps'
LANGUAGE 'c' IMMUTABLE STRICT PARALLEL SAFE
COST 1;

CREATE OPERATOR &&& (
    LEFTARG = gtype, RIGHTARG = gtype, PROCEDURE = geometry_overlaps_nd,
    COMMUTATOR = '&&&'--,
    --RESTRICT = gserialized_gist_sel_nd,
    --JOIN = gserialized_gist_joinsel_nd
);

CREATE OR REPLACE FUNCTION geometry_contains_nd(gtype, gtype)
RETURNS boolean
AS 'MODULE_PATHNAME' ,'gtype_gserialized_contains'
LANGUAGE 'c' IMMUTABLE STRICT PARALLEL SAFE
COST 1;

CREATE OPERATOR ~~ (
        LEFTARG = gtype, RIGHTARG = gtype, PROCEDURE = geometry_contains_nd--,
        --COMMUTATOR = '@@',
        --RESTRICT = gserialized_gist_sel_nd,
        --JOIN = gserialized_gist_joinsel_nd
);

CREATE OR REPLACE FUNCTION geometry_within_nd(gtype, gtype)
RETURNS boolean
AS 'MODULE_PATHNAME' ,'gtype_gserialized_within'
LANGUAGE 'c' IMMUTABLE STRICT PARALLEL SAFE
COST 1;

CREATE OPERATOR @@ (
    LEFTARG = gtype, RIGHTARG = gtype, PROCEDURE = geometry_within_nd,
    COMMUTATOR = '~~'/*,
    RESTRICT = gserialized_gist_sel_nd,
    JOIN = gserialized_gist_joinsel_nd*/
);

CREATE OR REPLACE FUNCTION geometry_same_nd(gtype, gtype)
RETURNS boolean
AS 'MODULE_PATHNAME' ,'gtype_gserialized_same'
LANGUAGE 'c' IMMUTABLE STRICT PARALLEL SAFE
COST 1;

CREATE OPERATOR ~~= (
        LEFTARG = gtype, RIGHTARG = gtype, PROCEDURE = geometry_same_nd,
        COMMUTATOR = '~~='/*,
        RESTRICT = gserialized_gist_sel_nd,
        JOIN = gserialized_gist_joinsel_nd*/
);

CREATE FUNCTION ST_DistanceCPA (gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 5000
AS 'MODULE_PATHNAME', 'gtype_ST_DistanceCPA';

CREATE OPERATOR |=| (
    LEFTARG = gtype, RIGHTARG = gtype,
    PROCEDURE = ST_DistanceCPA,
    COMMUTATOR = '|=|'
);

-- Availability: 2.2.0
CREATE OR REPLACE FUNCTION gtype_distance_centroid_nd(gtype, gtype)
RETURNS gtype
IMMUTABLE 
STRICT 
PARALLEL SAFE
LANGUAGE 'c' 
AS 'MODULE_PATHNAME', 'gtype_distance_nd';

-- Availability: 2.2.0
CREATE OPERATOR <<->> (
	LEFTARG = gtype,
    RIGHTARG = gtype,
	PROCEDURE = gtype_distance_centroid_nd,
	COMMUTATOR = '<<->>'
);


-- Availability: 2.2.0
CREATE OR REPLACE FUNCTION gtype_gist_distance_nd(internal,gtype,integer)
RETURNS gtype
PARALLEL SAFE
LANGUAGE 'c' 
AS 'MODULE_PATHNAME', 'gtype_gist_distance';

-- Availability: 2.0.0
CREATE OR REPLACE FUNCTION gtype_gist_consistent_nd(internal,gtype,integer)
RETURNS gtype
PARALLEL SAFE
LANGUAGE 'c' 
AS 'MODULE_PATHNAME', 'gtype_gist_consistent';

-- Availability: 2.0.0
CREATE OR REPLACE FUNCTION gtype_gist_compress_nd(internal)
RETURNS internal
PARALLEL SAFE
LANGUAGE 'c' 
AS 'MODULE_PATHNAME','gtype_gist_compress';

-- Availability: 2.0.0
CREATE OR REPLACE FUNCTION gtype_gist_penalty_nd(internal,internal,internal)
RETURNS internal
PARALLEL SAFE
LANGUAGE 'c' 
AS 'MODULE_PATHNAME' ,'gtype_gist_penalty';

-- Availability: 2.0.0
CREATE OR REPLACE FUNCTION gtype_gist_picksplit_nd(internal, internal)
RETURNS internal
PARALLEL SAFE
LANGUAGE 'c' 
AS 'MODULE_PATHNAME' ,'gtype_gist_picksplit';

-- Availability: 2.0.0
CREATE OR REPLACE FUNCTION gtype_gist_union_nd(bytea, internal)
RETURNS internal
PARALLEL SAFE
LANGUAGE 'c' 
AS 'MODULE_PATHNAME' ,'gtype_gist_union';

-- Availability: 2.0.0
CREATE OR REPLACE FUNCTION gtype_gist_same_nd(gtype, gtype, internal)
RETURNS internal
PARALLEL SAFE
LANGUAGE 'c' 
AS 'MODULE_PATHNAME' ,'gtype_gist_same';

-- Availability: 2.0.0
CREATE OR REPLACE FUNCTION gtype_gist_decompress_nd(internal)
RETURNS internal
PARALLEL SAFE
LANGUAGE 'c' 
AS 'MODULE_PATHNAME' ,'gtype_gist_decompress';


-- Availability: 2.0.0
CREATE OPERATOR CLASS gist_geometry_ops_nd
	FOR TYPE gtype USING GIST AS
	STORAGE 	public.gidx,
	OPERATOR        3        &&&	,
	-- Availability: 3.0.0
	OPERATOR        6        ~~=	,
	-- Availability: 3.0.0
	OPERATOR        7        ~~	,
	-- Availability: 3.0.0
	OPERATOR        8        @@	,
	-- Availability: 2.2.0
	OPERATOR        13       <<->> FOR ORDER BY pg_catalog.float_ops,
	-- Availability: 2.2.0
	OPERATOR        20       |=| FOR ORDER BY pg_catalog.float_ops,
	-- Availability: 2.2.0
	FUNCTION        8        gtype_gist_distance_nd (internal, gtype, integer),
	FUNCTION        1        gtype_gist_consistent_nd (internal, gtype, integer),
	FUNCTION        2        gtype_gist_union_nd (bytea, internal),
	FUNCTION        3        gtype_gist_compress_nd (internal),
	FUNCTION        4        gtype_gist_decompress_nd (internal),
	FUNCTION        5        gtype_gist_penalty_nd (internal, internal, internal),
	FUNCTION        6        gtype_gist_picksplit_nd (internal, internal),
	FUNCTION        7        gtype_gist_same_nd (gtype, gtype, internal);

--
-- PostGIS Functions
--

--
-- PostGIS Constructors
--
CREATE FUNCTION ST_MakePoint (gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_makepoint';

CREATE FUNCTION ST_MakePoint (gtype, gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_makepoint';

CREATE FUNCTION ST_MakePoint (gtype, gtype, gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_makepoint';

CREATE FUNCTION ST_MakePointM (gtype, gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_makepoint3dm';


CREATE FUNCTION ST_asEWKT (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 250
AS 'MODULE_PATHNAME', 'gtype_asEWKT';

CREATE FUNCTION asEWKT (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 250
AS 'MODULE_PATHNAME', 'gtype_asEWKT';

CREATE FUNCTION addbbox (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_addBBOX';

CREATE FUNCTION postgis_addbbox (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_addBBOX';

CREATE FUNCTION dropbbox (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_dropBBOX';

CREATE FUNCTION postgis_dropbbox (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_dropBBOX';


CREATE FUNCTION ST_X (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_x_point';

CREATE FUNCTION X (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_x_point';

CREATE FUNCTION ST_Y (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_y_point';

CREATE FUNCTION Y (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_y_point';

CREATE FUNCTION ST_Z (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_z_point';

CREATE FUNCTION Z (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_z_point';

CREATE FUNCTION ST_M (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_m_point';

CREATE FUNCTION M (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 1
AS 'MODULE_PATHNAME', 'gtype_m_point';

--
-- PostGIS Functions - Measures
--
CREATE FUNCTION ST_3DLength(gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT 
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_length_linestring';

CREATE FUNCTION ST_Length2d(gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_length2d_linestring';

CREATE FUNCTION ST_Length(gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_length2d_linestring';

CREATE FUNCTION ST_LengthSpheroid(gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 250
AS 'MODULE_PATHNAME', 'gtype_length_ellipsoid_linestring';



CREATE FUNCTION ST_IsPolygonCW (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_ST_IsPolygonCW';

CREATE FUNCTION ST_IsPolygonCCW (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_ST_IsPolygonCCW';

CREATE FUNCTION ST_DistanceSpheroid (gtype, gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 250
AS 'MODULE_PATHNAME', 'gtype_distance_ellipsoid';

CREATE FUNCTION ST_DistanceSpheroid (gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 250
AS 'MODULE_PATHNAME', 'gtype_distance_ellipsoid';

CREATE FUNCTION ST_azimuth (gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_azimuth';


--
-- PostGIS Functions - Temporal
--
CREATE FUNCTION ST_ClosestPointOfApproach (gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 5000
AS 'MODULE_PATHNAME', 'gtype_ST_ClosestPointOfApproach';



CREATE FUNCTION ST_IsValidTrajectory (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 5000
AS 'MODULE_PATHNAME', 'gtype_st_isvalidtrajectory';

CREATE FUNCTION ST_CPAWithin (gtype, gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 5000
AS 'MODULE_PATHNAME', 'gtype_ST_CPAWithin';


--
-- PostGIS Functions - GEOS
--
CREATE FUNCTION ST_Intersection (gtype, gtype, gtype DEFAULT '-1.0'::gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 5000
AS 'MODULE_PATHNAME', 'gtype_st_intersection';

CREATE FUNCTION ST_GeneratePoints (area gtype, npoints gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 250
AS 'MODULE_PATHNAME', 'gtype_st_generatepoints';

CREATE FUNCTION ST_GeneratePoints (area gtype, npoints gtype, seed gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 250
AS 'MODULE_PATHNAME', 'gtype_st_generatepoints';

CREATE FUNCTION ST_ConvexHull (gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 5000
AS 'MODULE_PATHNAME', 'gtype_convexhull';

CREATE FUNCTION ST_SymDifference (gtype, gtype, gtype DEFAULT '-1.0'::gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 5000
AS 'MODULE_PATHNAME', 'gtype_st_symdifference';

CREATE FUNCTION ST_SymDifference (gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 5000
AS 'MODULE_PATHNAME', 'gtype_st_symdifference';

CREATE FUNCTION ST_HausdorffDistance (gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 5000
AS 'MODULE_PATHNAME', 'gtype_hausdorffdistance';

CREATE FUNCTION ST_HausdorffDistance (gtype, gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 5000
AS 'MODULE_PATHNAME', 'gtype_hausdorffdistancedensify';

CREATE FUNCTION ST_FrechetDistance (gtype, gtype, gtype default '-1.0'::gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 5000
AS 'MODULE_PATHNAME', 'gtype_st_frechet_distance';

--
-- Algorithms
--
CREATE FUNCTION ST_Simplify (gtype, gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_ST_Simplify';

CREATE FUNCTION ST_Simplify (gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_ST_Simplify';

CREATE OR REPLACE FUNCTION ST_SnapToGrid(gtype, gtype, gtype, gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_snaptogrid';

CREATE OR REPLACE FUNCTION ST_SnapToGrid(gtype, gtype, gtype)
RETURNS gtype
AS $$ SELECT ST_SnapToGrid($1, '0'::gtype, '0'::gtype, $2, $3) $$
LANGUAGE sql 
IMMUTABLE 
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50;

CREATE OR REPLACE FUNCTION ST_SnapToGrid(gtype, gtype)
RETURNS gtype
AS $$ SELECT ST_SnapToGrid($1, '0'::gtype, '0'::gtype, $2, $2) $$
LANGUAGE sql
IMMUTABLE 
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50;

CREATE OR REPLACE FUNCTION ST_SnapToGrid(gtype, gtype, gtype, gtype, gtype, gtype)
RETURNS gtype
LANGUAGE c 
IMMUTABLE 
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_snaptogrid_pointoff';


--
-- PostGIS Affine Functions
--
CREATE FUNCTION ST_Affine (gtype, gtype, gtype, gtype, gtype, gtype, gtype, gtype, gtype, gtype, gtype, gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_affine';

CREATE OR REPLACE FUNCTION ST_Affine(gtype, gtype, gtype, gtype, gtype, gtype, gtype)
RETURNS gtype
AS 'SELECT ST_Affine($1,  $2, $3, 0,  $4, $5, 0,  0, 0, 1,  $6, $7, 0)'
LANGUAGE 'sql'
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50;

CREATE OR REPLACE FUNCTION ST_Rotate(gtype, gtype)
RETURNS gtype
AS 'SELECT ST_Affine($1,  cos($2), -sin($2), 0,  sin($2), cos($2), 0,  0, 0, 1,  0, 0, 0)'
LANGUAGE 'sql'
IMMUTABLE 
RETURNS NULL ON NULL INPUT 
PARALLEL SAFE
COST 50;

CREATE OR REPLACE FUNCTION ST_Rotate(gtype, gtype, gtype, gtype)
RETURNS gtype
AS 'SELECT ST_Affine($1,  cos($2), -sin($2), 0,  sin($2),  cos($2), 0, 0, 0, 1, $3 - cos($2) * $3 + sin($2) * $4, $4 - sin($2) * $3 - cos($2) * $4, 0)'
LANGUAGE 'sql' 
IMMUTABLE 
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50;

CREATE OR REPLACE FUNCTION ST_Rotate(gtype, gtype, gtype)
RETURNS gtype
AS 'SELECT ST_Affine($1,  cos($2), -sin($2), 0,  sin($2),  cos($2), 0, 0, 0, 1, ST_X($3) - cos($2) * ST_X($3) + sin($2) * ST_Y($3), ST_Y($3) - sin($2) * ST_X($3) - cos($2) * ST_Y($3), 0)'
LANGUAGE 'sql'
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50;

CREATE OR REPLACE FUNCTION ST_RotateZ(gtype, gtype)
RETURNS gtype
AS 'SELECT ST_Rotate($1, $2)'
LANGUAGE 'sql'
IMMUTABLE 
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50;

CREATE OR REPLACE FUNCTION ST_RotateX(gtype, gtype)
RETURNS gtype
AS 'SELECT ST_Affine($1, 1, 0, 0, 0, cos($2), -sin($2), 0, sin($2), cos($2), 0, 0, 0)'
LANGUAGE 'sql'
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50;

CREATE OR REPLACE FUNCTION ST_RotateY(gtype, gtype)
RETURNS gtype
AS 'SELECT ST_Affine($1,  cos($2), 0, sin($2),  0, 1, 0,  -sin($2), 0, cos($2), 0,  0, 0)'
LANGUAGE 'sql'
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50;

CREATE OR REPLACE FUNCTION ST_Translate(gtype, gtype, gtype, gtype)
RETURNS gtype
AS $$ SELECT ST_Affine($1, '1'::gtype, '0'::gtype, '0'::gtype, '0'::gtype, '1'::gtype, '0'::gtype, '0'::gtype, '0'::gtype, '1'::gtype, $2, $3, $4) $$
LANGUAGE 'sql'
IMMUTABLE 
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50;

CREATE OR REPLACE FUNCTION ST_Translate(gtype, gtype, gtype)
RETURNS gtype
AS 'SELECT ST_Translate($1, $2, $3, 0)'
LANGUAGE 'sql'
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50;

CREATE OR REPLACE FUNCTION ST_Scale(gtype, gtype)
RETURNS gtype
LANGUAGE c 
IMMUTABLE 
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_st_scale_2_args';

CREATE OR REPLACE FUNCTION ST_Scale(gtype, gtype, gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_st_scale';

CREATE OR REPLACE FUNCTION ST_Scale(gtype, gtype, gtype, gtype)
RETURNS gtype
AS 'SELECT ST_Scale($1, ST_MakePoint($2, $3, $4))'
LANGUAGE sql
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50;
/* XXX: Gonna Have to Get Creative With This One
CREATE OR REPLACE FUNCTION ST_Scale(gtype, gtype, gtype)
RETURNS gtype
AS 'SELECT ST_Scale($1, $2, $3, 1.0::float8::gtype)'
LANGUAGE sql
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50;
*/
--
-- MISC
--

CREATE OR REPLACE FUNCTION ST_Force2D(gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_force_2d';

CREATE OR REPLACE FUNCTION ST_Force3DZ(gtype, gtype default '0.0')
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_force_3dz';

CREATE OR REPLACE FUNCTION ST_Force3DM(gtype, gtype default '0.0')
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_force_3dm';

CREATE OR REPLACE FUNCTION ST_Force4D(gtype, gtype default '0.0', gtype default '0.0')
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_force_4d';

CREATE OR REPLACE FUNCTION ST_ForceCollection(gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_force_collection';

CREATE OR REPLACE FUNCTION ST_Multi(gtype)
RETURNS gtype
LANGUAGE c
IMMUTABLE
RETURNS NULL ON NULL INPUT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_force_multi';

---------------------------------------------------------------
-- 3D-functions
---------------------------------------------------------------

CREATE OR REPLACE FUNCTION ST_3DDistance(geom1 gtype, geom2 gtype)
RETURNS gtype
LANGUAGE 'c'
IMMUTABLE
STRICT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_3DDistance';

CREATE OR REPLACE FUNCTION ST_3DMaxDistance(geom1 gtype, geom2 gtype)
RETURNS gtype
LANGUAGE 'c'
IMMUTABLE
STRICT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_maxdistance3d';

CREATE OR REPLACE FUNCTION ST_3DClosestPoint(geom1 gtype, geom2 gtype)
RETURNS gtype
LANGUAGE 'c'
IMMUTABLE 
STRICT 
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_closestpoint3d';

CREATE OR REPLACE FUNCTION ST_3DShortestLine(geom1 gtype, geom2 gtype)
RETURNS gtype
LANGUAGE 'c'
IMMUTABLE
STRICT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_shortestline3d';

CREATE OR REPLACE FUNCTION ST_3DLongestLine(geom1 gtype, geom2 gtype)
RETURNS gtype
LANGUAGE 'c'
IMMUTABLE 
STRICT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME', 'gtype_longestline3d';


-- PostGIS equivalent function: AsText(geometry)
CREATE OR REPLACE FUNCTION ST_AsText(gtype)
RETURNS gtype
LANGUAGE 'c'
IMMUTABLE 
STRICT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME','gtype_asText';

-- Availability: 2.5.0
-- PostGIS equivalent function: AsText(gtype, gtype)
CREATE OR REPLACE FUNCTION ST_AsText(gtype, gtype)
RETURNS gtype
LANGUAGE 'c'
IMMUTABLE 
STRICT
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME','gtype_asText';


CREATE OR REPLACE FUNCTION ST_3DDWithin(geom1 gtype, geom2 gtype,gtype)
RETURNS gtype
LANGUAGE 'c' 
IMMUTABLE 
STRICT 
PARALLEL SAFE
COST 500
AS 'MODULE_PATHNAME', 'gtype_dwithin3d'
;


CREATE OR REPLACE FUNCTION ST_3DDFullyWithin(geom1 gtype, geom2 gtype,gtype)
RETURNS gtype
LANGUAGE 'c' 
IMMUTABLE 
STRICT 
PARALLEL SAFE
COST 500
AS 'MODULE_PATHNAME', 'gtype_dfullywithin3d'
;


CREATE OR REPLACE FUNCTION ST_Dimension(geom gtype)
RETURNS gtype
LANGUAGE 'c' 
IMMUTABLE 
STRICT 
PARALLEL SAFE
COST 10
AS 'MODULE_PATHNAME', 'gtype_dimension'
;


CREATE OR REPLACE FUNCTION ST_GeometryType(geom gtype)
RETURNS gtype
LANGUAGE 'c' 
IMMUTABLE 
STRICT 
PARALLEL SAFE
COST 10
AS 'MODULE_PATHNAME', 'gtype_geometrytype';


-- PostGIS equivalent function: NumGeometries(geometry)
CREATE OR REPLACE FUNCTION ST_NumGeometries(geom gtype)
RETURNS gtype
LANGUAGE 'c' 
IMMUTABLE 
STRICT 
PARALLEL SAFE
COST 10
AS 'MODULE_PATHNAME', 'gtype_numgeometries_collection';

-- PostGIS equivalent function: ST_GeometryN(geometry)
CREATE OR REPLACE FUNCTION ST_GeometryN(geom gtype, gtype)
RETURNS gtype
LANGUAGE 'c' 
IMMUTABLE 
STRICT 
PARALLEL SAFE
COST 10
AS 'MODULE_PATHNAME', 'gtype_geometryn_collection';


-- Availability: 2.0.0
CREATE OR REPLACE FUNCTION ST_NumPatches(gtype)
RETURNS gtype
LANGUAGE 'sql' 
IMMUTABLE
STRICT 
PARALLEL SAFE
COST 10
AS $$
	SELECT CASE WHEN @extschema@.ST_GeometryType($1) = ('"ST_PolyhedralSurface"'::gtype)
	THEN @extschema@.ST_NumGeometries($1)
	ELSE NULL END
$$;


-- Availability: 2.0.0
CREATE OR REPLACE FUNCTION ST_PatchN(gtype, gtype)
RETURNS gtype
LANGUAGE 'sql' 
IMMUTABLE
STRICT 
PARALLEL SAFE
COST 10
AS $$
	SELECT CASE WHEN @extschema@.ST_GeometryType($1) = ('"ST_PolyhedralSurface"'::gtype)
	THEN @extschema@.ST_GeometryN($1, $2)
	ELSE NULL END
$$;

-- PostGIS equivalent function: NumGeometries(geometry)
CREATE OR REPLACE FUNCTION ST_Reverse(geom gtype)
RETURNS gtype
LANGUAGE 'c' 
IMMUTABLE 
STRICT 
PARALLEL SAFE
COST 10
AS 'MODULE_PATHNAME', 'gtype_st_reverse';

CREATE OR REPLACE FUNCTION ST_GeomFromText(gtype)
RETURNS gtype
LANGUAGE 'c' 
IMMUTABLE 
STRICT 
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME','gtype_from_text';

-- PostGIS equivalent function: ST_GeometryFromText(gtype, gtype)
CREATE OR REPLACE FUNCTION ST_GeomFromText(gtype, gtype)
RETURNS gtype
LANGUAGE 'c' 
IMMUTABLE 
STRICT 
PARALLEL SAFE
COST 50
AS 'MODULE_PATHNAME','gtype_from_text';