#! /bin/sh

set -v

DB=${1:-geo}
REL=${2:-top_polys}
RESOL=${3:-1}
MAXFN=${4:-150000}

monitor -c "destroy pol_tmp" $DB

monitor -c "retrieve into pol_tmp (ap.all) from ap in all_plines*" $DB

monitor -c "define index pol_tmp_index on pol_tmp using btree (feat_num int4_ops)" $DB
monitor -c "define index pol_top_index on poly_topology using btree (centroid_fn int4_ops)" $DB

monitor -c "destroy $REL" $DB
monitor -c "create $REL (centroid_fn= int4, geo_bbox= box, pgn= POLYGON2)" $DB

monitor -c "vacuum" $DB # Allow better Query Optimisation

monitor -TN -c "
  retrieve (poly_topology.centroid_fn, pol_tmp.pln) where
    poly_topology.e_or_i != 'i'
    and poly_topology.edge_fn = pol_tmp.feat_num
    sort by centroid_fn" $DB |
$GEOHOME/conversion/topo.awk.copy.ge $DB $REL $RESOL

monitor -c "destroy pol_tmp" $DB
monitor -c "remove index pol_top_index" $DB