Geografic-Information-System/vector/v.voronoi/sw_output.c

#include <stdio.h>
#include <math.h>
#include <float.h>
#include <grass/vector.h>
#include "sw_defs.h"
#include "write.h"

int openpl(void)
{
    return 0;
}
int line(int a, int b, int c, int d)
{
    return 0;
}
int circle(int a, int b, int c)
{
    return 0;
}
int range(int a, int b, int c, int d)
{
    return 0;
}

double pxmin, pxmax, pymin, pymax, cradius;

int out_bisector(struct Edge *e)
{
    if (triangulate & plot & !debug) {
	line(e->reg[0]->coord.x, e->reg[0]->coord.y, e->reg[1]->coord.x,
	     e->reg[1]->coord.y);
    }

    /*
       if(!triangulate & !plot &!debug) {
       fprintf (stdout,"l %.*g %.*g %.*g \n", DBL_DIG, e->a, DBL_DIG, e->b, DBL_DIG, e->c);
       }

       if(debug) {
       fprintf (stdout,"line(%d) %gx+%gy=%g, bisecting %d %d\n", e->edgenbr,
       e->a, e->b, e->c, e->reg[le]->sitenbr, e->reg[re]->sitenbr);
       }
     */

    return 0;
}

int out_ep(struct Edge *e)
{
    static struct line_pnts *Points = NULL;
    static struct line_cats *Cats = NULL;

    if (!Points) {
	Points = Vect_new_line_struct();
	Cats = Vect_new_cats_struct();
    }

    if (!triangulate & plot)
	clip_line(e);
    if (!triangulate & !plot) {
	/*
	   fprintf (stdout,"e %d", e->edgenbr);
	   fprintf (stdout," %d ", e->ep[le] != (struct Site *)NULL ? e->ep[le]->sitenbr : -1);
	   fprintf (stdout,"%d\n", e->ep[re] != (struct Site *)NULL ? e->ep[re]->sitenbr : -1);
	 */

	write_ep(e);
    }

    return 0;
}

int out_vertex(struct Site *v)
{
    /*
       if(!triangulate & !plot &!debug)
       fprintf (stdout,"v %.*g %.*g\n", DBL_DIG, v->coord.x, DBL_DIG, v->coord.y);
       if(debug)
       fprintf (stdout,"vertex(%d) at %f %f\n", v->sitenbr, v->coord.x, v->coord.y);
     */

    return 0;
}


int out_site(struct Site *s)
{
    if (!triangulate & plot & !debug)
	circle(s->coord.x, s->coord.y, cradius);

    /*
       if(!triangulate & !plot & !debug)
       fprintf (stdout,"s %f %f\n", s->coord.x, s->coord.y);
       if(debug)
       fprintf (stdout,"site (%d) at %f %f\n", s->sitenbr, s->coord.x, s->coord.y);
     */

    return 0;
}

int out_triple(struct Site *s1, struct Site *s2, struct Site *s3)
{
    if (triangulate & !plot & !debug) {
	/* fprintf (stdout,"%d %d %d\n", s1->sitenbr, s2->sitenbr, s3->sitenbr); */

	write_triple(s1, s2, s3);
    }

    /*
       if(debug)
       fprintf (stdout,"circle through left=%d right=%d bottom=%d\n", 
       s1->sitenbr, s2->sitenbr, s3->sitenbr);
     */

    return 0;
}

int plotinit(void)
{
    double dx, dy, d;

    dy = ymax - ymin;
    dx = xmax - xmin;
    d = (dx > dy ? dx : dy) * 1.1;
    pxmin = xmin - (d - dx) / 2.0;
    pxmax = xmax + (d - dx) / 2.0;
    pymin = ymin - (d - dy) / 2.0;
    pymax = ymax + (d - dy) / 2.0;
    cradius = (pxmax - pxmin) / 350.0;
    openpl();
    range(pxmin, pymin, pxmax, pymax);

    return 0;
}


int clip_line(struct Edge *e)
{
    struct Site *s1, *s2;
    double x1, x2, y1, y2;

    if (e->a == 1.0 && e->b >= 0.0) {
	s1 = e->ep[1];
	s2 = e->ep[0];
    }
    else {
	s1 = e->ep[0];
	s2 = e->ep[1];
    }

    if (e->a == 1.0) {
	y1 = pymin;
	if (s1 != (struct Site *)NULL && s1->coord.y > pymin)
	    y1 = s1->coord.y;
	if (y1 > pymax)
	    return -1;
	x1 = e->c - e->b * y1;
	y2 = pymax;
	if (s2 != (struct Site *)NULL && s2->coord.y < pymax)
	    y2 = s2->coord.y;
	if (y2 < pymin)
	    return (0);
	x2 = e->c - e->b * y2;
	if (((x1 > pxmax) & (x2 > pxmax)) | ((x1 < pxmin) & (x2 < pxmin)))
	    return -1;

	if (x1 > pxmax) {
	    x1 = pxmax;
	    y1 = (e->c - x1) / e->b;
	}
	if (x1 < pxmin) {
	    x1 = pxmin;
	    y1 = (e->c - x1) / e->b;
	}
	if (x2 > pxmax) {
	    x2 = pxmax;
	    y2 = (e->c - x2) / e->b;
	}
	if (x2 < pxmin) {
	    x2 = pxmin;
	    y2 = (e->c - x2) / e->b;
	}
    }
    else {
	x1 = pxmin;
	if (s1 != (struct Site *)NULL && s1->coord.x > pxmin)
	    x1 = s1->coord.x;
	if (x1 > pxmax)
	    return (0);
	y1 = e->c - e->a * x1;
	x2 = pxmax;
	if (s2 != (struct Site *)NULL && s2->coord.x < pxmax)
	    x2 = s2->coord.x;
	if (x2 < pxmin)
	    return (0);
	y2 = e->c - e->a * x2;
	if (((y1 > pymax) & (y2 > pymax)) | ((y1 < pymin) & (y2 < pymin)))
	    return 0;

	if (y1 > pymax) {
	    y1 = pymax;
	    x1 = (e->c - y1) / e->a;
	}
	if (y1 < pymin) {
	    y1 = pymin;
	    x1 = (e->c - y1) / e->a;
	}
	if (y2 > pymax) {
	    y2 = pymax;
	    x2 = (e->c - y2) / e->a;
	}
	if (y2 < pymin) {
	    y2 = pymin;
	    x2 = (e->c - y2) / e->a;
	}
    }

    line(x1, y1, x2, y2);

    return 0;
}