static GtsSurface * happrox_list (GSList * points,
gboolean keep_enclosing,
gboolean closed,
CostFunc cost_func,
gpointer cost_data,
GtsStopFunc stop_func,
gpointer stop_data)
{
GtsSurface * s = gts_surface_new (gts_surface_class (),
GTS_FACE_CLASS (list_face_class ()),
gts_edge_class (),
gts_vertex_class ());
GtsTriangle * t;
GtsVertex * w1, * w2, * w3;
GtsListFace * f;
/* creates enclosing triangle */
t = gts_triangle_enclosing (gts_triangle_class (), points, 10.);
gts_triangle_vertices (t, &w1, &w2, &w3);
GTS_POINT (w1)->z = GTS_POINT (w2)->z = GTS_POINT (w3)->z =
keep_enclosing ? -10. : -1e30;
f = GTS_LIST_FACE (gts_face_new (s->face_class, t->e1, t->e2, t->e3));
gts_surface_add_face (s, GTS_FACE (f));
f->points = points;
/* refine surface */
surface_hf_refine (s, cost_func, cost_data, stop_func, stop_data);
/* destroy unused vertices */
gts_surface_foreach_face (s, (GtsFunc) destroy_unused, NULL);
/* destroy enclosing triangle */
if (!keep_enclosing) {
gts_allow_floating_vertices = TRUE;
gts_object_destroy (GTS_OBJECT (w1));
gts_object_destroy (GTS_OBJECT (w2));
gts_object_destroy (GTS_OBJECT (w3));
gts_allow_floating_vertices = FALSE;
}
else if (closed) {
GSList * l = gts_surface_boundary (s);
GtsFace * f;
g_assert (g_slist_length (l) == 3);
f = gts_face_new (s->face_class, l->data, l->next->data, l->next->next->data);
gts_surface_add_face (s, f);
if (!gts_face_is_compatible (f, s))
gts_triangle_revert (GTS_TRIANGLE (f));
g_slist_free (l);
gts_object_destroy (GTS_OBJECT (t));
}
else
gts_object_destroy (GTS_OBJECT (t));
return s;
}
static PyObject*
is_compatible(PygtsFace *self, PyObject *args)
{
PyObject *o1_=NULL;
GtsEdge *e=NULL;
PygtsTriangle *t=NULL;
PygtsSurface *s=NULL;
SELF_CHECK
/* Parse the args */
if(! PyArg_ParseTuple(args, "O", &o1_) )
return NULL;
/* Convert to PygtsObjects */
if( pygts_triangle_check(o1_) ) {
t = PYGTS_TRIANGLE(o1_);
}
else {
if( pygts_surface_check(o1_) ) {
s = PYGTS_SURFACE(o1_);
}
else {
PyErr_SetString(PyExc_TypeError, "expected a Triangle or Surface");
return NULL;
}
}
if(t!=NULL) {
if( (e = gts_triangles_common_edge(PYGTS_TRIANGLE_AS_GTS_TRIANGLE(self),
PYGTS_TRIANGLE_AS_GTS_TRIANGLE(t)))
== NULL ) {
PyErr_SetString(PyExc_RuntimeError, "Faces do not share common edge");
return NULL;
}
if(gts_triangles_are_compatible(PYGTS_TRIANGLE_AS_GTS_TRIANGLE(self),
PYGTS_TRIANGLE_AS_GTS_TRIANGLE(t),
e)==TRUE) {
Py_INCREF(Py_True);
return Py_True;
}
}
else {
if(gts_face_is_compatible(PYGTS_FACE_AS_GTS_FACE(self),
PYGTS_SURFACE_AS_GTS_SURFACE(s))==TRUE) {
Py_INCREF(Py_True);
return Py_True;
}
}
Py_INCREF(Py_False);
return Py_False;
}
